Purpose

LA_PPSV computes the solution to a linear system of equations $A\,X = 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. The Cholesky decomposition is used to factor $A$ as

$$A = U^H U \mbox{ if {\bf UPLO} = 'U', or } A = L\:L^H \mbox{ if {\bf UPLO} = 'L'}$$

where $U$ is an upper triangular matrix and $L$ is a lower triangular matrix ($L=U^H$). The factored form of $A$ is then used to solve the above system.