Bug report 1/2 from Lawrence Mulholland (NAG)
See: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4970
Posted on Wed May 25th
From Lawrence
=============
In dsbgvd.f (ssbgvd.f)
with the case WANTZ = .False.
LWMIN = 2*N (line 282)
But, (line 140)
Call dsbgst(.....,work(INDWRK), iinfo)
where INDWRK = N+1 and dsbgst requires a workspace of 2*N
Therefore either
i. line 282 should be LWMIN = 3*N; or,
ii. line 140 should use WORK(INDE)
From Julien
===========
Let us go with fix (ii) then.
Note: I simply used "WORK", not "WORK(INDE)"