[petsc-users] DMCreateRestriction with DMDA

[zakaryah .]

Hi everyone,

I am wondering about how to use PETSc to create a restriction on a
structured mesh. The documentation for DMCreateRestriction says:

For DMDA objects this only works for "uniform refinement", that is the
> refined mesh was obtained DMRefine() or the coarse mesh was obtained by
> DMCoarsen().
>

If the fine DMDA is dmf, shouldn't I be able to do the following?

DMCoarsen(dmf,comm,&dmc);
> DMCreateRestriction(dmc,dmf,&restriction);
>

When I try this, I get the error

> [0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type
> [0]PETSC ERROR: DM type da does not implement DMCreateRestriction
>

It seems to me it shouldn't be too hard to get the restriction for a
uniform coarsening, so I'm sure I'm missing something. Thanks in advance
for your help!

Cheers, Zak
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200515/1eb0d9ea/attachment.html>