Dear Yao Tao,
Indeed, those extra lines are not needed.
We just removed them.
The fix will be included in the next LAPACK release.
On Apr 23, 2012, at 8:28 PM, Tao Yao wrote:
i found a subtle difference between dsyevd.f and ssyevd.f when using lapack
dsyevd updates LOPT after calling dsytrd and dlacpy. but those codes are not
visible in ssyevd.f, cheevd.f and zheevd.f.
pls refer to line 329, 344 in dsyevd.f.
i wonder whether those codes are necessary, because dsyevd has precalculated
LOPT with at least 1+6*N+2*N**2.
if those codes must be there, why not for ssyevd?
thanks in advance.
GCG Systems & Technology Lab
Building 28, ZhongGuanCun SoftPark
No.8 Dong Bei Wang West Road
Haidian District Beijing P.R.China 100193
Lapack mailing list