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

Send all phases for surface volume and saturation.

Browse Source 
Previously, we only did send the information for the
first phase which led to computing wrong values in
a simulation. With this patch we now send all the data
for all the phases.
This commit is contained in:
Markus Blatt
2015-05-05 19:14:27 +02:00
parent db06fa48f5
commit a31aa46744
1 changed files with 10 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
opm/autodiff/RedistributeDataHandles.hpp
View File

@@ -73,8 +73,10 @@ public:
buffer.write(sendState_.rv()[e.index()]); buffer.write(sendState_.rv()[e.index()]);
buffer.write(sendState_.pressure()[e.index()]); buffer.write(sendState_.pressure()[e.index()]);
buffer.write(sendState_.temperature()[e.index()]); buffer.write(sendState_.temperature()[e.index()]);
buffer.write(sendState_.saturation()[e.index()]); for ( int i=0; i<sendState_.numPhases(); ++i )
{
buffer.write(sendState_.saturation()[e.index()*sendState_.numPhases()+i]);
}
for ( int i=0; i<sendGrid_.numCellFaces(e.index()); ++i ) for ( int i=0; i<sendGrid_.numCellFaces(e.index()); ++i )
{ {
buffer.write(sendState_.facepressure()[sendGrid_.cellFace(e.index(), i)]); buffer.write(sendState_.facepressure()[sendGrid_.cellFace(e.index(), i)]);
@@ -95,7 +97,7 @@ public:
for ( int i=0; i<recvState_.numPhases(); ++i ) for ( int i=0; i<recvState_.numPhases(); ++i )
{ {
buffer.read(val); buffer.read(val);
recvState_.surfacevol()[e.index()]=val; recvState_.surfacevol()[e.index()*sendState_.numPhases()+i]=val;
} }
buffer.read(val); buffer.read(val);
recvState_.gasoilratio()[e.index()]=val; recvState_.gasoilratio()[e.index()]=val;
@@ -105,9 +107,11 @@ public:
recvState_.pressure()[e.index()]=val; recvState_.pressure()[e.index()]=val;
buffer.read(val); buffer.read(val);
recvState_.temperature()[e.index()]=val; recvState_.temperature()[e.index()]=val;
buffer.read(val); for ( int i=0; i<recvState_.numPhases(); ++i )
recvState_.saturation()[e.index()]=val; {
buffer.read(val);
recvState_.saturation()[e.index()*sendState_.numPhases()+i]=val;
}
for ( int i=0; i<recvGrid_.numCellFaces(e.index()); ++i ) for ( int i=0; i<recvGrid_.numCellFaces(e.index()); ++i )
{ {
buffer.read(val); buffer.read(val);