OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-01 13:29:08 -06:00
Code Issues Packages Projects Releases Wiki Activity

Enforce row-major ordering of the pressure Jacobian

Browse Source 
Otherwise, we will use the transposed Jacobian as a coefficient
matrix in the 'linsolver_'.  This is wrong when solving compressible
problems for which the Jacobian matrix is not symmetric.
This commit is contained in:
Bård Skaflestad 2013-05-14 11:53:09 +02:00
parent 3f0f0363c3
commit 55165bd68d
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
ImpesTPFAAD.hpp
View File

@ -278,7 +278,7 @@ namespace Opm {
assemble(dt, state, well_state);
const int nc = grid_.number_of_cells;
M matr = cell_residual_.derivative()[0];
Eigen::SparseMatrix<double, Eigen::RowMajor> matr = cell_residual_.derivative()[0];
#if HACK_INCOMPRESSIBLE_GRAVITY
matr.coeffRef(0, 0) *= 2;