[petsc-dev] DMCreateSubDM error

Tue, 05 Dec 2017 16:48:10 -0800 

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

Attachment: EXOTest12.c
Description: EXOTest12.c

Attachment: OneHex8.exo
Description: OneHex8.exo

Attachment: twobytwo.exo
Description: twobytwo.exo

Attachment: TwoQuads.exo
Description: TwoQuads.exo

Reply via email to