20 October 2023 11:06:32.127 AM linplus_test(): FORTRAN77 version: Test linplus(). TEST01 C83_CR_FA factors a complex tridiagonal matrix; C83_CR_SLS solves 1 or more factored systems. Matrix order N = 10 Solution: Col: 1 Row --- 1: 1.00 10.0 2: 2.00 20.0 3: 3.00 30.0 4: 4.00 40.0 5: 5.00 50.0 6: 6.00 60.0 7: 7.00 70.0 8: 8.00 80.0 9: 9.00 90.0 10: 10.0 100. TEST0196 For a general matrix, R8VEC_TO_R8GE converts a real vector to an R8GE matrix. R8GE_TO_R8VEC converts an R8GE matrix to a real vector. Matrix rows M = 4 Matrix columns N = 6 The R8GE indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 21.0000 22.0000 23.0000 24.0000 25.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 41.0000 42.0000 43.0000 44.0000 45.0000 Col: 6 Row --- 1 16.0000 2 26.0000 3 36.0000 4 46.0000 1 1 1 11.0000 2 1 2 21.0000 3 1 3 31.0000 4 1 4 41.0000 1 2 5 12.0000 2 2 6 22.0000 3 2 7 32.0000 4 2 8 42.0000 1 3 9 13.0000 2 3 10 23.0000 3 3 11 33.0000 4 3 12 43.0000 1 4 13 14.0000 2 4 14 24.0000 3 4 15 34.0000 4 4 16 44.0000 1 5 17 15.0000 2 5 18 25.0000 3 5 19 35.0000 4 5 20 45.0000 1 6 21 16.0000 2 6 22 26.0000 3 6 23 36.0000 4 6 24 46.0000 The recovered R8GE indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 21.0000 22.0000 23.0000 24.0000 25.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 41.0000 42.0000 43.0000 44.0000 45.0000 Col: 6 Row --- 1 16.0000 2 26.0000 3 36.0000 4 46.0000 R8BTO_SL_TEST R8BTO_SL solves a block Toeplitz system. Block order M = 2 Block number L = 3 Matrix order N = 6 The block Toeplitz matrix: Col: 1 2 3 4 5 Row --- 1 9. 1. 3. 4. 5. 2 2. 8. 6. 6. 7. 3 7. 8. 9. 1. 3. 4 8. 8. 2. 8. 6. 5 9. 0. 7. 8. 9. 6 9. 9. 8. 8. 2. Col: 6 Row --- 1 6. 2 7. 3 4. 4 6. 5 1. 6 8. The right hand side B: 1 97.000000 2 137.00000 3 93.000000 4 128.00000 5 113.00000 6 141.00000 The computed solution X: 1 1.0000000 2 2.0000000 3 3.0000000 4 4.0000000 5 5.0000000 6 6.0000000 R8GE_DET_TEST R8GE_DET computes the determinant of an R8GE matrix. Matrix order N = 4 R8GE_DET computes the determinant = 112.000 Exact determinant = 112.000 TEST295 For a matrix in general storage, R8GE_DILU returns the DILU factors. Matrix rows M = 9 Matrix columns N = 9 Matrix A: Col: 1 2 3 4 5 Row --- 1 4.00000 -1.00000 0.00000 -1.00000 0.00000 2 -1.00000 4.00000 -1.00000 0.00000 -1.00000 3 0.00000 -1.00000 4.00000 -1.00000 0.00000 4 -1.00000 0.00000 -1.00000 4.00000 -1.00000 5 0.00000 -1.00000 0.00000 -1.00000 4.00000 6 0.00000 0.00000 -1.00000 0.00000 -1.00000 7 0.00000 0.00000 0.00000 -1.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 -1.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 Row --- 1 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 3 -1.00000 0.00000 0.00000 0.00000 4 0.00000 -1.00000 0.00000 0.00000 5 -1.00000 0.00000 -1.00000 0.00000 6 4.00000 -1.00000 0.00000 -1.00000 7 -1.00000 4.00000 -1.00000 0.00000 8 0.00000 -1.00000 4.00000 -1.00000 9 -1.00000 0.00000 -1.00000 4.00000 DILU factor: 1 0.25000000 2 0.26666667 3 0.26785714 4 0.28717949 5 0.29017857 6 0.29053178 7 0.29220211 8 0.29260134 9 0.29266578 TEST31 For a matrix in general storage, R8GE_FA computes the LU factors, R8GE_SL solves a factored system. Matrix order N = 5 The matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.661187E-01 0.617272E-01 0.183837E-02 0.859097 2 0.956318 0.257578 0.449539 0.897504 0.840847 3 0.829509 0.109957 0.401306 0.350752 0.123104 4 0.561695 0.438290E-01 0.754673 0.945448E-01 0.751236E-02 5 0.415307 0.633966 0.797287 0.136169E-01 0.260303 The compressed LU factors: Col 1 2 3 4 5 Row 1 0.956318 0.257578 0.449539 0.897504 0.840847 2 -0.228395 0.522106 0.602062 -0.376149 -0.104858 3 -0.867399 0.217324 0.614552 -0.510026 -0.507943 4 -0.587352 0.205820 -0.231419 -0.391459 -0.511487 5 -0.434277 -0.139612E-01 0.803036E-01 -0.610161 0.939815 The pivot vector P: 1 2 2 5 3 4 4 4 5 5 Solution: 1 1.0000000 2 2.0000000 3 3.0000000 4 4.0000000 5 5.0000000 TEST315 For a matrix in general storage, R8GE_ILU returns the ILU factors. Matrix rows M = 9 Matrix columns N = 9 Matrix A: Col: 1 2 3 4 5 Row --- 1 4.00000 -1.00000 0.00000 -1.00000 0.00000 2 -1.00000 4.00000 -1.00000 0.00000 -1.00000 3 0.00000 -1.00000 4.00000 -1.00000 0.00000 4 -1.00000 0.00000 -1.00000 4.00000 -1.00000 5 0.00000 -1.00000 0.00000 -1.00000 4.00000 6 0.00000 0.00000 -1.00000 0.00000 -1.00000 7 0.00000 0.00000 0.00000 -1.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 -1.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 Row --- 1 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 3 -1.00000 0.00000 0.00000 0.00000 4 0.00000 -1.00000 0.00000 0.00000 5 -1.00000 0.00000 -1.00000 0.00000 6 4.00000 -1.00000 0.00000 -1.00000 7 -1.00000 4.00000 -1.00000 0.00000 8 0.00000 -1.00000 4.00000 -1.00000 9 -1.00000 0.00000 -1.00000 4.00000 Factor L: Col: 1 2 3 4 5 Row --- 1 1.00000 0.00000 0.00000 0.00000 0.00000 2 -0.250000 1.00000 0.00000 0.00000 0.00000 3 0.00000 -0.266667 1.00000 0.00000 0.00000 4 -0.250000 0.00000 -0.267857 1.00000 0.00000 5 0.00000 -0.266667 0.00000 -0.287179 1.00000 6 0.00000 0.00000 -0.267857 0.00000 -0.290179 7 0.00000 0.00000 0.00000 -0.287179 0.00000 8 0.00000 0.00000 0.00000 0.00000 -0.290179 9 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 Row --- 1 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 0.00000 6 1.00000 0.00000 0.00000 0.00000 7 -0.290532 1.00000 0.00000 0.00000 8 0.00000 -0.292202 1.00000 0.00000 9 -0.290532 0.00000 -0.292601 1.00000 Factor U: Col: 1 2 3 4 5 Row --- 1 4.00000 -1.00000 0.00000 -1.00000 0.00000 2 0.00000 3.75000 -1.00000 0.00000 -1.00000 3 0.00000 0.00000 3.73333 -1.00000 0.00000 4 0.00000 0.00000 0.00000 3.48214 -1.00000 5 0.00000 0.00000 0.00000 0.00000 3.44615 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 Col: 6 7 8 9 Row --- 1 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 3 -1.00000 0.00000 0.00000 0.00000 4 0.00000 -1.00000 0.00000 0.00000 5 -1.00000 0.00000 -1.00000 0.00000 6 3.44196 -1.00000 0.00000 -1.00000 7 0.00000 3.42229 -1.00000 0.00000 8 0.00000 0.00000 3.41762 -1.00000 9 0.00000 0.00000 0.00000 3.41687 Product L*U: Col: 1 2 3 4 5 Row --- 1 4.00000 -1.00000 0.00000 -1.00000 0.00000 2 -1.00000 4.00000 -1.00000 0.250000 -1.00000 3 0.00000 -1.00000 4.00000 -1.00000 0.266667 4 -1.00000 0.250000 -1.00000 4.00000 -1.00000 5 0.00000 -1.00000 0.266667 -1.00000 4.00000 6 0.00000 0.00000 -1.00000 0.267857 -1.00000 7 0.00000 0.00000 0.00000 -1.00000 0.287179 8 0.00000 0.00000 0.00000 0.00000 -1.00000 9 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 Row --- 1 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 3 -1.00000 0.00000 0.00000 0.00000 4 0.267857 -1.00000 0.00000 0.00000 5 -1.00000 0.287179 -1.00000 0.00000 6 4.00000 -1.00000 0.290179 -1.00000 7 -1.00000 4.00000 -1.00000 0.290532 8 0.290179 -1.00000 4.00000 -1.00000 9 -1.00000 0.290532 -1.00000 4.00000 TEST317 For a matrix in general storage, R8GE_INDICATOR sets up the indicator matrix. Matrix rows M = 7 Matrix columns N = 5 The R8GE indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 21.0000 22.0000 23.0000 24.0000 25.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 41.0000 42.0000 43.0000 44.0000 45.0000 5 51.0000 52.0000 53.0000 54.0000 55.0000 6 61.0000 62.0000 63.0000 64.0000 65.0000 7 71.0000 72.0000 73.0000 74.0000 75.0000 TEST34 For a matrix in general storage, R8GE_FS factors and solves a linear system. Matrix order N = 10 Solution: 1 1.0000000 2 2.0000000 3 3.0000000 4 4.0000000 5 5.0000000 6 6.0000000 7 7.0000000 8 8.0000000 9 9.0000000 10 10.000000 TEST345 For a matrix in general storage, R8GE_FSS factors and solves multiple linear system. Matrix order N = 10 Solutions: Col 1 2 3 Row 1 1.00000 1.00000 1.00000 2 1.00000 2.00000 2.00000 3 1.00000 3.00000 3.00000 4 1.00000 4.00000 1.00000 5 1.00000 5.00000 2.00000 6 1.00000 6.00000 3.00000 7 1.00000 7.00000 1.00000 8 1.00000 8.00000 2.00000 9 1.00000 9.00000 3.00000 10 1.00000 10.0000 1.00000 linplus_test(): Normal end of execution. 20 October 2023 11:06:32.129 AM