Add an inverse function for symmetric positive definite matrices

This includes both unit tests and a reference manual entry.