[petsc-users] DMUninterpolate of periodic mesh

[Matthew Knepley]

On Wed, Feb 5, 2025 at 12:27 PM Matteo Semplice via petsc-users <
petsc-users at mcs.anl.gov> wrote:

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

This looks like a bug. I will fix it.

Note that you can get what you want by passing interpolate = PETSC_FALSE

  Thanks,

     Matt


> [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
>
>
>

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://urldefense.us/v3/__https://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!YZMEed4K2IwyEQhBwJUypWI3aKs5Wmwawwo7BynU-0Vn_H-W4_qFcSC3h-JFkofKZsXlP63vBUj4owuXbY5T$  <https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!YZMEed4K2IwyEQhBwJUypWI3aKs5Wmwawwo7BynU-0Vn_H-W4_qFcSC3h-JFkofKZsXlP63vBUj4o7DdLz_A$ >
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20250205/6750b10f/attachment.html>