09 May 2025 9:29:09.330 PM r8bto_test(): FORTRAN90 version: Test r8bto(). R8BTO_DIF2_TEST R8BTO_DIF2 sets up an R8BTO version of the second difference matrix (assuming M = 1 ). Block order M = 1 Block number L = 5 Matrix order N = 5 The R8BTO second difference matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 0.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 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 R8BTO_INDICATOR_TEST R8BTO_INDICATOR sets up an indicator matrix for a real block Toeplitz matrix, Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 21.0000 22.0000 23.0000 24.0000 25.0000 3 31.0000 32.0000 11.0000 12.0000 13.0000 4 41.0000 42.0000 21.0000 22.0000 23.0000 5 51.0000 52.0000 31.0000 32.0000 11.0000 6 61.0000 62.0000 41.0000 42.0000 21.0000 Col: 6 Row --- 1 16.0000 2 26.0000 3 14.0000 4 24.0000 5 12.0000 6 22.0000 R8BTO_MTV_TEST R8BTO_MTV computes A'* x for a block Toeplitz matrix. Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 5.00000 5.00000 6.00000 6.00000 7.00000 3 7.00000 8.00000 1.00000 2.00000 3.00000 4 8.00000 8.00000 5.00000 5.00000 6.00000 5 9.00000 0.00000 7.00000 8.00000 1.00000 6 9.00000 9.00000 8.00000 8.00000 5.00000 Col: 6 Row --- 1 6.00000 2 7.00000 3 4.00000 4 6.00000 5 2.00000 6 5.00000 The "vector" x: Col: 1 2 3 Row --- 1 11.0000 12.0000 13.0000 2 21.0000 22.0000 23.0000 The product A'*x: Col: 1 2 3 Row --- 1 700.000 556.000 498.000 2 606.000 592.000 534.000 R8BTO_MV_TEST R8BTO_MV computes A * x for a block Toeplitz matrix. Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 5.00000 5.00000 6.00000 6.00000 7.00000 3 7.00000 8.00000 1.00000 2.00000 3.00000 4 8.00000 8.00000 5.00000 5.00000 6.00000 5 9.00000 0.00000 7.00000 8.00000 1.00000 6 9.00000 9.00000 8.00000 8.00000 5.00000 Col: 6 Row --- 1 6.00000 2 7.00000 3 4.00000 4 6.00000 5 2.00000 6 5.00000 The "vector" x: Col: 1 2 3 Row --- 1 11.0000 12.0000 13.0000 2 21.0000 22.0000 23.0000 The product A*x: Col: 1 2 3 Row --- 1 380.000 432.000 418.000 2 616.000 642.000 740.000 R8BTO_PRINT_TEST R8BTO_PRINT prints an R8BTO matrix, Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 21.0000 22.0000 23.0000 24.0000 25.0000 3 31.0000 32.0000 11.0000 12.0000 13.0000 4 41.0000 42.0000 21.0000 22.0000 23.0000 5 51.0000 52.0000 31.0000 32.0000 11.0000 6 61.0000 62.0000 41.0000 42.0000 21.0000 Col: 6 Row --- 1 16.0000 2 26.0000 3 14.0000 4 24.0000 5 12.0000 6 22.0000 R8BTO_PRINT_SOME_TEST R8BTO_PRINT_SOME prints some of an R8BTO matrix, Block order M = 2 Block number L = 3 Matrix order N = 6 Row (1:6), Cols (3:4): Col: 3 4 Row --- 1 13.0000 14.0000 2 23.0000 24.0000 3 11.0000 12.0000 4 21.0000 22.0000 5 31.0000 32.0000 6 41.0000 42.0000 R8BTO_RANDOM_TEST R8BTO_RANDOM returns a random R8BTO matrix. Block order M = 2 Block number L = 3 Matrix order N = 6 The random R8BTO matrix: Col: 1 2 3 4 5 Row --- 1 0.223053 0.288186 0.909550 0.413598 0.969371 2 0.966337 0.729953 0.771159 0.693302 0.282265 3 0.412330 0.745685E-01 0.223053 0.288186 0.909550 4 0.981536 0.936296 0.966337 0.729953 0.771159 5 0.457951 0.133943 0.412330 0.745685E-01 0.223053 6 0.826796 0.877936 0.981536 0.936296 0.966337 Col: 6 Row --- 1 0.637015 2 0.691158 3 0.413598 4 0.693302 5 0.288186 6 0.729953 R8BTO_SL_TEST R8BTO_SL solves a block Toeplitz system. Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 9.00000 1.00000 3.00000 4.00000 5.00000 2 2.00000 8.00000 6.00000 6.00000 7.00000 3 7.00000 8.00000 9.00000 1.00000 3.00000 4 8.00000 8.00000 2.00000 8.00000 6.00000 5 9.00000 0.00000 7.00000 8.00000 9.00000 6 9.00000 9.00000 8.00000 8.00000 2.00000 Col: 6 Row --- 1 6.00000 2 7.00000 3 4.00000 4 6.00000 5 1.00000 6 8.00000 The right hand side B: 1 97.0000 2 137.000 3 93.0000 4 128.000 5 113.000 6 141.000 The computed solution X: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 R8BTO_TO_R8GE_TEST R8BTO_TO_R8GE converts an R8BTO matrix to R8GE format. Block order M = 2 Block number L = 3 Matrix order N = 6 The R8BTO matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 21.0000 22.0000 23.0000 24.0000 25.0000 3 31.0000 32.0000 11.0000 12.0000 13.0000 4 41.0000 42.0000 21.0000 22.0000 23.0000 5 51.0000 52.0000 31.0000 32.0000 11.0000 6 61.0000 62.0000 41.0000 42.0000 21.0000 Col: 6 Row --- 1 16.0000 2 26.0000 3 14.0000 4 24.0000 5 12.0000 6 22.0000 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 21.0000 22.0000 23.0000 24.0000 25.0000 3 31.0000 32.0000 11.0000 12.0000 13.0000 4 41.0000 42.0000 21.0000 22.0000 23.0000 5 51.0000 52.0000 31.0000 32.0000 11.0000 6 61.0000 62.0000 41.0000 42.0000 21.0000 Col: 6 Row --- 1 16.0000 2 26.0000 3 14.0000 4 24.0000 5 12.0000 6 22.0000 R8BTO_ZEROS_TEST R8BTO_ZEROS zeros an R8BTO matrix. Block order M = 2 Block number L = 3 Matrix order N = 6 The zero R8BTO matrix: Col: 1 2 3 4 5 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 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 Row --- 1 0.00000 2 0.00000 3 0.00000 4 0.00000 5 0.00000 6 0.00000 r8bto_test(): Normal end of execution. 09 May 2025 9:29:09.331 PM