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-ma...@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
