[dune-pdelab] Getting linear solver reduction from newton method

[René Heß]

Hi Amir,

Newton doesn't store the intermediate defects. You would need to create
a vector in the NewtonResult class and make sure that the intermediate
defects are stored in this vector (adjust the apply method in
NewtonSolver).

Side note: The old Newton is deprecated and might be removed at some
point so it might make sense to use the new class called
NewtonMethod.


Best regards,
René

Amir Peiraviminaei <apminaei at gmail.com> writes:

> Hi everyone,
> I'm using Newton method to solve a system of PDEs,
> using PDESolver = Dune::PDELab::Newton<IGO, LS, U>;
> Normally I get the first defect and final defect form
> PDESolver.result().first_defect and PDESolver.result().defect, but I don't
> know how to get defects in between.
> For example, if newton method converges in 5 iterations then how one can
> get defect_1, defect_2, .. , defect_4?
>
> Best wishes,
> Amir
> _______________________________________________
> dune-pdelab mailing list
> dune-pdelab at lists.dune-project.org
> https://lists.dune-project.org/mailman/listinfo/dune-pdelab
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 832 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20210217/235d3bba/attachment.sig>