Fri Aug 19 17:25:07 2022 r8bto_test(): Python version: 3.6.9 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: 0 1 2 3 4 Row --- 0 : 2 -1 0 0 0 1 : -1 2 -1 0 0 2 : 0 -1 2 -1 0 3 : 0 0 -1 2 -1 4 : 0 0 0 -1 2 R8BTO_INDICATOR_TEST R8BTO_INDICATOR sets up an R8BTO indicator matrix Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 0 1 2 3 4 Row --- 0 : 11 12 13 14 15 1 : 21 22 23 24 25 2 : 40 41 11 12 13 3 : 50 51 21 22 23 4 : 60 61 40 41 11 5 : 70 71 50 51 21 Col: 5 Row --- 0 : 16 1 : 26 2 : 14 3 : 24 4 : 12 5 : 22 R8BTO_MTV_TEST R8BTO_MTV computes A'* x for an R8BTO matrix. Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 0 1 2 3 4 Row --- 0 : 1 2 3 4 5 1 : 5 5 6 6 7 2 : 7 8 1 2 3 3 : 8 8 5 5 6 4 : 9 0 7 8 1 5 : 9 9 8 8 5 Col: 5 Row --- 0 : 6 1 : 7 2 : 4 3 : 6 4 : 2 5 : 5 The matrix x: [[11. 12. 13. 14. 15. 16.] [21. 22. 23. 24. 25. 26.]] The product A'*x: [[430. 448. 364.] [756. 650. 608.]] R8BTO_MV_TEST R8BTO_MV computes A * x. Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 0 1 2 3 4 Row --- 0 : 1 2 3 4 5 1 : 5 5 6 6 7 2 : 7 8 1 2 3 3 : 8 8 5 5 6 4 : 9 0 7 8 1 5 : 9 9 8 8 5 Col: 5 Row --- 0 : 6 1 : 7 2 : 4 3 : 6 4 : 2 5 : 5 The matrix x: [[11. 12. 13. 14. 15. 16.] [21. 22. 23. 24. 25. 26.]] The product A*x: [[430. 448. 364.] [756. 650. 608.]] R8BTO_PRINT_TEST R8BTO_PRINT prints an R8BTO matrix Block order M = 2 Block number L = 3 Matrix order N = 6 The R8BTO matrix: Col: 0 1 2 3 4 Row --- 0 : 11 12 13 14 15 1 : 21 22 23 24 25 2 : 40 41 11 12 13 3 : 50 51 21 22 23 4 : 60 61 40 41 11 5 : 70 71 50 51 21 Col: 5 Row --- 0 : 16 1 : 26 2 : 14 3 : 24 4 : 12 5 : 22 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 : 24 25 2 : 12 13 3 : 22 23 4 : 41 11 5 : 51 21 R8BTO_RANDOM_TEST R8BTO_RANDOM returns a random R8BTO matrix Block order M = 2 Block number L = 3 Matrix order N = 6 The random RBTO matrix: Col: 0 1 2 3 4 Row --- 0 : 0.509711 0.378267 0.488665 0.605465 0.813057 1 : 0.0670433 0.376888 0.269608 0.23133 0.240986 2 : 0.385218 0.559602 0.509711 0.378267 0.488665 3 : 0.379579 0.20132 0.0670433 0.376888 0.269608 4 : 0.76544 0.928624 0.385218 0.559602 0.509711 5 : 0.666475 0.123407 0.379579 0.20132 0.0670433 Col: 5 Row --- 0 : 0.943707 1 : 0.123308 2 : 0.605465 3 : 0.23133 4 : 0.378267 5 : 0.376888 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: 0 1 2 3 4 Row --- 0 : 11 12 13 14 15 1 : 21 22 23 24 25 2 : 40 41 11 12 13 3 : 50 51 21 22 23 4 : 60 61 40 41 11 5 : 70 71 50 51 21 Col: 5 Row --- 0 : 16 1 : 26 2 : 14 3 : 24 4 : 12 5 : 22 The R8GE matrix: [[11. 12. 13. 14. 15. 16.] [21. 22. 23. 24. 25. 26.] [40. 41. 11. 12. 13. 14.] [50. 51. 21. 22. 23. 24.] [60. 61. 40. 41. 11. 12.] [70. 71. 50. 51. 21. 22.]] R8BTO_ZEROS_TEST R8BTO_ZEROS zeros an R8BTO matrix. Block order M = 3 Block number L = 2 Matrix order N = 6 The zero R8BTO matrix: Col: 0 1 2 3 4 Row --- 0 : 0 0 0 0 0 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 Col: 5 Row --- 0 : 0 1 : 0 2 : 0 3 : 0 4 : 0 5 : 0 r8bto_test(): Normal end of execution. Fri Aug 19 17:25:07 2022