11 April 2023 8:59:42.101 AM r8cb_test(): FORTRAN90 version: Test r8cb(). R8CB_DET_TEST R8CB_DET computes the determinant of a matrix factored by R8CB_NP_FA. Matrix order N = 10 Lower bandwidth ML = 2 Upper bandwidth MU = 3 The compact band matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.561695 0.109957 0.401306 2 0.956318 0.415307 0.438290E-01 0.754673 0.945448E-01 3 0.829509 0.661187E-01 0.633966 0.797287 0.136169E-01 4 0.257578 0.617272E-01 0.183837E-02 0.859097 5 0.449539 0.897504 0.840847 6 0.350752 0.123104 7 0.751236E-02 Col: 6 7 8 9 10 Row --- 3 0.260303 4 0.912484 0.692066 5 0.113664 0.561662 0.188955 6 0.351629 0.861216 0.761492 0.617205 7 0.822887 0.453794 0.396988 0.361529 0.299329 8 0.267132 0.911977 0.185314 0.212930 0.825003 9 0.597917 0.574366 0.714471 0.824660 10 0.367027 0.117707 0.618618E-01 R8CB_DET computes the determinant = 0.183588E-01 R8GE_DET computes the determinant = 0.183588E-01 R8CB_DIF2_TEST R8CB_DIF2 computes an R8CB second difference matrix Matrix rows M = 8 Matrix columns N = 10 Lower bandwidth ML = 2 Upper bandwidth MU = 3 The R8CB second difference matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 -1.00000 2.00000 6 0.00000 -1.00000 7 0.00000 Col: 6 7 8 9 10 Row --- 3 0.00000 4 0.00000 0.00000 5 -1.00000 0.00000 0.00000 6 2.00000 -1.00000 0.00000 0.00000 7 -1.00000 2.00000 -1.00000 0.00000 0.00000 8 0.00000 -1.00000 2.00000 -1.00000 0.00000 R8CB_INDICATOR_TEST R8CB_INDICATOR computes an R8CB indicator matrix Matrix rows M = 8 Matrix columns N = 10 Lower bandwidth ML = 2 Upper bandwidth MU = 3 The R8CB indicator matrix: Col: 1 2 3 4 5 Row --- 1 101.000 102.000 103.000 104.000 2 201.000 202.000 203.000 204.000 205.000 3 301.000 302.000 303.000 304.000 305.000 4 402.000 403.000 404.000 405.000 5 503.000 504.000 505.000 6 604.000 605.000 7 705.000 Col: 6 7 8 9 10 Row --- 3 306.000 4 406.000 407.000 5 506.000 507.000 508.000 6 606.000 607.000 608.000 609.000 7 706.000 707.000 708.000 709.000 710.000 8 806.000 807.000 808.000 809.000 810.000 R8CB_ML_TEST R8CB_ML computes A*x or A'*x for a compact band matrix A which has been factored by R8CB_FA. Matrix order N = 10 Lower bandwidth ML = 1 Upper bandwidth MU = 2 A*x and PLU*x 1 2.07579 2.07579 2 5.38830 5.38830 3 5.18076 5.18076 4 8.02059 8.02059 5 8.06762 8.06762 6 12.7308 12.7308 7 8.52067 8.52067 8 15.4017 15.4017 9 20.7578 20.7578 10 13.7301 13.7301 A'*x and (PLU)'*x 1 1.71194 1.71194 2 2.49071 2.49071 3 3.53774 3.53774 4 7.10439 7.10439 5 6.80503 6.80503 6 14.9862 14.9862 7 18.7469 18.7469 8 7.88328 7.88328 9 12.7733 12.7733 10 2.65418 2.65418 R8CB_MTV_TEST R8CB_MTV computes b=A'*x, where A is an R8CB matrix. Matrix rows M = 8 Matrix columns N = 8 Lower bandwidth ML = 2 Upper bandwidth MU = 1 The R8CB matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 31.0000 32.0000 33.0000 34.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 6 64.0000 65.0000 7 75.0000 Col: 6 7 8 Row --- 5 56.0000 6 66.0000 67.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 The vector x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 The product b=A'*x: 1 146.000 2 320.000 3 582.000 4 932.000 5 1370.00 6 1896.00 7 1637.00 8 1250.00 R8CB_MV_TEST R8CB_MV computes b=A*x, where A is an R8CB matrix. Matrix rows M = 8 Matrix columns N = 8 Lower bandwidth ML = 2 Upper bandwidth MU = 1 The R8CB matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 31.0000 32.0000 33.0000 34.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 6 64.0000 65.0000 7 75.0000 Col: 6 7 8 Row --- 5 56.0000 6 66.0000 67.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 The vector x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 The product b=A*x: 1 35.0000 2 134.000 3 330.000 4 614.000 5 986.000 6 1446.00 7 1994.00 8 1829.00 R8CB_NP_FA_TEST R8CB_NP_FA factors an R8CB matrix, no pivoting: Matrix order N = 10 Lower bandwidth ML = 1 Upper bandwidth MU = 2 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution to transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8CB_NP_SL_TEST R8CB_NP_SL solves a linear system factored by R8CB_NP_FA. Matrix order N = 10 Lower bandwidth ML = 1 Upper bandwidth MU = 2 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution to transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8CB_PRINT_TEST R8CB_PRINT prints an R8CB matrix Matrix rows M = 8 Matrix columns N = 10 Lower bandwidth ML = 2 Upper bandwidth MU = 3 The R8CB matrix: Col: 1 2 3 4 5 Row --- 1 101.000 102.000 103.000 104.000 2 201.000 202.000 203.000 204.000 205.000 3 301.000 302.000 303.000 304.000 305.000 4 402.000 403.000 404.000 405.000 5 503.000 504.000 505.000 6 604.000 605.000 7 705.000 Col: 6 7 8 9 10 Row --- 3 306.000 4 406.000 407.000 5 506.000 507.000 508.000 6 606.000 607.000 608.000 609.000 7 706.000 707.000 708.000 709.000 710.000 8 806.000 807.000 808.000 809.000 810.000 R8CB_PRINT_SOME_TEST R8CB_PRINT_SOME prints some of an R8CB matrix Matrix rows M = 8 Matrix columns N = 10 Lower bandwidth ML = 2 Upper bandwidth MU = 3 Rows 3-6, Cols 3-6: Col: 3 4 5 6 Row --- 3 303.000 304.000 305.000 306.000 4 403.000 404.000 405.000 406.000 5 503.000 504.000 505.000 506.000 6 604.000 605.000 606.000 R8CB_RANDOM_TEST R8CB_RANDOM randomizes an R8CB matrix. Matrix rows M = 8 Matrix columns N = 10 Lower bandwidth ML = 2 Upper bandwidth MU = 3 The random R8CB matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.561695 0.109957 0.401306 2 0.956318 0.415307 0.438290E-01 0.754673 0.945448E-01 3 0.829509 0.661187E-01 0.633966 0.797287 0.136169E-01 4 0.257578 0.617272E-01 0.183837E-02 0.859097 5 0.449539 0.897504 0.840847 6 0.350752 0.123104 7 0.751236E-02 Col: 6 7 8 9 10 Row --- 3 0.260303 4 0.912484 0.692066 5 0.113664 0.561662 0.597917 6 0.351629 0.861216 0.188955 0.185314 7 0.822887 0.453794 0.761492 0.574366 0.617205 8 0.267132 0.911977 0.396988 0.367027 0.361529 R8CB_TO_R8GE_TEST R8CB_TO_R8GE converts an R8CB matrix to R8GE format. Matrix rows M = 5 Matrix columns N = 8 Lower bandwidth ML = 2 Upper bandwidth MU = 1 The R8CB matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 31.0000 32.0000 33.0000 34.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 Col: 6 7 8 Row --- 5 56.0000 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 0.00000 2 21.0000 22.0000 23.0000 0.00000 0.00000 3 31.0000 32.0000 33.0000 34.0000 0.00000 4 0.00000 42.0000 43.0000 44.0000 45.0000 5 0.00000 0.00000 53.0000 54.0000 55.0000 Col: 6 7 8 Row --- 1 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 5 56.0000 0.00000 0.00000 R8CB_TO_R8VEC_TEST R8CB_TO_R8VEC converts an R8CB matrix to an R8VEC. Matrix rows M = 5 Matrix columns N = 8 Lower bandwidth ML = 2 Upper bandwidth MU = 1 The R8CB matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 31.0000 32.0000 33.0000 34.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 Col: 6 7 8 Row --- 5 56.0000 1 1 1 0.00000 2 1 2 11.0000 3 1 3 21.0000 4 1 4 31.0000 1 2 5 12.0000 2 2 6 22.0000 3 2 7 32.0000 4 2 8 42.0000 1 3 9 23.0000 2 3 10 33.0000 3 3 11 43.0000 4 3 12 53.0000 1 4 13 34.0000 2 4 14 44.0000 3 4 15 54.0000 4 4 16 0.00000 1 5 17 45.0000 2 5 18 55.0000 3 5 19 0.00000 4 5 20 0.00000 1 6 21 56.0000 2 6 22 0.00000 3 6 23 0.00000 4 6 24 0.00000 1 7 25 0.00000 2 7 26 0.00000 3 7 27 0.00000 4 7 28 0.00000 1 8 29 0.00000 2 8 30 0.00000 3 8 31 0.00000 4 8 32 0.00000 The recovered R8CB matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 31.0000 32.0000 33.0000 34.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 Col: 6 7 8 Row --- 5 56.0000 R8CB_ZEROS_TEST R8CB_ZEROS zeros an R8CB matrix. Matrix rows M = 8 Matrix columns N = 10 Lower bandwidth ML = 2 Upper bandwidth MU = 3 The R8CB zero matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 6 0.00000 0.00000 7 0.00000 Col: 6 7 8 9 10 Row --- 3 0.00000 4 0.00000 0.00000 5 0.00000 0.00000 0.00000 6 0.00000 0.00000 0.00000 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 R8VEC_TO_R8CB_TEST R8VEC_TO_R8CB converts an R8VEC to an R8CB matrix. Matrix rows M = 5 Matrix columns N = 8 Lower bandwidth ML = 2 Upper bandwidth MU = 1 The R8CB indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 31.0000 32.0000 33.0000 34.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 Col: 6 7 8 Row --- 5 56.0000 1 1 1 0.00000 2 1 2 11.0000 3 1 3 21.0000 4 1 4 31.0000 1 2 5 12.0000 2 2 6 22.0000 3 2 7 32.0000 4 2 8 42.0000 1 3 9 23.0000 2 3 10 33.0000 3 3 11 43.0000 4 3 12 53.0000 1 4 13 34.0000 2 4 14 44.0000 3 4 15 54.0000 4 4 16 0.00000 1 5 17 45.0000 2 5 18 55.0000 3 5 19 0.00000 4 5 20 0.00000 1 6 21 56.0000 2 6 22 0.00000 3 6 23 0.00000 4 6 24 0.00000 1 7 25 0.00000 2 7 26 0.00000 3 7 27 0.00000 4 7 28 0.00000 1 8 29 0.00000 2 8 30 0.00000 3 8 31 0.00000 4 8 32 0.00000 The recovered R8CB indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 31.0000 32.0000 33.0000 34.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 Col: 6 7 8 Row --- 5 56.0000 r8cb_test(): Normal end of execution. 11 April 2023 8:59:42.102 AM