OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-28 18:21:00 -06:00
Code Issues Packages Projects Releases Wiki Activity

Add dune ILU when using only one process

Browse Source
This commit is contained in:
Tobias Meyer Andersen 2023-11-27 16:49:07 +01:00
parent 006aa14fbd
commit 8b5ab973e2
1 changed files with 16 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
opm/simulators/linalg/PreconditionerFactory_impl.hpp
View File

@ -167,6 +167,10 @@ struct StandardPreconditioners
using M = typename F::Matrix; using M = typename F::Matrix;
using V = typename F::Vector; using V = typename F::Vector;
using P = PropertyTree; using P = PropertyTree;
// TODO: add DuneILU with correct communicator
// F::addCreator("DuneILU", [](const O& op, const P& prm, const std::function<V()>&, std::size_t, const C& comm) {
// return wrapBlockPreconditioner<Dune::Preconditioner::SeqILU<M, V, V>>(comm, op.getmat());
// });
F::addCreator("ILU0", [](const O& op, const P& prm, const std::function<V()>&, std::size_t, const C& comm) { F::addCreator("ILU0", [](const O& op, const P& prm, const std::function<V()>&, std::size_t, const C& comm) {
return createParILU(op, prm, comm, 0); return createParILU(op, prm, comm, 0);
}); });
@ -176,6 +180,12 @@ struct StandardPreconditioners
F::addCreator("ILUn", [](const O& op, const P& prm, const std::function<V()>&, std::size_t, const C& comm) { F::addCreator("ILUn", [](const O& op, const P& prm, const std::function<V()>&, std::size_t, const C& comm) {
return createParILU(op, prm, comm, prm.get<int>("ilulevel", 0)); return createParILU(op, prm, comm, prm.get<int>("ilulevel", 0));
}); });
// F::addCreator("DuneILU", [](const O& op, const P& prm, const std::function<V()>&, std::size_t) {
// const double w = prm.get<double>("relaxation", 1.0);
// const int n = prm.get<int>("ilulevel", 0);
// const bool resort = prm.get<bool>("resort", false);
// return wrapBlockPreconditioner<Dune::SeqILU<M, V, V>>(op.getmat(), n, w, resort);
// });
F::addCreator("DILU", [](const O& op, const P& prm, const std::function<V()>&, std::size_t, const C& comm) { F::addCreator("DILU", [](const O& op, const P& prm, const std::function<V()>&, std::size_t, const C& comm) {
DUNE_UNUSED_PARAMETER(prm); DUNE_UNUSED_PARAMETER(prm);
return wrapBlockPreconditioner<MultithreadDILU<M, V, V>>(comm, op.getmat()); return wrapBlockPreconditioner<MultithreadDILU<M, V, V>>(comm, op.getmat());
@ -360,6 +370,12 @@ struct StandardPreconditioners<Operator,Dune::Amg::SequentialInformation>
return std::make_shared<Opm::ParallelOverlappingILU0<M, V, V, C>>( return std::make_shared<Opm::ParallelOverlappingILU0<M, V, V, C>>(
op.getmat(), 0, w, Opm::MILU_VARIANT::ILU); op.getmat(), 0, w, Opm::MILU_VARIANT::ILU);
}); });
F::addCreator("DuneILU", [](const O& op, const P& prm, const std::function<V()>&, std::size_t) {
const double w = prm.get<double>("relaxation", 1.0);
const int n = prm.get<int>("ilulevel", 0);
const bool resort = prm.get<bool>("resort", false);
return wrapPreconditioner<Dune::SeqILU<M, V, V>>(op.getmat(), n, w, resort);
});
F::addCreator("ParOverILU0", [](const O& op, const P& prm, const std::function<V()>&, std::size_t) { F::addCreator("ParOverILU0", [](const O& op, const P& prm, const std::function<V()>&, std::size_t) {
const double w = prm.get<double>("relaxation", 1.0); const double w = prm.get<double>("relaxation", 1.0);
const int n = prm.get<int>("ilulevel", 0); const int n = prm.get<int>("ilulevel", 0);