From 1bbdcbedc9f9b482824de3c47df961178ab9a80a Mon Sep 17 00:00:00 2001 From: Julius Witte Date: Wed, 25 Sep 2019 12:45:12 -0400 Subject: [PATCH] Store 1D quadrature --- include/deal.II/matrix_free/shape_info.h | 6 ++++++ include/deal.II/matrix_free/shape_info.templates.h | 1 + 2 files changed, 7 insertions(+) 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, -- 2.39.5