def test_geodesic_log_euclidean():
old_dist = geodesic_logeuclid(m1.numpy_array, m2.numpy_array, 0.5)
m1.reset_fields()
m2.reset_fields()
new_dist = Geodesic.log_euclidean(m1, m2, 0.5)
return _get_state(old_dist, new_dist, "geodesic riemannian")
def compute(matrice_order, loop_number):
covmats = CovMats.random(1, matrice_order)
for i in range(0, loop_number):
mean = Mean.log_euclidean(covmats)
covmats.add_all(
[Geodesic.log_euclidean(covmat, mean) for covmat in covmats])
def compute(matrice_order, loop_number):
covmats = CovMats.random(1, matrice_order)
for i in range(0, loop_number):
mean = Mean.log_euclidean(covmats)
covmats.add_all([Geodesic.log_euclidean(covmat, mean) for covmat in covmats])
