> To be honest, I don't even see the point of having 9 integration points for
> a 4-noded quadrilateral or am I overseeing something?
Yes, you are missing some background information:
Quadrature formulas are exact for polynomials up to a certain degree.
For Gauss quadrature you need degree n to in
Dear S. A. Mohseni,
the Gauss points are assigned in that program as always in deal.II:
QGauss quadrature_formula(degree + 2);
and in the argument, you can change to whatever number you want:
QGauss quadrature_formula(n);
with n=1 or n=2 or n = 3 for instance.
Best regards,
Thomas W.
--
+
Dear Timo Heister and Thomas Wick,
I am trying to run your phase-field crack propagation example, namely the
single edge notched tension test, with less than 9 Gauss points.
Hence, the fe variable is initialized using FE_Q(degree) where degree
cannot be 0. This means degree should always be >1