dsyev eigenvalue solver

Open forum for general discussions relating to PLASMA.

Re: dsyev eigenvalue solver

Postby nchilton » Mon Apr 16, 2012 8:36 pm

Right....I've got pretty limited experience with c/c++, so would it be easiest to install the netlib lapacke and use that? Can I use the netlib lapacke in conjunction with the mkl blas, or do I need to use the entire netlib lapack/lapacke/blas suite solely?

Nick
nchilton
 
Posts: 16
Joined: Wed Apr 04, 2012 9:04 pm

Re: dsyev eigenvalue solver

Postby mateo70 » Tue Apr 17, 2012 10:06 am

Hello Nick,

No problem to use lapacke from Netlib in conjunction of MKL BLAS. That's what the plasma installer is doing when you configure it with MKL. To be sure to use the one from netlib and not the functions included in MKL, you have to make sure the -llapacke is placed just before the list of -lmkl_... that you need.

Mathieu
mateo70
 
Posts: 92
Joined: Fri May 07, 2010 3:48 pm

Re: dsyev eigenvalue solver

Postby nchilton » Wed Apr 18, 2012 2:20 am

Mathieu,

Fantastic, compilation achieved!

Now to debug! First question - the PLASMA_dsyev routine needs integers for jobz and uplo, compared to character for lapack/mkl etc. What are the integer versions for 'N' or 'V' and 'U' or 'L', respectively?

Cheers!
Nick
nchilton
 
Posts: 16
Joined: Wed Apr 04, 2012 9:04 pm

Re: dsyev eigenvalue solver

Postby nchilton » Fri May 11, 2012 2:10 am

Anyone?
nchilton
 
Posts: 16
Joined: Wed Apr 04, 2012 9:04 pm

Re: dsyev eigenvalue solver

Postby mateo70 » Fri May 11, 2012 9:42 am

Hello,

sorry for the late answer, I havn't seen the post before.
You have all those infomations in PLASMA headers files in include directory.
In your case, it will be PlasmaNoVec/PlasmaVec and PlasmaUpper/PlasmaLower

Mathieu
mateo70
 
Posts: 92
Joined: Fri May 07, 2010 3:48 pm

Re: dsyev eigenvalue solver

Postby nchilton » Sat May 12, 2012 2:21 am

Mathieu, thank you very much, all sorted! PLASMA is certainly quick, I'm getting good speed-ups! It seems irrespective of my matrix size, DYSEV works optimally with 8 cores. Any ideas why?

Cheers!
Nick
nchilton
 
Posts: 16
Joined: Wed Apr 04, 2012 9:04 pm

Re: dsyev eigenvalue solver

Postby haidar » Sat May 12, 2012 10:17 am

nchilton wrote:Mathieu, thank you very much, all sorted! PLASMA is certainly quick, I'm getting good speed-ups! It seems irrespective of my matrix size,
DYSEV works optimally with 8 cores. Any ideas why?
Cheers!
Nick


Hi Nick,
could you specify what did you mean by irrespective?
It should scale up to 96 cores that was what we already test it. However your matrix need to be large enough.
what is your matrix size? and what kind of system are you running on?
Thanks
Azzam
haidar
 
Posts: 13
Joined: Tue Sep 07, 2010 12:01 pm

Re: dsyev eigenvalue solver

Postby nchilton » Tue May 15, 2012 2:34 am

Hi Azzam,

I was timing the diagonalization of various matrices ranging in size from 2000 to 10,000 square. All of the results showed essentially flat results (only very slight increases) past 8 cores. The machine has 2 Quad-core Intel Xeon L5520 @ 2.27 GHz with a total of 48 GB DDR3 @ 1066 MHz and now that I think about it, it might be due to hyper threading... Do you use MPI to access large numbers of cores?

Cheers,
Nick
nchilton
 
Posts: 16
Joined: Wed Apr 04, 2012 9:04 pm

Re: dsyev eigenvalue solver

Postby admin » Tue May 15, 2012 8:07 am

PLASMA only uses Pthreads.
No MPI.
If your machine is a dual quad-core, then you only have eight physical cores.
You will not get any speedup from using more threads than physical cores.
In fact I am surprised that you are not seeing a slowdown.
Jakub
admin
Site Admin
 
Posts: 79
Joined: Wed May 13, 2009 1:27 pm

Re: dsyev eigenvalue solver

Postby mateo70 » Tue May 15, 2012 9:57 am

Hello Nick,

Can you post the result of the following command on your system, please ?
pkg-config --libs hwloc

Tio check that your threads are correctly binded, and also what value to you give to PLASMA_Init ?
If you give something bigger than 8, your threads won't be binded to the cores and the results shouldn't be as good as with only 8 threads as mentionned by Jakub.

Mathieu
mateo70
 
Posts: 92
Joined: Fri May 07, 2010 3:48 pm

PreviousNext

Return to User discussion

Who is online

Users browsing this forum: Bing [Bot] and 1 guest

cron