11 April 2023 10:37:00.863 AM r83t_test(): FORTRAN90 version: Test r83t(). R83T_CG_TEST R83T_CG solves an R83T linear system using the conjugate gradient method. M = 5 N = 5 The R83T matrix A: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 2 -1.00000 2.00000 -1.00000 3 -1.00000 2.00000 -1.00000 4 -1.00000 2.00000 -1.00000 5 -1.00000 2.00000 The right hand side B: 1: 0.0000000 2: 0.0000000 3: 0.0000000 4: 0.0000000 5: 6.0000000 The solution X: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 R83T_DIF2_TEST R83T_DIF2 sets up an R83T second difference matrix. M = 5 N = 5 The R83T second difference matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 2 -1.00000 2.00000 -1.00000 3 -1.00000 2.00000 -1.00000 4 -1.00000 2.00000 -1.00000 5 -1.00000 2.00000 R83T_GS_SL_TEST R83T_GS_SL solves a linear system using Gauss-Seidel iteration, with R83T matrix storage. Matrix order N = 10 Iterations per call = 25 Current solution estimate: 1: 0.63460551 2: 1.3265045 3: 2.0956645 4: 2.9545086 5: 3.9074095 6: 4.9508817 7: 6.0744243 8: 7.2619014 9: 8.4932953 10: 9.7466477 Current solution estimate: 1: 0.95356655 2: 1.9145041 3: 2.8853282 4: 3.8675698 5: 4.8617324 6: 5.8673331 7: 6.8830194 8: 7.9067459 9: 8.9359910 10: 9.9679955 Current solution estimate: 1: 0.99412614 2: 1.9891847 3: 2.9854940 4: 3.9832476 5: 4.9825092 6: 5.9832177 7: 6.9852020 8: 7.9882034 9: 8.9919029 10: 9.9959514 R83T_INDICATOR_TEST R83T_INDICATOR sets up an R83T indicator matrix. M = 5 N = 5 The R83T indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 32.0000 33.0000 34.0000 4 43.0000 44.0000 45.0000 5 54.0000 55.0000 R83T_JAC_SL_TEST R83T_JAC_SL solves a linear system using Jacobi iteration, with R83T matrix storage. Matrix order N = 10 Iterations per call = 25 Current solution estimate: 1: 0.31517148 2: 0.72779730 3: 1.1404231 4: 1.8275799 5: 2.5147367 6: 3.5904667 7: 4.6661967 8: 6.1282006 9: 7.5902045 10: 9.2951023 Current solution estimate: 1: 0.75754530 2: 1.5150906 3: 2.3493609 4: 3.1836311 5: 4.1472871 6: 5.1109431 7: 6.2158120 8: 7.3206809 9: 8.5336643 10: 9.7466477 Current solution estimate: 1: 0.91002076 2: 1.8343234 3: 2.7586261 4: 3.7212423 5: 4.6838586 6: 5.6966602 7: 6.7094617 8: 7.7683880 9: 8.8273142 10: 9.9136571 R83T_MTV_TEST R83T_MTV multiplies an R83T matrix transposed times a vector. M = 5 N = 6 The R83T matrix A: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 32.0000 33.0000 34.0000 4 43.0000 44.0000 45.0000 5 54.0000 55.0000 Col: 6 Row --- 5 56.0000 The vector x: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 The product b = A'*x: 1: 53.000000 2: 152.00000 3: 317.00000 4: 548.00000 5: 455.00000 6: 280.00000 R83T_MV_TEST R83T_MV multiplies an R83T matrix times a vector. M = 5 N = 6 The R83T matrix A: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 32.0000 33.0000 34.0000 4 43.0000 44.0000 45.0000 5 54.0000 55.0000 Col: 6 Row --- 5 56.0000 The vector x: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 6: 6.0000000 The product b = A*x: 1: 35.000000 2: 134.00000 3: 299.00000 4: 530.00000 5: 827.00000 R83T_PRINT_TEST R83T_PRINT prints an R83T matrix. M = 5 N = 5 The R83T matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 32.0000 33.0000 34.0000 4 43.0000 44.0000 45.0000 5 54.0000 55.0000 R83T_PRINT_SOME_TEST R83T_PRINT_SOME prints some of an R83T matrix. M = 9 N = 9 Rows 3:6, Cols 5:8: Col: 5 6 7 8 Row --- 4 45.0000 5 55.0000 56.0000 6 65.0000 66.0000 67.0000 R83T_RANDOM_TEST R83T_RANDOM sets up an R83T random matrix. M = 5 N = 5 The R83T random matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 2 0.829509 0.561695 0.415307 3 0.661187E-01 0.257578 0.109957 4 0.438290E-01 0.633966 0.617272E-01 5 0.449539 0.401306 R83T_RES_TEST R83T_RES evaluates the residual given an approximate solution of a linear system A*x=b. M = 5 N = 5 The R83T matrix A: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 2 -1.00000 2.00000 -1.00000 3 -1.00000 2.00000 -1.00000 4 -1.00000 2.00000 -1.00000 5 -1.00000 2.00000 The right hand side B: 1: 0.0000000 2: 0.0000000 3: 0.0000000 4: 0.0000000 5: 6.0000000 The solution X: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 The residual b-A*x: 1: 0.22204460E-15 2: -0.44408921E-15 3: 0.88817842E-15 4: 0.0000000 5: 0.0000000 R83T_TO_R8GE_TEST R83T_TO_R8GE converts an R83T matrix to R8GE format. M = 5 N = 5 The R83T indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 2 21.0000 22.0000 23.0000 3 32.0000 33.0000 34.0000 4 43.0000 44.0000 45.0000 5 54.0000 55.0000 The R8GE format matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 0.00000 2 21.0000 22.0000 23.0000 0.00000 0.00000 3 0.00000 32.0000 33.0000 34.0000 0.00000 4 0.00000 0.00000 43.0000 44.0000 45.0000 5 0.00000 0.00000 0.00000 54.0000 55.0000 R83T_ZEROS_TEST R83T_ZEROS sets up an R83T zero matrix. M = 5 N = 5 The R83T zero matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 0.00000 2 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 5 0.00000 0.00000 r83t_test(): Normal end of execution. 11 April 2023 10:37:00.863 AM