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