I am implementing QR decomposition (in Fortran) for a complex-valued matrix, using Classical Gram-Schmidt and Modified Gram-Schmidt (and Householder). I was expecting that the Classical Gram-Schmidt would yield slightly different results than the Modified, due to its numerical instability (round-off errors). Instead, my obtained R matrices are the same for both occasions. Does anyone know why is this happening?
I'm using double precision.