Purpose

LA_PPSVX computes the solution to a linear system of equations $AX = B$, where $A$ is real symmetric or complex Hermitian and, in either case, positive definite, and where $X$ and $B$ are rectangular matrices or vectors. $A$ is stored in packed format.
LA_PPSVX can also optionally equilibrate the system if $A$ is poorly scaled, estimate the condition number of (the equilibrated) $A$, and compute error bounds.