On Mon, Mar 27, 2023 at 12:30 AM Adrian Croucher <a.crouc...@auckland.ac.nz> wrote:
> hi > > I am testing my code against the latest PETSc release branch, having > built it mostly with PETSc 3.15 until now. > > An error is coming up when I read in a DMPlex, distribute it, add extra > cells etc. to it (to represent dual porosity subsurface flow) and then > re-distribute it for load balancing: > > [0]PETSC ERROR: Petsc has generated inconsistent data > [0]PETSC ERROR: Point SF contains 25 which is a cell > [0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting. > [0]PETSC ERROR: Petsc Release Version 3.18.5, unknown > [0]PETSC ERROR: > /home/acro018/software/geothermal/waiwera/code/petsc-3.19-debug/mesh_test > on a v3.19-debug named en-354401 by acro018 Mon Mar 27 14:28:07 2023 > [0]PETSC ERROR: Configure options --with-x --download-cmake > --download-hdf5 --download-zlib --download-netcdf --download-pnetcdf > --download-exodusii --download-triangle --download-ptscotch > --download-chaco --download-hypre > [0]PETSC ERROR: #1 DMPlexCheckPointSF() at > /home/acro018/software/PETSc/code/src/dm/impls/plex/plex.c:8875 > [0]PETSC ERROR: #2 DMPlexCreateNumbering_Plex() at > /home/acro018/software/PETSc/code/src/dm/impls/plex/plex.c:8147 > [0]PETSC ERROR: #3 DMPlexCreatePartitionerGraph_Native() at > /home/acro018/software/PETSc/code/src/dm/impls/plex/plexpartition.c:153 > [0]PETSC ERROR: #4 DMPlexCreatePartitionerGraph() at > /home/acro018/software/PETSc/code/src/dm/impls/plex/plexpartition.c:514 > [0]PETSC ERROR: #5 PetscPartitionerDMPlexPartition() at > /home/acro018/software/PETSc/code/src/dm/impls/plex/plexpartition.c:781 > [0]PETSC ERROR: #6 DMPlexDistribute() at > /home/acro018/software/PETSc/code/src/dm/impls/plex/plexdistribute.c:1820 > [0]PETSC ERROR: #7 ../src/mesh.F90:3618 > > The code still runs and appears to produce correct results. Has some new > checking been added which might be showing up a problem that was already > there? Any other clues as to how to approach it? > Yes, it is a new check. I do not expect cells that are shared until we create an overlap. We can just put in a flag to shut it off, since otherwise it is a great check for finding bugs. Can you show me how you add cells, so we can make a small example that tests this to go in PETSc? Thanks, Matt > - Adrian > > -- > Dr Adrian Croucher > Senior Research Fellow > Department of Engineering Science > University of Auckland, New Zealand > email: a.crouc...@auckland.ac.nz > tel: +64 (0)9 923 4611 > > -- 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 https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>