19 July 2016 8:03:44.616 AM R8GD_PRB FORTRAN90 version: Test the R8GD library. 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.218418 0.956318 0.00000 0.00000 0.00000 2 0.00000 0.561695 0.415307 0.00000 0.00000 3 0.661187E-01 0.00000 0.257578 0.109957 0.00000 4 0.00000 0.438290E-01 0.00000 0.633966 0.617272E-01 5 0.00000 0.00000 0.449539 0.00000 0.401306 6 0.00000 0.00000 0.00000 0.797287 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.350752 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.829509 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.754673 0.00000 0.00000 0.00000 0.00000 6 0.183837E-02 0.897504 0.00000 0.00000 0.00000 7 0.00000 0.945448E-01 0.136169E-01 0.00000 0.00000 8 0.859097 0.00000 0.840847 0.123104 0.00000 9 0.00000 0.751236E-02 0.00000 0.260303 0.912484 10 0.00000 0.00000 0.113664 0.00000 0.351629 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 0.416775 2 2.25502 3 3.85104 4 7.64945 5 4.70871 6 10.6572 7 6.11445 8 7.95874 9 3.32756 10 12.5581 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.218418 0.956318 0.00000 0.00000 0.00000 2 0.00000 0.561695 0.415307 0.00000 0.00000 3 0.661187E-01 0.00000 0.257578 0.109957 0.00000 4 0.00000 0.438290E-01 0.00000 0.633966 0.617272E-01 5 0.00000 0.00000 0.449539 0.00000 0.401306 6 0.00000 0.00000 0.00000 0.797287 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.350752 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.829509 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.754673 0.00000 0.00000 0.00000 0.00000 6 0.183837E-02 0.897504 0.00000 0.00000 0.00000 7 0.00000 0.945448E-01 0.136169E-01 0.00000 0.00000 8 0.859097 0.00000 0.840847 0.123104 0.00000 9 0.00000 0.751236E-02 0.00000 0.260303 0.912484 10 0.00000 0.00000 0.113664 0.00000 0.351629 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 10.4261 2 2.36931 3 1.27868 4 2.93216 5 7.88319 6 9.48271 7 2.52451 8 12.9893 9 11.5202 10 4.42560 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.218418 0.956318 0.00000 0.00000 0.00000 2 0.00000 0.561695 0.415307 0.00000 0.00000 3 0.661187E-01 0.00000 0.257578 0.109957 0.00000 4 0.00000 0.438290E-01 0.00000 0.633966 0.617272E-01 5 0.00000 0.00000 0.449539 0.00000 0.401306 6 0.00000 0.00000 0.00000 0.797287 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.350752 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.829509 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.754673 0.00000 0.00000 0.00000 0.00000 6 0.183837E-02 0.897504 0.00000 0.00000 0.00000 7 0.00000 0.945448E-01 0.136169E-01 0.00000 0.00000 8 0.859097 0.00000 0.840847 0.123104 0.00000 9 0.00000 0.751236E-02 0.00000 0.260303 0.912484 10 0.00000 0.00000 0.113664 0.00000 0.351629 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.218418 0.956318 0.00000 0.00000 0.00000 2 0.00000 0.561695 0.415307 0.00000 0.00000 3 0.661187E-01 0.00000 0.257578 0.109957 0.00000 4 0.00000 0.438290E-01 0.00000 0.633966 0.617272E-01 5 0.00000 0.00000 0.449539 0.00000 0.401306 6 0.00000 0.00000 0.00000 0.797287 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.350752 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.829509 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.754673 0.00000 0.00000 0.00000 0.00000 6 0.183837E-02 0.897504 0.00000 0.00000 0.00000 7 0.00000 0.945448E-01 0.136169E-01 0.00000 0.00000 8 0.859097 0.00000 0.840847 0.123104 0.00000 9 0.00000 0.751236E-02 0.00000 0.260303 0.912484 10 0.00000 0.00000 0.113664 0.00000 0.351629 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_PRB Normal end of execution. 19 July 2016 8:03:44.619 AM