[petsc-dev] DMPlex redistribution

Adrian Croucher a.croucher at auckland.ac.nz
Sun Mar 26 23:29:46 CDT 2023


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
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 
[0]PETSC ERROR: #2 DMPlexCreateNumbering_Plex() at 
[0]PETSC ERROR: #3 DMPlexCreatePartitionerGraph_Native() at 
[0]PETSC ERROR: #4 DMPlexCreatePartitionerGraph() at 
[0]PETSC ERROR: #5 PetscPartitionerDMPlexPartition() at 
[0]PETSC ERROR: #6 DMPlexDistribute() at 
[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

More information about the petsc-dev mailing list