22 August 2022 01:20:04 PM r8gd_test(): C version Test r8gd(). R8GD_DIF2_TEST R8GD_DIF2 sets up an R8GD second difference matrix; Matrix order N = 5 Number of diagonals NDIAG = 3 The offset vector: 0: -1 1: 0 2: 1 The R8GD second difference matrix: Col: 1 2 3 4 5 Row --- 1 2 -1 0 0 0 2 -1 2 -1 0 0 3 0 -1 2 -1 0 4 0 0 -1 2 -1 5 0 0 0 -1 2 R8GD_INDICATOR_TEST R8GD_INDICATOR sets up an R8GD indicator matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 0: -2 1: 0 2: 1 3: 3 The R8GD indicator matrix: Col: 1 2 3 4 5 Row --- 1 101 102 0 104 0 2 0 202 203 0 205 3 301 0 303 304 0 4 0 402 0 404 405 5 0 0 503 0 505 6 0 0 0 604 0 7 0 0 0 0 705 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 306 0 0 0 0 4 0 407 0 0 0 5 506 0 508 0 0 6 606 607 0 609 0 7 0 707 708 0 710 8 806 0 808 809 0 9 0 907 0 909 910 10 0 0 1008 0 1010 R8GD_MTV_TEST R8GD_MTV computes A'*x where A is an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 0: -2 1: 0 2: 1 3: 9 The raw matrix: Col: 1 2 3 4 Row --- 1 0 0.218418 0.956318 0.829509 2 0 0.561695 0.415307 0 3 0.0661187 0.257578 0.109957 0 4 0.043829 0.633966 0.0617272 0 5 0.449539 0.401306 0.754673 0 6 0.797287 0.00183837 0.897504 0 7 0.350752 0.0945448 0.0136169 0 8 0.859097 0.840847 0.123104 0 9 0.00751236 0.260303 0.912484 0 10 0.113664 0.351629 0 0 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0 0 0 2 0 0.561695 0.415307 0 0 3 0.0661187 0 0.257578 0.109957 0 4 0 0.043829 0 0.633966 0.0617272 5 0 0 0.449539 0 0.401306 6 0 0 0 0.797287 0 7 0 0 0 0 0.350752 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0.829509 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0.754673 0 0 0 0 6 0.00183837 0.897504 0 0 0 7 0 0.0945448 0.0136169 0 0 8 0.859097 0 0.840847 0.123104 0 9 0 0.00751236 0 0.260303 0.912484 10 0 0 0.113664 0 0.351629 A' * x: 0 0.416775 1 2.255024 2 3.851042 3 7.649455 4 4.708707 5 10.657173 6 6.114449 7 7.958738 8 3.327558 9 12.558149 R8GD_MV_TEST R8GD_MV computes A * x where A is an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 0: -2 1: 0 2: 1 3: 9 The raw matrix: Col: 1 2 3 4 Row --- 1 0 0.218418 0.956318 0.829509 2 0 0.561695 0.415307 0 3 0.0661187 0.257578 0.109957 0 4 0.043829 0.633966 0.0617272 0 5 0.449539 0.401306 0.754673 0 6 0.797287 0.00183837 0.897504 0 7 0.350752 0.0945448 0.0136169 0 8 0.859097 0.840847 0.123104 0 9 0.00751236 0.260303 0.912484 0 10 0.113664 0.351629 0 0 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0 0 0 2 0 0.561695 0.415307 0 0 3 0.0661187 0 0.257578 0.109957 0 4 0 0.043829 0 0.633966 0.0617272 5 0 0 0.449539 0 0.401306 6 0 0 0 0.797287 0 7 0 0 0 0 0.350752 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0.829509 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0.754673 0 0 0 0 6 0.00183837 0.897504 0 0 0 7 0 0.0945448 0.0136169 0 0 8 0.859097 0 0.840847 0.123104 0 9 0 0.00751236 0 0.260303 0.912484 10 0 0 0.113664 0 0.351629 A * x: 0 10.426146 1 2.369312 2 1.278679 3 2.932157 4 7.883189 5 9.482706 6 2.524510 7 12.989296 8 11.520151 9 4.425599 R8GD_PRINT_TEST R8GD_PRINT prints an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 0: -2 1: 0 2: 1 3: 3 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0 0.829509 0 2 0 0.561695 0.415307 0 0.0661187 3 0.257578 0 0.109957 0.043829 0 4 0 0.0617272 0 0.449539 0.401306 5 0 0 0.797287 0 0.00183837 6 0 0 0 0.0945448 0 7 0 0 0 0 0.123104 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 0.633966 0 0 0 0 4 0 0.754673 0 0 0 5 0.897504 0 0.350752 0 0 6 0.0136169 0.859097 0 0.840847 0 7 0 0.00751236 0.260303 0 0.912484 8 0.113664 0 0.351629 0.822887 0 9 0 0.267132 0 0.692066 0.561662 10 0 0 0.861216 0 0.453794 R8GD_PRINT_SOME_TEST R8GD_PRINT_SOME prints some of an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 0: -2 1: 0 2: 1 3: 3 Rows 3-6, Cols 3-6: Col: 3 4 5 6 Row --- 3 0.109957 0.043829 0 0.633966 4 0 0.449539 0.401306 0 5 0.797287 0 0.00183837 0.897504 6 0 0.0945448 0 0.0136169 R8GD_RANDOM_TEST R8GD_RANDOM randomly generates an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 0: -2 1: 0 2: 1 3: 9 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0 0 0 2 0 0.561695 0.415307 0 0 3 0.0661187 0 0.257578 0.109957 0 4 0 0.043829 0 0.633966 0.0617272 5 0 0 0.449539 0 0.401306 6 0 0 0 0.797287 0 7 0 0 0 0 0.350752 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0.829509 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0.754673 0 0 0 0 6 0.00183837 0.897504 0 0 0 7 0 0.0945448 0.0136169 0 0 8 0.859097 0 0.840847 0.123104 0 9 0 0.00751236 0 0.260303 0.912484 10 0 0 0.113664 0 0.351629 R8GD_TO_R8GE_TEST R8GD_TO_R8GE converts an R8GD matrix to R8GE format; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 0: -2 1: 0 2: 1 3: 3 The R8GD matrix: Col: 1 2 3 4 5 Row --- 1 101 102 0 104 0 2 0 202 203 0 205 3 301 0 303 304 0 4 0 402 0 404 405 5 0 0 503 0 505 6 0 0 0 604 0 7 0 0 0 0 705 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 306 0 0 0 0 4 0 407 0 0 0 5 506 0 508 0 0 6 606 607 0 609 0 7 0 707 708 0 710 8 806 0 808 809 0 9 0 907 0 909 910 10 0 0 1008 0 1010 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 101 102 0 104 0 2 0 202 203 0 205 3 301 0 303 304 0 4 0 402 0 404 405 5 0 0 503 0 505 6 0 0 0 604 0 7 0 0 0 0 705 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 306 0 0 0 0 4 0 407 0 0 0 5 506 0 508 0 0 6 606 607 0 609 0 7 0 707 708 0 710 8 806 0 808 809 0 9 0 907 0 909 910 10 0 0 1008 0 1010 R8GD_ZEROS_TEST R8GD_ZEROS zeros an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 0: -2 1: 0 2: 1 3: 3 The R8GD zero matrix: Col: 1 2 3 4 5 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0 0 0 0 0 6 0 0 0 0 0 7 0 0 0 0 0 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0 0 0 0 0 6 0 0 0 0 0 7 0 0 0 0 0 8 0 0 0 0 0 9 0 0 0 0 0 10 0 0 0 0 0 r8gd_test(): Normal end of execution. 22 August 2022 01:20:04 PM