[petsc-users] Element to local dof map using dmplex

[Morten Nobel-Jørgensen]

Dear PETSc developers / Matt

Thanks for your suggestions regarding our use of dmplex in a FEM context. However, Matt's advise on  using the PetscFE is not sufficient for our needs (our end goal is a topology optimization framework - not just FEM) and we must honestly admit that we do not see how we can use the MATIS and the MatSetValuesClosure or DMPlexMatSetClosure to solve our current issues as Stefano has suggested.

We have therefore created a more representative, yet heavily oversimplified, code example that demonstrates our problem. That is, the dof handling is only correct on a single process and goes wrong on np>1.

We hope very much that you can help us to overcome our problem.

Thank you for an excellent toolkit
Morten and Niels
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161012/b8778d17/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dmplex_test.tar
Type: application/x-tar
Size: 54272 bytes
Desc: dmplex_test.tar
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161012/b8778d17/attachment-0001.tar>