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

simplify: use matrix-vector product in DRS for matrix

Browse Source
This commit is contained in:
Arne Morten Kvarving 2019-03-14 11:25:45 +01:00
parent e12b959def
commit f98be9d43b
1 changed files with 3 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/autodiff/BlackoilAmg.hpp
View File

@ -118,13 +118,9 @@ scaleMatrixDRS(const Operator& op, const Communication& comm,
for (auto j = (*i).begin(); j != endj; ++j) {
BlockVector bvec(0.0);
Block& block = *j;
for (std::size_t ii = 0; ii < Block::rows; ii++) {
for (std::size_t jj = 0; jj < Block::cols; jj++) {
// should introduce limmits which also change the weights
bvec[jj] += bw[ii]*block[ii][jj];
}
}
block[pressureIndex] = bvec;
BlockVector& bvec = block[pressureIndex];
// should introduce limits which also change the weights
block.mtv(bw, bvec);
}
}
}