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

Implemented case det(jacobian)=0 in Newton method.

Browse Source
This commit is contained in:
Xavier Raynaud 2012-03-06 09:09:37 +01:00
parent 02296478b1
commit 801f975b5b
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/polymer/TransportModelPolymer.cpp
View File

@ -492,6 +492,10 @@ namespace Opm
direction[0] *= -1.0; direction[0] *= -1.0;
direction[1] *= -1.0; direction[1] *= -1.0;
} }
if ((std::abs(direction[0]) + std::abs(direction[0])) == 0) {
direction[0] = end_point[0]-x[0];
direction[1] = end_point[1]-x[1];
}
bool t0_exists = true; bool t0_exists = true;
double t0; double t0;
if (direction[0] > 0) { if (direction[0] > 0) {
@ -613,7 +617,7 @@ namespace Opm
// a specified piecewise linear curve. At each iteration, we use a robust 1d solver. // a specified piecewise linear curve. At each iteration, we use a robust 1d solver.
void TransportModelPolymer::solveSingleCellNewton(int cell) void TransportModelPolymer::solveSingleCellNewton(int cell)
{ {
const int max_iters_falsi = 20; const int max_iters_falsi = 200;
const double tol = 1e-7; const double tol = 1e-7;
// the tolerance for 1d solver is set as a function of the residual // the tolerance for 1d solver is set as a function of the residual
// The tolerance falsi_tol is improved by (reduc_factor_falsi_tol * "previous residual") at each step // The tolerance falsi_tol is improved by (reduc_factor_falsi_tol * "previous residual") at each step