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.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611