[petsc-users] DMUninterpolate of periodic mesh

[Matteo Semplice]

Dear all

     I have updated a code of mine to Petsc3.22, when trying to 
uninterpolate a periodic mesh and get the following error

[0]PETSC ERROR: --------------------- Error Message 
--------------------------------------------------------------
[0]PETSC ERROR: No support for this operation for this object type
[0]PETSC ERROR: Missing local coordinate vector
[0]PETSC ERROR: WARNING! There are unused option(s) set! Could be the 
program crashed before usage or a spelling mistake, etc!
[0]PETSC ERROR:   Option left: name:-dm_view (no value) source: command line
[0]PETSC ERROR: See https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!faxBLw4JxiiCgvgXrNR1BSgniTuJMn0UxlQ7HDbFsTvh1NfjdD1Od-Ac3koxw3kaxHCmCf6MByHS5ty49s-1W9xW81DQIhYr8dp4Zw$  for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.22.0, unknown
[0]PETSC ERROR: ./testDM with 1 MPI process(es) and PETSC_ARCH  on 
signalkuppe by matteo Wed Feb  5 18:21:24 2025
[0]PETSC ERROR: Configure options: 
--prefix=/home/matteo/software/petscsaved/3.22-opt/ 
PETSC_DIR=/home/matteo/software/petsc --PETSC_ARCH=opt 
--with-debugging=0 --COPTFLAGS="-O3 -march=native -mtune=native -mavx2" 
--CXXOPTFLAGS="-O3 -march=native -mtune=native -mavx2" --FOPTFLAGS="-O3 
-march=native -mtune=native -mavx2" --with-strict-petscerrorcode 
--download-hdf5 --download-ml --with-metis --with-parmetis --with-gmsh 
--with-triangle --with-zlib --with-p4est-dir=~/software/p4est/local/
[0]PETSC ERROR: #1 DMLocalizeCoordinates() at 
/home/matteo/software/petsc/src/dm/interface/dmperiodicity.c:368
[0]PETSC ERROR: #2 DMPlexCopy_Internal() at 
/home/matteo/software/petsc/src/dm/impls/plex/plexcreate.c:37
[0]PETSC ERROR: #3 DMPlexUninterpolate() at 
/home/matteo/software/petsc/src/dm/impls/plex/plexinterpolate.c:1892
[0]PETSC ERROR: #4 main() at ../src/testDM.cpp:32
[0]PETSC ERROR: PETSc Option Table entries:
[0]PETSC ERROR: -dm_view (source: command line)
[0]PETSC ERROR: -dm_view_early (source: command line)
[0]PETSC ERROR: ----------------End of Error Message -------send entire 
error message to petsc-maint at mcs.anl.gov----------

Basically I am creating a quad mesh on [0,1]x[0,1] with the command

   PetscCall( DMPlexCreateBoxMesh(MPI_COMM_WORLD, dim, simplex, faces, 
lower, upper, periodicity, interpolate, 0, PETSC_TRUE, &dmMesh) );

and then calling

   PetscCall( DMPlexUninterpolate(dmMesh, &dmMeshUnint) );

raises the error. This seems independent of the values that I pass for 
the new parameters localizationHeight and sparseLocalize.

Do I need to change something else in addition to adding the new 
parameters in the DMPlexCreateBoxMesh call?

Thanks in advance

     Matteo