Help needed: Rum time error from lapack function calls.

Open discussion regarding features, bugs, issues, vendors, etc.

Help needed: Rum time error from lapack function calls.

Postby vg3104 » Mon Sep 13, 2010 5:53 pm

Hi,
Can anyone please provide some advice on how to deal with a run-time error from a lapack function (used to calculate eigen values). The codes can be compiled but I got the following run-time error:
MKL ERROR: Parameter 5 was incorrect on entry to DSYEV

Some background information:
1. I'm trying to use dsyev_ function for calculating eigenvalue/vectors by doing the following:
dsyev_("V","U", &n, &mat[0][0], &lda, &vec[0], work, &lwork, &info);
2. Environment is Linux and this is my LIB_DIRS:
LIB_DIRS = -L/xenv/ICC/X/11.0_081/mkl/lib/em64t /xenv/ICC/X/11.0_081/mkl/lib/em64t/libmkl_solver_ilp 64_sequential.a -Wl,--start-group /xenv/ICC/X/11.0_081/mkl/lib/em64t/libmkl_intel_ilp64.a /xenv/ICC/X /11.0_081/mkl/lib/em64t/libmkl_sequential.a /xenv/ICC/X/11.0_081/mkl/lib/em64t/libmkl_core.a -Wl,--end-group -lpthread -llapack
I have been struggling with this error for a long time, can someone please let me know what might be causing this?
Any information will be highly appreciated.
Thanks,
vg3104
 
Posts: 1
Joined: Mon Sep 13, 2010 5:48 pm

Re: Help needed: Rum time error from lapack function calls.

Postby jteg68 » Sun Sep 26, 2010 12:12 pm

Hi,

the run time error tells you what is wrong, parameter 5 is incorrect.
Check the description in http://www.netlib.org/lapack/double/dsyev.f or some other documentation.

In this case lda was not large enough.
If A is NxN then lda needs to be at least N. If you stored A as a contigous 1-dim array, then lda=N.

/Jörgen
jteg68
 
Posts: 9
Joined: Wed Sep 02, 2009 8:18 am


Return to User Discussion

Who is online

Users browsing this forum: No registered users and 1 guest

cron