SYEVX Example in fortran

General CULA Dense (LAPACK & BLAS) support and troubleshooting. Use this forum if you are having a general problem or have encountered a bug.

SYEVX Example in fortran

Postby quinioue » Tue Jun 08, 2010 6:53 am

I am looking for an fortran example using SYEVX.
I have some problem with the function arguments, and I can't obtain eigenvalues.
Many thanks for response.
quinioue
CULA Academic
 
Posts: 9
Joined: Tue Feb 23, 2010 3:27 am

Re:SYEVX Example in fortran

Postby kyle » Tue Jun 08, 2010 7:02 am

Can you post a snippet of your current code or give some more details into what you are trying to do? Also, the SYEVX is the "expert" driver that gives more options. Be sure you actually need those additional options.
kyle
Administrator
 
Posts: 301
Joined: Fri Jun 12, 2009 7:47 pm

Re:SYEVX Example in fortran

Postby quinioue » Tue Jun 08, 2010 8:17 am

You could find in attached file the source of my test program. When I use syev function the eigenvalues are Ok. Perhaps, I make some errors in the SYEVX parameters, but I don't know where.

I am using cula premium 2.0 on linux 64 bit. I would like calculate the i lowest eigenvalues/eigenvectors of a NxN symetric matrice.

Many thanks for your help.

PS: I tested this function with a C program, and it works fine in double and single precision. Perhaps there is a bug in the fortran calling ?

[file name=syevx.txt size=3228]http://www.culatools.com/images/fbfiles/files/syevx.txt[/file]
quinioue
CULA Academic
 
Posts: 9
Joined: Tue Feb 23, 2010 3:27 am

Re:SYEVX Example in fortran

Postby kyle » Wed Jun 09, 2010 10:03 am

It seems there is a small bug in our FORTRAN interface for this function. This will be fixed in the next release; there is nothing wrong with your calling style. In the mean time, you can use the C/C++ interface if you'd like to test out some code.

Sorry for the inconvenience, CULA 2.0 should be available shortly with this fixed.
kyle
Administrator
 
Posts: 301
Joined: Fri Jun 12, 2009 7:47 pm

Re:SYEVX Example in fortran

Postby quinioue » Thu Jun 10, 2010 1:52 am

Many thanks for your response, I look forward the new release.
quinioue
CULA Academic
 
Posts: 9
Joined: Tue Feb 23, 2010 3:27 am

Re:SYEVX Example in fortran

Postby kyle » Thu Jun 10, 2010 9:19 am

CULA 2.0 should be out very shortly, we are just waiting on NVIDIA to release CUDA 3.1.

Here is the expected output we are now seeing:
Initializing CULA
Calling CULA_DSYEV
M = 5
W(1) = 0.71001198910815744
W(2) = 0.82403591500104556
W(3) = 6.5766620655512185
W(4) = 7.5254814352281558
W(5) = 7.9438085784221197
Shutting down CULA

The only change you'll need to your code is that N cannot be used as input and output. Use another variable, such as M, to denote the number of eigenvalues found.

STATUS = CULA_DSYEVX('V', 'I', 'U', N, A, N, VL, VU, IL,
$ IU, ABSTOL, M, W, Z, N, IFAIL)
kyle
Administrator
 
Posts: 301
Joined: Fri Jun 12, 2009 7:47 pm

Re: SYEVX Example in fortran

Postby ink » Tue Jun 29, 2010 8:02 am

SYEVX in CULA 2.0 preview does not work.
SYEV does.
ink
CULA Academic
 
Posts: 5
Joined: Wed May 19, 2010 8:14 am

Re: SYEVX Example in fortran

Postby john » Tue Jun 29, 2010 9:05 am

Both should be fixed and functional in CULA 2.0 Final which is available as of yesterday.

Please do note that SYEVX returns a value in the integer "m" and as such that variable should not be the same variable as the other integer arguments that are input only. Kyle's post (two above this one) details this further.
john
Administrator
 
Posts: 587
Joined: Thu Jul 23, 2009 2:31 pm

Re: SYEVX Example in fortran

Postby ink » Tue Jun 29, 2010 9:10 am

yes, I noticed CULA 2.0 has been released. I'd like to test but I can't login into download page anymore.
ink
CULA Academic
 
Posts: 5
Joined: Wed May 19, 2010 8:14 am

Re: SYEVX Example in fortran

Postby quinioue » Thu Jul 01, 2010 5:32 am

DSYEVX still do not work in the new release !
quinioue
CULA Academic
 
Posts: 9
Joined: Tue Feb 23, 2010 3:27 am

Re: SYEVX Example in fortran

Postby ink » Thu Jul 01, 2010 5:35 am

I've tested it and it works fine although not as quickly as I was hoping for.
ink
CULA Academic
 
Posts: 5
Joined: Wed May 19, 2010 8:14 am

Re: SYEVX Example in fortran

Postby quinioue » Thu Jul 01, 2010 5:44 am

Have you tested in fortran interface ?
quinioue
CULA Academic
 
Posts: 9
Joined: Tue Feb 23, 2010 3:27 am

Re: SYEVX Example in fortran

Postby ink » Thu Jul 01, 2010 6:01 am

yes
ink
CULA Academic
 
Posts: 5
Joined: Wed May 19, 2010 8:14 am

Re: SYEVX Example in fortran

Postby quinioue » Thu Jul 01, 2010 8:00 am

Sorry, I have tested by error with the old library, but now I need to upgrade the NVIDIA driver.
I'll test again later.
quinioue
CULA Academic
 
Posts: 9
Joined: Tue Feb 23, 2010 3:27 am


Return to CULA Dense Support

Who is online

Users browsing this forum: No registered users and 1 guest

cron