OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge pull request #2130 from bska/report-size-mismatch

Browse Source 
Report Element Counts for Solution Vector Size Mismatch
This commit is contained in:
Joakim Hove
2020-11-24 14:05:26 +01:00
committed by GitHub
parent b5012898d8 1c726b85c9
commit 52df226da9
1 changed files with 29 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

47
src/opm/output/eclipse/RestartIO.cpp
View File

@@ -50,18 +50,21 @@
#include <algorithm>
#include <cassert>
#include <cmath>
#include <cstddef>
#include <fstream>
#include <initializer_list>
#include <iomanip>
#include <iterator>
#include <sstream>
#include <stdexcept>
#include <string>
#include <unordered_set>
#include <utility>
#include <vector>
#include <fmt/format.h>
#include <fmt/chrono.h>
namespace Opm { namespace RestartIO {
namespace {
@@ -190,9 +193,13 @@ namespace {
const RestartValue& restart_value,
const EclipseGrid& grid)
{
for (const auto& elm: restart_value.solution)
if (elm.second.data.size() != grid.getNumActive())
throw std::runtime_error("Wrong size on solution vector: " + elm.first);
for (const auto& [name, vector] : restart_value.solution)
if (vector.data.size() != grid.getNumActive()) {
const auto msg = fmt::format("Incorrectly sized solution vector {}. "
"Expected {} elements, but got {}.", name,
grid.getNumActive(), vector.data.size());
throw std::runtime_error(msg);
}
if (es.getSimulationConfig().getThresholdPressure().size() > 0) {
// If the the THPRES option is active the restart_value should have a
@@ -621,28 +628,32 @@ namespace {
}
}
int numChar(const std::size_t num_reports)
{
return static_cast<int>(
1 + std::floor(std::log10(static_cast<double>(num_reports))));
}
void logRestartOutput(const int report_step,
const std::size_t num_reports,
const std::vector<int>& inteHD)
{
using namespace fmt::literals;
using Ix = ::Opm::RestartIO::Helpers::VectorItems::intehead;
std::ostringstream logmsg;
auto timepoint = std::tm{};
timepoint.tm_year = inteHD[Ix::YEAR] - 1900;
timepoint.tm_mon = inteHD[Ix::MONTH] - 1;
timepoint.tm_mday = inteHD[Ix::DAY];
logmsg << "Restart file written for report step: "
<< std::setw(numChar(num_reports)) << report_step << '/'
<< std::setw(0) << num_reports << ". Date: "
<< std::setw(4) << inteHD[Ix::YEAR] << '/'
<< std::setw(2) << std::setfill('0') << inteHD[Ix::MONTH] << '/'
<< std::setw(2) << std::setfill('0') << inteHD[Ix::DAY];
timepoint.tm_hour = inteHD[Ix::IHOURZ];
timepoint.tm_min = inteHD[Ix::IMINTS];
timepoint.tm_sec = inteHD[Ix::ISECND] / (1000 * 1000);
::Opm::OpmLog::info(logmsg.str());
const auto msg =
fmt::format("Restart file written for report step "
"{report_step:>{width}}/{num_reports}, "
"date = {timepoint:%d-%b-%Y %H:%M:%S}",
"width"_a = fmt::formatted_size("{}", num_reports),
"report_step"_a = report_step,
"num_reports"_a = num_reports,
"timepoint"_a = timepoint);
::Opm::OpmLog::info(msg);
}
} // Anonymous namespace