changes to lu.cpp