Commit 1464f051 authored by Dominic Kempf's avatar Dominic Kempf

Disable line search in those case where we already are at the solution

This potentially fixes #117.
parent 5c670ca6
Pipeline #15997 passed with stage
in 16 minutes and 11 seconds
......@@ -641,7 +641,7 @@ namespace Dune
virtual void line_search(TrialVector& z, TestVector& r) override
{
if (strategy_ == noLineSearch)
if ((strategy_ == noLineSearch) || (this->res_.defect < this->abs_limit_))
{
this->u_->axpy(-1.0, z);
this->defect(r);
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment