def test_geodesic_riemann(): old_dist = geodesic_riemann(m1.numpy_array, m2.numpy_array, 0.5) m1.reset_fields() m2.reset_fields() new_dist = Geodesic.riemannian(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.riemannian(covmats) covmats.add_all([Geodesic.riemannian(covmat, mean) for covmat in covmats])