[petsc-users] segfault in MatAssemblyEnd() when using large matrices on multi-core MAC OS-X

[Ronald M. Caplan]

Hello,

I am running a simple test code which takes a sparse AIJ matrix in PETSc
and multiplies it by a vector.

The matrix is defined as an AIJ MPI matrix.

When I run the program on a single core, it runs fine.

When I run it using MPI with multiple threads (I am on a 4-core, 8-thread
MAC) I can get the code to run correctly for matrices under a certain size
(2880 X 2880), but when the matrix is set to be larger, the code crashes
with a segfault and the error says it was in the MatAssemblyEnd().
Sometimes it works with -n 2, but typically it always crashes when using
multi-core.

Any ideas on what it could be?

Thanks,

Ron Caplan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120727/4906d396/attachment.html>