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