C++ lufp gives correct solution (i.e. Ax=b) for reduced system when Nr <= Nc