Matrix Q = I;
// check the expected properties of Q
assertTrue(Q.epsilonEquals(Q.getTranspose(),1e-6));
assertTrue(Q.epsilonEquals(Q.inverse(),1e-6));
// SimpleMatrix result = Q.mult(A.extractMatrix(w,width,w,width));
AStridedMatrix result = Multiplications.multiply(Q, A.subMatrix(w,width-w,w,width-w));
for( int i = 1; i < width-w; i++ ) {