def test_geodesic_logeuclid_middle():
"""Test log euclidean geodesic when alpha = 0.5"""
A = 0.5*np.eye(3)
B = 2*np.eye(3)
Ctrue = 1*np.eye(3)
assert_array_almost_equal(geodesic_logeuclid(A,B,0.5),Ctrue)
def test_geodesic_logeuclid_1():
"""TTest log euclidean geodesic when alpha = 1"""
A = 0.5*np.eye(3)
B = 2*np.eye(3)
assert_array_almost_equal(geodesic_logeuclid(A,B,1),B)
