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

Initialize the communication object in initialize().

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2024-01-17 12:39:52 +01:00
parent 63208c38f6
commit a52bcf51b1
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -272,6 +272,13 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
f.interiorCellNum_ = interiorCellNum_;
}
#if HAVE_MPI
if (isParallel()) {
const std::size_t size = simulator_.vanguard().grid().leafGridView().size(0);
detail::copyParValues(parallelInformation_, size, *comm_);
}
#endif
// Print parameters to PRT/DBG logs.
if (on_io_rank && parameters_[activeSolverNum_].linear_solver_print_json_definition_) {
std::ostringstream os;
@ -308,12 +315,6 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
void initPrepare(const Matrix& M, Vector& b)
{
const bool firstcall = (matrix_ == nullptr);
#if HAVE_MPI
if (firstcall && isParallel()) {
const std::size_t size = M.N();
detail::copyParValues(parallelInformation_, size, *comm_);
}
#endif
// update matrix entries for solvers.
if (firstcall) {