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