OilfieldToolsNet/opm-simulators
4
0 Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-10-26 07:30:08 -05:00
Code Issues Packages Projects Releases Wiki Activity

Add more output data to NLDD simulation summary

Browse Source
This commit is contained in:
jakobtorben committed 2025-03-31 10:58:29 +02:00
1 parent 91d3d5bb2e
commit c968518f7c
3 files changed
+57 -3
Download Patch File Download Diff File Expand all files Collapse all files

No files matched your search

opm/simulators/flow/BlackoilModelNldd.hpp
+10
View File
@@ -451,6 +451,13 @@ private:
++overlap_cells;
}
}
// Store data for summary output
local_reports_accumulated_.success.num_wells = num_wells;
local_reports_accumulated_.success.num_domains = num_domains;
local_reports_accumulated_.success.num_overlap_cells = overlap_cells;
local_reports_accumulated_.success.num_owned_cells = owned_cells;
// Gather data from all ranks
std::vector<int> all_owned(comm.size());
std::vector<int> all_overlap(comm.size());
@@ -1019,6 +1026,7 @@ private:
const Scalar acceptable_local_cnv_sum = 1.0;
if (mb_sum < acceptable_local_mb_sum && cnv_sum < acceptable_local_cnv_sum) {
local_report.converged = true;
local_report.accepted_unconverged_domains += 1;
logger.debug(fmt::format("Accepting solution in unconverged domain {} on rank {}.", domain.index, rank_));
logger.debug(fmt::format("Value of mb_sum: {} cnv_sum: {}", mb_sum, cnv_sum));
} else {
@@ -1032,9 +1040,11 @@ private:
}
}
if (local_report.converged) {
local_report.converged_domains += 1;
auto local_solution = Details::extractVector(solution, domain.cells);
Details::setGlobal(local_solution, domain.cells, locally_solved);
} else {
local_report.unconverged_domains += 1;
wellModel_.setPrimaryVarsDomain(domain.index, initial_local_well_primary_vars);
Details::setGlobal(initial_local_solution, domain.cells, solution);
model_.simulator().model().invalidateAndUpdateIntensiveQuantities(/*timeIdx=*/0, domain);
opm/simulators/timestepping/SimulatorReport.cpp
+30 -3
View File
@@ -35,7 +35,8 @@ namespace Opm
return SimulatorReportSingle{1.0, 2.0, 3.0, 4.0, 5.0, 6.0,
7.0, 8.0, 9.0, 10.0, 11.0, 12.0,
13, 14, 15, 16, 17, 18,
true, false, 19, 20.0, 21.0};
true, false, 19, 20.0, 21.0,
22, 23, 24, 25, 26, 27, 28};
}
bool SimulatorReportSingle::operator==(const SimulatorReportSingle& rhs) const
@@ -62,7 +63,14 @@ namespace Opm
this->well_group_control_changed == rhs.well_group_control_changed &&
this->exit_status == rhs.exit_status &&
this->global_time == rhs.global_time &&
this->timestep_length == rhs.timestep_length;
this->timestep_length == rhs.timestep_length &&
this->num_domains == rhs.num_domains &&
this->num_wells == rhs.num_wells &&
this->num_overlap_cells == rhs.num_overlap_cells &&
this->num_owned_cells == rhs.num_owned_cells &&
this->converged_domains == rhs.converged_domains &&
this->unconverged_domains == rhs.unconverged_domains &&
this->accepted_unconverged_domains == rhs.accepted_unconverged_domains;
}
void SimulatorReportSingle::operator+=(const SimulatorReportSingle& sr)
@@ -88,6 +96,10 @@ namespace Opm
}
max_linear_iterations = std::max(max_linear_iterations, sr.total_linear_iterations);
converged_domains += sr.converged_domains;
unconverged_domains += sr.unconverged_domains;
accepted_unconverged_domains += sr.accepted_unconverged_domains;
// It makes no sense adding time points. Therefore, do not
// overwrite the value of global_time which gets set in
// NonlinearSolver.hpp by the line:
@@ -228,11 +240,15 @@ namespace Opm
void SimulatorReportSingle::reportNLDD(std::ostream& os, const SimulatorReportSingle* failureReport) const
{
os << fmt::format("Owned + overlap cells: {:7}\n", num_owned_cells + num_overlap_cells);
os << fmt::format("Number of wells: {:7}\n", num_wells);
os << fmt::format("Number of domains: {:7}\n", num_domains);
os << fmt::format("-------------------------------------------------------\n");
double t = total_time + (failureReport ? failureReport->total_time : 0.0);
os << fmt::format("Total time: {:9.2f} s\n", t);
t = pre_post_time + (failureReport ? failureReport->pre_post_time : 0.0);
os << fmt::format(" Pre/post time: {:7.2f} s\n", t);
os << fmt::format(" Pre/post/wait time: {:7.2f} s\n", t);
t = solver_time + (failureReport ? failureReport->solver_time : 0.0);
os << fmt::format(" Solver time: {:9.2f} s", t);
@@ -314,6 +330,17 @@ namespace Opm
100.0*failureReport->total_linear_iterations/noZero(n));
}
os << std::endl;
os << fmt::format("-------------------------------------------------------\n");
n = converged_domains + (failureReport ? failureReport->converged_domains : 0);
os << fmt::format("Converged domain solves: {:7}", n);
os << std::endl;
n = accepted_unconverged_domains + (failureReport ? failureReport->accepted_unconverged_domains : 0);
os << fmt::format(" Accepted with relaxed tol: {:7}", n);
os << std::endl;
n = unconverged_domains + (failureReport ? failureReport->unconverged_domains : 0);
os << fmt::format("Unconverged domain solves: {:7}", n);
os << std::endl << std::endl;
}
SimulatorReport SimulatorReport::serializationTestObject()
opm/simulators/timestepping/SimulatorReport.hpp
+17
View File
@@ -59,6 +59,16 @@ namespace Opm
double global_time = 0.0;
double timestep_length = 0.0;
// NLDD specific data
int num_domains = 0;
int num_wells = 0;
int num_overlap_cells = 0;
int num_owned_cells = 0;
int converged_domains = 0;
int unconverged_domains = 0;
int accepted_unconverged_domains = 0;
static SimulatorReportSingle serializationTestObject();
bool operator==(const SimulatorReportSingle&) const;
@@ -95,6 +105,13 @@ namespace Opm
serializer(exit_status);
serializer(global_time);
serializer(timestep_length);
serializer(num_domains);
serializer(num_wells);
serializer(num_overlap_cells);
serializer(num_owned_cells);
serializer(converged_domains);
serializer(unconverged_domains);
serializer(accepted_unconverged_domains);
}
};