U.set(0,0,1); // this is not explicity set and is assumed to be 1
Matrix I = Matrix.createIdentity(width-w);
// SimpleMatrix Q = I.minus(U.mult(U.transpose()).scale(qr.getGamma()));
Matrix temp1 = Multiplications.multiply(U, U.getTranspose());
temp1.scale(qr.getGamma());
I.sub(temp1);
Matrix Q = I;
// check the expected properties of Q