lusolve(a,b,ps,zero=0.0)
Instance Public methods
Solves a*x = b for x, using LU decomposition.
a is a matrix, b is a constant vector, x is the solution vector.
ps is the pivot, a vector which indicates the permutation of rows performed
during LU decomposition.