def setUp(self): # Funny setup to ensure that the 32 bit memory region is # exhausted. Then we do some (grid-based) pointer # dereferencing in the 64 bit memory region - this has been # problematic in the past. vecList = [] for i in range(12500): vec = DoubleVector() vec[81920] = 0 vecList.append(vec) self.grid1 = ecl.EclGrid(case) self.grid2 = ecl.EclGrid(case)
# FITNESS FOR A PARTICULAR PURPOSE. # # See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html> # for more details. import random import copy from ert.util.tvector import DoubleVector from ert.util.tvector import IntVector from ert.util.tvector import BoolVector from ert.util.stringlist import StringList import ert.util.stat as stat from ert.util.util_func import * v = DoubleVector(0) print v v[0] = 77.25 v[1] = 123.25 v[2] = 66.25 v[3] = 56.25 v[4] = 111.25 v[5] = 99.25 v[12] = 12 v.printf() v.clear() for i in range(100000): r = random.random() v.append(r)
zcorn = EclKW.read_grdecl(fileH, "ZCORN") coord = EclKW.read_grdecl(fileH, "COORD") actnum = EclKW.read_grdecl(fileH, "ACTNUM", ecl_type=ecl.ECL_INT_TYPE) grid = ecl.EclGrid.create(specgrid, zcorn, coord, actnum) return grid init_file = EclFile("data/eclipse/case/ECLIPSE.INIT") egrid_file = "data/eclipse/case/ECLIPSE.EGRID" grid_file = "data/eclipse/case/ECLIPSE.GRID" grdecl_file = "data/eclipse/case/include/example_grid_sim.GRDECL" grid = load_grdecl(grdecl_file) grid = load_grid(grid_file) grid = load_egrid(egrid_file) print "Thickness(10,11,12): %g" % grid.cell_dz(ijk=(10, 11, 12)) permx_column = DoubleVector(-999) grid.load_column(init_file.iget_named_kw("PERMX", 0), 5, 5, permx_column) permx_column.printf() print "top2 : %g depth(10,10,0) : %g " % (grid.top( 10, 10), grid.depth(ijk=(10, 10, 0))) print "bottom2 : %g depth(10,10,nz-1) : %g " % (grid.bottom( 10, 10), grid.depth(ijk=(10, 10, grid.nz - 1))) kw_list = init_file[1:7] print kw_list
fileH = open( grdecl_file , "r") specgrid = EclKW.read_grdecl( fileH , "SPECGRID" , ecl_type = ecl.ECL_INT_TYPE , strict = False) zcorn = EclKW.read_grdecl( fileH , "ZCORN" ) coord = EclKW.read_grdecl( fileH , "COORD" ) actnum = EclKW.read_grdecl( fileH , "ACTNUM" , ecl_type = ecl.ECL_INT_TYPE ) grid = ecl.EclGrid.create( specgrid , zcorn , coord , actnum ) return grid init_file = EclFile( "data/eclipse/case/ECLIPSE.INIT" ) egrid_file = "data/eclipse/case/ECLIPSE.EGRID" grid_file = "data/eclipse/case/ECLIPSE.GRID" grdecl_file = "data/eclipse/case/include/example_grid_sim.GRDECL" grid = load_grdecl( grdecl_file ) grid = load_grid( grid_file ) grid = load_egrid( egrid_file ) print "Thickness(10,11,12): %g" % grid.cell_dz( ijk=(10,11,12) ) permx_column = DoubleVector( -999 ) grid.load_column( init_file.iget_named_kw( "PERMX" , 0 ) , 5 , 5 , permx_column) permx_column.printf() print "top2 : %g depth(10,10,0) : %g " % (grid.top( 10, 10) , grid.depth( ijk=(10,10,0))) print "bottom2 : %g depth(10,10,nz-1) : %g " % (grid.bottom( 10 , 10 ) , grid.depth( ijk=(10,10,grid.nz - 1))) kw_list = init_file[1:7] print kw_list
# # See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html> # for more details. import random import copy from ert.util.tvector import DoubleVector from ert.util.tvector import IntVector from ert.util.tvector import BoolVector from ert.util.stringlist import StringList import ert.util.stat as stat from ert.util.util_func import * v = DoubleVector( 0 ) print v v[0] = 77.25 v[1] = 123.25 v[2] = 66.25 v[3] = 56.25 v[4] = 111.25 v[5] = 99.25 v[12] = 12 v.printf( ) v.clear() for i in range( 100000 ): r = random.random() v.append( r )