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

Use diagonal matrix block pointers for well source terms.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen
2023-01-03 10:48:19 +01:00
parent e705dc41c2
commit 2964bd409e
2 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
opm/simulators/wells/BlackoilWellModel.hpp
View File

@@ -273,7 +273,7 @@ namespace Opm {
// add source from wells to the reservoir matrix
void addReservoirSourceTerms(GlobalEqVector& residual,
SparseMatrixAdapter& jacobian) const;
std::vector<typename SparseMatrixAdapter::MatrixBlock*>& diagMatAddress) const;
// called at the beginning of a report step
void beginReportStep(const int time_step);

8
opm/simulators/wells/BlackoilWellModel_impl.hpp
View File

@@ -1250,7 +1250,8 @@ namespace Opm {
template <typename TypeTag>
void BlackoilWellModel<TypeTag>::
addReservoirSourceTerms(GlobalEqVector& residual, SparseMatrixAdapter& jacobian) const
addReservoirSourceTerms(GlobalEqVector& residual,
std::vector<typename SparseMatrixAdapter::MatrixBlock*>& diagMatAddress) const
{
// NB this loop may write multiple times to the same element
// if a cell is perforated by more than one well, so it should
@@ -1264,14 +1265,13 @@ namespace Opm {
for (unsigned perfIdx = 0; perfIdx < rates.size(); ++perfIdx) {
unsigned cellIdx = cells[perfIdx];
auto rate = rates[perfIdx];
// Scalar volume = ebosSimulator_.problem().volume(cellIdx,0);
rate *= -1.0;
VectorBlockType res(0.0);
using MatrixBlockType = Opm::MatrixBlock<Scalar, numEq, numEq >;
using MatrixBlockType = typename SparseMatrixAdapter::MatrixBlock;
MatrixBlockType bMat(0.0);
ebosSimulator_.model().linearizer().setResAndJacobi(res, bMat, rate);
residual[cellIdx] += res;
jacobian.addToBlock(cellIdx, cellIdx, bMat);
*diagMatAddress[cellIdx] += bMat;
}
}
}