Function lu 

pub fn lu<T: Numeric<T>>(
    a: &Sprs<T>,
    s: &mut Symb,
    tol: T,
) -> Result<Nmrc<T>, Error>

(L,U,Pinv) = lu(A, (Q lnz unz)). lnz and unz can be guess

See: sqr(...)