def test_add_edge(self):
createRegion = RegionDefinition( 1 )
with self.assertRaises(TypeError):
createRegion.addEdge("Poly1")
with self.assertRaises(TypeError):
createRegion.addEdge(112)
def test_create(self):
with self.assertRaises(TypeError):
regionDef = RegionDefinition( "GG" )
with self.assertRaises(ValueError):
regionDef = RegionDefinition.create( 1 , self.faults , ["Test" , 1])
regionDef = RegionDefinition.create( 1 , self.faults , ["DF28_LC" , self.poly_file1 , self.poly_file2 ])
def test_create(self):
with self.assertRaises(TypeError):
regionDef = RegionDefinition("GG")
with self.assertRaises(ValueError):
regionDef = RegionDefinition.create(1, self.faults, ["Test", 1])
regionDef = RegionDefinition.create(
1, self.faults, ["DF28_LC", self.poly_file1, self.poly_file2])
def test_create(self):
defRegion = RegionDefinition.create( 5 , self.faults , ["DF4_MC","DF15_MC","DF43_MC","DF25_LC","DF21_C","DF1_C",self.poly_file6, "DF26_MC","DF34_MCS","DF32_MC",self.poly_file5])
region_kw = EclKW.create( "REGIONS" , self.grid.getGlobalSize() , EclTypeEnum.ECL_INT_TYPE )
region_kw.assign( 0 )
for k in range(self.grid.getNZ()):
with self.assertRaises(NotImplementedError):
block_list = defRegion.findInternalBlocks( self.grid , self.fault_blocks[k] )
def test_add_edge(self):
createRegion = RegionDefinition(1)
with self.assertRaises(TypeError):
createRegion.addEdge("Poly1")
with self.assertRaises(TypeError):
createRegion.addEdge(112)
def test_create(self):
defRegion = RegionDefinition.create(5, self.faults, [
"DF4_MC", "DF15_MC", "DF43_MC", "DF25_LC", "DF21_C", "DF1_C",
self.poly_file6, "DF26_MC", "DF34_MCS", "DF32_MC", self.poly_file5
])
region_kw = EclKW.create("REGIONS", self.grid.getGlobalSize(),
EclTypeEnum.ECL_INT_TYPE)
region_kw.assign(0)
for k in range(self.grid.getNZ()):
with self.assertRaises(NotImplementedError):
block_list = defRegion.findInternalBlocks(
self.grid, self.fault_blocks[k])
def test_update(self):
reg1 = RegionDefinition( 1 )
reg2 = RegionDefinition( 2 )
self.assertTrue( not reg1.hasPolygon() )
reg1.addEdge( self.poly1 )
reg2.addEdge( self.poly2 )
self.assertTrue( reg1.hasPolygon() )
region_kw = EclKW.create( "REGIONS" , self.grid.getGlobalSize() , EclTypeEnum.ECL_INT_TYPE )
region_kw.assign( 0 )
block_list1 = []
block_list2 = []
with self.assertRaises(NotImplementedError):
block_list1 = reg1.findInternalBlocks(self.grid , reg1.splitFaultBlocks(self.grid , self.fault_blocks[0]))
with self.assertRaises(NotImplementedError):
block_list2 = reg2.findInternalBlocks(self.grid , reg2.splitFaultBlocks(self.grid , self.fault_blocks[0]))
if block_list1:
for block in block_list1:
region_id = reg1.getRegionID()
block.assignToRegion( region_id )
for g in block.getGlobalIndexList():
region_kw[g] = region_id
if block_list2:
for block in block_list2:
region_id = reg2.getRegionID()
block.assignToRegion( region_id )
for g in block.getGlobalIndexList():
region_kw[g] = region_id
for j in range(4):
for i in range(4):
g1 = i + j*self.grid.getNX()
g2 = i + 12 + (j + 12)* self.grid.getNX()
def test_update(self):
reg1 = RegionDefinition(1)
reg2 = RegionDefinition(2)
self.assertTrue(not reg1.hasPolygon())
reg1.addEdge(self.poly1)
reg2.addEdge(self.poly2)
self.assertTrue(reg1.hasPolygon())
region_kw = EclKW.create("REGIONS", self.grid.getGlobalSize(),
EclTypeEnum.ECL_INT_TYPE)
region_kw.assign(0)
block_list1 = []
block_list2 = []
with self.assertRaises(NotImplementedError):
block_list1 = reg1.findInternalBlocks(
self.grid,
reg1.splitFaultBlocks(self.grid, self.fault_blocks[0]))
with self.assertRaises(NotImplementedError):
block_list2 = reg2.findInternalBlocks(
self.grid,
reg2.splitFaultBlocks(self.grid, self.fault_blocks[0]))
if block_list1:
for block in block_list1:
region_id = reg1.getRegionID()
block.assignToRegion(region_id)
for g in block.getGlobalIndexList():
region_kw[g] = region_id
if block_list2:
for block in block_list2:
region_id = reg2.getRegionID()
block.assignToRegion(region_id)
for g in block.getGlobalIndexList():
region_kw[g] = region_id
for j in range(4):
for i in range(4):
g1 = i + j * self.grid.getNX()
g2 = i + 12 + (j + 12) * self.grid.getNX()
