Mon Aug 22 16:54:43 2022 r8gd_test(): Python version: 3.6.9 Test r8gd(). R8GD_DIF2_TEST R8GD_DIF2 sets up an R8GD second difference matrix. Matrix order N = 5 Number of diagonals NDIAG = 3 The offset vector: [-1 0 1] The R8GD 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 R8GD_INDICATOR_TEST R8GD_INDICATOR sets up an R8GD indicator matrix. Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: [-2 0 1 9] The general diagonal matrix: Col: 0 1 2 3 4 Row 0 : 101 102 0 0 0 1 : 0 202 203 0 0 2 : 301 0 303 304 0 3 : 0 402 0 404 405 4 : 0 0 503 0 505 5 : 0 0 0 604 0 6 : 0 0 0 0 705 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 110 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 506 0 0 0 0 5 : 606 607 0 0 0 6 : 0 707 708 0 0 7 : 806 0 808 809 0 8 : 0 907 0 909 910 9 : 0 0 1008 0 1010 R8GD_MTV_TEST R8GD_MTV computes A'*x, where A is an R8GD matrix Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: [-2 0 1 9] The general diagonal matrix: Col: 0 1 2 3 4 Row 0 : 0.374975 0.439676 0 0 0 1 : 0 0.731458 0.799019 0 0 2 : 0.772089 0 0.58984 0.0846239 0 3 : 0 0.18238 0 0.0314934 0.694309 4 : 0 0 0.364944 0 0.18902 5 : 0 0 0 0.921829 0 6 : 0 0 0 0 0.981032 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.433831 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.0978926 0 0 0 0 5 : 0.243231 0.00873172 0 0 0 6 : 0 0.328384 0.0649966 0 0 7 : 0.777328 0 0.0176033 0.707108 0 8 : 0 0.362127 0 0.891784 0.882074 9 : 0 0 0.577512 0 0.778072 A' * x: 0: 2.69124 1: 2.63211 2: 5.19228 3: 5.91082 4: 10.5896 5: 8.16747 6: 5.61022 7: 6.37092 8: 13.6829 9: 16.1532 R8GD_MV_TEST R8GD_MV computes A * x, where A is an R8GD matrix Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: [-2 0 1 9] The general diagonal matrix: Col: 0 1 2 3 4 Row 0 : 0.290032 0.960899 0 0 0 1 : 0 0.138293 0.115868 0 0 2 : 0.255814 0 0.817914 0.028612 0 3 : 0 0.901269 0 0.48794 0.294213 4 : 0 0 0.835257 0 0.139335 5 : 0 0 0 0.959494 0 6 : 0 0 0 0 0.695622 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.944539 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.990033 0 0 0 0 5 : 0.482711 0.822529 0 0 0 6 : 0 0.317735 0.102322 0 0 7 : 0.473873 0 0.403863 0.290535 0 8 : 0 0.652735 0 0.23912 0.841714 9 : 0 0 0.102746 0 0.449047 A * x: 0: 11.6572 1: 0.624188 2: 2.824 3: 5.22536 4: 9.14264 5: 12.4919 6: 6.52083 7: 8.68895 8: 15.1384 9: 5.31244 R8GD_PRINT_TEST R8GD_PRINT prints an R8GD matrix. Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: [-2 0 1 9] The R8GD matrix: Col: 0 1 2 3 4 Row 0 : 0.853072 0.264666 0 0 0 1 : 0 0.0501588 0.317042 0 0 2 : 0.155777 0 0.338486 0.784137 0 3 : 0 0.0654758 0 0.109948 0.773401 4 : 0 0 0.123364 0 0.544954 5 : 0 0 0 0.462741 0 6 : 0 0 0 0 0.979941 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.627961 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.0151172 0 0 0 0 5 : 0.983608 0.402053 0 0 0 6 : 0 0.19769 0.118285 0 0 7 : 0.404366 0 0.946367 0.749132 0 8 : 0 0.88496 0 0.864922 0.20413 9 : 0 0 0.264813 0 0.608911 R8GD_PRINT_SOME_TEST R8GD_PRINT_SOME prints some of an R8GD matrix. Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: [-2 0 1 9] Rows 3-6, Cols 3-6: Col: 3 4 5 6 Row 3 : 404 405 0 0 4 : 0 505 506 0 5 : 604 0 606 607 6 : 0 705 0 707 R8GD_RANDOM_TEST R8GD_RANDOM generates a random R8GD matrix. Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: [-2 0 1 9] The general diagonal matrix: Col: 0 1 2 3 4 Row 0 : 0.859286 0.676871 0 0 0 1 : 0 0.0384812 0.524393 0 0 2 : 0.62913 0 0.604685 0.16748 0 3 : 0 0.785607 0 0.00566174 0.926915 4 : 0 0 0.899258 0 0.895193 5 : 0 0 0 0.399245 0 6 : 0 0 0 0 0.938226 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.543417 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.771648 0 0 0 0 5 : 0.32459 0.376284 0 0 0 6 : 0 0.345027 0.807534 0 0 7 : 0.251565 0 0.553309 0.422735 0 8 : 0 0.401584 0 0.356304 0.114771 9 : 0 0 0.0429702 0 0.418288 R8GD_TO_R8GE_TEST R8GD_TO_R8GE converts an R8GD matrix to R8GE format. Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: [-2 0 1 9] The R8GD matrix: Col: 0 1 2 3 4 Row 0 : 101 102 0 0 0 1 : 0 202 203 0 0 2 : 301 0 303 304 0 3 : 0 402 0 404 405 4 : 0 0 503 0 505 5 : 0 0 0 604 0 6 : 0 0 0 0 705 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 110 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 506 0 0 0 0 5 : 606 607 0 0 0 6 : 0 707 708 0 0 7 : 806 0 808 809 0 8 : 0 907 0 909 910 9 : 0 0 1008 0 1010 The R8GE matrix: [[ 101. 102. 0. 0. 0. 0. 0. 0. 0. 110.] [ 0. 202. 203. 0. 0. 0. 0. 0. 0. 0.] [ 301. 0. 303. 304. 0. 0. 0. 0. 0. 0.] [ 0. 402. 0. 404. 405. 0. 0. 0. 0. 0.] [ 0. 0. 503. 0. 505. 506. 0. 0. 0. 0.] [ 0. 0. 0. 604. 0. 606. 607. 0. 0. 0.] [ 0. 0. 0. 0. 705. 0. 707. 708. 0. 0.] [ 0. 0. 0. 0. 0. 806. 0. 808. 809. 0.] [ 0. 0. 0. 0. 0. 0. 907. 0. 909. 910.] [ 0. 0. 0. 0. 0. 0. 0. 1008. 0. 1010.]] R8GD_ZEROS_TEST R8GD_ZEROS zeros an R8GD matrix. Matrix order N = 5 Number of diagonals NDIAG = 3 The offset vector: [-1 0 1] The zero R8GD 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 r8gd_test(): Normal end of execution. Mon Aug 22 16:54:43 2022