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)