mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
943d7fc2ce
Enable master and slaves to start at any report step. In the previous commits, only first report step was supported.
79 lines
2.9 KiB
C++
79 lines
2.9 KiB
C++
/*
|
|
Copyright 2024 Equinor AS
|
|
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#include <config.h>
|
|
|
|
#include <opm/simulators/flow/ReservoirCoupling.hpp>
|
|
#include <opm/simulators/utils/ParallelCommunication.hpp>
|
|
#include <opm/common/OpmLog/OpmLog.hpp>
|
|
|
|
#include <fmt/format.h>
|
|
|
|
namespace Opm {
|
|
namespace ReservoirCoupling {
|
|
|
|
void custom_error_handler_(MPI_Comm* comm, int* err, const std::string &msg)
|
|
{
|
|
// It can be useful to have a custom error handler for debugging purposes.
|
|
// If not, MPI will use the default error handler which aborts the program, and
|
|
// it can be difficult to determine exactly where the error occurred. With a custom
|
|
// error handler you can at least set a breakpoint here to get a backtrace.
|
|
int rank;
|
|
MPI_Comm_rank(*comm, &rank);
|
|
char err_string[MPI_MAX_ERROR_STRING];
|
|
int len;
|
|
MPI_Error_string(*err, err_string, &len);
|
|
const std::string error_msg = fmt::format(
|
|
"Reservoir coupling MPI error handler {} rank {}: {}", msg, rank, err_string);
|
|
// NOTE: The output to terminal vie stderr or stdout has been redirected to files for
|
|
// the slaves, see Main.cpp. So the following output will not be visible in the terminal
|
|
// if we are called from a slave process.
|
|
// std::cerr << error_msg << std::endl;
|
|
OpmLog::error(error_msg); // Output to log file, also shows the message on stdout for master.
|
|
MPI_Abort(*comm, *err);
|
|
}
|
|
|
|
void custom_error_handler_slave_(MPI_Comm* comm, int* err, ...)
|
|
{
|
|
custom_error_handler_(comm, err, "slave");
|
|
}
|
|
|
|
void custom_error_handler_master_(MPI_Comm* comm, int* err, ...)
|
|
{
|
|
custom_error_handler_(comm, err, "master");
|
|
}
|
|
|
|
void setErrhandler(MPI_Comm comm, bool is_master)
|
|
{
|
|
MPI_Errhandler errhandler;
|
|
// NOTE: Lambdas with captures cannot be used as C function pointers, also
|
|
// converting lambdas that use ellipsis "..." as last argument to a C function pointer
|
|
// is not currently possible, so we need to use static functions instead.
|
|
if (is_master) {
|
|
MPI_Comm_create_errhandler(custom_error_handler_master_, &errhandler);
|
|
}
|
|
else {
|
|
MPI_Comm_create_errhandler(custom_error_handler_slave_, &errhandler);
|
|
}
|
|
MPI_Comm_set_errhandler(comm, errhandler);
|
|
}
|
|
|
|
} // namespace ReservoirCoupling
|
|
} // namespace Opm
|