Setting Element Order Independently

[Shane Sawyer]

Hello everyone,

I had a question about the possibility of setting directly the polynomial order of elements. In my problem, I have a priori knowledge about the order of the elements in my (1d) mesh. I would like to be able to loop over the active cells and then assign the order directly, rather than using an estimator to perform the p-adaptivity. I have looked at step-27, but this didn't appear to provide a mechanism for per element modification. Is it possible to make such modifications directly to the elements?

Thank you in advance,

Shane 

[Bruno Turcksin]

Shane,

Yes, it is possible. There are two different functions depending on what you want to do. If you want to set the order during setup you can use this function It is being used in step-46. If you want to control the p-refinement yourself, you can use this function. The order of the finite element will be changed during the refinement.

Best,

Bruno

[Shane Sawyer]

Bruno,

Thank you for the quick reply. These look like exactly what I need!