templates_test(): Fortran77 version Test templates(). 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 = 5.96E-08 CONVERGENCE TEST TOLERANCE = 1.00E-07 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 5.56E-08 1.66E-02 6 Chebyshev 8.79E-08 3.73E-02 97 SOR 8.66E-08 1.45E-02 72 BiCG 5.56E-08 1.66E-02 6 CGS 2.16E-05 1.22E+00 122 BiCGSTAB 1.42E-08 1.86E-02 6 GMRESm 3.47E+05 2.32E+05 144 1 QMR 9.93E-08 2.07E-02 6 Jacobi 9.21E-08 2.90E-02 136 WARNING: MIN EIGENVALUE < 0: SET TO ZERO Order 36 SPD 2-d Poisson matrix (Jacobi preconditioning) CG 5.56E-08 1.66E-02 6 Chebyshev 1.20E-06 1.16E-01 144 1 SOR 8.66E-08 1.45E-02 72 BiCG 5.56E-08 1.66E-02 6 CGS 2.16E-05 1.22E+00 122 BiCGSTAB 1.42E-08 1.86E-02 6 GMRESm 1.53E+00 6.10E+04 144 1 QMR 9.93E-08 2.07E-02 6 Order 21 SPD Wathen matrix (no preconditioning) CG 8.47E-08 6.64E-03 23 Chebyshev 1.52E-02 2.32E+02 84 1 SOR 9.65E-08 4.15E-03 56 BiCG 8.47E-08 6.64E-03 23 CGS 1.79E-08 3.17E-01 48 BiCGSTAB 4.07E-08 1.16E-02 21 GMRESm NaN NaN 84 1 QMR 2.99E-07 4.98E-03 84 1 WARNING: MIN EIGENVALUE < 0: SET TO ZERO Order 21 SPD Wathen matrix (Jacobi preconditioning) CG 4.36E-08 5.81E-03 15 Chebyshev 5.58E-01 5.28E+03 84 1 SOR 9.65E-08 4.15E-03 56 BiCG 4.36E-08 5.81E-03 15 CGS 3.73E+10 7.76E+04 84 1 BiCGSTAB 1.14E-08 3.95E+04 51 X GMRESm 2.37E+01 1.88E+04 84 1 QMR 1.06E-07 9.96E-03 84 1 Order 27 SPD 3-d Poisson matrix (no preconditioning) CG 1.16E-08 8.94E-03 4 Chebyshev 7.01E-08 2.68E-02 30 SOR 9.87E-08 1.79E-02 12 BiCG 1.16E-08 8.94E-03 4 CGS 5.67E-09 1.79E-02 4 BiCGSTAB 6.69E-10 1.12E-02 4 GMRESm NaN NaN 108 1 QMR 4.08E-08 1.79E-02 4 Jacobi 6.94E-08 1.12E-02 47 WARNING: MIN EIGENVALUE < 0: SET TO ZERO Order 27 SPD 3-d Poisson matrix (Jacobi preconditioning) CG 8.26E-08 1.12E-02 4 Chebyshev 8.94E-08 1.01E-01 51 SOR 9.87E-08 1.79E-02 12 BiCG 8.26E-08 1.12E-02 4 CGS 1.87E-08 1.79E-02 4 BiCGSTAB 4.67E-10 2.01E-02 4 GMRESm 9.80E-01 5.56E+04 108 1 QMR 1.42E-07 2.01E-02 108 1 Order 36 PDE4 nonsymmetric matrix (no preconditioning) BiCG 5.31E-08 3.36E-02 48 CGS 5.49E-01 1.86E+03 144 1 BiCGSTAB 3.05E+00 2.56E+04 103 -10 GMRESm NaN NaN 144 1 QMR 7.59E-06 4.11E-02 144 1 Order 36 PDE4 nonsymmetric matrix (Jacobi preconditioning) BiCG 5.31E-08 3.36E-02 48 CGS 5.49E-01 1.86E+03 144 1 BiCGSTAB 3.05E+00 2.56E+04 103 -10 GMRESm NaN NaN 144 1 QMR 7.59E-06 4.11E-02 144 1 SYMMETRIC POSITIVE DEFINITE MATRICES: ( 6 TESTS ) THERE ARE 1 CRITICAL ERRORS AND 13 SUSPICIOUS RESULTS. SEE FILE test.results FOR ADDITIONAL INFORMATION NONSYMMETRIC MATRICES: ( 2 TESTS ) THERE ARE 0 CRITICAL ERRORS AND 8 SUSPICIOUS RESULTS. SEE FILE test.results FOR ADDITIONAL INFORMATION templates_test(): Normal end of execution.