]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Store 1D quadrature
authorJulius Witte <julius.witte@iwr.uni-heidelberg.de>
Wed, 25 Sep 2019 16:45:12 +0000 (12:45 -0400)
committerJulius Witte <julius.witte@iwr.uni-heidelberg.de>
Thu, 26 Sep 2019 13:59:50 +0000 (09:59 -0400)
include/deal.II/matrix_free/shape_info.h
include/deal.II/matrix_free/shape_info.templates.h

index 79804d776ca8edf70776ff22fd99e45d14280355..6442772347ee1ab1dc466f76d6032a156b8c86d4 100644 (file)
@@ -232,6 +232,12 @@ namespace internal
        */
       AlignedVector<Number> 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
index 959e2ebecaa4bc5c1878dde9f70d12368d6043bf..998b10cdfb20e4fe57cbd9f2fc624f624a74faf6 100644 (file)
@@ -75,6 +75,7 @@ namespace internal
                               const FiniteElement<dim> &fe_in,
                               const unsigned int        base_element_number)
     {
+      this->quadrature             = quad;
       const FiniteElement<dim> *fe = &fe_in.base_element(base_element_number);
 
       Assert(fe->n_components() == 1,

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.