Well I could try to reproduce this example with older version of deal.II 
but this also requires a specific version of Petsc. I forgot to mention 
above that deal.II is compiled with petsc-3.15.0. A hint on the timeline, 
e.g. the day the results are produced, could be useful.

Best
Giang

On Wednesday, June 2, 2021 at 9:45:50 PM UTC+2 Wolfgang Bangerth wrote:

> On 6/2/21 3:07 AM, [email protected] wrote: 
> > 
> > compile but does not run correctly with deal.II 9.2.0 . The analysis can 
> run 
> > until step 204 and then diverge. Would it be possible changes in 
> > preconditioner destroy the convergence? 
> > 
> > 
> ************************************************************************************************
>  
>
> > Time step = 201, at t = 2.010000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 4.052763e-10 
> > 
> ************************************************************************************************
>  
>
> > Time step = 202, at t = 2.020000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 7.194398e-10 
> > 
> ************************************************************************************************
>  
>
> > Time step = 203, at t = 2.030000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 1.361073e-09 
> > 
> ************************************************************************************************
>  
>
> > Time step = 204, at t = 2.040000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 4.896128e-09 
> > 
> ************************************************************************************************
>  
>
> > Time step = 205, at t = 2.050000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 3.220123e-08 
> > 
> ************************************************************************************************
>  
>
> > Time step = 206, at t = 2.060000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 7.647033e-07 
> > 
> ************************************************************************************************
>  
>
> > Time step = 207, at t = 2.070000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 3.355695e-04 
> > 
> ************************************************************************************************
>  
>
> > Time step = 208, at t = 2.080000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 5.007715e+01 
> > 
> ************************************************************************************************
>  
>
> > Time step = 209, at t = 2.090000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 1.262003e+12 
> > 
> ************************************************************************************************
>  
>
> > Time step = 210, at t = 2.100000e-01 
> >  GMRES_ITR = 8   GMRES_RES = 1.279221e+33 
>
> It looks like the GMRES solver converged just fine, but that the solution 
> diverges anyway. 
>
> Could you try to run the program with previous versions of deal.II? Did it 
> work correctly back then? What version is the last one that worked? We may 
> also want to rope in the author of the program, Jie Cheng (CC:ed), maybe 
> he knows! 
>
> Best 
> W. 
>
> -- 
> ------------------------------------------------------------------------ 
> Wolfgang Bangerth email: [email protected] 
> www: http://www.math.colostate.edu/~bangerth/ 
>
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/c50b2468-94b5-4ad5-8064-36c08de38ee4n%40googlegroups.com.

Reply via email to