There are some inaccuracies in the documentation for the matrix Z in the DLAQR? family of routines.
According to the documentation:
- DHSEQR takes Z of dimension (LDZ, N), LDZ >= N. It passes it to DLAQR0.
-- DLAQR0 takes Z of dimension (LDZ, IHI), LDZ >= IHIZ (IHIZ <= N). It passes it to DLAQR3, DLAQR5, DLAHQR.
--- DLAQR3 takes Z of dimension (LDZ, N), specification of LDZ is referred to the calling routine.
--- DLAQR5 takes Z of dimension (LDZ, IHI), LDZ >= N. However, DLAQR5 does not have IHI parameter. Also, DLAQR5 says that "...transformation is accumulated into Z(ILOZ:IHIZ,ILO:IHI)..." but there is no ILO parameter either.
--- DLAHQR takes Z of dimension (LDZ, N), LDZ >= N.

