E.g., this can be used to solve constraint equations of the form J.I.Jt.X = B where J is a Jacobian (Jt its transpose), I is an Inverse mas matrix, and B the residual. In this case, I is symmetric positive definite, and therefore so is JIJt.