def test_rx(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
rx18 = surf18.get_rx_distance(mesh)[0]
rx19 = surf19.get_rx_distance(mesh)[0]
rx20 = surf20.get_rx_distance(mesh)[0]
aac([rx18, rx19, rx20], [51.610675, 54.441119, -60.205692])
surfa = MultiSurface(surf18.surfaces + surf19.surfaces)
surfb = MultiSurface(surf19.surfaces + surf20.surfaces)
rxa = surfa.get_rx_distance(mesh)[0]
rxb = surfb.get_rx_distance(mesh)[0]
aac([rxa, rxb], [53.034889, -56.064366])
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])