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 forceSerial parameter to ISTLSolverEbos.

Browse Source 
Used for single-domain solves.
This commit is contained in:
Atgeirr Flø Rasmussen 2023-08-22 16:33:18 +02:00
parent 622103a367
commit e1dd2bf148
3 changed files with 11 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/simulators/flow/BlackoilModelEbosNldd.hpp
View File

@ -168,7 +168,8 @@ public:
loc_param.linear_solver_reduction_ = 1e-2;
}
loc_param.linear_solver_print_json_definition_ = false;
domain_linsolvers_.emplace_back(model_.ebosSimulator(), loc_param);
const bool force_serial = true;
domain_linsolvers_.emplace_back(model_.ebosSimulator(), loc_param, force_serial);
}
assert(int(domains_.size()) == num_domains);

3
opm/simulators/linalg/ISTLSolverEbos.cpp
View File

@ -124,11 +124,12 @@ void FlexibleSolverInfo<Matrix,Vector,Comm>::create(const Matrix& matrix,
const PropertyTree& prm,
std::size_t pressureIndex,
std::function<Vector()> trueFunc,
const bool forceSerial,
[[maybe_unused]] Comm& comm)
{
// Write sizes of linear systems on all ranks to debug log.
{
if (!forceSerial) {
#if HAVE_MPI
auto basic_comm = comm.communicator();
#else

10
opm/simulators/linalg/ISTLSolverEbos.hpp
View File

@ -103,6 +103,7 @@ struct FlexibleSolverInfo
const PropertyTree& prm,
std::size_t pressureIndex,
std::function<Vector()> trueFunc,
const bool forceSerial,
Comm& comm);
std::unique_ptr<AbstractSolverType> solver_;
@ -177,13 +178,14 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
/// \param[in] simulator The opm-models simulator object
/// \param[in] parameters Explicit parameters for solver setup, do not
/// read them from command line parameters.
ISTLSolverEbos(const Simulator& simulator, const FlowLinearSolverParameters& parameters)
ISTLSolverEbos(const Simulator& simulator, const FlowLinearSolverParameters& parameters, bool forceSerial = false)
: simulator_(simulator),
iterations_( 0 ),
calls_( 0 ),
converged_(false),
matrix_(nullptr),
parameters_(parameters)
parameters_(parameters),
forceSerial_(forceSerial)
{
initialize();
}
@ -377,7 +379,7 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
bool isParallel() const {
#if HAVE_MPI
return comm_->communicator().size() > 1;
return !forceSerial_ && comm_->communicator().size() > 1;
#else
return false;
#endif
@ -403,6 +405,7 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
prm_,
pressureIndex,
trueFunc,
forceSerial_,
*comm_);
}
else
@ -502,6 +505,7 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
bool useWellConn_;
FlowLinearSolverParameters parameters_;
bool forceSerial_ = false;
PropertyTree prm_;
std::shared_ptr< CommunicationType > comm_;