Hi Kyle,

I have a quick question. Should I be expecting FE_Simplex of cubic order to work in an FESystem? In case the answer is "yes", I'm attaching a minimal example that I would expect to work, but it fails some assert in the library. I tried checking, but could not really understand why this is happening.

It's an out of range index checks that triggers. I believe that the code in FESystem that builds hanging node constraints does not know how to deal with triangles yet.

Would you mind opening a bug report on github for this, and attach the minimal test case? I can see that I fix this in the next week or so.

Thanks for reporting! Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/


--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/729433d5-f627-4f74-9be1-9414e19bf722%40colostate.edu.

Reply via email to