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

Bugfix: must always initialize gpress_omegaweighted_, even if no gravity.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-02-20 17:05:28 +01:00
parent 2aac2255b5
commit 22bf83dda7
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/core/pressure/IncompTpfa.cpp
View File

@ -47,13 +47,15 @@ namespace Opm
tpfa_htrans_compute(gg, permeability, &htrans_[0]); tpfa_htrans_compute(gg, permeability, &htrans_[0]);
if (gravity) { if (gravity) {
gpress_.resize(g.cell_facepos[ g.number_of_cells ], 0.0); gpress_.resize(g.cell_facepos[ g.number_of_cells ], 0.0);
gpress_omegaweighted_.resize(g.cell_facepos[ g.number_of_cells ], 0.0);
mim_ip_compute_gpress(gg->number_of_cells, gg->dimensions, gravity, mim_ip_compute_gpress(gg->number_of_cells, gg->dimensions, gravity,
gg->cell_facepos, gg->cell_faces, gg->cell_facepos, gg->cell_faces,
gg->face_centroids, gg->cell_centroids, gg->face_centroids, gg->cell_centroids,
&gpress_[0]); &gpress_[0]);
} }
// gpress_omegaweighted_ is sent to assembler always, and it dislikes
// getting a zero pointer.
gpress_omegaweighted_.resize(g.cell_facepos[ g.number_of_cells ], 0.0);
h_ = ifs_tpfa_construct(gg); h_ = ifs_tpfa_construct(gg);
} }