On 4/12/19 1:55 PM, luca.heltai wrote:
> Wolfgang, is that true also for mass matrices? I’d agree with you for
> stiffness matrices, but I’d surprised this worked ok for mass
> matrices as well.
I'm pretty sure. The theory goes like this: instead of computing the
matrix and rhs using the
Wolfgang, is that true also for mass matrices? I’d agree with you for stiffness
matrices, but I’d surprised this worked ok for mass matrices as well.
If so, I’ve always been over integrating in my life…
:)
L.
> On 12 Apr 2019, at 21:15, Wolfgang Bangerth wrote:
>
> On 4/12/19 8:41 AM,
On 4/12/19 8:41 AM, Robert Spartus wrote:
>
> That is some fascinating information! It seems like step-44, for
> instance, does not follow this recommendation, as there the polynomial
> degree is 2, while the quadrature degree is 3
Actually, Gauss quadrature with degree+1 points in each
Hello!
I have been writing a scalar advection code in DG using the Meshworker
framework. The next step is to incorporate p-adativity for which I'll need
hp elements. Is there a way of using Meshworker based code and modifying it
to accommodate hp elements?
Thank you for your effort.
--
Apurva
Dear Luca,
That is some fascinating information! It seems like step-44, for instance,
does not follow this recommendation, as there the polynomial degree is 2,
while the quadrature degree is 3, instead of the recommended 5 (
If you plan to use any domain that is not a square (or an affine
transformation), you have to make sure you integrate exactly the product of two
polynomials of order degree and of the determinant of the Jacobian. This last
term is constant only for simple meshes, but it is the square root of a