Compute the Bures-Wasserstein Exponential

This function computes the Riemannian exponential map using the Bures-Wasserstein metric for symmetric positive-definite matrices. The map operates by solving a Lyapunov equation and then constructing the exponential.

bures_wasserstein_exp(sigma, v)

Arguments

sigma

A symmetric positive-definite matrix of class dppMatrix, representing the reference point.

v

A symmetric matrix of class dspMatrix, representing the tangent vector to be mapped.

Value

A symmetric positive-definite matrix of class dppMatrix, representing the point on the manifold.