KSPSkipConverged

2007-08-29 Thread Todd Munson
ED_ITS- iters: 5, ||r|| = 1.232695e-03, > > > x0,N-1=0.88,9.999315 > > > lcd- DIVERGED_ITS - iters: 5, ||r|| = 5.629995e-02, > > > x0,N-1=1.52,9.771470 > > > preonly- CONVERGED_ITS- iters: 1, ||r|| =

KSPSkipConverged

2007-08-28 Thread Lisandro Dalcin
> lcd- DIVERGED_ITS - iters: 5, ||r|| = 5.629995e-02, > > x0,N-1=1.52,9.771470 > > preonly- CONVERGED_ITS- iters: 1, ||r|| = 0.00e+00, > > x0,N-1=1.00,1.000000 > > gmres - DIVERGED_ITS - iters

KSPSkipConverged

2007-08-28 Thread Lisandro Dalcin
onverged), so perhaps there are a problem with the loop index. > > Okay, should be fixed in dev. > > > Matt, from your previous mail, something is not clear to me: > > > > Can I modify KSPSkipConverged for release-2.3.3 ? > > Yes, that seems fine. Have you made a

KSPSkipConverged

2007-08-28 Thread Lisandro Dalcin
em are special purpose. However, GLTR does not seems to stopping because of happy breakdown, it always stop at maxit-1 (with my skip converged), so perhaps there are a problem with the loop index. Matt, from your previous mail, something is not clear to me: Can I modify KSPSkipConverged for

KSPSkipConverged

2007-08-28 Thread Lisandro Dalcin
eration 4 and not 5, as it should). 2- Modify KSPSkipConverged and push on petsc-dev. Or perhpas can we also push this on release-2.3.3? The previous way is rather buggy, especially in conjunction with KSP_NORM_NO. Below the result (petsc4py is a nice tool for test/debug, isn't it?) tfqmr

KSPSkipConverged

2007-08-28 Thread Lisandro Dalcin
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

KSPSkipConverged

2007-08-28 Thread Matthew Knepley
in dev. > Matt, from your previous mail, something is not clear to me: > > Can I modify KSPSkipConverged for release-2.3.3 ? Yes, that seems fine. Have you made a patch to release before? Thanks, Matt > > > > Matt > > > > On 8/28/07, Matthew Knepley wr

KSPSkipConverged

2007-08-28 Thread Matthew Knepley
s to do > > this). It also seems that I have to review KSP type GLTR (it stopped > > at iteration 4 and not 5, as it should). > > > > 2- Modify KSPSkipConverged and push on petsc-dev. Or perhpas can we > > also push this on release-2.3.3? The previous way is rather bugg

KSPSkipConverged

2007-08-28 Thread Matthew Knepley
and TCQMR seems to do > this). It also seems that I have to review KSP type GLTR (it stopped > at iteration 4 and not 5, as it should). > > 2- Modify KSPSkipConverged and push on petsc-dev. Or perhpas can we > also push this on release-2.3.3? The previous way is rather buggy, > esp

KSPSkipConverged

2007-08-28 Thread Matthew Knepley
Yes, definitely. Go ahead and push it. Matt On 8/28/07, Lisandro Dalcin 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 conv