Iterate: factorisation: changed from QR to LU