11 April 2023 10:19:18.872 AM r8blt_test(): FORTRAN90 version: Test r8blt(). R8BLT_DET_TEST R8BLT_DET computes the determinant of an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 2 0.101333-319 0.00000 3 0.00000 0.00000 0.467823-309 4 0.00000 0.00000 0.692679-309 0.494066-323 5 0.00000 0.00000 0.467822-309 0.00000 6 0.467823-309 0.695331-309 0.692679-309 7 0.00000 0.695331-309 8 0.692679-309 Col: 6 7 8 9 10 Row --- 6 0.00000 7 0.138524-308 0.692679-309 8 0.205470E-80 0.695331-309 0.695331-309 9 -0.280776E-49 0.00000 0.692679-309 0.695331-309 10 0.467822-309 0.741098-322 0.695331-309 0.695331-309 Determinant = 0.00000 R8BLT_INDICATOR_TEST R8BLT_INDICATOR sets up an R8BLT indicator matrix. Matrix order N = 6 Lower bandwidth ML = 2 The R8BLT indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 2 21.0000 22.0000 3 31.0000 32.0000 33.0000 4 42.0000 43.0000 44.0000 5 53.0000 54.0000 55.0000 6 64.0000 65.0000 Col: 6 Row --- 6 66.0000 R8BLT_MTV_TEST R8BLT_MTV computes b=A'*x, where A is an R8BLT matrix. Matrix order N = 5 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 2 0.138524-308 31.0000 3 11.0000 22.0000 33.0000 4 21.0000 32.0000 43.0000 54.0000 5 42.0000 53.0000 64.0000 0.00000 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b=A'*x: 1 117.000 2 466.000 3 536.000 4 536.000 5 0.00000 R8BLT_MV_TEST R8BLT_MV computes b=A*x, where A is an R8BLT matrix. Matrix order N = 5 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 2 0.138524-308 31.0000 3 11.0000 22.0000 33.0000 4 21.0000 32.0000 43.0000 54.0000 5 42.0000 53.0000 64.0000 0.00000 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b=A*x: 1 0.00000 2 62.0000 3 154.000 4 430.000 5 499.000 R8BLT_PRINT_TEST R8BLT_PRINT prints an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 2 0.00000 0.467823-309 3 0.00000 0.00000 1.00000 4 0.00000 0.00000 2.00000 5.00000 5 0.00000 3.00000 0.00000 154.000 6 4.00000 0.00000 430.000 7 62.0000 499.000 8 0.692679-309 Col: 6 7 8 9 10 Row --- 6 0.00000 7 0.138524-308 31.0000 8 11.0000 22.0000 33.0000 9 21.0000 32.0000 43.0000 54.0000 10 42.0000 53.0000 64.0000 0.00000 R8BLT_PRINT_SOME_TEST R8BLT_PRINT_SOME prints some of an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 Rows 2:5, Cols 3:5: Col: 3 4 5 Row --- 3 303.000 4 403.000 404.000 5 503.000 504.000 505.000 R8BLT_RANDOM_TEST R8BLT_RANDOM randomizes an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 101.000 2 201.000 202.000 3 301.000 302.000 303.000 4 401.000 402.000 403.000 404.000 5 502.000 503.000 504.000 505.000 6 603.000 604.000 605.000 7 704.000 705.000 8 805.000 Col: 6 7 8 9 10 Row --- 6 606.000 7 706.000 707.000 8 806.000 807.000 808.000 9 906.000 907.000 908.000 909.000 10 1007.00 1008.00 1009.00 1010.00 R8BLT_SL_TEST R8BLT_SL solves A*x=b where A is an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 101.000 2 201.000 202.000 3 301.000 302.000 303.000 4 401.000 402.000 403.000 404.000 5 502.000 503.000 504.000 505.000 6 603.000 604.000 605.000 7 704.000 705.000 8 805.000 Col: 6 7 8 9 10 Row --- 6 606.000 7 706.000 707.000 8 806.000 807.000 808.000 9 906.000 907.000 908.000 909.000 10 1007.00 1008.00 1009.00 1010.00 The right hand side: 1 101.000 2 605.000 3 1814.00 4 4030.00 5 7054.00 6 10886.0 7 15526.0 8 20974.0 9 27230.0 10 34294.0 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 R8BLT_SLT_TEST R8BLT_SLT solves A'*x=b where A is an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 101.000 2 201.000 202.000 3 301.000 302.000 303.000 4 401.000 402.000 403.000 404.000 5 502.000 503.000 504.000 505.000 6 603.000 604.000 605.000 7 704.000 705.000 8 805.000 Col: 6 7 8 9 10 Row --- 6 606.000 7 706.000 707.000 8 806.000 807.000 808.000 9 906.000 907.000 908.000 909.000 10 1007.00 1008.00 1009.00 1010.00 The right hand side: 1 3010.00 2 5428.00 3 8654.00 4 12688.0 5 17530.0 6 23180.0 7 29638.0 8 24716.0 9 18271.0 10 10100.0 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 R8BLT_TO_R8GE_TEST R8BLT_TO_R8GE converts a matrix from R8BLT to R8GE format. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 101.000 2 201.000 202.000 3 301.000 302.000 303.000 4 401.000 402.000 403.000 404.000 5 502.000 503.000 504.000 505.000 6 603.000 604.000 605.000 7 704.000 705.000 8 805.000 Col: 6 7 8 9 10 Row --- 6 606.000 7 706.000 707.000 8 806.000 807.000 808.000 9 906.000 907.000 908.000 909.000 10 1007.00 1008.00 1009.00 1010.00 The R8GE matrix Col: 1 2 3 4 5 Row --- 1 101.000 0.00000 0.00000 0.00000 0.00000 2 201.000 202.000 0.00000 0.00000 0.00000 3 301.000 302.000 303.000 0.00000 0.00000 4 401.000 402.000 403.000 404.000 0.00000 5 0.00000 502.000 503.000 504.000 505.000 6 0.00000 0.00000 603.000 604.000 605.000 7 0.00000 0.00000 0.00000 704.000 705.000 8 0.00000 0.00000 0.00000 0.00000 805.000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 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 0.00000 5 0.00000 0.00000 0.00000 0.00000 0.00000 6 606.000 0.00000 0.00000 0.00000 0.00000 7 706.000 707.000 0.00000 0.00000 0.00000 8 806.000 807.000 808.000 0.00000 0.00000 9 906.000 907.000 908.000 909.000 0.00000 10 0.00000 1007.00 1008.00 1009.00 1010.00 R8BLT_ZEROS_TEST R8BLT_ZEROS zeros an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 2 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 0.00000 6 0.00000 0.00000 0.00000 7 0.00000 0.00000 8 0.00000 Col: 6 7 8 9 10 Row --- 6 0.00000 7 0.00000 0.00000 8 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 r8blt_test(): Normal end of execution. 11 April 2023 10:19:18.873 AM