<div dir="ltr">Hello Christian,<div>Thank you for your response! Having a non-conforming FEM is not an issue for me at this point, so I have been looking into the variableqkdgfem.hh implementation as per your suggestion, and think that this is exactly what I was looking for.</div><div><br></div><div>However, I had a follow-up question I was hoping to get feedback on: on line 24 of variableqkdgfem.hh, "Dune::QkDGLocalFiniteElement" is called, and my compiler complains that this does not name a template type. Should this be defined for me somewhere in dune/localfunctions/, where classes like "QkLocalFiniteElement" and "PkLocalFiniteElement" are defined? I don't see this anywhere in my code (I am using v2.5.1), so I was curious if I am missing some file or if I have to define the class "QkDGLocalFiniteElement" myself. I appreciate any feedback you or anyone else may have!</div><div><br></div><div>Thank you for your time,</div><div>Bryan</div><div><br></div><div class="gmail_extra"><div class="gmail_quote">On Thu, Mar 22, 2018 at 2:15 PM, Christian Engwer <span dir="ltr"><<a href="mailto:christian.engwer@uni-muenster.de" target="_blank">christian.engwer@uni-<wbr>muenster.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Bryan,<br>
<br>
there is a limited amount of support for p-refinement. There are<br>
different variable*fem.hh implementations of DG bases with variable<br>
polynomial degree. You can specify the degree via a non-standard<br>
method<br>
      void setOrder (const EntityType& e, unsigned int p)<br>
<br>
This should work quite straight forward. If you want to use conforming<br>
FEM things are a bit more complicated. You could implement something<br>
very similar to the current implementation, but how do you ensure<br>
conformity. This has to be done either by the finite-element-map by<br>
providing a specific finite-element which has a lower degree on a<br>
face, or by introducing appropriate constraints. Currently<br>
dune-localfunctions is not able to create such mixed degree elements<br>
and I don't think this is a viable apporach. So you want the<br>
constraints... PDELab has a couple of constraints implementations, but<br>
non that will work for this case.<br>
<br>
So first question is ... what exactly do you need?<br>
<br>
Best<br>
<span class="gmail-m_-5141631560027768196HOEnZb"><font color="#888888">Christian<br>
</font></span><div class="gmail-m_-5141631560027768196HOEnZb"><div class="gmail-m_-5141631560027768196h5"><br>
<br>
> Hello,<br>
> I am attempting to implement the ability to locally refine the polynomial<br>
> degree of the function space within my DUNE code, using PDELab. My first<br>
> thought was to use DUNE's "CompositeGridFunction" functionality, but I do<br>
> not believe this allows for modification of the grid space degree across<br>
> the domain.<br>
><br>
> To do this, would I instead have to code my own finite-element map whose<br>
> degree changes across the domain, and make a grid function space out of<br>
> that? Are there any tutorials available where such refinement is done? I<br>
> don't have any ideas of a starting point, so any thoughts or pointers would<br>
> be greatly appreciated!<br>
><br>
> Thank you,<br>
> Bryan<br>
</div></div></blockquote></div><br></div></div>