[petsc-users] seg fault with VecGetArray
Matthew Knepley
knepley at gmail.com
Sat May 12 05:16:04 CDT 2012
On Sat, May 12, 2012 at 4:33 AM, Mohammad Mirzadeh <mirzadeh at gmail.com>wrote:
> Hi guys,
>
> I'm having a really weird issue here! My code seg faults for certain
> problem size and after using gdb I have been able to pinpoint the problem
> to a VecGetArray call. Here's a series of things I have tried so far
>
> 1) -on_error_attach_debugger -----> unsuccessful; does not launch debugger
> 2) -start_in_debugger -------> unsuccessful; does not start debugger
>
Run with -log_summary. It will tell you what options the program got. Also,
are there errors relating to X? Send
all output to petsc-maint at mcs.anl.gov
> 3) attaching debugger myself -----> code runs in debugger and seg faults
> when calling VecGetArray
>
Is this a debug build? What dereference is causing the SEGV? Is the Vec a
valid object? It sounds like
it has been corrupted.
> 4) using ierr=VecGetArray;CHKERRQ(ierr) ------> PETSc does not produce
> error messages; the code simply seg faults and terminates
> 5) checking the values of ierr inside the debugger ---------> They are all
> 0 up untill the code terminates; I think this means petsc does not generate
> error?
> 6) checking for memory leak with valgrind -----------> All I get are leaks
> from OpenMPI and PetscInitialize and PetscFinalize; I think these are
> just routine and safe?
>
How can we say anything without the valgrind output?
Matt
>
> What else can I try to find the problem? Any recommendation is really
> appreciated!
>
> Thanks,
> Mohammad
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120512/c0fbe9c3/attachment.htm>
More information about the petsc-users
mailing list