Added LU decomposition to linear solver