On Fri, 1 Jul 2011, Truman Ellis wrote:
> Won't FEMContext only give my the FEType of my current element?
The FEType is constant for any given variable. Even in the case of p
refinement the FEType merely encodes the minimum-p case's type, and
reinit functions then examine the p refinement level
> FEMContext::elem->neighbor(FEMContext::side) is the neighboring
> element.
So this will give me an Elem * pointer, which looks like the geometry of
the element, but I can't figure out how to get the appropriate FEType to
build a FEBase. The Elem::type() function only returns the inherited
ElemTy
On Thu, 30 Jun 2011, Truman Ellis wrote:
> How do I access neighbor information using
> the DiffSystem framework,
FEMContext::elem->neighbor(FEMContext::side) is the neighboring
element.
> specifically to calculate jumps and averages for calculating fluxes?
Currently that's only possible by do
I am trying to implement discontinuous Galerkin to solve the Laplace
equation using DiffSystem. Example 21 solves the same problem via the
older equation_system method. How do I access neighbor information using
the DiffSystem framework, specifically to calculate jumps and averages
for calculating