def test_distances_rotations(self):
for axis, angle in axis_angle_sequence():
s = random_direction()
R = rotation_from_axis_angle(axis, angle)
s2 = np.dot(R, s)
dist = geodesic_distance_on_sphere(s, s2)
# Note: this is == only if axis is orthogonal to s
assert dist <= angle
def test_distances_rotations(self):
for axis, angle in axis_angle_sequence():
s = random_direction()
R = rotation_from_axis_angle(axis, angle)
s2 = np.dot(R, s)
dist = geodesic_distance_on_sphere(s, s2)
# Note: this is == only if axis is orthogonal to s
assert dist <= angle
def any_distant_direction_test():
for s in directions_sequence():
z = any_distant_direction(s)
d = geodesic_distance_on_sphere(z, s)
assert d > np.pi / 6
def any_distant_direction_test():
for s in directions_sequence():
z = any_distant_direction(s)
d = geodesic_distance_on_sphere(z, s)
assert d > np.pi / 6
