ZHBGST rwork ref to unitialized elements

Open discussion regarding features, bugs, issues, vendors, etc.

ZHBGST rwork ref to unitialized elements

Postby LawrenceMulholland » Wed May 06, 2009 7:33 am

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.
Posts: 3
Joined: Fri Mar 23, 2007 5:59 am
Location: NAG Ltd

Re: ZHBGST rwork ref to unitialized elements

Postby Julien Langou » Sun May 10, 2009 2:29 pm

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.
Julien Langou
Posts: 824
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA

Return to User Discussion

Who is online

Users browsing this forum: wclodius and 4 guests