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

test_props: uncommented tests.

Browse Source
This commit is contained in:
Steinar Foss 2019-09-13 15:38:00 +02:00
parent 41a05a7a98
commit c826e24607
3 changed files with 19 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
python/cxx/eclipse_grid.cpp
View File

@ -52,8 +52,8 @@ void python::common::export_EclipseGrid(py::module& module) {
.def_property_readonly( "cartesianSize", &getCartesianSize )
.def( "globalIndex", &getGlobalIndex )
.def( "getIJK", &getIJK )
.def( "_cellVolume1G", &cellVolume1G)
.def( "_cellVolume3", &cellVolume3)
.def( "getCellVolume", &cellVolume1G)
.def( "getCellVolume", &cellVolume3)
;
}

2
python/cxx/eclipse_state.cpp
View File

@ -92,7 +92,7 @@ void python::common::export_EclipseState(py::module& module) {
py::class_< EclipseState >( module, "EclipseState" )
.def(py::init<const Deck&>())
.def_property_readonly( "title", &EclipseState::getTitle )
.def( "_props", &EclipseState::get3DProperties, ref_internal)
.def( "props", &EclipseState::get3DProperties, ref_internal)
.def( "grid", &EclipseState::getInputGrid, ref_internal)
.def( "_cfg", &EclipseState::cfg, ref_internal)
.def( "tables", &EclipseState::getTableManager, ref_internal)

31
python/tests/test_props.py
View File

@ -1,6 +1,10 @@
import unittest
import opm.io
from opm.io.parser import Parser
from opm.io.ecl_state import EclipseState
class TestProps(unittest.TestCase):
def assertClose(self, expected, observed, epsilon=1e-08):
@ -9,12 +13,11 @@ class TestProps(unittest.TestCase):
self.assertTrue(diff <= epsilon, msg=err_msg)
def setUp(self):
self.spe3 = opm.io.parse('tests/spe3/SPE3CASE1.DATA').state
"""
parser = Parser()
deck = parser.parse('tests/spe3/SPE3CASE1.DATA')
self.spe3 = EclipseState(deck)
self.props = self.spe3.props()
"""
"""
def test_contains(self):
p = self.props
self.assertTrue('PORO' in p)
@ -42,29 +45,27 @@ class TestProps(unittest.TestCase):
print('set(PERMX) = %s' % set(permx))
# 130mD, 40mD, 20mD, and 150mD, respectively, top to bottom
darcys = {0:md2si(130), 1:md2si(40), 2:md2si(20), 3:md2si(150)}
for i in range(grid.getNX()):
for j in range(grid.getNY()):
for k in range(grid.getNZ()):
for i in range(grid.nx):
for j in range(grid.ny):
for k in range(grid.nz):
g_idx = grid.globalIndex(i,j,k)
perm = permx[g_idx]
darcy = darcys[k]
self.assertClose(darcy, perm)
def test_volume(self):
e3dp = self.props
grid = self.spe3.grid()
for i in range(grid.getNX()):
for j in range(grid.getNY()):
for k in range(grid.getNZ()):
for i in range(grid.nx):
for j in range(grid.ny):
for k in range(grid.nz):
g_idx = grid.globalIndex(i,j,k)
exp = 293.3 * 293.3 * 30 # cubicfeet = 73 078.6084 cubic meter
exp *= (12*0.0254)**3 # cubic feet to cubic meter
if k == 0:
self.assertClose(exp, grid.getCellVolume(g_idx))
self.assertEqual(grid.getCellVolume(g_idx), grid.getCellVolume(None, i, j, k))
"""
self.assertEqual(grid.getCellVolume(g_idx), grid.getCellVolume(i, j, k))
if __name__ == "__main__":
unittest.main()