def test_add_edge(self):
createRegion = RegionDefinition( 1 )
with self.assertRaises(TypeError):
createRegion.addEdge("Poly1")
with self.assertRaises(TypeError):
createRegion.addEdge(112)
def test_add_edge(self):
createRegion = RegionDefinition(1)
with self.assertRaises(TypeError):
createRegion.addEdge("Poly1")
with self.assertRaises(TypeError):
createRegion.addEdge(112)
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()
