def angles(v):
"""Print the angles from u to the x-axis and target respectively (in
degrees).
"""
return ('x-axis angle: %s, target angle: %s' %
(degrees_between_unsigned(v, x_axis),
degrees_between_unsigned(v, target)))
"""Print the angles from u to the x-axis and target respectively (in
degrees).
"""
return ('x-axis angle: %s, target angle: %s' %
(degrees_between_unsigned(v, x_axis),
degrees_between_unsigned(v, target)))
def prob(u, v):
"""Compute a success probability from sampled and target vectors."""
angle = math.acos(abs(numpy.dot(u, v)))
return math.exp(-PROB_CONCENTRATION * angle ** 2)
x1 = normalized(numpy.array([1, 0]))
x2 = normalized(numpy.array([3, 1]))
x3 = normalized(numpy.array([0, 1]))
print(degrees_between_unsigned(x1, x1), prob(x1, x1))
print(degrees_between_unsigned(x1, x2), prob(x1, x2))
print(degrees_between_unsigned(x1, x3), prob(x1, x3))
rand = random.Random(SEED)
target = normalized(UNNORMED_TARGET)
x_axis = numpy.zeros(DIMS)
x_axis[0] = 1.0
print('Target:', target)
print('Angle from x-axis to target:', degrees_between_unsigned(x_axis, target))
model = make_basic_bingham_wishart_model(DIMS, KAPPA, DOF, rand)
print(model.inv_scale())
