Skip to content
Snippets Groups Projects
Commit d7803af7 authored by Timo Koch's avatar Timo Koch
Browse files

[patches] Update istl patch to 2.6

parent db599111
No related branches found
No related tags found
1 merge request!1454[patches] Update istl patch to 2.6
- If the AMGBackend should be used without SuperLU as coarse grid solver, it can
be benefitial to decrease the corresponding tolerance. To do so, apply the
patch istl-2.4.1.patch in your directory containing
dune-istl 2.4.1, e.g.:
patch -p1 <../dumux/patches/istl-2.4.1.patch
If the AMGBackend should be used without SuperLU as coarse grid solver, it can
be benefitial to decrease the corresponding tolerance. To do so, apply the patch
- `istl-2.6.patch` in your directory containing dune-istl 2.6
```bash
patch -p1 <../dumux/patches/istl-2.6.patch
```
diff --git a/dune/istl/paamg/amg.hh b/dune/istl/paamg/amg.hh diff --git a/dune/istl/paamg/amg.hh b/dune/istl/paamg/amg.hh
index 56b1d8d..d672756 100644 index b4b90a97..5bd66cda 100644
--- a/dune/istl/paamg/amg.hh --- a/dune/istl/paamg/amg.hh
+++ b/dune/istl/paamg/amg.hh +++ b/dune/istl/paamg/amg.hh
@@ -526,13 +526,13 @@ namespace Dune @@ -523,14 +523,14 @@ namespace Dune
// We are still participating on this level // in order to ensure that we fulfill the alignement requirements
solver_.reset(new BiCGSTABSolver<X>(const_cast<M&>(matrices_->matrices().coarsest().getRedistributed()), solver_.reset(new BiCGSTABSolver<X>(const_cast<M&>(matrices_->matrices().coarsest().getRedistributed()),
*scalarProduct_, *scalarProduct_,
- *coarseSmoother_, 1E-2, 1000, 0)); - *coarseSmoother_, 1E-2, 1000, 0));
...@@ -11,10 +11,11 @@ index 56b1d8d..d672756 100644 ...@@ -11,10 +11,11 @@ index 56b1d8d..d672756 100644
else else
solver_.reset(); solver_.reset();
}else }else
solver_.reset(new BiCGSTABSolver<X>(const_cast<M&>(*matrices_->matrices().coarsest()), {
*scalarProduct_, solver_.reset(new BiCGSTABSolver<X>(const_cast<M&>(*matrices_->matrices().coarsest()),
- *coarseSmoother_, 1E-2, 1000, 0)); *scalarProduct_,
+ *coarseSmoother_, 1E-12, 1000, 0)); - *coarseSmoother_, 1E-2, 1000, 0));
} + *coarseSmoother_, 1E-12, 1000, 0));
} // // we have to allocate these types using the rebound allocator
// // in order to ensure that we fulfill the alignement requirements
// using Alloc = typename A::template rebind<BiCGSTABSolver<X>>::other;
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment