We follow
[1]}, which relies on reformulating the
log-determinant in terms of the Jacobian expressed via an orthonormal
basis of the tangent space. The Jacobian determinant term is similar
to the Euclidean case suggesting that the high dimensional case can
reuse techniques from [2]}.
tab:rezenderuntime shows our running times are comparable to
[1]}.
