def test_check_orthogonality(): azi1 = 1 azi2 = 91 assert rotate.check_orthogonality(azi1, azi2) == "left-hand" assert rotate.check_orthogonality(azi2, azi1) == "right-hand" azi1 = 315 azi2 = 45 assert rotate.check_orthogonality(azi1, azi2) == "left-hand" assert rotate.check_orthogonality(azi2, azi1) == "right-hand" azi1 = 405 azi2 = 495 assert rotate.check_orthogonality(azi1, azi2) == "left-hand" assert rotate.check_orthogonality(azi2, azi1) == "right-hand" azi1 = 46 azi2 = 137 assert not rotate.check_orthogonality(azi1, azi2) assert not rotate.check_orthogonality(azi2, azi1) azi1 = 46 azi2 = 314 assert not rotate.check_orthogonality(azi1, azi2) assert not rotate.check_orthogonality(azi2, azi1)
def test_rotate_ne_12(): n = np.array([1.0, 0.0]) e = np.array([0.0, 1.0]) dnew1, dnew2 = rotate.rotate_ne_12(n, e, 30, 120) assert rotate.check_orthogonality(30, 120) == "left-hand" dnew1_true = np.array([np.sqrt(3) / 2.0, 0.5]) dnew2_true = np.array([-0.5, np.sqrt(3) / 2.0]) npt.assert_allclose(dnew1, dnew1_true) npt.assert_allclose(dnew2, dnew2_true)
def test_rotate_ne_12(): n = np.array([1.0, 0.0]) e = np.array([0.0, 1.0]) dnew1, dnew2 = rotate.rotate_ne_12(n, e, 30, 120) assert rotate.check_orthogonality(30, 120) == "left-hand" dnew1_true = np.array([np.sqrt(3)/2.0, 0.5]) dnew2_true = np.array([-0.5, np.sqrt(3)/2.0]) npt.assert_allclose(dnew1, dnew1_true) npt.assert_allclose(dnew2, dnew2_true)
def test_check_orthogonality_2(): azi1 = -180 azi2 = -90 assert rotate.check_orthogonality(azi1, azi2) == "left-hand" assert rotate.check_orthogonality(azi2, azi1) == "right-hand" azi1 = -315 azi2 = 135 assert rotate.check_orthogonality(azi1, azi2) == "left-hand" assert rotate.check_orthogonality(azi2, azi1) == "right-hand" azi1 = -181 azi2 = -90 assert not rotate.check_orthogonality(azi1, azi2) assert not rotate.check_orthogonality(azi2, azi1)