Performs an incomplete LU decomposition on a sparse matrix
[L[,U]] = sp_luinc(M[,alpha])
the real sparse matrix on which we will perform the LU decomposition
a coefficient used in the partial pivoting (default value: 1).
partial pivoting and Markowitz criterion:
|a[p][k]| >= alpha * max_i |a[i][k]|
the L sparse matrix from the LU decomposition
the U sparse matrix from the LU decomposition
Performs an incomplete LU decomposition on a sparse matrix.
At the end of the decomposition, the number of non zeros elements in the L and U matrixes is less or equal to the number of non zeros elements in the M matrix.