[petsc-users] Problem with interpolating projected field of quadratic order

[Sander Arens]

Hello all,

Recently I've been trying to run a fem problem using DMPlex, PetscFE, etc.
using a mesh created with gmsh. I used DMPlexReverseCell to get rid of the
negative jacobian determinants. However, I noticed another problem: when I
projected the coordinates on a quadratic fem-field and interpolated them
they didn't match the interpolated values of the original coordinate field
anymore.

I attached a simple test to reproduce this. With linear interpolation I see
no problem, it's only when I start using quadratic interpolation.

I have no idea why this is giving problems. Am I missing something simple
here or is this a bug?

Thanks,
Sander
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160618/11db7a91/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.c
Type: text/x-csrc
Size: 9033 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160618/11db7a91/attachment-0001.c>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: makefile
Type: application/octet-stream
Size: 2749 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160618/11db7a91/attachment-0001.obj>