On 8/28/07, Matthew Knepley <knepley at gmail.com> wrote: > I will ask Todd about gltr since it might be supposed to > do something > funny. You really do not want to look at it.
Now I understand way you said that!! Well, let me say that actually GLTR is almost OK, the only two things to remark are: 1- KSPMonitor is called allways with it=0; perhaps this is OK, but perhaps no. You get always output like '0 KSP Residual norm' 2 - At some point, it converges and returns without incrementing the 'ksp->its' counter. That's the reason I though there could be a problem. However, as the implementation is a bit contrived and 'ksp->its' is used everywere, I think I'll leave GLTR for review by the author, please FWD this mail to him... Regards, > > Matt > > > Below the result (petsc4py is a nice tool for test/debug, isn't it?) > > > > > > tfqmr - CONVERGED_ITS - iters: 5, ||r|| = 4.881889e-03, > > x0,N-1=0.999996,9.978373 > > minres - CONVERGED_ITS - iters: 5, ||r|| = 5.356222e-02, > > x0,N-1=1.000009,9.701171 > > fgmres - DIVERGED_ITS - iters: 5, ||r|| = 5.356222e-02, > > x0,N-1=1.000009,9.701171 > > stcg - CONVERGED_ITS - iters: 5, ||r|| = 5.629995e-02, > > x0,N-1=1.000052,9.771470 > > qcg - CONVERGED_ITS - iters: 5, ||r|| = 5.629995e-02, > > x0,N-1=-1.000052,-9.771470 > > cg - CONVERGED_ITS - iters: 5, ||r|| = 5.629995e-02, > > x0,N-1=1.000052,9.771470 > > lgmres - DIVERGED_ITS - iters: 5, ||r|| = 5.356222e-02, > > x0,N-1=1.000009,9.701171 > > cgne - CONVERGED_ITS - iters: 5, ||r|| = 7.192229e-02, > > x0,N-1=1.000000,7.116166 > > chebychev - DIVERGED_ITS - iters: 5, ||r|| = 2.591834e+00, > > x0,N-1=0.636559,0.708271 > > cgs - CONVERGED_ITS - iters: 5, ||r|| = 1.457830e-03, > > x0,N-1=1.000000,9.994777 > > bicg - CONVERGED_ITS - iters: 5, ||r|| = 5.629995e-02, > > x0,N-1=1.000052,9.771470 > > lsqr - CONVERGED_ITS - iters: 5, ||r|| = 4.612376e-01, > > x0,N-1=1.000000,7.116166 > > gltr - CONVERGED_ITS - iters: 4, ||r|| = 5.629995e-02, > > x0,N-1=1.000052,9.771470 > > tcqmr - DIVERGED_ITS - iters: 5, ||r|| = 0.000000e+00, > > x0,N-1=1.000009,9.701171 > > bcgs - CONVERGED_ITS - iters: 5, ||r|| = 2.566301e-03, > > x0,N-1=0.999703,9.982615 > > cr - CONVERGED_ITS - iters: 5, ||r|| = 5.356222e-02, > > x0,N-1=1.000009,9.701171 > > symmlq - CONVERGED_ITS - iters: 5, ||r|| = 5.629995e-02, > > x0,N-1=1.000052,9.771470 > > bcgsl - CONVERGED_ITS - iters: 5, ||r|| = 1.232695e-03, > > x0,N-1=0.999988,9.999315 > > lcd - DIVERGED_ITS - iters: 5, ||r|| = 5.629995e-02, > > x0,N-1=1.000052,9.771470 > > preonly - CONVERGED_ITS - iters: 1, ||r|| = 0.000000e+00, > > x0,N-1=1.000000,1.000000 > > gmres - DIVERGED_ITS - iters: 5, ||r|| = 5.356222e-02, > > x0,N-1=1.000009,9.701171 > > richardson - DIVERGED_ITS - iters: 5, ||r|| = 1.215430e+00, > > x0,N-1=1.000000,4.095100 > > > > > > > > > > > > > > On 8/28/07, Lisandro Dalcin <dalcinl at gmail.com> wrote: > > > > Does it make sense to change KSPSkipConverged to return > > > > KSP_CONVERGED_ITS if iternum==maxit ? > > > > > > > > KSP_DIVERGED_ITS means convergence failure, but IMHO, KSPSkipConverged > > > > should not imply convergence failure (this has implications in SNES). > > > > > > > > > > > > -- > > > > Lisandro Dalc?n > > > > --------------- > > > > Centro Internacional de M?todos Computacionales en Ingenier?a (CIMEC) > > > > Instituto de Desarrollo Tecnol?gico para la Industria Qu?mica (INTEC) > > > > Consejo Nacional de Investigaciones Cient?ficas y T?cnicas (CONICET) > > > > PTLC - G?emes 3450, (3000) Santa Fe, Argentina > > > > Tel/Fax: +54-(0)342-451.1594 > > > > > > > > > > > > > > > > > -- > > > What most experimenters take for granted before they begin their > > > experiments is infinitely more interesting than any results to which > > > their experiments lead. > > > -- Norbert Wiener > > > > > > > > > > > > -- > > Lisandro Dalc?n > > --------------- > > Centro Internacional de M?todos Computacionales en Ingenier?a (CIMEC) > > Instituto de Desarrollo Tecnol?gico para la Industria Qu?mica (INTEC) > > Consejo Nacional de Investigaciones Cient?ficas y T?cnicas (CONICET) > > PTLC - G?emes 3450, (3000) Santa Fe, Argentina > > Tel/Fax: +54-(0)342-451.1594 > > > > > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which > their experiments lead. > -- Norbert Wiener > > -- Lisandro Dalc?n --------------- Centro Internacional de M?todos Computacionales en Ingenier?a (CIMEC) Instituto de Desarrollo Tecnol?gico para la Industria Qu?mica (INTEC) Consejo Nacional de Investigaciones Cient?ficas y T?cnicas (CONICET) PTLC - G?emes 3450, (3000) Santa Fe, Argentina Tel/Fax: +54-(0)342-451.1594