ScaLAPACK and already distributed LHS matrix

Post here if you have a question about LAPACK or ScaLAPACK algorithm or data format

ScaLAPACK and already distributed LHS matrix

Postby listboss » Tue Jan 04, 2011 2:55 pm


In the current version of my code I build the LHS matrix 'A' on different processors of our cluster
since it's too big to fit on a single machine. (A is a dense matrix)
The way I build the matrix on different processors is compatible with "two-dimensional block cyclic decomposition"
and I have all the information about which blocks of matrix are on which machines.
So I was wondering if I can use ScaLAPACK to solve A*X=B equation and at the end just gather the solution onto master machine.

any input and/or hint is appreciated.

thanks in advance
Posts: 1
Joined: Tue Jan 04, 2011 2:15 pm

Re: ScaLAPACK and already distributed LHS matrix

Postby Powergon » Sun Jan 27, 2013 2:23 pm

I am a new user of ScaLapack and have a similar question. I have a matrix that has to be computed over distributed nodes (each entry is obtained by calling some subprogram). But I cannot gather the entries back to form a single matrix ('cuz it's too large). How can I then call, say, PDSYEVR routine to compute its eigenvalues ? In an eigen-solver routine in Scalapck, the matrix variable "A" is the actual whole matrix or what ?

I really appreciate anyone that can give a hint.
Posts: 13
Joined: Sat Sep 22, 2012 12:42 pm

Return to Algorithm / Data

Who is online

Users browsing this forum: No registered users and 2 guests