From: Julius Witte Date: Wed, 25 Sep 2019 16:45:12 +0000 (-0400) Subject: Store 1D quadrature X-Git-Tag: v9.2.0-rc1~1010^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=1bbdcbedc9f9b482824de3c47df961178ab9a80a;p=dealii.git Store 1D quadrature --- diff --git a/include/deal.II/matrix_free/shape_info.h b/include/deal.II/matrix_free/shape_info.h index 79804d776c..6442772347 100644 --- a/include/deal.II/matrix_free/shape_info.h +++ b/include/deal.II/matrix_free/shape_info.h @@ -232,6 +232,12 @@ namespace internal */ AlignedVector hessians_within_subface[2]; + /** + * We store a copy of the one-dimensional quadrature formula + * used for initialization. + */ + Quadrature<1> quadrature; + /** * Renumbering from deal.II's numbering of cell degrees of freedom to * lexicographic numbering used inside the FEEvaluation schemes of the diff --git a/include/deal.II/matrix_free/shape_info.templates.h b/include/deal.II/matrix_free/shape_info.templates.h index 959e2ebeca..998b10cdfb 100644 --- a/include/deal.II/matrix_free/shape_info.templates.h +++ b/include/deal.II/matrix_free/shape_info.templates.h @@ -75,6 +75,7 @@ namespace internal const FiniteElement &fe_in, const unsigned int base_element_number) { + this->quadrature = quad; const FiniteElement *fe = &fe_in.base_element(base_element_number); Assert(fe->n_components() == 1,