<div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small"><span style="background-color:rgb(255,255,255)"><font color="#444444">Hi everyone,<br></font></span></div><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small"><span style="background-color:rgb(255,255,255)"><font color="#444444">I'm using </font></span></div><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small"><div style="font-size:14px;line-height:19px;white-space:pre"><div style=""><font style="background-color:rgb(255,255,255)" color="#444444">Dune::PDELab::ISTLBackend_BCGS_AMG_SSOR<IGO> LS</font></div><div style="">as linear solver within newton method to slove a system of PDEs discretized with DG method.</div><div style="">the problem is when it build a hierarchy of at least 2 levels, it does not converge. </div><div style="">this is the message I get:</div><div style=""><br></div><div style="">When AMG preconditioner sets for 2 levels:<br><br>TIME STEP [implicit Euler]      1 time (from):   0.0000e+00 dt:   3.1450e+01 time (to):   3.1450e+01<br>STAGE 1 time (to):   3.1450e+01.<br>  Initial defect:   1.2878e+01<br>Using a direct coarse solver (UMFPack)<br>Building hierarchy of 2 levels (including coarse solver) took 1.351 seconds.<br>=== Dune::BiCGSTABSolver<br>=== rate=0.9541, T=72.91, TIT=0.7291, IT=100<br>Catched Error, Dune reported error: NewtonLinearSolverError [linearSolve:../dune-pdelab/dune/pdelab/newton/newton.hh:244]: NewtonSolver::linearSolve(): Linear solver did not converge in 100 iterations<br><br><br><br>When AMG preconditioner sets for 1 level:<br><br>TIME STEP [implicit Euler]      1 time (from):   0.0000e+00 dt:   3.1450e+01 time (to):   3.1450e+01<br>STAGE 1 time (to):   3.1450e+01.<br>  Initial defect:   1.7495e+01<br>Using a direct coarse solver (UMFPack)<br>Building hierarchy of 1 levels (including coarse solver) took 12.65 seconds.<br>=== Dune::BiCGSTABSolver<br>=== rate=7.55e-29, T=0.6516, TIT=inf, IT=0<br>  Newton iteration  1.  New defect:   4.8970e-03.  Reduction (this):   2.7991e-04.  Reduction (total):   2.7991e-04<br></div><div style=""><br></div><div style="">Is it possible that this implementation of Linear solver(LS) does not work with DG method?</div><div style="">How to use linear solvers from dune-solvers module for newton method from dune-pdelab?</div><div style=""><br></div><div style="">Thank you.</div><div style="">Amir</div></div></div></div>