Skip to content

Draft: Implementation of non-negative Newton damping method.

Chaofan Chen requested to merge ChaofanChen/ogs:Non-negative_Newton_constraint into master

In this MR, a constraint of Non-negative Newton damping method is added to avoid negative solutions during Newton iterations. The constraint is similar to the Line search method, the damping factor formulation used in this method is

    x_{n+1} = x_n + \delta \frac{-f(x_n)}{f'(x_n)}

where

     \delta = \frac{1}{max(1.0, -\frac{\Delta x}{\alpha \cdot x})}

This work is part of the wellbore simulator process !4481

  1. Feature description was added to the changelog
  2. Tests covering your feature were added?
  3. Any new feature or behavior change was documented?
Edited by Dmitry Yu. Naumov

Merge request reports