Add LU decomposition for determinant and inversion support