spd_learn.functional.transport_logchol#

spd_learn.functional.transport_logchol(X, A, B)[source]#

Parallel transport for log-Cholesky metric.

The parallel transport of matrices \(\mathbf{X}\) in tangent space from an initial SPD/HPD matrix \(\mathbf{A}\) to a final SPD/HPD matrix \(\mathbf{B}\) for log-Cholesky metric is given in Proposition 7 of [1].

Warning: this function must be applied to matrices \(\mathbf{X}\) already projected in tangent space with a logarithmic map at \(\mathbf{A}\), not to SPD/HPD matrices in manifold.

Parameters:
  • X (ndarray, shape (..., n, n)) – Symmetric/Hermitian matrices in tangent space at A.

  • A (ndarray, shape (n, n)) – Initial SPD/HPD matrix.

  • B (ndarray, shape (n, n)) – Final SPD/HPD matrix.

Returns:

X_new – Matrices in tangent space transported from A to B.

Return type:

ndarray, shape (…, n, n)

Notes

Added in version 0.10.

Changed in version 0.11: Correct formula for HPD matrices.

Changed in version 0.12: Add support for NumPy and PyTorch.

See also

transport

References

[1]

Riemannian geometry of symmetric positive definite matrices via Cholesky decomposition Z. Lin. SIAM J Matrix Anal Appl, 2019, 40(4), pp. 1353-1370.