granat wrote: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?
Julien Langou wrote:Hello Robert,
thanks for the patch. Is there an easy to check this by simply changing the HRD.dat file for example?
Users browsing this forum: No registered users and 5 guests