20 January 2017 11:43:24 AM R8TO_PRB C version Test the R8TO library. R8TO_DIF2_TEST R8TO_DIF2 sets the second difference as an R8TO matrix. Matrix order N = 5 The R8TO matrix: Col: 1 2 3 4 5 Row --- 1 2 -1 0 0 0 2 -1 2 -1 0 0 3 0 -1 2 -1 0 4 0 0 -1 2 -1 5 0 0 0 -1 2 R8TO_INDICATOR_TEST R8TO_INDICATOR sets up an R8TO indicator matrix. Matrix order N = 5 The R8TO indicator matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 14 15 2 21 11 12 13 14 3 31 21 11 12 13 4 41 31 21 11 12 5 51 41 31 21 11 R8TO_MTV_TEST R8TO_MTV computes b=A'*x=b, where A is an R8TO matrix. Matrix order N = 5 The Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 14 15 2 21 11 12 13 14 3 31 21 11 12 13 4 41 31 21 11 12 5 51 41 31 21 11 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 b=A'*x: 0 565.000000 1 426.000000 2 309.000000 3 225.000000 4 185.000000 R8TO_MV_TEST R8TO_MV computes b=A*x=b, where A is an R8TO matrix. Matrix order N = 5 The Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 14 15 2 21 11 12 13 14 3 31 21 11 12 13 4 41 31 21 11 12 5 51 41 31 21 11 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 b=A*x: 0 205.000000 1 201.000000 2 219.000000 3 270.000000 4 365.000000 R8TO_PRINT_TEST R8TO_PRINT prints an R8TO matrix. Matrix order N = 5 The R8TO indicator matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 14 15 2 21 11 12 13 14 3 31 21 11 12 13 4 41 31 21 11 12 5 51 41 31 21 11 R8TO_PRINT_SOME_TEST R8TO_PRINT_SOME prints some of an R8TO matrix. Matrix order N = 5 Rows 1:4, Cols 0:2 Col: 0 1 2 Row --- 1 21 11 12 2 31 21 11 3 41 31 21 4 51 41 31 R8TO_RANDOM_TEST R8TO_RANDOM randomizes an R8TO matrix. Matrix order N = 5 The Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.829509 0.561695 0.415307 2 0.0661187 0.218418 0.956318 0.829509 0.561695 3 0.257578 0.0661187 0.218418 0.956318 0.829509 4 0.109957 0.257578 0.0661187 0.218418 0.956318 5 0.043829 0.109957 0.257578 0.0661187 0.218418 R8TO_SL_TEST R8TO_SL solves a Toeplitz system. Matrix order N = 5 The Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.829509 0.561695 0.415307 2 0.0661187 0.218418 0.956318 0.829509 0.561695 3 0.257578 0.0661187 0.218418 0.956318 0.829509 4 0.109957 0.257578 0.0661187 0.218418 0.956318 5 0.043829 0.109957 0.257578 0.0661187 0.218418 Solution: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 R8TO_SLT_TEST R8TO_SLT solves A'*x=b, where A is an R8TO matrix. Matrix order N = 5 The Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.829509 0.561695 0.415307 2 0.0661187 0.218418 0.956318 0.829509 0.561695 3 0.257578 0.0661187 0.218418 0.956318 0.829509 4 0.109957 0.257578 0.0661187 0.218418 0.956318 5 0.043829 0.109957 0.257578 0.0661187 0.218418 Solution to transposed system: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 R8TO_TO_R8GE_TEST R8TO_TO_R8GE converts an R8TO matrix to R8GE format. Matrix order N = 5 The R8TO matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.829509 0.561695 0.415307 2 0.0661187 0.218418 0.956318 0.829509 0.561695 3 0.257578 0.0661187 0.218418 0.956318 0.829509 4 0.109957 0.257578 0.0661187 0.218418 0.956318 5 0.043829 0.109957 0.257578 0.0661187 0.218418 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.829509 0.561695 0.415307 2 0.0661187 0.218418 0.956318 0.829509 0.561695 3 0.257578 0.0661187 0.218418 0.956318 0.829509 4 0.109957 0.257578 0.0661187 0.218418 0.956318 5 0.043829 0.109957 0.257578 0.0661187 0.218418 R8TO_ZEROS_TEST R8TO_ZEROS zeros an R8TO matrix. Matrix order N = 5 The R8TO matrix: Col: 1 2 3 4 5 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0 0 0 0 0 R8TO_PRB Normal end of execution. 20 January 2017 11:43:24 AM