06 October 2025 6:33:32.280 PM r8gd_test(): Fortran90 version: Test r8gd(). R8GD_DIF2_TEST R8GD_DIF2 sets up an R8GD second difference matrix. Matrix order N = 5 Number of diagonals NDIAG = 3 The R8GS second difference matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 0.00000 0.00000 0.00000 0.00000 2 0.00000 2.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 2.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 2.00000 0.00000 5 0.00000 0.00000 0.00000 0.00000 2.00000 R8GD_INDICATOR_TEST R8GD_INDICATOR sets up an indicator matrix for a general diagonal matrix. Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 9 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 101.000 102.000 0.00000 0.00000 0.00000 2 0.00000 202.000 203.000 0.00000 0.00000 3 301.000 0.00000 303.000 304.000 0.00000 4 0.00000 402.000 0.00000 404.000 405.000 5 0.00000 0.00000 503.000 0.00000 505.000 6 0.00000 0.00000 0.00000 604.000 0.00000 7 0.00000 0.00000 0.00000 0.00000 705.000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.00000 110.000 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 506.000 0.00000 0.00000 0.00000 0.00000 6 606.000 607.000 0.00000 0.00000 0.00000 7 0.00000 707.000 708.000 0.00000 0.00000 8 806.000 0.00000 808.000 809.000 0.00000 9 0.00000 907.000 0.00000 909.000 910.000 10 0.00000 0.00000 1008.00 0.00000 1010.00 R8GD_MTV_TEST R8GD_MTV computes A' * x for a general diagonal matrix. Matrix order N = 10 Number of diagonals NDIAG = 4 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.804585 0.990787 0.00000 0.00000 0.00000 2 0.00000 0.351237 0.464576 0.00000 0.00000 3 0.423226 0.00000 0.878318 0.625311 0.00000 4 0.00000 0.784440E-01 0.00000 0.226459 0.725259 5 0.00000 0.00000 0.407916 0.00000 0.156185 6 0.00000 0.00000 0.00000 0.779539E-01 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.877506 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.00000 0.862259 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.148022 0.00000 0.00000 0.00000 0.00000 6 0.345861 0.457927 0.00000 0.00000 0.00000 7 0.00000 0.534903 0.149943 0.00000 0.00000 8 0.701202 0.00000 0.555350 0.652882 0.00000 9 0.00000 0.847291 0.00000 0.532063 0.512769 10 0.00000 0.00000 0.798979E-02 0.00000 0.645536E-01 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 b=A'*x: 1 2.07426 2 2.00704 3 5.60368 4 3.24949 5 9.82451 6 8.42490 7 14.1175 8 5.57230 9 10.0116 10 6.12272 R8GD_MV_TEST R8GD_MV computes A * x for a general diagonal matrix. Matrix order N = 10 Number of diagonals NDIAG = 4 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.623107 0.668428 0.00000 0.00000 0.00000 2 0.00000 0.738786 0.134894 0.00000 0.00000 3 0.381966 0.00000 0.594773 0.582221 0.00000 4 0.00000 0.396907 0.00000 0.396695 0.970086 5 0.00000 0.00000 0.165531 0.00000 0.752276 6 0.00000 0.00000 0.00000 0.627484 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.403373 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.00000 0.269901 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.539986 0.00000 0.00000 0.00000 0.00000 6 0.601536 0.647212E-01 0.00000 0.00000 0.00000 7 0.00000 0.382269 0.627317 0.00000 0.00000 8 0.739051 0.00000 0.358554 0.889375 0.00000 9 0.00000 0.590386 0.00000 0.300398E-01 0.784051 10 0.00000 0.00000 0.318234 0.00000 0.534779 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 b=A*x: 1 4.65897 2 1.88225 3 4.49517 4 7.23102 5 7.49789 6 6.57220 7 9.71129 8 15.3071 9 12.2436 10 7.89366 R8GD_PRINT_TEST R8GD_PRINT prints a general diagonal matrix Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 9 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.271966 0.648129E-01 0.00000 0.00000 0.00000 2 0.00000 0.669628 0.321098 0.00000 0.00000 3 0.410130 0.00000 0.447899 0.487697 0.00000 4 0.00000 0.233630 0.00000 0.751123 0.254449 5 0.00000 0.00000 0.713925 0.00000 0.629098 6 0.00000 0.00000 0.00000 0.598130 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.841181E-01 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.00000 0.860458 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.682251 0.00000 0.00000 0.00000 0.00000 6 0.887165 0.279552 0.00000 0.00000 0.00000 7 0.00000 0.671076 0.138966 0.00000 0.00000 8 0.433189 0.00000 0.694712 0.791448 0.00000 9 0.00000 0.360211 0.00000 0.603855E-01 0.176990 10 0.00000 0.00000 0.203805 0.00000 0.668492 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: 1 -2 2 0 3 1 4 9 Rows 3-6, Cols 3-6: Col: 3 4 5 6 Row --- 3 303.000 304.000 0.00000 0.00000 4 0.00000 404.000 405.000 0.00000 5 503.000 0.00000 505.000 506.000 6 0.00000 604.000 0.00000 606.000 R8GD_RANDOM_TEST R8GD_RANDOM randomly generates a general diagonal matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 9 The general diagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.771333 0.322932 0.00000 0.00000 0.00000 2 0.00000 0.145134 0.201864 0.00000 0.00000 3 0.322825 0.00000 0.265406 0.504257 0.00000 4 0.00000 0.936679E-01 0.00000 0.996473 0.256140 5 0.00000 0.00000 0.640704 0.00000 0.589237 6 0.00000 0.00000 0.00000 0.607395 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.527701 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.00000 0.844720 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.880457 0.00000 0.00000 0.00000 0.00000 6 0.322586 0.111062 0.00000 0.00000 0.00000 7 0.00000 0.440504 0.460737 0.00000 0.00000 8 0.704766 0.00000 0.614819 0.139358E-01 0.00000 9 0.00000 0.974830 0.00000 0.756183 0.472475 10 0.00000 0.00000 0.296105 0.00000 0.883195E-01 R8GD_TO_R8GE_TEST R8GD_TO_R8GE converts an R8GD matrix to R8GD format. Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 9 The R8GD matrix: Col: 1 2 3 4 5 Row --- 1 101.000 102.000 0.00000 0.00000 0.00000 2 0.00000 202.000 203.000 0.00000 0.00000 3 301.000 0.00000 303.000 304.000 0.00000 4 0.00000 402.000 0.00000 404.000 405.000 5 0.00000 0.00000 503.000 0.00000 505.000 6 0.00000 0.00000 0.00000 604.000 0.00000 7 0.00000 0.00000 0.00000 0.00000 705.000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.00000 110.000 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 506.000 0.00000 0.00000 0.00000 0.00000 6 606.000 607.000 0.00000 0.00000 0.00000 7 0.00000 707.000 708.000 0.00000 0.00000 8 806.000 0.00000 808.000 809.000 0.00000 9 0.00000 907.000 0.00000 909.000 910.000 10 0.00000 0.00000 1008.00 0.00000 1010.00 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 101.000 102.000 0.00000 0.00000 0.00000 2 0.00000 202.000 203.000 0.00000 0.00000 3 301.000 0.00000 303.000 304.000 0.00000 4 0.00000 402.000 0.00000 404.000 405.000 5 0.00000 0.00000 503.000 0.00000 505.000 6 0.00000 0.00000 0.00000 604.000 0.00000 7 0.00000 0.00000 0.00000 0.00000 705.000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.00000 110.000 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 506.000 0.00000 0.00000 0.00000 0.00000 6 606.000 607.000 0.00000 0.00000 0.00000 7 0.00000 707.000 708.000 0.00000 0.00000 8 806.000 0.00000 808.000 809.000 0.00000 9 0.00000 907.000 0.00000 909.000 910.000 10 0.00000 0.00000 1008.00 0.00000 1010.00 R8GD_ZEROS_TEST R8GD_ZEROS returns a zero R8GD matrix; Matrix order N = 10 Number of diagonals NDIAG = 4 The offset vector: 1 -2 2 0 3 1 4 9 The zero R8GD 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 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 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 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 r8gd_test(): Normal end of execution. 06 October 2025 6:33:32.281 PM