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 #761 from atgeirr/fix-multflt-early-bug

Browse Source 
Allow updateStoredTransmissibilities() before createMatrix_().
This commit is contained in:
Atgeirr Flø Rasmussen 2022-12-01 08:14:05 +01:00 committed by GitHub
commit 9fabd9b49d
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
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_();
Stencil stencil(gridView_(), model_().dofMapper());
@ -505,6 +510,12 @@ private:
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();
#ifdef _OPENMP
#pragma omp parallel for