Should cubic FE_SimplexP work in an FESystem on master?

[Kyle Schwiebert]

Hello,

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.

Unfortunately, I wouldn't have any comment on which specific commit might have caused the issue (if it is an issue) as I just switched to my own build of deal.II and was previously using 9.5.1. I'm not sure if this new instalation could be the issue, so I'll describe what I did:

I compiled deal.II with

cmake -DCMAKE_INSTALL_PREFIX=/usr/local -DDEAL_II_WITH_MPI=ON -DCMAKE_EXPORT_COMPILE_COMMANDS=ON -DDEAL_II_WITH_TRILINOS=ON -DDEAL_II_WITH_PETSC=ON ..

followed by 

sudo make -j 6 install

I did not get any notable errors or warnings.

I'm compiling and running the example with

cmake ..

make example

./example

Thanks for any help you can offer.

Regards,
Kyle

[Wolfgang Bangerth]

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: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/

[Kyle Schwiebert]

Wolfgang,

Thank you very much for your willingness to look into the issue. I made a github issue. I couldn't upload a .cpp source file, so I included the example as a code block.

Regards,

Kyle

--
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 a topic in the Google Groups "deal.II User Group" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/dealii/QmcWWwdUnFI/unsubscribe.
To unsubscribe from this group and all its topics, send an email to dealii+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/729433d5-f627-4f74-9be1-9414e19bf722%40colostate.edu.