def test_rjb(self):
mesh = Mesh(numpy.array([-118.]), numpy.array([33])) # 1 point
surf18 = MultiSurface.from_csv(cd / 'msurface18.csv') # 2 planes
surf19 = MultiSurface.from_csv(cd / 'msurface19.csv') # 2 planes
surf20 = MultiSurface.from_csv(cd / 'msurface20.csv') # 1 plane
rjb18 = surf18.get_joyner_boore_distance(mesh)[0]
rjb19 = surf19.get_joyner_boore_distance(mesh)[0]
rjb20 = surf20.get_joyner_boore_distance(mesh)[0]
aac([rjb18, rjb19, rjb20], [85.676294, 89.225542, 92.937021])
surfa = MultiSurface(surf18.surfaces + surf19.surfaces)
surfb = MultiSurface(surf19.surfaces + surf20.surfaces)
rjba = surfa.get_joyner_boore_distance(mesh)[0]
rjbb = surfb.get_joyner_boore_distance(mesh)[0]
aac([rjba, rjbb], [85.676294, 89.225542])
def test_joyner_boore_distance_mesh1D(self):
surf = MultiSurface(self.surfaces_mesh1D)
numpy.testing.assert_equal(surf.get_joyner_boore_distance(self.mesh1D),
numpy.array([-1., 2., 2.]))
def test_joyner_boore_distance_mesh2D(self):
surf = MultiSurface(self.surfaces_mesh2D)
numpy.testing.assert_equal(
surf.get_joyner_boore_distance(self.mesh2D), numpy.array([[-1.0, 2.0, 2.0], [4.0, 4.0, 5.0]])
)
def test_joyner_boore_distance_mesh1D(self):
surf = MultiSurface(self.surfaces_mesh1D)
numpy.testing.assert_equal(surf.get_joyner_boore_distance(self.mesh1D),
numpy.array([-1., 2., 2.]))