Scroll to navigation

damped_newton(4rheolef) rheolef-7.0 damped_newton(4rheolef)

NAME

damped_newton -- damped Newton nonlinear algorithm

DESCRIPTION

Nonlinear damped Newton algorithm for the resolution of the following problem:



       F(u) = 0
A simple call to the algorithm writes: 



    my_problem P;


    field uh (Vh);


    damped_newton (P, uh, tol, max_iter);
In addition to the members required for the Newton algorithm (see newton(4)), the space_norm is required for the damped Newton line search algorithm: 



    class my_problem {


    public:


      ...


      value_type derivative_trans_mult (const value_type& mrh) const;


      Float space_norm (const value_type& uh) const;


    };

IMPLEMENTATION

template <class Problem, class Field, class Real, class Size>
int damped_newton (const Problem& F, Field& u, Real& tol, Size& max_iter, odiststream* p_derr=0)

SEE ALSO

newton(4)

COPYRIGHT

Copyright (C) 2000-2018 Pierre Saramito <Pierre.Saramito@imag.fr> GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>. This is free software: you are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law.

rheolef-7.0 rheolef-7.0