Page 1 of 1

ZHBGST rwork ref to unitialized elements

PostPosted: Wed May 06, 2009 7:33 am
by LawrenceMulholland
The DO 240 and DO 470 loops of the code have a lower limit of
I2 + KA. However I find that the smallest index element of RWORK
that is set is J2-M meaning that the lower limit should be such that
J-KA-M = J2-M => J=J2+KA.

I have substituted J2+KA for I2+KA in the two loops and find that the routine
passes all our tests.

Is it possible for someone to confirm this.

Re: ZHBGST rwork ref to unitialized elements

PostPosted: Sun May 10, 2009 2:29 pm
by Julien Langou
Hello Lawrence,
thanks for the patch,
I have corrected our svn repos with it (rev 656). Both CHBGST and ZHBGST.
See: http://www.netlib.org/lapack/Errata/.
Julien Langou.