OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-15 02:43:25 -06:00
Code Issues Packages Projects Releases Wiki Activity

- made heuristic scaling of trueimpes

Browse Source
This commit is contained in:
hnil 2022-10-21 12:51:05 +02:00
parent cf6423bfce
commit 8ede0e67b2
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/simulators/linalg/getQuasiImpesWeights.hpp
View File

@ -122,7 +122,11 @@ namespace Amg
VectorBlockType bweights; VectorBlockType bweights;
MatrixBlockType block_transpose = Details::transposeDenseMatrix(block); MatrixBlockType block_transpose = Details::transposeDenseMatrix(block);
block_transpose.solve(bweights, rhs); block_transpose.solve(bweights, rhs);
bweights /= 1000.0; // given normal densities this scales weights to about 1. double abs_max = *std::max_element(
bweights.begin(), bweights.end(), [](double a, double b) { return std::fabs(a) < std::fabs(b); });
// probably a scaling which could give approximately total compressibility would be better
bweights /= std::fabs(abs_max); // given normal densities this scales weights to about 1.
weights[index] = bweights; weights[index] = bweights;
++index; ++index;
} }