debugging LU factorization code