def test_join_faults(self): grid = EclGrid.createRectangular((100, 100, 10), (1, 1, 1)) # Fault1 Fault4 # | | # | | # | | # | ------- Fault2 | # | | # | | # # -------- Fault3 # fault1 = Fault(grid, "Fault1") fault2 = Fault(grid, "Fault2") fault3 = Fault(grid, "Fault3") fault4 = Fault(grid, "Fault4") fault1.addRecord(1, 1, 10, grid.getNY() - 1, 0, 0, "X") fault2.addRecord(5, 10, 15, 15, 0, 0, "Y") fault3.addRecord(5, 10, 5, 5, 0, 0, "Y") fault4.addRecord(20, 20, 10, grid.getNY() - 1, 0, 0, "X") rays = fault1.getEndRays(0) self.assertEqual(rays[0], [(2, 10), (0, -1)]) self.assertEqual(rays[1], [(2, 100), (0, 1)]) extra = Fault.joinFaults(fault1, fault3, 0) self.assertEqual(extra, [(2, 10), (2, 6), (5, 6)])
def test_join_faults(self): grid = EclGrid.createRectangular( (100,100,10) , (1,1,1)) # Fault1 Fault4 # | | # | | # | | # | ------- Fault2 | # | | # | | # # -------- Fault3 # fault1 = Fault(grid , "Fault1") fault2 = Fault(grid , "Fault2") fault3 = Fault(grid , "Fault3") fault4 = Fault(grid , "Fault4") fault1.addRecord(1 , 1 , 10 , grid.getNY() - 1 , 0 , 0 , "X") fault2.addRecord(5 , 10 , 15 , 15 , 0 , 0 , "Y") fault3.addRecord(5 , 10 , 5 , 5 , 0 , 0 , "Y") fault4.addRecord(20 , 20 , 10 , grid.getNY() - 1 , 0 , 0 , "X") rays = fault1.getEndRays(0) self.assertEqual( rays[0] , [(2,10) , (0,-1)]) self.assertEqual( rays[1] , [(2,100) , (0,1)]) extra = Fault.joinFaults( fault1 , fault3 , 0) self.assertEqual( extra , [(2,10) , (2,6) , (5,6)] )
def join_faults(self, fault1, fault2): if not fault1.intersectsFault(fault2, self.getK()): layer = self.getGeoLayer() try: layer.addIJBarrier(Fault.joinFaults(fault1, fault2, self.getK())) except ValueError: err = 'Failed to join faults %s and %s' names = (fault1.getName(), fault2.getName()) print(err % names) raise ValueError(err % names)
def join_faults(self, fault1, fault2): if not fault1.intersectsFault(fault2, self.getK()): layer = self.getGeoLayer() try: layer.addIJBarrier(Fault.joinFaults(fault1, fault2, self.getK())) except ValueError: err = 'Failed to join faults %s and %s' names = (fault1.getName(), fault2.getName()) print(err % names) raise ValueError(err % names)