The do loop
DO 10 J = JL, ILO+JA+NB-IOFF-1, -NB
JB = MIN( JA+IHI-J, NB )
in pcgehdrv makes JB such that the following call to
CALL PCLARFT( 'Forward', 'Columnwise', IHI-K, JB, A, I+1, J,
$ DESCA, TAU, WORK( IPT ), WORK( IPW ) )
adresses TAU at position N, even though TAU only consist of N-1
initiated/computed values from PCGEHRD.
This is a possible bug - anyone who knows who wrote this code?