11 April 2023 10:15:32.126 AM 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.218418 0.661187E-01 0.956318 0.257578 0.829509 2 0.617272E-01 0.183837E-02 0.449539 0.897504 0.401306 3 0.561695 0.438290E-01 0.218418 0.661187E-01 0.956318 4 0.754673 0.945448E-01 0.617272E-01 0.183837E-02 0.449539 5 0.415307 0.633966 0.561695 0.438290E-01 0.218418 6 0.797287 0.136169E-01 0.754673 0.945448E-01 0.617272E-01 Col: 6 Row --- 1 0.109957 2 0.350752 3 0.257578 4 0.897504 5 0.661187E-01 6 0.183837E-02 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. 11 April 2023 10:15:32.126 AM