********** TUTORIAL MULTIGRID LABORATORY ******************** ib1 ib2 ib3 ib4 ib5 ib6 ib7 ib8 ib9 ib10 ib11 ib12 ib13 1 1 2 2 0 2 4 1 0 3 0 3 4 icase igal iout ism levels nc 2 0 4 0 0 4 2 4 2 nmg ns(1) ns(2) 10 1 0 aa bb cc ss f0 f1 0.100E-01 0.100E+01 0.000E+00 0.000E+00-0.100E+01 0.200E+01 om spar delta tol 0.500E+00 0.100E+01 0.100E-01 0.100E-06 cycle = V grid pointer array 1 1 2 3 3 7 4 15 matrix at grid 1 1 0.000E+00 0.106E+01 -0.200E-01 2 -0.102E+01 0.106E+01 0.000E+00 matrix at grid 2 1 0.000E+00 0.112E+01 -0.400E-01 2 -0.104E+01 0.108E+01 -0.400E-01 3 -0.104E+01 0.108E+01 -0.400E-01 4 -0.104E+01 0.112E+01 0.000E+00 matrix at grid 3 1 0.000E+00 0.124E+01 -0.800E-01 2 -0.108E+01 0.116E+01 -0.800E-01 3 -0.108E+01 0.116E+01 -0.800E-01 4 -0.108E+01 0.116E+01 -0.800E-01 5 -0.108E+01 0.116E+01 -0.800E-01 6 -0.108E+01 0.116E+01 -0.800E-01 7 -0.108E+01 0.116E+01 -0.800E-01 8 -0.108E+01 0.124E+01 0.000E+00 matrix at grid 4 1 0.000E+00 0.148E+01 -0.160E+00 2 -0.116E+01 0.132E+01 -0.160E+00 3 -0.116E+01 0.132E+01 -0.160E+00 4 -0.116E+01 0.132E+01 -0.160E+00 5 -0.116E+01 0.132E+01 -0.160E+00 6 -0.116E+01 0.132E+01 -0.160E+00 7 -0.116E+01 0.132E+01 -0.160E+00 8 -0.116E+01 0.132E+01 -0.160E+00 9 -0.116E+01 0.132E+01 -0.160E+00 10 -0.116E+01 0.132E+01 -0.160E+00 11 -0.116E+01 0.132E+01 -0.160E+00 12 -0.116E+01 0.132E+01 -0.160E+00 13 -0.116E+01 0.132E+01 -0.160E+00 14 -0.116E+01 0.132E+01 -0.160E+00 15 -0.116E+01 0.132E+01 -0.160E+00 16 -0.116E+01 0.148E+01 0.000E+00 Right hand side vector at grid 4 1 -0.132E+01 2 0.000E+00 3 0.000E+00 4 0.000E+00 5 0.000E+00 6 0.000E+00 7 0.000E+00 8 0.000E+00 9 0.000E+00 10 0.000E+00 11 0.000E+00 12 0.000E+00 13 0.000E+00 14 0.000E+00 15 0.000E+00 16 0.640E+00 matrix at grid 1 1 0.000E+00 0.000E+00 0.000E+00 2 0.000E+00 0.000E+00 0.000E+00 matrix at grid 2 1 0.000E+00 0.000E+00 0.000E+00 2 0.000E+00 0.000E+00 0.000E+00 3 0.000E+00 0.000E+00 0.000E+00 4 0.000E+00 0.000E+00 0.000E+00 matrix at grid 3 1 0.000E+00 0.000E+00 0.000E+00 2 0.000E+00 0.000E+00 0.000E+00 3 0.000E+00 0.000E+00 0.000E+00 4 0.000E+00 0.000E+00 0.000E+00 5 0.000E+00 0.000E+00 0.000E+00 6 0.000E+00 0.000E+00 0.000E+00 7 0.000E+00 0.000E+00 0.000E+00 8 0.000E+00 0.000E+00 0.000E+00 matrix at grid 4 1 0.000E+00 0.980E+00 0.340E+00 2 -0.660E+00 0.320E+00 0.340E+00 3 -0.660E+00 0.320E+00 0.340E+00 4 -0.660E+00 0.320E+00 0.340E+00 5 -0.660E+00 0.320E+00 0.340E+00 6 -0.660E+00 0.320E+00 0.340E+00 7 -0.660E+00 0.320E+00 0.340E+00 8 -0.660E+00 0.320E+00 0.340E+00 9 -0.660E+00 0.320E+00 0.340E+00 10 -0.660E+00 0.320E+00 0.340E+00 11 -0.660E+00 0.320E+00 0.340E+00 12 -0.660E+00 0.320E+00 0.340E+00 13 -0.660E+00 0.320E+00 0.340E+00 14 -0.660E+00 0.320E+00 0.340E+00 15 -0.660E+00 0.320E+00 0.340E+00 16 -0.660E+00 -0.200E-01 0.000E+00 Right hand side vector at grid 4 1 -0.132E+01 2 0.000E+00 3 0.000E+00 4 0.000E+00 5 0.000E+00 6 0.000E+00 7 0.000E+00 8 0.000E+00 9 0.000E+00 10 0.000E+00 11 0.000E+00 12 0.000E+00 13 0.000E+00 14 0.000E+00 15 0.000E+00 16 -0.136E+01 Exact discrete solution on finest grid vector at grid 4 1 -0.100E+01 2 -0.100E+01 3 -0.999E+00 4 -0.100E+01 5 -0.996E+00 6 -0.101E+01 7 -0.984E+00 8 -0.103E+01 9 -0.939E+00 10 -0.112E+01 11 -0.769E+00 12 -0.145E+01 13 -0.129E+00 14 -0.269E+01 15 0.228E+01 16 -0.738E+01 l2-norm of initial residual 0.367E+00 iteration number 1 reduction factor 0.263E+00 l2-norm of residual 0.966E-01 l2-norm of res. when leaving mg 0.966E-01 total iterations 1 average reduction factor 0.263E+00 Right hand side for defect correction vector at grid 4 1 -0.137E+01 2 0.797E-01 3 -0.674E-01 4 0.906E-01 5 -0.811E-01 6 0.368E-01 7 -0.294E-01 8 0.626E-01 9 -0.569E-01 10 0.126E-01 11 -0.823E-02 12 0.679E-02 13 -0.340E-02 14 -0.456E-01 15 -0.175E+00 16 0.531E+00 l2-norm of initial residual 0.752E-01 iteration number 1 reduction factor 0.213E+00 l2-norm of residual 0.160E-01 l2-norm of res. when leaving mg 0.160E-01 total iterations 1 average reduction factor 0.213E+00 Right hand side for defect correction vector at grid 4 1 -0.135E+01 2 0.516E-01 3 -0.544E-01 4 0.645E-01 5 -0.476E-01 6 0.240E-01 7 -0.315E-01 8 0.497E-01 9 -0.340E-01 10 0.131E-01 11 -0.129E-01 12 0.250E-01 13 -0.122E-01 14 0.273E-01 15 -0.327E+00 16 0.401E+00 l2-norm of initial residual 0.510E-01 iteration number 1 reduction factor 0.116E+00 l2-norm of residual 0.594E-02 l2-norm of res. when leaving mg 0.594E-02 total iterations 1 average reduction factor 0.116E+00 Right hand side for defect correction vector at grid 4 1 -0.134E+01 2 0.386E-01 3 -0.434E-01 4 0.434E-01 5 -0.282E-01 6 0.190E-01 7 -0.283E-01 8 0.312E-01 9 -0.170E-01 10 0.618E-02 11 -0.939E-02 12 0.190E-01 13 -0.152E-01 14 0.107E+00 15 -0.367E+00 16 0.231E+00 l2-norm of initial residual 0.471E-01 iteration number 1 reduction factor 0.425E-01 l2-norm of residual 0.200E-02 l2-norm of res. when leaving mg 0.200E-02 total iterations 1 average reduction factor 0.425E-01 Multigrid solution on finest grid vector at grid 4 1 -0.101E+01 2 -0.983E+00 3 -0.101E+01 4 -0.981E+00 5 -0.100E+01 6 -0.990E+00 7 -0.101E+01 8 -0.986E+00 9 -0.100E+01 10 -0.995E+00 11 -0.998E+00 12 -0.985E+00 13 -0.994E+00 14 -0.938E+00 15 -0.119E+01 16 -0.781E+00 Local error on finest grid vector at grid 4 1 -0.117E-01 2 0.177E-01 3 -0.157E-01 4 0.213E-01 5 -0.807E-02 6 0.183E-01 7 -0.269E-01 8 0.460E-01 9 -0.613E-01 10 0.124E+00 11 -0.230E+00 12 0.464E+00 13 -0.865E+00 14 0.175E+01 15 -0.348E+01 16 0.659E+01 Maximum and L2 norm of the error = 0.65947E+01 0.19316E+01 Work unit distribution over the grids 1 0.50E+00 2 0.10E+01 3 0.20E+01 4 0.40E+01 .