19 December 2023 2:40:43.437 PM SPARSEKIT_PRB13 Fortran77 version Test the preconditioners and iterative solvers in SPARSKIT. Methods 1-5 1. Try ILU(0) Preconditioner Precon set-up returned with ierr 0 its = 0 res. norm = 0.498082E+38 its = 1 res. norm = 0.497784E+38 its = 2 res. norm = 0.497556E+38 its = 3 res. norm = 0.497556E+38 its = 4 res. norm = 0.497542E+38 its = 5 res. norm = 0.400401E+38 its = 6 res. norm = 0.645173E+37 its = 7 res. norm = 0.565808E+37 its = 8 res. norm = 0.539531E+37 its = 9 res. norm = 0.217382E+37 its = 10 res. norm = 0.198316E+37 its = 11 res. norm = 0.199024E+37 its = 12 res. norm = 0.199023E+37 its = 13 res. norm = 0.199000E+37 its = 14 res. norm = 0.199000E+37 its = 15 res. norm = 0.168135E+37 its = 16 res. norm = 0.925905E+36 its = 17 res. norm = 0.693906E+36 its = 18 res. norm = 0.398088E+36 its = 19 res. norm = 0.375771E+36 its = 20 res. norm = 0.369362E+36 its = 21 res. norm = 0.389562E+36 its = 22 res. norm = 0.389562E+36 its = 23 res. norm = 0.319131E+36 its = 24 res. norm = 0.319131E+36 its = 25 res. norm = 0.253024E+36 its = 26 res. norm = 0.180636E+36 its = 27 res. norm = 0.160191E+36 its = 28 res. norm = 0.160126E+36 its = 29 res. norm = 0.140310E+36 its = 30 res. norm = 0.129970E+36 its = 31 res. norm = 0.146438E+36 its = 32 res. norm = 0.146438E+36 its = 33 res. norm = 0.115091E+36 its = 34 res. norm = 0.115091E+36 its = 35 res. norm = 0.105517E+36 its = 36 res. norm = 0.994489E+35 its = 37 res. norm = 0.667177E+35 its = 38 res. norm = 0.543198E+35 its = 39 res. norm = 0.480009E+35 its = 40 res. norm = 0.393956E+35 its = 41 res. norm = 0.420539E+35 its = 42 res. norm = 0.420539E+35 its = 43 res. norm = 0.420539E+35 its = 44 res. norm = 0.418945E+35 its = 45 res. norm = 0.296956E+35 its = 46 res. norm = 0.159840E+35 its = 47 res. norm = 0.158832E+35 its = 48 res. norm = 0.140984E+35 its = 49 res. norm = 0.832903E+34 its = 50 res. norm = 0.787796E+34 its = 51 res. norm = 0.814545E+34 its = 52 res. norm = 0.814545E+34 its = 53 res. norm = 0.814545E+34 its = 54 res. norm = 0.814233E+34 its = 55 res. norm = 0.572987E+34 its = 56 res. norm = 0.444219E+34 its = 57 res. norm = 0.442372E+34 its = 58 res. norm = 0.210749E+34 its = 59 res. norm = 0.191258E+34 its = 60 res. norm = 0.144830E+34 its = 61 res. norm = 0.149379E+34 its = 62 res. norm = 0.149379E+34 its = 63 res. norm = 0.149379E+34 its = 64 res. norm = 0.149072E+34 its = 65 res. norm = 0.105057E+34 its = 66 res. norm = 0.984781E+33 its = 67 res. norm = 0.960949E+33 its = 68 res. norm = 0.430126E+33 its = 69 res. norm = 0.429758E+33 its = 70 res. norm = 0.357290E+33 its = 71 res. norm = 0.377145E+33 its = 72 res. norm = 0.377144E+33 its = 73 res. norm = 0.264634E+33 its = 74 res. norm = 0.263698E+33 its = 75 res. norm = 0.258595E+33 its = 76 res. norm = 0.258595E+33 its = 77 res. norm = 0.256972E+33 its = 78 res. norm = 0.205181E+33 its = 79 res. norm = 0.205181E+33 its = 80 res. norm = 0.187204E+33 its = 81 res. norm = 0.228690E+33 its = 82 res. norm = 0.228690E+33 its = 83 res. norm = 0.228690E+33 its = 84 res. norm = 0.228684E+33 its = 85 res. norm = 0.161944E+33 its = 86 res. norm = 0.161692E+33 its = 87 res. norm = 0.148004E+33 its = 88 res. norm = 0.141089E+33 its = 89 res. norm = 0.122467E+33 its = 90 res. norm = 0.121905E+33 its = 91 res. norm = 0.152640E+33 its = 92 res. norm = 0.152640E+33 its = 93 res. norm = 0.152640E+33 its = 94 res. norm = 0.152639E+33 its = 95 res. norm = 0.104157E+33 its = 96 res. norm = 0.103523E+33 its = 97 res. norm = 0.103475E+33 its = 98 res. norm = 0.735820E+32 its = 99 res. norm = 0.577321E+32 its = 100 res. norm = 0.488792E+32 its = 101 res. norm = 0.538952E+32 its = 102 res. norm = 0.538951E+32 its = 103 res. norm = 0.538951E+32 its = 104 res. norm = 0.538951E+32 its = 105 res. norm = 0.365719E+32 its = 106 res. norm = 0.308828E+32 its = 107 res. norm = 0.304638E+32 its = 108 res. norm = 0.291957E+32 its = 109 res. norm = 0.216723E+32 its = 110 res. norm = 0.188111E+32 pgmres returned with ierr = 1 Maximum error in solution = 8.3507732176055747E+018 2. Try MILU(0) Preconditioner Precon set-up returned with ierr 0 its = 0 res. norm = 0.498082E+38 its = 1 res. norm = 0.469474E+38 its = 2 res. norm = 0.469250E+38 its = 3 res. norm = 0.358837E+38 its = 4 res. norm = 0.341645E+38 its = 5 res. norm = 0.334659E+38 its = 6 res. norm = 0.311493E+38 its = 7 res. norm = 0.311434E+38 its = 8 res. norm = 0.681917E+37 its = 9 res. norm = 0.671801E+37 its = 10 res. norm = 0.109219E+37 its = 11 res. norm = 0.109219E+37 its = 12 res. norm = 0.109219E+37 its = 13 res. norm = 0.109219E+37 its = 14 res. norm = 0.109219E+37 its = 15 res. norm = 0.109218E+37 its = 16 res. norm = 0.109218E+37 its = 17 res. norm = 0.109156E+37 its = 18 res. norm = 0.108687E+37 its = 19 res. norm = 0.984310E+36 its = 20 res. norm = 0.324861E+36 its = 21 res. norm = 0.324861E+36 its = 22 res. norm = 0.324861E+36 its = 23 res. norm = 0.324861E+36 its = 24 res. norm = 0.324861E+36 its = 25 res. norm = 0.324861E+36 its = 26 res. norm = 0.324861E+36 its = 27 res. norm = 0.324851E+36 its = 28 res. norm = 0.324766E+36 its = 29 res. norm = 0.324384E+36 its = 30 res. norm = 0.319289E+36 its = 31 res. norm = 0.319289E+36 its = 32 res. norm = 0.319289E+36 its = 33 res. norm = 0.319289E+36 its = 34 res. norm = 0.319289E+36 its = 35 res. norm = 0.319289E+36 its = 36 res. norm = 0.319289E+36 its = 37 res. norm = 0.319289E+36 its = 38 res. norm = 0.319289E+36 its = 39 res. norm = 0.319268E+36 its = 40 res. norm = 0.319268E+36 its = 41 res. norm = 0.319268E+36 its = 42 res. norm = 0.319268E+36 its = 43 res. norm = 0.319268E+36 its = 44 res. norm = 0.319268E+36 its = 45 res. norm = 0.319268E+36 its = 46 res. norm = 0.319268E+36 its = 47 res. norm = 0.319268E+36 its = 48 res. norm = 0.319268E+36 its = 49 res. norm = 0.319268E+36 its = 50 res. norm = 0.319268E+36 its = 51 res. norm = 0.319268E+36 its = 52 res. norm = 0.319268E+36 its = 53 res. norm = 0.319268E+36 its = 54 res. norm = 0.319268E+36 its = 55 res. norm = 0.319268E+36 its = 56 res. norm = 0.319268E+36 its = 57 res. norm = 0.319268E+36 its = 58 res. norm = 0.319268E+36 its = 59 res. norm = 0.319268E+36 its = 60 res. norm = 0.319268E+36 its = 61 res. norm = 0.319268E+36 its = 62 res. norm = 0.319268E+36 its = 63 res. norm = 0.319268E+36 its = 64 res. norm = 0.319268E+36 its = 65 res. norm = 0.319268E+36 its = 66 res. norm = 0.319268E+36 its = 67 res. norm = 0.319268E+36 its = 68 res. norm = 0.319268E+36 its = 69 res. norm = 0.319267E+36 its = 70 res. norm = 0.319267E+36 its = 71 res. norm = 0.319267E+36 its = 72 res. norm = 0.319267E+36 its = 73 res. norm = 0.319267E+36 its = 74 res. norm = 0.319267E+36 its = 75 res. norm = 0.319267E+36 its = 76 res. norm = 0.319267E+36 its = 77 res. norm = 0.319267E+36 its = 78 res. norm = 0.319267E+36 its = 79 res. norm = 0.319267E+36 its = 80 res. norm = 0.319267E+36 its = 81 res. norm = 0.319267E+36 its = 82 res. norm = 0.319267E+36 its = 83 res. norm = 0.319267E+36 its = 84 res. norm = 0.319267E+36 its = 85 res. norm = 0.319267E+36 its = 86 res. norm = 0.319267E+36 its = 87 res. norm = 0.319267E+36 its = 88 res. norm = 0.319267E+36 its = 89 res. norm = 0.319267E+36 its = 90 res. norm = 0.319267E+36 its = 91 res. norm = 0.319267E+36 its = 92 res. norm = 0.319267E+36 its = 93 res. norm = 0.319267E+36 its = 94 res. norm = 0.319267E+36 its = 95 res. norm = 0.319267E+36 its = 96 res. norm = 0.319267E+36 its = 97 res. norm = 0.319267E+36 its = 98 res. norm = 0.319267E+36 its = 99 res. norm = 0.319266E+36 its = 100 res. norm = 0.319266E+36 its = 101 res. norm = 0.319266E+36 its = 102 res. norm = 0.319266E+36 its = 103 res. norm = 0.319266E+36 its = 104 res. norm = 0.319266E+36 its = 105 res. norm = 0.319266E+36 its = 106 res. norm = 0.319266E+36 its = 107 res. norm = 0.319266E+36 its = 108 res. norm = 0.319266E+36 its = 109 res. norm = 0.319266E+36 its = 110 res. norm = 0.319266E+36 pgmres returned with ierr = 1 Maximum error in solution = 8810424091560.4453 3. Try ILUT Preconditioner with tol = 0.001, lfil=1. Precon set-up returned with ierr 5 its = 0 res. norm = 0.498082E+38 its = 1 res. norm = 0.493748E+38 its = 2 res. norm = 0.489952E+38 its = 3 res. norm = 0.489951E+38 its = 4 res. norm = 0.242584E+38 its = 5 res. norm = 0.731282E+37 its = 6 res. norm = 0.723348E+37 its = 7 res. norm = 0.671894E+37 its = 8 res. norm = 0.669702E+37 its = 9 res. norm = 0.229901E+37 its = 10 res. norm = 0.139316E+37 its = 11 res. norm = 0.139399E+37 its = 12 res. norm = 0.139399E+37 its = 13 res. norm = 0.139399E+37 its = 14 res. norm = 0.619590E+36 its = 15 res. norm = 0.517880E+36 its = 16 res. norm = 0.465893E+36 its = 17 res. norm = 0.158987E+36 its = 18 res. norm = 0.125791E+36 its = 19 res. norm = 0.479010E+35 its = 20 res. norm = 0.343122E+35 its = 21 res. norm = 0.343341E+35 its = 22 res. norm = 0.343341E+35 its = 23 res. norm = 0.272868E+35 its = 24 res. norm = 0.267890E+35 its = 25 res. norm = 0.144296E+35 its = 26 res. norm = 0.108508E+35 its = 27 res. norm = 0.789167E+34 its = 28 res. norm = 0.566902E+34 its = 29 res. norm = 0.533668E+34 its = 30 res. norm = 0.489220E+34 its = 31 res. norm = 0.492658E+34 its = 32 res. norm = 0.492658E+34 its = 33 res. norm = 0.347052E+34 its = 34 res. norm = 0.346406E+34 its = 35 res. norm = 0.284200E+34 its = 36 res. norm = 0.280501E+34 its = 37 res. norm = 0.261473E+34 its = 38 res. norm = 0.241069E+34 its = 39 res. norm = 0.211971E+34 its = 40 res. norm = 0.698312E+33 its = 41 res. norm = 0.710605E+33 its = 42 res. norm = 0.710605E+33 its = 43 res. norm = 0.710605E+33 its = 44 res. norm = 0.386085E+33 its = 45 res. norm = 0.381319E+33 its = 46 res. norm = 0.312342E+33 its = 47 res. norm = 0.310614E+33 its = 48 res. norm = 0.285111E+33 its = 49 res. norm = 0.178166E+33 its = 50 res. norm = 0.177849E+33 its = 51 res. norm = 0.187294E+33 its = 52 res. norm = 0.187294E+33 its = 53 res. norm = 0.124642E+33 its = 54 res. norm = 0.980517E+32 its = 55 res. norm = 0.936919E+32 its = 56 res. norm = 0.932034E+32 its = 57 res. norm = 0.924707E+32 its = 58 res. norm = 0.758981E+32 its = 59 res. norm = 0.593626E+32 its = 60 res. norm = 0.521601E+32 its = 61 res. norm = 0.554573E+32 its = 62 res. norm = 0.554573E+32 its = 63 res. norm = 0.505556E+32 its = 64 res. norm = 0.364200E+32 its = 65 res. norm = 0.317523E+32 its = 66 res. norm = 0.235923E+32 its = 67 res. norm = 0.188534E+32 its = 68 res. norm = 0.188219E+32 its = 69 res. norm = 0.123948E+32 its = 70 res. norm = 0.123663E+32 its = 71 res. norm = 0.132854E+32 its = 72 res. norm = 0.132596E+32 its = 73 res. norm = 0.985483E+31 its = 74 res. norm = 0.938901E+31 its = 75 res. norm = 0.849346E+31 its = 76 res. norm = 0.566786E+31 its = 77 res. norm = 0.565627E+31 its = 78 res. norm = 0.379073E+31 pgmres returned with ierr = 0 Maximum error in solution = 47183636107.799881 4. Try ILUT Preconditioner with tol = 0.001, lfil=5. Precon set-up returned with ierr 5 its = 0 res. norm = 0.498082E+38 its = 1 res. norm = 0.493748E+38 its = 2 res. norm = 0.489952E+38 its = 3 res. norm = 0.489951E+38 its = 4 res. norm = 0.242584E+38 its = 5 res. norm = 0.731282E+37 its = 6 res. norm = 0.723348E+37 its = 7 res. norm = 0.671894E+37 its = 8 res. norm = 0.669702E+37 its = 9 res. norm = 0.229901E+37 its = 10 res. norm = 0.139316E+37 its = 11 res. norm = 0.139399E+37 its = 12 res. norm = 0.139399E+37 its = 13 res. norm = 0.139399E+37 its = 14 res. norm = 0.619590E+36 its = 15 res. norm = 0.517880E+36 its = 16 res. norm = 0.465893E+36 its = 17 res. norm = 0.158987E+36 its = 18 res. norm = 0.125791E+36 its = 19 res. norm = 0.479010E+35 its = 20 res. norm = 0.343122E+35 its = 21 res. norm = 0.343341E+35 its = 22 res. norm = 0.343341E+35 its = 23 res. norm = 0.272868E+35 its = 24 res. norm = 0.267890E+35 its = 25 res. norm = 0.144296E+35 its = 26 res. norm = 0.108508E+35 its = 27 res. norm = 0.789167E+34 its = 28 res. norm = 0.566902E+34 its = 29 res. norm = 0.533668E+34 its = 30 res. norm = 0.489220E+34 its = 31 res. norm = 0.492658E+34 its = 32 res. norm = 0.492658E+34 its = 33 res. norm = 0.347052E+34 its = 34 res. norm = 0.346406E+34 its = 35 res. norm = 0.284200E+34 its = 36 res. norm = 0.280501E+34 its = 37 res. norm = 0.261473E+34 its = 38 res. norm = 0.241069E+34 its = 39 res. norm = 0.211971E+34 its = 40 res. norm = 0.698312E+33 its = 41 res. norm = 0.710605E+33 its = 42 res. norm = 0.710605E+33 its = 43 res. norm = 0.710605E+33 its = 44 res. norm = 0.386085E+33 its = 45 res. norm = 0.381319E+33 its = 46 res. norm = 0.312342E+33 its = 47 res. norm = 0.310614E+33 its = 48 res. norm = 0.285111E+33 its = 49 res. norm = 0.178166E+33 its = 50 res. norm = 0.177849E+33 its = 51 res. norm = 0.187294E+33 its = 52 res. norm = 0.187294E+33 its = 53 res. norm = 0.124642E+33 its = 54 res. norm = 0.980517E+32 its = 55 res. norm = 0.936919E+32 its = 56 res. norm = 0.932034E+32 its = 57 res. norm = 0.924707E+32 its = 58 res. norm = 0.758981E+32 its = 59 res. norm = 0.593626E+32 its = 60 res. norm = 0.521601E+32 its = 61 res. norm = 0.554573E+32 its = 62 res. norm = 0.554573E+32 its = 63 res. norm = 0.505556E+32 its = 64 res. norm = 0.364200E+32 its = 65 res. norm = 0.317523E+32 its = 66 res. norm = 0.235923E+32 its = 67 res. norm = 0.188534E+32 its = 68 res. norm = 0.188219E+32 its = 69 res. norm = 0.123948E+32 its = 70 res. norm = 0.123663E+32 its = 71 res. norm = 0.132854E+32 its = 72 res. norm = 0.132596E+32 its = 73 res. norm = 0.985483E+31 its = 74 res. norm = 0.938901E+31 its = 75 res. norm = 0.849346E+31 its = 76 res. norm = 0.566786E+31 its = 77 res. norm = 0.565627E+31 its = 78 res. norm = 0.379073E+31 pgmres returned with ierr = 0 Maximum error in solution = 47183636107.799881 5. Try ILUT Preconditioner with tol = .0001, lfil=7. Precon set-up returned with ierr 5 its = 0 res. norm = 0.498082E+38 its = 1 res. norm = 0.493748E+38 its = 2 res. norm = 0.465118E+38 its = 3 res. norm = 0.354091E+38 its = 4 res. norm = 0.337484E+38 its = 5 res. norm = 0.332807E+38 its = 6 res. norm = 0.312545E+38 its = 7 res. norm = 0.312507E+38 its = 8 res. norm = 0.145458E+38 its = 9 res. norm = 0.121935E+38 its = 10 res. norm = 0.201088E+37 its = 11 res. norm = 0.193035E+37 its = 12 res. norm = 0.193034E+37 its = 13 res. norm = 0.193005E+37 its = 14 res. norm = 0.193004E+37 its = 15 res. norm = 0.193000E+37 its = 16 res. norm = 0.193000E+37 its = 17 res. norm = 0.193000E+37 its = 18 res. norm = 0.192995E+37 its = 19 res. norm = 0.183465E+37 its = 20 res. norm = 0.175781E+37 its = 21 res. norm = 0.171270E+37 its = 22 res. norm = 0.171270E+37 its = 23 res. norm = 0.171270E+37 its = 24 res. norm = 0.171270E+37 its = 25 res. norm = 0.171270E+37 its = 26 res. norm = 0.171270E+37 its = 27 res. norm = 0.171270E+37 its = 28 res. norm = 0.171270E+37 its = 29 res. norm = 0.169124E+37 its = 30 res. norm = 0.154069E+37 its = 31 res. norm = 0.150238E+37 its = 32 res. norm = 0.150238E+37 its = 33 res. norm = 0.150238E+37 its = 34 res. norm = 0.150238E+37 its = 35 res. norm = 0.150238E+37 its = 36 res. norm = 0.150238E+37 its = 37 res. norm = 0.150238E+37 its = 38 res. norm = 0.150237E+37 its = 39 res. norm = 0.149905E+37 its = 40 res. norm = 0.143735E+37 its = 41 res. norm = 0.136313E+37 its = 42 res. norm = 0.136313E+37 its = 43 res. norm = 0.136313E+37 its = 44 res. norm = 0.136313E+37 its = 45 res. norm = 0.136313E+37 its = 46 res. norm = 0.136313E+37 its = 47 res. norm = 0.136313E+37 its = 48 res. norm = 0.136305E+37 its = 49 res. norm = 0.133038E+37 its = 50 res. norm = 0.130433E+37 its = 51 res. norm = 0.129957E+37 its = 52 res. norm = 0.129957E+37 its = 53 res. norm = 0.129957E+37 its = 54 res. norm = 0.129957E+37 its = 55 res. norm = 0.129957E+37 its = 56 res. norm = 0.129957E+37 its = 57 res. norm = 0.129957E+37 its = 58 res. norm = 0.129948E+37 its = 59 res. norm = 0.129806E+37 its = 60 res. norm = 0.125236E+37 its = 61 res. norm = 0.122293E+37 its = 62 res. norm = 0.122293E+37 its = 63 res. norm = 0.122293E+37 its = 64 res. norm = 0.122293E+37 its = 65 res. norm = 0.122293E+37 its = 66 res. norm = 0.122291E+37 its = 67 res. norm = 0.122289E+37 its = 68 res. norm = 0.122288E+37 its = 69 res. norm = 0.117724E+37 its = 70 res. norm = 0.105511E+37 its = 71 res. norm = 0.115509E+37 its = 72 res. norm = 0.115509E+37 its = 73 res. norm = 0.115509E+37 its = 74 res. norm = 0.115509E+37 its = 75 res. norm = 0.115509E+37 its = 76 res. norm = 0.115503E+37 its = 77 res. norm = 0.115500E+37 its = 78 res. norm = 0.114906E+37 its = 79 res. norm = 0.113300E+37 its = 80 res. norm = 0.107926E+37 its = 81 res. norm = 0.112450E+37 its = 82 res. norm = 0.112450E+37 its = 83 res. norm = 0.112450E+37 its = 84 res. norm = 0.112450E+37 its = 85 res. norm = 0.112450E+37 its = 86 res. norm = 0.112450E+37 its = 87 res. norm = 0.112449E+37 its = 88 res. norm = 0.112443E+37 its = 89 res. norm = 0.111159E+37 its = 90 res. norm = 0.108905E+37 its = 91 res. norm = 0.111859E+37 its = 92 res. norm = 0.111859E+37 its = 93 res. norm = 0.111859E+37 its = 94 res. norm = 0.111859E+37 its = 95 res. norm = 0.111859E+37 its = 96 res. norm = 0.111859E+37 its = 97 res. norm = 0.111859E+37 its = 98 res. norm = 0.111845E+37 its = 99 res. norm = 0.111814E+37 its = 100 res. norm = 0.106520E+37 its = 101 res. norm = 0.109150E+37 its = 102 res. norm = 0.109150E+37 its = 103 res. norm = 0.109150E+37 its = 104 res. norm = 0.109150E+37 its = 105 res. norm = 0.109150E+37 its = 106 res. norm = 0.109150E+37 its = 107 res. norm = 0.109150E+37 its = 108 res. norm = 0.109150E+37 its = 109 res. norm = 0.109143E+37 its = 110 res. norm = 0.109140E+37 pgmres returned with ierr = 1 Maximum error in solution = 15334309366194.100 SPARSEKIT_PRB13 Normal end of execution. 19 December 2023 2:40:43.442 PM