OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Now the relative volume discrepancy takes std::fabs(), so max_element() will be ok.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2011-02-02 10:12:12 +01:00
parent 4f78892dc4
commit c6835bb9a7
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
dune/porsol/blackoil/BlackoilFluid.hpp
View File

@ -136,7 +136,7 @@ namespace Opm
totphasevol[cell] = state.total_phase_volume_; totphasevol[cell] = state.total_phase_volume_;
double pv = rock.porosity(cell)*grid.cellVolume(cell); double pv = rock.porosity(cell)*grid.cellVolume(cell);
voldiscr[cell] = (totphasevol[cell] - pv)/dt; voldiscr[cell] = (totphasevol[cell] - pv)/dt;
relvoldiscr[cell] = (totphasevol[cell] - pv)/pv; relvoldiscr[cell] = std::fabs(totphasevol[cell] - pv)/pv;
saturation[cell] = state.saturation_; saturation[cell] = state.saturation_;
rel_perm[cell] = state.relperm_; rel_perm[cell] = state.relperm_;
viscosity[cell] = state.viscosity_; viscosity[cell] = state.viscosity_;