Dear S. A. Mohseni,

the Gauss points are assigned in that program as always in deal.II:

 QGauss<dim> quadrature_formula(degree + 2);

and in the argument, you can change to whatever number you want:

QGauss<dim> quadrature_formula(n);

with n=1 or n=2 or n = 3 for instance.

Best regards,

Thomas W.

--
++--------------------------------------------++
Dr. Thomas Wick
Maitre de conferences / Assistant Professor

Centre de Mathematiques Appliquees (CMAP)
Ecole Polytechnique
91128 Palaiseau cedex, France

Email:  [email protected]
Phone:  0033 1 69 33 4579
www:    http://www.cmap.polytechnique.fr/~wick/
++--------------------------------------------++
--

On 03/30/2017 03:13 PM, 'Seyed Ali Mohseni' via deal.II User Group wrote:
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<dim>(degree) where degree cannot be 0. This means degree should always be >1 here in order to make things work, but how can one compute the examples with 1 or 4 Gauss points? Is the code not suited for what I am trying to achieve or where are the parts I have to change in order to make it work for less integration points. 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?

Thank you in advance.

Kind regards,
S. A. Mohseni
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected] <mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout.

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
For more options, visit https://groups.google.com/d/optout.

Reply via email to