On Thu, Aug 25, 2022 at 7:12 PM Mike Michell <mi.mike1...@gmail.com> wrote:
> Hi, this is a duplication of > https://lists.mcs.anl.gov/pipermail/petsc-users/2022-August/046746.html > for in-depth question. > > I wrote a short code as attached to test interpolation between two DMPlex > objects. Goal is to map solution field defined on cell-centroid(DMPlex-1) > into vertex(DMPlex-2) field or vice versa. > > Basically, DMCreateInterpolation() fails in the example code and it was > not allowed to get useful error messages. The DMPlex is created by loading > a 2D square box in Gmsh. Can I get some comments on that? > Sure, I am looking at it. Thanks, Matt > Thanks, > -- 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://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>