Hi,I am using Scalapack example sample_pssyev_call.f to solve the eigenproblem for big matrices. So far I have used the code successfully to solve matrix of size 23400x23400. The code fails when I try to solve a larger matrix of size 39600x39600, however. According to the message below the problem seems to be with LWORK. I am defining LWORK as LWORK=MAXN*MAXN*8 and it has been working perfectly fine for smaller matrices. I tried several things such as LWORK query, increasing the number of processors, changed the value of NB, etc, but none of those solved the problem. I thought that maybe something was going wrong during definition of the input matrix A (for example elements in the term diag([...can take on very small values that would make PSSYEV to crash later on maybe because of some division by zero or a maximum interation reached, etc). Then I gave a try using my own data but the error is exactly the same, the problem with LWORK persists. As you can see in the message below, the code executes completely, however, the computation is incorrect. I am only printing the first 10 largest eigenvalues. I am not exhausting the memory, I have plenty (24 GB) for this type of calculation. Does anyone one has an idea on why the code is crashing..??.

Hi Julie,THANKS a lot for your response. Sorry for the delay but I was traveling. As I said in my previous post I did try the workspace query but it didn't solve the problem, I got the same problem. You can find the code below. I would like to know if you can reproduce the same problem if you use the code below. The code posted uses a matrix of 23400x23400 and for that case the workspace query gave me LWORK=137053791. When I use small matrices there is not problems. For example the following are the results for N=50, 5000 and 10000:

Now, for very large matrices (that still fit my memory of course) the code crashes after a few hours running and got the following error: mpirun noticed that process rank 0 with PID 8457 on node arxt27 exited on signal 9 (Killed). In all the experiments I always kept the parameter MAXMAXN with a value larger than the workspace. Do you have any idea where is the problem ?. THANKS a lot for your help.JOVI