[petsc-users] Segfault in DMMeshSetChart using Fortran

[Chris Eldred]

I am getting a segfault in DMMeshSetChart and I was wondering if someone
could shed a little light on it.

PETSC (v3.3 p2) was compiled with gcc/gfortran/g++ (4.4.3) using the
options: --with-clanguage=C++ --with-sieve --download-f-blas-lapack
--download-parmetis --download-hdf5 --download-boost --download-metis.

Running under gdb, the following error was obtained:
Program received signal SIGSEGV, Segmentation fault.
0x083352aa in ALE::Obj<ALE::IFSieve<int, ALE::malloc_allocator<int> >,
ALE::malloc_allocator<ALE::IFSieve<int, ALE::malloc_allocator<int> > >
>::operator-> (
    this=0x1c)
    at /home/celdred/Desktop/NTM/petsc-3.3-p2/include/sieve/ALE_mem.hh:675
675         X*   operator->() const {return objPtr;};

The relevant code snippet is:

call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
CHKERRQ(ierr)
call MPI_Comm_rank(PETSC_COMM_WORLD,rank,ierr)
CHKERRQ(ierr)
call DMMeshCreate(PETSC_COMM_WORLD,model_mesh,ierr)
CHKERRQ(ierr)

ncells = nx*ny
nvertices = nx*ny
nedges = 2*nx*ny

start_index = 0
end_index = ncells + nvertices + nedges - 1
call DMMeshSetChart(model_mesh, start_index, end_index , ierr)
CHKERRQ(ierr)

Any ideas?

-- 
Chris Eldred
DOE Computational Science Graduate Fellow
Graduate Student, Atmospheric Science, Colorado State University
B.S. Applied Computational Physics, Carnegie Mellon University, 2009
chris.eldred at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120725/8caee05b/attachment.html>