def test_compute_transmissibility(self):
problem = TwoDimReservoir(self.inputs)
np.testing.assert_allclose(problem.compute_transmissibility(0, 0),
4000.0)
return
def compute_transmissibility(self, i, j):
if self.p[i] > self.p[j]:
kr = self.kr[i]
else:
kr = self.kr[j]
return kr * TwoDimReservoir.compute_transmissibility(self, i, j)
