LA_SPRFS / LA_HPRFS improve the computed solution to a system
of linear equations when the coefficient matrix is a real symmetric / complex
symmetric / complex Hermitian, indefinite and packed, and provides
error bounds and backward error estimates for the solution.
References: See [1] and [9,21,20].
-----------------------------------