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

[Jed Brown]

1. Check for memory leaks using Valgrind.

2. Be sure to run --with-debugging=1 (the default) when trying to find the
error.

3. Send the full error message and the relevant bit of code.

On Fri, Jul 27, 2012 at 3:17 PM, Ronald M. Caplan <caplanr at predsci.com>wrote:

> 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/a6f74854/attachment.html>