Hello,

I am writing a package for manipulating positive definite matrices on their Riemannian manifold.

Those matrix may be real or complex and the code should be able to treat either case alike.

My code involves standard linear algebra operations such as, for example, Cholesky and eigenvalue-eigenvector decomposition.

For performance concern, should i take care to flag the matrices as Symmetric is they are real and as Hermitian if they are complex, or should I flag them as Hermitian in both cases?