OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#3383 Make transmissibility condenser behave well with zero internal degrees of freedom.

Browse Source
This commit is contained in:
Gaute Lindkvist 2018-09-19 13:41:49 +02:00
parent 477b7dd018
commit 38d3b666e6
1 changed files with 13 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
ApplicationCode/ReservoirDataModel/Completions/RigTransmissibilityCondenser.cpp
View File

@ -415,12 +415,21 @@ void RigTransmissibilityCondenser::calculateCondensedTransmissibilities()
int externalEquationCount = totalEquationCount - internalEquationCount; int externalEquationCount = totalEquationCount - internalEquationCount;
MatrixXd condensedSystem;
MatrixXd Tee = totalSystem.bottomRightCorner(externalEquationCount, externalEquationCount); MatrixXd Tee = totalSystem.bottomRightCorner(externalEquationCount, externalEquationCount);
MatrixXd Tei = totalSystem.bottomLeftCorner(externalEquationCount, internalEquationCount);
m_TiiInv = totalSystem.topLeftCorner(internalEquationCount, internalEquationCount).inverse();
m_Tie = totalSystem.topRightCorner(internalEquationCount, externalEquationCount);
MatrixXd condensedSystem = Tee - Tei * m_TiiInv * m_Tie; if (internalEquationCount == 0)
{
condensedSystem = Tee;
}
else
{
MatrixXd Tei = totalSystem.bottomLeftCorner(externalEquationCount, internalEquationCount);
m_TiiInv = totalSystem.topLeftCorner(internalEquationCount, internalEquationCount).inverse();
m_Tie = totalSystem.topRightCorner(internalEquationCount, externalEquationCount);
condensedSystem = Tee - Tei * m_TiiInv * m_Tie;
}
// std::cout << "Te = " << std::endl << condensedSystem << std::endl << std::endl; // std::cout << "Te = " << std::endl << condensedSystem << std::endl << std::endl;