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_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])
