SimpleMatrix Q = I.minus(U.mult(U.transpose()).scale(qr.getGamma()));
// check the expected properties of Q
assertTrue(Q.isIdentical(Q.transpose(),1e-6));
assertTrue(Q.isIdentical(Q.invert(),1e-6));
SimpleMatrix result = Q.mult(A.extractMatrix(w,width,w,width));
for( int i = 1; i < width-w; i++ ) {
assertEquals(0,result.get(i,0),1e-5);