This function solves a system of linear equation with symmetric positive definite (SPD) matrix along with multiple right-hand side vector
\[Ax=B\]
where \(A\) is a dense SPD triangular matrix of size \(m \times m\), \(x\) is a vector need to be computed, and \(B\) is input vector. The maximum matrix size supported in FPGA is templated by NMAX.