20 January 2017 10:17:35 AM R85_PRB C version Test the R85 library. R85_DIF2_TEST R85_DIF2 sets up an R85 second difference matrix. Matrix order N = 5 The R85 second difference matrix: Col: 1 2 3 4 5 Row --- 1 2.000000 -1.000000 0.000000 2 -1.000000 2.000000 -1.000000 0.000000 3 0.000000 -1.000000 2.000000 -1.000000 0.000000 4 0.000000 -1.000000 2.000000 -1.000000 5 0.000000 -1.000000 2.000000 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.000000 12.000000 13.000000 2 21.000000 22.000000 23.000000 24.000000 3 31.000000 32.000000 33.000000 34.000000 35.000000 4 42.000000 43.000000 44.000000 45.000000 5 53.000000 54.000000 55.000000 R85_MTV_TEST R85_MTV computes b=A'*x, where A is an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 1 2 3 4 5 Row --- 1 11.000000 12.000000 13.000000 2 21.000000 22.000000 23.000000 24.000000 3 31.000000 32.000000 33.000000 34.000000 35.000000 4 42.000000 43.000000 44.000000 45.000000 5 53.000000 54.000000 55.000000 The vector x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 The product b=A'*x: 0 11.000000 1 112.000000 2 395.000000 3 508.000000 4 275.000000 R85_MV_TEST R85_MV computes b=A*x, where A is an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 1 2 3 4 5 Row --- 1 11.000000 12.000000 13.000000 2 21.000000 22.000000 23.000000 24.000000 3 31.000000 32.000000 33.000000 34.000000 35.000000 4 42.000000 43.000000 44.000000 45.000000 5 53.000000 54.000000 55.000000 The vector x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 The product b=A*x: 0 146.000000 1 320.000000 2 595.000000 3 596.000000 4 560.000000 R85_NP_FS_TEST R85_NP_FS factors and solves a pentadiagonal linear system, with no pivoting. Matrix order N = 10 The pentadiagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.350752 0.043829 0.218418 2 0.351629 0.094545 0.633966 0.956318 3 0.188955 0.822887 0.013617 0.061727 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.066119 7 0.861216 0.007512 0.797287 0.257578 8 0.367027 0.453794 0.260303 0.001838 0.109957 9 0.617205 0.911977 0.912484 0.897504 10 0.361529 0.597917 0.113664 Right hand side: 0 1.620874 1 5.747549 2 4.580672 3 10.106421 4 15.881457 5 13.956661 6 15.842357 7 19.883225 8 16.009274 9 10.093831 Solution to A*x=b: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 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.000000 12.000000 13.000000 2 21.000000 22.000000 23.000000 24.000000 3 31.000000 32.000000 33.000000 34.000000 35.000000 4 42.000000 43.000000 44.000000 45.000000 5 53.000000 54.000000 55.000000 R85_PRINT_SOME_TEST R85_PRINT_SOME prints some of an R85 matrix. Matrix order N = 8 Rows 1-5, Cols 2-4: Col: 2 3 4 Row --- 1 12.000000 13.000000 2 22.000000 23.000000 24.000000 3 32.000000 33.000000 34.000000 4 42.000000 43.000000 44.000000 5 53.000000 54.000000 R85_RANDOM_TEST R85_RANDOM randomizes an R85 matrix. Matrix order N = 5 The R85 random matrix: Col: 1 2 3 4 5 Row --- 1 0.109957 0.561695 0.218418 2 0.401306 0.043829 0.415307 0.956318 3 0.897504 0.754673 0.633966 0.066119 0.829509 4 0.350752 0.797287 0.061727 0.257578 5 0.094545 0.001838 0.449539 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.000000 12.000000 13.000000 2 21.000000 22.000000 23.000000 24.000000 3 31.000000 32.000000 33.000000 34.000000 35.000000 4 42.000000 43.000000 44.000000 45.000000 5 53.000000 54.000000 55.000000 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11 21 31 0 0 2 12 22 32 42 0 3 13 23 33 43 53 4 0 24 34 44 54 5 0 0 35 45 55 R85_ZEROS_TEST R85_ZEROS zeros an R85 matrix. Matrix order N = 5 The R85 zero matrix: Col: 1 2 3 4 5 Row --- 1 0.000000 0.000000 0.000000 2 0.000000 0.000000 0.000000 0.000000 3 0.000000 0.000000 0.000000 0.000000 0.000000 4 0.000000 0.000000 0.000000 0.000000 5 0.000000 0.000000 0.000000 R85_PRB Normal end of execution. 20 January 2017 10:17:35 AM