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

Change the number of integration points in the initialization

Browse Source 
The number is changed from 100->2000 to increase accuracy.
This commit is contained in:
Tor Harald Sandve 2015-10-30 08:14:10 +01:00
parent 3dea26ece4
commit 7542aca827
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/core/simulator/initStateEquil_impl.hpp
View File

@ -366,9 +366,9 @@ namespace Opm
typedef Details::RK4IVP<ODE> WPress;
std::array<WPress,2> wpress = {
{
WPress(drho, up , p0, 100)
WPress(drho, up , p0, 2000)
,
WPress(drho, down, p0, 100)
WPress(drho, down, p0, 2000)
}
};
@ -426,9 +426,9 @@ namespace Opm
typedef Details::RK4IVP<ODE> OPress;
std::array<OPress,2> opress = {
{
OPress(drho, up , p0, 100)
OPress(drho, up , p0, 2000)
,
OPress(drho, down, p0, 100)
OPress(drho, down, p0, 2000)
}
};
@ -488,9 +488,9 @@ namespace Opm
typedef Details::RK4IVP<ODE> GPress;
std::array<GPress,2> gpress = {
{
GPress(drho, up , p0, 100)
GPress(drho, up , p0, 2000)
,
GPress(drho, down, p0, 100)
GPress(drho, down, p0, 2000)
}
};