19 July 2016 09:21:18 PM R8GD_PRB C++ version Test the R8GD library. R8GD_DIF2_TEST R8GD_DIF2 sets up an R8GD second difference matrix; Matrix order N = 5 Number of diagonals NDIAG = 3 The offset vector: 1 -1 2 0 3 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: 1 -2 2 0 3 1 4 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 for an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 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: 1 0.416775 2 2.25502 3 3.85104 4 7.64945 5 4.70871 6 10.6572 7 6.11445 8 7.95874 9 3.32756 10 12.5581 R8GD_MV_TEST R8GD_MV computes A * x for an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 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: 1 10.4261 2 2.36931 3 1.27868 4 2.93216 5 7.88319 6 9.48271 7 2.52451 8 12.9893 9 11.5202 10 4.4256 R8GD_PRINT_TEST R8GD_PRINT prints an R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 3 The R8GD 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: 1 -2 2 0 3 1 4 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; R8GD_MTV computes A'*x; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 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 = 5 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 3 The R8GD matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 14 0 2 0 22 23 0 25 3 31 0 33 34 0 4 0 42 0 44 45 5 0 0 53 0 55 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 14 0 2 0 22 23 0 25 3 31 0 33 34 0 4 0 42 0 44 45 5 0 0 53 0 55 R8GD_ZEROS_TEST R8GD_ZEROS zeros an R8GD matrix; Matrix order N = 5 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 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 R8GD_PRB Normal end of execution. 19 July 2016 09:21:18 PM