Wed Oct 8 08:50:30 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.543499 0.708358 0 0 0 1 : 0 0.679207 0.00831241 0 0 2 : 0.931732 0 0.0783891 0.0248894 0 3 : 0 0.817921 0 0.35215 0.786162 4 : 0 0 0.363797 0 0.566969 5 : 0 0 0 0.533684 0 6 : 0 0 0 0 0.631339 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.352611 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.442345 0 0 0 0 5 : 0.119798 0.416513 0 0 0 6 : 0 0.700492 0.329899 0 0 7 : 0.215792 0 0.929882 0.996492 0 8 : 0 0.561812 0 0.284769 0.509431 9 : 0 0 0.825482 0 0.321513 A' * x: 0: 3.3387 1: 5.33846 2: 2.07078 3: 4.68537 4: 10.3989 5: 4.65685 6: 12.4588 7: 18.0032 8: 10.5349 9: 8.15262 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.714723 0.857571 0 0 0 1 : 0 0.001271 0.206396 0 0 2 : 0.502935 0 0.9272 0.118196 0 3 : 0 0.869225 0 0.677096 0.164834 4 : 0 0 0.985713 0 0.0101878 5 : 0 0 0 0.788936 0 6 : 0 0 0 0 0.0122765 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.428599 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.613038 0 0 0 0 5 : 0.081504 0.480853 0 0 0 6 : 0 0.837871 0.733599 0 0 7 : 0.190865 0 0.559291 0.389104 0 8 : 0 0.401944 0 0.0170173 0.782517 9 : 0 0 0.390319 0 0.098574 A * x: 0: 6.71586 1: 0.62173 2: 3.75732 3: 5.271 4: 6.68631 5: 7.01074 6: 11.7953 7: 9.12146 8: 10.7919 9: 4.10829 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.892978 0.540709 0 0 0 1 : 0 0.673254 0.246262 0 0 2 : 0.757503 0 0.812532 0.697651 0 3 : 0 0.965871 0 0.571482 0.166674 4 : 0 0 0.162947 0 0.522754 5 : 0 0 0 0.164197 0 6 : 0 0 0 0 0.842435 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.913583 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.762102 0 0 0 0 5 : 0.516993 0.427082 0 0 0 6 : 0 0.96687 0.0179927 0 0 7 : 0.844449 0 0.515661 0.698794 0 8 : 0 0.241198 0 0.679276 0.170225 9 : 0 0 0.357308 0 0.517724 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.978786 0.545347 0 0 0 1 : 0 0.587165 0.5964 0 0 2 : 0.677439 0 0.406662 0.736802 0 3 : 0 0.17401 0 0.317332 0.818833 4 : 0 0 0.215166 0 0.734666 5 : 0 0 0 0.540454 0 6 : 0 0 0 0 0.760234 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.23034 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0.200715 0 0 0 0 5 : 0.678998 0.761366 0 0 0 6 : 0 0.197786 0.890302 0 0 7 : 0.142911 0 0.296962 0.57955 0 8 : 0 0.0431161 0 0.997246 0.384991 9 : 0 0 0.666037 0 0.0989704 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. Wed Oct 8 08:50:31 2025