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

print warning when off-diagonal element is found in Schur complement instead of OPM_THROW.

Browse Source
This commit is contained in:
Robert Kloefkorn 2014-10-07 10:00:38 +02:00
parent ce996c2a6c
commit 62cefb3a3e
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/autodiff/NewtonIterationBlackoilCPR.cpp
View File

@ -241,7 +241,8 @@ namespace Opm
// std::cout << "++++++++++++++++++++++++++++++++++++++++++++\n" // std::cout << "++++++++++++++++++++++++++++++++++++++++++++\n"
// << D // << D
// << "++++++++++++++++++++++++++++++++++++++++++++\n" << std::endl; // << "++++++++++++++++++++++++++++++++++++++++++++\n" << std::endl;
OPM_THROW(std::logic_error, "Cannot do Schur complement with respect to non-diagonal block."); std::cerr << "WARNING (ignored): Cannot do Schur complement with respect to non-diagonal block." << std::endl;
//OPM_THROW(std::logic_error, "Cannot do Schur complement with respect to non-diagonal block.");
} }
V diag = D.diagonal(); V diag = D.diagonal();
Eigen::DiagonalMatrix<double, Eigen::Dynamic> invD = (1.0 / diag).matrix().asDiagonal(); Eigen::DiagonalMatrix<double, Eigen::Dynamic> invD = (1.0 / diag).matrix().asDiagonal();
@ -296,7 +297,8 @@ namespace Opm
// Find inv(D). // Find inv(D).
const M& D = equation.derivative()[n]; const M& D = equation.derivative()[n];
if (!isDiagonal(D)) { if (!isDiagonal(D)) {
OPM_THROW(std::logic_error, "Cannot do Schur complement with respect to non-diagonal block."); std::cerr << "WARNING (ignored): Cannot do Schur complement with respect to non-diagonal block." << std::endl;
//OPM_THROW(std::logic_error, "Cannot do Schur complement with respect to non-diagonal block.");
} }
V diag = D.diagonal(); V diag = D.diagonal();
Eigen::DiagonalMatrix<double, Eigen::Dynamic> invD = (1.0 / diag).matrix().asDiagonal(); Eigen::DiagonalMatrix<double, Eigen::Dynamic> invD = (1.0 / diag).matrix().asDiagonal();