[petsc-dev] DMPlex redistribution

[Adrian Croucher]

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?

- Adrian

-- 
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.croucher at auckland.ac.nz
tel: +64 (0)9 923 4611