Tue May 20 22:27:34 2025 r8gd_test(): python version: 3.10.12 numpy version: 1.26.4 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.415165 0.243361 0 0 0 1 : 0 0.764844 0.515545 0 0 2 : 0.226735 0 0.838621 0.104074 0 3 : 0 0.776823 0 0.722692 0.571084 4 : 0 0 0.387471 0 0.559132 5 : 0 0 0 0.0442874 0 6 : 0 0 0 0 0.587651 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.379442 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.436501 0 0 0 0 5 : 0.776444 0.0997956 0 0 0 6 : 0 0.224209 0.966575 0 0 7 : 0.402654 0 0.521791 0.294316 0 8 : 0 0.673401 0 0.624081 0.992616 9 : 0 0 0.565932 0 0.0409563 A' * x: 0: 1.09537 1: 4.88034 2: 5.48431 3: 3.46871 4: 9.19356 5: 10.0624 6: 8.22885 7: 16.5997 8: 7.97126 9: 9.72255 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.498096 0.843247 0 0 0 1 : 0 0.439995 0.121731 0 0 2 : 0.00675643 0 0.295072 0.906265 0 3 : 0 0.733284 0 0.664596 0.667991 4 : 0 0 0.668451 0 0.0905497 5 : 0 0 0 0.325539 0 6 : 0 0 0 0 0.0353862 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.935537 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.387329 0 0 0 0 5 : 0.965327 0.95131 0 0 0 6 : 0 0.440378 0.0548459 0 0 7 : 0.732351 0 0.799204 0.221937 0 8 : 0 0.461089 0 0.914885 0.734439 9 : 0 0 0.173443 0 0.289981 A * x: 0: 11.54 1: 1.24518 2: 4.51703 3: 7.46491 4: 4.78208 5: 13.7533 6: 3.69835 7: 12.7852 8: 18.806 9: 4.28735 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.243838 0.638669 0 0 0 1 : 0 0.0237142 0.216276 0 0 2 : 0.0495138 0 0.169679 0.318093 0 3 : 0 0.969297 0 0.926937 0.947378 4 : 0 0 0.47111 0 0.769818 5 : 0 0 0 0.670496 0 6 : 0 0 0 0 0.687188 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.229997 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.581727 0 0 0 0 5 : 0.939484 0.589998 0 0 0 6 : 0 0.999215 0.965955 0 0 7 : 0.173532 0 0.303333 0.296118 0 8 : 0 0.521474 0 0.407443 0.597637 9 : 0 0 0.164159 0 0.915823 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.903787 0.513042 0 0 0 1 : 0 0.970544 0.709769 0 0 2 : 0.95326 0 0.931052 0.721521 0 3 : 0 0.0108551 0 0.844887 0.752138 4 : 0 0 0.978168 0 0.67905 5 : 0 0 0 0.370292 0 6 : 0 0 0 0 0.705515 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.126141 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.332093 0 0 0 0 5 : 0.408493 0.207221 0 0 0 6 : 0 0.0913218 0.9317 0 0 7 : 0.207624 0 0.0884228 0.528466 0 8 : 0 0.856123 0 0.186042 0.361665 9 : 0 0 0.901139 0 0.613808 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. Tue May 20 22:27:34 2025