Completion of LU decomposition related code