1. 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)
2. In zhbgvd.f (chbgvd.f)
Line 312: LRWMIN = N
Line 374: Call ZHBGST(....,RWORK(INDWRK), IINFO)
INDWRK = N+1 and ZHBGST requires RWORK(N)
Therefore either
i. line 312 should be LWMIN = 2*N; or,
ii. line 374 should use RWORK(INDE)