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

Merge pull request #4734 from atgeirr/add-linear-system-size-printout

Browse Source 
Add output of linear system sizes to DBG file.
This commit is contained in:
Atgeirr Flø Rasmussen 2023-07-25 10:51:29 +02:00 committed by GitHub
commit 0d2d8dfe21
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 20 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -31,6 +31,8 @@
#include <opm/simulators/linalg/ParallelIstlInformation.hpp>
#include <opm/simulators/utils/ParallelCommunication.hpp>
#include <fmt/format.h>
#if COMPILE_BDA_BRIDGE
#include <opm/simulators/linalg/bda/BdaBridge.hpp>
#include <opm/simulators/linalg/bda/WellContributions.hpp>
@ -123,6 +125,24 @@ void FlexibleSolverInfo<Matrix,Vector,Comm>::create(const Matrix& matrix,
[[maybe_unused]] Comm& comm)
{
// Write sizes of linear systems on all ranks to debug log.
{
auto basic_comm = comm.communicator();
std::ostringstream os;
os << "Linear system ";
if (basic_comm.size() > 1) {
os << fmt::format("on MPI rank: {} ", basic_comm.rank());
}
os << fmt::format("blocksize: {} size: {:7d} block nonzeroes: {:9d}",
Matrix::block_type::rows, matrix.N(), matrix.nonzeroes());
DeferredLogger local_logger;
local_logger.debug(os.str());
auto global_logger = gatherDeferredLogger(local_logger, basic_comm);
if (basic_comm.rank() == 0) {
global_logger.logMessages();
}
}
std::function<Vector()> weightsCalculator =
getWeightsCalculator<Vector>(prm, matrix, pressureIndex, trueFunc);