templated_test(): Fortran77 version Test templated(). DETAILS OF ITERATIVE TEMPLATES TEST: Univ. of Tennessee and Oak Ridge National Laboratory October 1, 1993 Details of these algorithms are described in "Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, Eijkhout, Pozo, Romine, and van der Vorst, SIAM Publications, 1993. (ftp netlib2.cs.utk.edu; cd linalg; get templates.ps). MACHINE PRECISION = 1.11E-16 CONVERGENCE TEST TOLERANCE = 1.00E-15 For a detailed description of the following information, ===================================================== CONVERGENCE NORMALIZED NUM METHOD CRITERION RESIDUAL ITER INFO FLAG ===================================================== Order 36 SPD 2-d Poisson matrix (no preconditioning) CG 1.28E-16 3.08E-03 6 Chebyshev 6.17E-11 3.29E+02 144 1 SOR 2.18E-14 2.39E-01 144 1 BiCG 1.28E-16 3.08E-03 6 CGS NaN NaN 144 1 BiCGSTAB 1.29E-17 3.85E-03 6 GMRESm NaN NaN 144 1 QMR 2.58E-16 5.40E-03 6 Jacobi 3.02E-08 6.16E+05 144 1 CHEBYSHEV WARNING: COMPUTED MIN EIGENVALUE <= 0: SET TO EPSILON Order 36 SPD 2-d Poisson matrix (Jacobi preconditioning) CG 1.28E-16 3.08E-03 6 Chebyshev 1.20E-06 6.41E+06 144 1 SOR 2.18E-14 2.39E-01 144 1 BiCG 1.28E-16 3.08E-03 6 CGS NaN NaN 144 1 BiCGSTAB 1.29E-17 3.85E-03 6 GMRESm 1.44E+00 9.98E+12 144 1 QMR 2.58E-16 5.40E-03 6 Order 21 SPD Wathen matrix (no preconditioning) CG 2.74E-18 1.24E-03 26 Chebyshev 1.52E-02 2.32E+10 84 1 SOR 8.41E-11 3.28E+02 84 1 BiCG 2.74E-18 1.24E-03 26 CGS 4.07E-16 6.80E-02 27 BiCGSTAB 1.24E-16 1.55E-03 27 GMRESm NaN NaN 84 1 QMR 3.36E-16 1.08E-03 26 CHEBYSHEV WARNING: COMPUTED MIN EIGENVALUE <= 0: SET TO EPSILON Order 21 SPD Wathen matrix (Jacobi preconditioning) CG 9.10E-22 6.18E-04 20 Chebyshev 5.58E-01 5.28E+11 84 1 SOR 8.41E-11 3.28E+02 84 1 BiCG 9.10E-22 6.18E-04 20 CGS NaN NaN 84 1 BiCGSTAB 6.78E-16 2.86E+12 74 X GMRESm 5.86E+00 1.91E+12 84 1 QMR 4.12E-16 1.24E-03 20 Order 27 SPD 3-d Poisson matrix (no preconditioning) CG 2.83E-17 3.33E-03 4 Chebyshev 9.50E-16 1.41E-02 68 SOR 9.92E-16 2.50E-03 24 BiCG 2.83E-17 3.33E-03 4 CGS 1.34E-17 2.70E-03 4 BiCGSTAB 1.42E-18 3.33E-03 4 GMRESm NaN NaN 108 1 QMR 4.15E-16 6.45E-03 4 Jacobi 7.35E-16 1.33E-02 98 CHEBYSHEV WARNING: COMPUTED MIN EIGENVALUE <= 0: SET TO EPSILON Order 27 SPD 3-d Poisson matrix (Jacobi preconditioning) CG 2.41E-17 4.58E-03 4 Chebyshev 6.75E-14 4.74E-01 108 1 SOR 9.92E-16 2.50E-03 24 BiCG 2.41E-17 4.58E-03 4 CGS 4.09E-17 3.33E-03 4 BiCGSTAB 4.29E-18 4.99E-03 4 GMRESm 1.07E+00 6.00E+12 108 1 QMR 2.50E-16 3.74E-03 4 Order 125 PDE1 nonsymmetric matrix (no preconditioning) BiCG 2.13E-16 1.35E-03 65 CGS 8.11E-16 6.68E-03 91 BiCGSTAB 6.90E-16 5.12E-03 100 GMRESm NaN NaN 500 1 QMR 1.44E-14 1.41E-02 500 1 Order 125 PDE1 nonsymmetric matrix (Jacobi preconditioning) BiCG 7.38E-16 7.58E-02 60 CGS 3.76E-16 2.29E-02 71 BiCGSTAB 2.96E-16 4.05E-03 84 GMRESm NaN NaN 500 1 QMR 9.12E-14 6.65E-02 500 1 Order 125 PDE2 nonsymmetric matrix (no preconditioning) BiCG 7.58E-16 3.55E-02 31 CGS 9.46E-16 1.88E+00 40 BiCGSTAB 9.61E-03 1.36E+11 248 -10 GMRESm NaN NaN 500 1 QMR 1.03E-14 9.08E-02 500 1 Order 125 PDE2 nonsymmetric matrix (Jacobi preconditioning) BiCG 8.97E-16 4.95E-02 31 CGS 1.68E-16 1.87E+00 37 BiCGSTAB 2.62E-09 2.93E+04 500 1 GMRESm NaN NaN 500 1 QMR 1.44E-14 9.14E-02 500 1 Order 125 PDE3 nonsymmetric matrix (no preconditioning) BiCG 1.18E-14 3.86E-03 500 1 CGS 8.28E+12 4.14E+12 500 1 BiCGSTAB 1.11E+02 5.48E+11 116 -10 GMRESm NaN NaN 500 1 QMR 6.97E-13 8.26E-03 500 1 Order 125 PDE3 nonsymmetric matrix (Jacobi preconditioning) BiCG 8.28E-16 2.34E-03 443 CGS 2.64E+08 3.61E+12 500 1 BiCGSTAB 4.69E+00 9.06E+11 106 -10 GMRESm NaN NaN 500 1 QMR 6.01E-14 2.22E-03 500 1 Order 36 PDE4 nonsymmetric matrix (no preconditioning) BiCG 1.10E-16 1.11E-02 42 CGS 9.41E-01 4.60E+11 144 1 BiCGSTAB 1.25E-01 6.84E+10 144 1 GMRESm NaN NaN 144 1 QMR 1.47E-14 8.35E-03 144 1 Order 36 PDE4 nonsymmetric matrix (Jacobi preconditioning) BiCG 1.10E-16 1.11E-02 42 CGS 9.41E-01 4.60E+11 144 1 BiCGSTAB 1.25E-01 6.84E+10 144 1 GMRESm NaN NaN 144 1 QMR 1.47E-14 8.35E-03 144 1 SYMMETRIC POSITIVE DEFINITE MATRICES: ( 6 TESTS ) THERE ARE 1 CRITICAL ERRORS AND 19 SUSPICIOUS RESULTS. SEE FILE test.results FOR ADDITIONAL INFORMATION NONSYMMETRIC MATRICES: ( 8 TESTS ) THERE ARE 0 CRITICAL ERRORS AND 27 SUSPICIOUS RESULTS. SEE FILE test.results FOR ADDITIONAL INFORMATION templated_test(): Normal end of execution.