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

Allow updateStoredTransmissibilities() before createMatrix_().

Browse Source 
A bug was triggered by a case placing MULTFLT at the start of the SCHEDULE.
Then the required structures would be empty and the simulator bombed. With
the change, we increase robustness and address this.
This commit is contained in:
Atgeirr Flø Rasmussen 2022-11-30 11:36:36 +01:00
parent 6cd200a254
commit 684332bb3f
1 changed files with 11 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
opm/models/discretization/common/tpfalinearizer.hh
View File

@ -310,6 +310,11 @@ private:
// Construct the BCRS matrix for the Jacobian of the residual function // Construct the BCRS matrix for the Jacobian of the residual function
void createMatrix_() void createMatrix_()
{ {
if (!neighborInfo_.empty()) {
// It is ok to call this function multiple times, but it
// should not do anything if already called.
return;
}
const auto& model = model_(); const auto& model = model_();
Stencil stencil(gridView_(), model_().dofMapper()); Stencil stencil(gridView_(), model_().dofMapper());
@ -505,6 +510,12 @@ private:
void updateStoredTransmissibilities() void updateStoredTransmissibilities()
{ {
if (neighborInfo_.empty()) {
// This function was called before createMatrix_() was called.
// We call initFirstIteration_(), not createMatrix_(), because
// that will also initialize the residual consistently.
initFirstIteration_();
}
unsigned numCells = model_().numTotalDof(); unsigned numCells = model_().numTotalDof();
#ifdef _OPENMP #ifdef _OPENMP
#pragma omp parallel for #pragma omp parallel for