Use a conjugate gradient for a normal equation to solve the system A.x = b
[xsol[,iter]] = sp_cgne(A,b,tol[,maxit[,M[,x0]]])
the A real sparse matrix of the system to be solved.
the b real vector of the system to be solved.
the tolerance for the resolution of the equation system.
the maximum number of iteration of the solver.
the preconditionning matrix (an identity matrix by default). The preconditionning matrix must be symmetric.
an initial starting solution.
the solution vector of the system.
the number of iterations performed to solve the system.
Use a conjugate gradient for a normal equation to solve the system A.x = b