OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

PreconditionerFactory: use Scalar type from operator

Browse Source
This commit is contained in:
Arne Morten Kvarving 2024-04-13 15:18:32 +02:00
parent b7bc7b7bf5
commit b9ee637d78
1 changed files with 27 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -50,11 +50,11 @@
#include <opm/simulators/linalg/PreconditionerFactoryGPUIncludeWrapper.hpp> #include <opm/simulators/linalg/PreconditionerFactoryGPUIncludeWrapper.hpp>
namespace Opm namespace Opm {
{
template <class Smoother> template <class Smoother>
struct AMGSmootherArgsHelper { struct AMGSmootherArgsHelper
{
static auto args(const PropertyTree& prm) static auto args(const PropertyTree& prm)
{ {
using SmootherArgs = typename Dune::Amg::SmootherTraits<Smoother>::Arguments; using SmootherArgs = typename Dune::Amg::SmootherTraits<Smoother>::Arguments;
@ -69,7 +69,8 @@ struct AMGSmootherArgsHelper {
}; };
template <class M, class V, class C> template <class M, class V, class C>
struct AMGSmootherArgsHelper<ParallelOverlappingILU0<M, V, V, C>> { struct AMGSmootherArgsHelper<ParallelOverlappingILU0<M, V, V, C>>
{
static auto args(const PropertyTree& prm) static auto args(const PropertyTree& prm)
{ {
using Smoother = ParallelOverlappingILU0<M, V, V, C>; using Smoother = ParallelOverlappingILU0<M, V, V, C>;
@ -88,7 +89,6 @@ struct AMGSmootherArgsHelper<ParallelOverlappingILU0<M, V, V, C>> {
} }
}; };
// trailing return type with decltype used for detecting existence of setUseFixedOrder member function by overloading the setUseFixedOrder function // trailing return type with decltype used for detecting existence of setUseFixedOrder member function by overloading the setUseFixedOrder function
template <typename C> template <typename C>
auto setUseFixedOrder(C criterion, bool booleanValue) -> decltype(criterion.setUseFixedOrder(booleanValue)) auto setUseFixedOrder(C criterion, bool booleanValue) -> decltype(criterion.setUseFixedOrder(booleanValue))
@ -279,7 +279,8 @@ struct StandardPreconditioners {
OPM_THROW(std::logic_error, OPM_THROW(std::logic_error,
"Pressure index out of bounds. It needs to specified for CPR"); "Pressure index out of bounds. It needs to specified for CPR");
} }
using LevelTransferPolicy = PressureTransferPolicy<O, Comm, double, false>; using Scalar = typename V::field_type;
using LevelTransferPolicy = PressureTransferPolicy<O, Comm, Scalar, false>;
return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy, Comm>>( return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy, Comm>>(
op, prm, weightsCalculator, pressureIndex, comm); op, prm, weightsCalculator, pressureIndex, comm);
}); });
@ -294,7 +295,8 @@ struct StandardPreconditioners {
OPM_THROW(std::logic_error, OPM_THROW(std::logic_error,
"Pressure index out of bounds. It needs to specified for CPR"); "Pressure index out of bounds. It needs to specified for CPR");
} }
using LevelTransferPolicy = PressureTransferPolicy<O, Comm, double, true>; using Scalar = typename V::field_type;
using LevelTransferPolicy = PressureTransferPolicy<O, Comm, Scalar, true>;
return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy, Comm>>( return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy, Comm>>(
op, prm, weightsCalculator, pressureIndex, comm); op, prm, weightsCalculator, pressureIndex, comm);
}); });
@ -311,7 +313,8 @@ struct StandardPreconditioners {
OPM_THROW(std::logic_error, OPM_THROW(std::logic_error,
"Pressure index out of bounds. It needs to specified for CPR"); "Pressure index out of bounds. It needs to specified for CPR");
} }
using LevelTransferPolicy = PressureBhpTransferPolicy<O, Comm, double, false>; using Scalar = typename V::field_type;
using LevelTransferPolicy = PressureBhpTransferPolicy<O, Comm, Scalar, false>;
return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy, Comm>>( return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy, Comm>>(
op, prm, weightsCalculator, pressureIndex, comm); op, prm, weightsCalculator, pressureIndex, comm);
}); });
@ -513,11 +516,16 @@ struct StandardPreconditioners<Operator, Dune::Amg::SequentialInformation> {
} }
}); });
F::addCreator("famg", [](const O& op, const P& prm, const std::function<V()>&, std::size_t) { F::addCreator("famg", [](const O& op, const P& prm, const std::function<V()>&, std::size_t) {
auto crit = AMGHelper<O, C, M, V>::criterion(prm); if constexpr (std::is_same_v<typename V::field_type, float>) {
Dune::Amg::Parameters parms; OPM_THROW(std::logic_error, "famg requires UMFPack which is not available for floats");
parms.setNoPreSmoothSteps(1); return nullptr;
parms.setNoPostSmoothSteps(1); } else {
return getRebuildOnUpdateWrapper<Dune::Amg::FastAMG<O, V>>(op, crit, parms); auto crit = AMGHelper<O, C, M, V>::criterion(prm);
Dune::Amg::Parameters parms;
parms.setNoPreSmoothSteps(1);
parms.setNoPostSmoothSteps(1);
return getRebuildOnUpdateWrapper<Dune::Amg::FastAMG<O, V>>(op, crit, parms);
}
}); });
} }
if constexpr (std::is_same_v<O, WellModelMatrixAdapter<M, V, V, false>>) { if constexpr (std::is_same_v<O, WellModelMatrixAdapter<M, V, V, false>>) {
@ -527,8 +535,9 @@ struct StandardPreconditioners<Operator, Dune::Amg::SequentialInformation> {
if (pressureIndex == std::numeric_limits<std::size_t>::max()) { if (pressureIndex == std::numeric_limits<std::size_t>::max()) {
OPM_THROW(std::logic_error, "Pressure index out of bounds. It needs to specified for CPR"); OPM_THROW(std::logic_error, "Pressure index out of bounds. It needs to specified for CPR");
} }
using Scalar = typename V::field_type;
using LevelTransferPolicy using LevelTransferPolicy
= PressureBhpTransferPolicy<O, Dune::Amg::SequentialInformation, double, false>; = PressureBhpTransferPolicy<O, Dune::Amg::SequentialInformation, Scalar, false>;
return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy>>( return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy>>(
op, prm, weightsCalculator, pressureIndex); op, prm, weightsCalculator, pressureIndex);
}); });
@ -540,7 +549,8 @@ struct StandardPreconditioners<Operator, Dune::Amg::SequentialInformation> {
if (pressureIndex == std::numeric_limits<std::size_t>::max()) { if (pressureIndex == std::numeric_limits<std::size_t>::max()) {
OPM_THROW(std::logic_error, "Pressure index out of bounds. It needs to specified for CPR"); OPM_THROW(std::logic_error, "Pressure index out of bounds. It needs to specified for CPR");
} }
using LevelTransferPolicy = PressureTransferPolicy<O, Dune::Amg::SequentialInformation, double, false>; using Scalar = typename V::field_type;
using LevelTransferPolicy = PressureTransferPolicy<O, Dune::Amg::SequentialInformation, Scalar, false>;
return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy>>( return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy>>(
op, prm, weightsCalculator, pressureIndex); op, prm, weightsCalculator, pressureIndex);
}); });
@ -550,7 +560,8 @@ struct StandardPreconditioners<Operator, Dune::Amg::SequentialInformation> {
if (pressureIndex == std::numeric_limits<std::size_t>::max()) { if (pressureIndex == std::numeric_limits<std::size_t>::max()) {
OPM_THROW(std::logic_error, "Pressure index out of bounds. It needs to specified for CPR"); OPM_THROW(std::logic_error, "Pressure index out of bounds. It needs to specified for CPR");
} }
using LevelTransferPolicy = PressureTransferPolicy<O, Dune::Amg::SequentialInformation, double, true>; using Scalar = typename V::field_type;
using LevelTransferPolicy = PressureTransferPolicy<O, Dune::Amg::SequentialInformation, Scalar, true>;
return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy>>( return std::make_shared<OwningTwoLevelPreconditioner<O, V, LevelTransferPolicy>>(
op, prm, weightsCalculator, pressureIndex); op, prm, weightsCalculator, pressureIndex);
}); });