SYMINV: An Algorithm for the Inversion of a Positive Definite Matrix by the Cholesky Decomposition Terry Seaks Econometrica 1972 40(5), 961