11 April 2023 10:29:37.593 AM r85_test(): FORTRAN90 version: Test r85(). R85_DIF2_TEST R85_DIF2 sets up an R85 second difference matrix. Matrix order N = 8 The R85 second difference matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 -1.00000 2.00000 6 0.00000 -1.00000 7 0.00000 Col: 6 7 8 Row --- 4 0.00000 5 -1.00000 0.00000 6 2.00000 -1.00000 0.00000 7 -1.00000 2.00000 -1.00000 8 0.00000 -1.00000 2.00000 R85_INDICATOR_TEST R85_INDICATOR sets up an R85 indicator matrix. Matrix order N = 5 The R85 indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 R85_MTV_TEST R85_MTV computes b=A'*x, where A is an R85 matrix. Matrix order N = 8 The R85 indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 6 64.0000 65.0000 7 75.0000 Col: 6 7 8 Row --- 4 46.0000 5 56.0000 57.0000 6 66.0000 67.0000 68.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 The vector 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 The product b=A'*x 1 74.0000 2 230.000 3 505.000 4 890.000 5 1385.00 6 1990.00 7 1994.00 8 1829.00 R85_MV_TEST R85_MV computes b=A*x, where A is an R85 matrix. Matrix order N = 8 The R85 indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 6 64.0000 65.0000 7 75.0000 Col: 6 7 8 Row --- 4 46.0000 5 56.0000 57.0000 6 66.0000 67.0000 68.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 The vector 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 The product b=A*x 1 146.000 2 320.000 3 595.000 4 980.000 5 1475.00 6 2080.00 7 1922.00 8 1658.00 R85_NP_FS_TEST R85_NP_FS factors and solves an R85 linear system, with no pivoting. Matrix order N = 10 The pentadiagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.350752 0.438290E-01 0.218418 2 0.351629 0.945448E-01 0.633966 0.956318 3 0.188955 0.822887 0.136169E-01 0.617272E-01 0.829509 4 0.761492 0.267132 0.859097 0.449539 5 0.396988 0.692066 0.840847 6 0.185314 0.561662 7 0.574366 Col: 6 7 8 9 10 Row --- 4 0.561695 5 0.401306 0.415307 6 0.123104 0.754673 0.661187E-01 7 0.861216 0.751236E-02 0.797287 0.257578 8 0.367027 0.453794 0.260303 0.183837E-02 0.109957 9 0.617205 0.911977 0.912484 0.897504 10 0.361529 0.597917 0.113664 The right hand side: 1 1.62087 2 5.74755 3 4.58067 4 10.1064 5 15.8815 6 13.9567 7 15.8424 8 19.8832 9 16.0093 10 10.0938 Solution: 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 R85_PRINT_TEST R85_PRINT prints an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 R85_PRINT_SOME_TEST R85_PRINT_SOME prints some of an R85 matrix. Matrix order N = 8 Rows 2-6, Cols 3-5: Col: 3 4 5 Row --- 2 23.0000 24.0000 3 33.0000 34.0000 35.0000 4 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 6 64.0000 65.0000 R85_RANDOM_TEST R85_RANDOM sets up a random R85 matrix. Matrix order N = 8 The random R85 matrix: Col: 1 2 3 4 5 Row --- 1 0.754673 0.257578 0.218418 2 0.840847 0.797287 0.109957 0.956318 3 0.822887 0.123104 0.183837E-02 0.438290E-01 0.829509 4 0.267132 0.751236E-02 0.897504 0.633966 5 0.692066 0.260303 0.350752 6 0.561662 0.912484 7 0.861216 Col: 6 7 8 Row --- 4 0.561695 5 0.617272E-01 0.415307 6 0.945448E-01 0.449539 0.661187E-01 7 0.113664 0.136169E-01 0.401306 8 0.453794 0.351629 0.859097 R85_TO_R8GE_TEST R85_TO_R8GE converts an R85 matrix to R8GE format. Matrix order N = 5 The R85 matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 42.0000 43.0000 44.0000 45.0000 5 53.0000 54.0000 55.0000 The R8GE format Col: 1 2 3 4 5 Row --- 1 11.0000 21.0000 31.0000 0.00000 0.00000 2 12.0000 22.0000 32.0000 42.0000 0.00000 3 13.0000 23.0000 33.0000 43.0000 53.0000 4 0.00000 24.0000 34.0000 44.0000 54.0000 5 0.00000 0.00000 35.0000 45.0000 55.0000 R85_ZEROS_TEST R85_ZEROS zeros an R85 matrix. Matrix order N = 8 The zero R85 matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 0.00000 0.00000 2 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 5 0.00000 0.00000 0.00000 6 0.00000 0.00000 7 0.00000 Col: 6 7 8 Row --- 4 0.00000 5 0.00000 0.00000 6 0.00000 0.00000 0.00000 7 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 r85_test(): Normal end of execution. 11 April 2023 10:29:37.594 AM