[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