Re: [Libmesh-users] Discontinuous Galerkin in DiffSystem

2011-07-01 Thread Roy Stogner
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

Re: [Libmesh-users] Discontinuous Galerkin in DiffSystem

2011-07-01 Thread Truman Ellis
> 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

Re: [Libmesh-users] Discontinuous Galerkin in DiffSystem

2011-06-30 Thread Roy Stogner
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

[Libmesh-users] Discontinuous Galerkin in DiffSystem

2011-06-30 Thread Truman Ellis
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