[petsc-dev] DMCreateSubDM error
Matthew Knepley
knepley at gmail.com
Thu Dec 7 19:01:24 CST 2017
On Tue, Dec 5, 2017 at 7:47 PM, Blaise A Bourdin <bourdin at lsu.edu> wrote:
> Hi,
>
> I am building torture test for DMPlex with somewhat complicated dof
> layout. In some situations, DMCreateSubDM fails.
>
> DMUAS has 3 fields corresponding to
> U: corresponding to vector-valued P1 or P2 Lagrange FE, Sigma,
> vector-valued cell centered, and Alpha, same element type as U but scalar.
> For Tri / Tet, P1 or P2, all is well, but with Quads / Hex and polynomial
> order 2, when trying to extract the DM for the (U,alpha) pair of fields, I
> get a "Nonconforming object sizes” error.
>
> Try EXOTest12 -i TwoQuads.exo -order 2 for instance
>
> I am attaching an code illustrating this behavior and three example meshes.
>
I have incorporated your test as
src/dm/impls/plex/examples/tests/ex26
in
knepley/feature-subdm-test
It runs now. Feel free to include other tests at the bottom of the source.
Thanks,
Matt
> Regards,
> Blaise
>
>
>
> --
> Department of Mathematics and Center for Computation & Technology
> Louisiana State University, Baton Rouge, LA 70803, USA
> Tel. +1 (225) 578 1612 <(225)%20578-1612>, Fax +1 (225) 578 4276
> <(225)%20578-4276> http://www.math.lsu.edu/~bourdin
>
>
>
>
>
>
>
>
--
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.caam.rice.edu/~mk51/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20171207/9b6fe5c2/attachment.html>
More information about the petsc-dev
mailing list