LAPACK Archives

[Lapack] does LOPT need update in dsyevd.f

Dear Yao Tao,
Indeed, those extra lines are not needed.
We just removed them.
The fix will be included in the next LAPACK release.
Thank you

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.

Yao Tao

GCG Systems & Technology Lab
Building 28, ZhongGuanCun SoftPark
No.8 Dong Bei Wang West Road
Haidian District Beijing P.R.China 100193
Tel: 86-10-82454864
Email: yaotao@Domain.Removed
Lapack mailing list

-------------- next part --------------
An HTML attachment was scrubbed...

<Prev in Thread] Current Thread [Next in Thread>

For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or