Re: [petsc-dev] DMCreateSubDM error

[Matthew Knepley]

On Tue, Dec 5, 2017 at 7:47 PM, Blaise A Bourdin  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/ 

[petsc-dev] DMCreateSubDM error

[Blaise A Bourdin]

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.

Regards,
Blaise



--
Department of Mathematics and Center for Computation & Technology
Louisiana State University, Baton Rouge, LA 70803, USA
Tel. +1 (225) 578 1612, Fax  +1 (225) 578 4276 http://www.math.lsu.edu/~bourdin









EXOTest12.c
Description: EXOTest12.c


OneHex8.exo
Description: OneHex8.exo


twobytwo.exo
Description: twobytwo.exo


TwoQuads.exo
Description: TwoQuads.exo