Error when using magmaf_dsyevd subroutine

Open discussion for MAGMA

Error when using magmaf_dsyevd subroutine

Postby ssk87 » Mon Oct 17, 2011 1:46 pm

hey,

I tried to replace the existing dsyevd call with magmaf_dysevd call

ie)
I changed the following call from


call dsyevd('V','L',n, dbl_mb(k_v),n,
$ evals, WORK, LWORK,IWORK,LIWORK, info)


to
call magmaf_dsyevd('V','L',n, dbl_mb(k_v),n,
$ evals, WORK, LWORK,IWORK,LIWORK, info)

I ran the executable with dsyevd call and i can run all my input files without any issue
After changing my call to magmaf_dysevd I can run some the executable for some inputs
but for some inputs i am not getting an error.
The error is due to eigen values and eigen vector calculated. The eigen values and vector calculated for both the functions are not matching for some iterations
For dsyevd:
ex)
Eigen values -41.5371014014666
Eigen vectors -0.999978541295458
info = 0

____________________________________________
Eigen values -1.14087720918238
Eigen vectors 1.348890343184306E-016
info = 0
___________________________________________

For magmaf_dsyevd:

Eigen values -417.329141005850
Eigen vectors -1.575938597845697E+088
info = 0
__________________________________________
Eigen values 1.04850447237086
Eigen vectors 7.700221792289171E-017
info = 0
____________________________________________

Do you know how does it happens?
Because magmaf_dsyevd works for some examples and not working for some!!
How can I fix this issue?

Regards
Selva
ssk87
 
Posts: 4
Joined: Tue Sep 06, 2011 12:51 pm

Return to User discussion

Who is online

Users browsing this forum: Bing [Bot], Google [Bot] and 3 guests

cron