**LA_SPSV** computes the solution to a linear system of
equations , where is a real or complex symmetric
matrix stored in packed format and and are
rectangular matrices or vectors.
A diagonal pivoting method is used to factor as

where (or ) is a product of permutation and unit upper (or lower) triangular matrices, and is a symmetric block diagonal matrix with and diagonal blocks. The factored form of is then used to solve the above system.

where (or ) is a product of permutation and unit upper (or lower) triangular matrices, and is a complex Hermitian block diagonal matrix with and diagonal blocks. The factored form of is then used to solve the above system.