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

Add property tree to constructor

Browse Source
This commit is contained in:
jakobtorben 2024-11-07 16:08:25 +01:00
parent 6d22b4f72a
commit 9d47aa605e
2 changed files with 10 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
opm/simulators/linalg/HyprePreconditioner.hpp
View File

@ -31,11 +31,11 @@
#include <vector>
#include <numeric>
namespace Dune {
namespace Hypre {
/// Wrapper for Hypre's BoomerAMG preconditioner
template<class M, class X, class Y>
class HyprePreconditioner : public PreconditionerWithUpdate<X,Y> {
class HyprePreconditioner : public Dune::PreconditionerWithUpdate<X,Y> {
public:
//! \brief The matrix type the preconditioner is for
using matrix_type = M;
@ -47,8 +47,8 @@ public:
using field_type = typename X::field_type;
// Constructor
HyprePreconditioner (const M& A)
: A_(A)
HyprePreconditioner (const M& A, const Opm::PropertyTree& prm)
: A_(A), prm_(prm)
{
OPM_TIMEBLOCK(prec_construct);
@ -129,8 +129,8 @@ public:
DUNE_UNUSED_PARAMETER(x);
}
SolverCategory::Category category() const override {
return SolverCategory::sequential;
Dune::SolverCategory::Category category() const override {
return Dune::SolverCategory::sequential;
}
bool hasPerfectUpdate() const override
@ -194,6 +194,8 @@ private:
}
const M& A_;
const Opm::PropertyTree& prm_;
HYPRE_Solver solver_ = nullptr;
HYPRE_IJMatrix A_hypre_ = nullptr;
HYPRE_ParCSRMatrix parcsr_A_ = nullptr;

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

@ -549,9 +549,8 @@ struct StandardPreconditioners<Operator, Dune::Amg::SequentialInformation> {
});
// Only add Hypre for scalar matrices
if constexpr (M::block_type::rows == 1 && M::block_type::cols == 1) {
F::addCreator("HypreBoomerAMG", [](const O& op, const P& prm, const std::function<V()>&, std::size_t) {
DUNE_UNUSED_PARAMETER(prm);
return std::make_shared<HyprePreconditioner<M, V, V>>(op.getmat());
F::addCreator("hypre", [](const O& op, const P& prm, const std::function<V()>&, std::size_t) {
return std::make_shared<Hypre::HyprePreconditioner<M, V, V>>(op.getmat(), prm);
});
}
}