Module SPOIR in SLATECSolve a positive definite symmetric system of linear equations. Iterative refinement is used to obtain an error estimate.
CALL SPOIR (A, LDA, N, V, ITASK, IND, WORK)| D2b1b | Solution of systems of linear equations, including matrix inverses, and LU and related factorizations (real symmetric positive definite matrices) |