[petsc-users] sieve-dev Segfault in DMMeshSetChart using Fortran

Matthew Knepley knepley at gmail.com
Wed Jul 25 12:35:36 CDT 2012


On Wed, Jul 25, 2012 at 10:31 AM, Chris Eldred <chris.eldred at gmail.com>wrote:

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

Please use DMComplex. It is mostly the same and MUCH easier to debug. I
will help with any conversion
you need.

  Thanks,

    Matt


> 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
>



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


More information about the petsc-users mailing list