From: Simon Sticko <simon@sticko.se>
Date: Fri, 13 May 2022 09:29:11 +0000 (+0200)
Subject: Set is_tensor_product_flag in Quadrature::initialize(..)
X-Git-Tag: v9.4.0-rc1~222^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b3407f138503cddbecf90f51dda30d8996e2f757;p=dealii.git

Set is_tensor_product_flag in Quadrature::initialize(..)

If you first create a quadrature which is a tensor product and then
call initialize with points and weights which does not correspond
to a tensor product, you can make is_tensor_product_flag have the
wrong value. Set the flag in intialize(..) to avoid this.
---

diff --git a/source/base/quadrature.cc b/source/base/quadrature.cc
index 1d99787305..8fc7271a43 100644
--- a/source/base/quadrature.cc
+++ b/source/base/quadrature.cc
@@ -53,8 +53,9 @@ Quadrature<dim>::initialize(const std::vector<Point<dim>> &p,
                             const std::vector<double> &    w)
 {
   AssertDimension(w.size(), p.size());
-  quadrature_points = p;
-  weights           = w;
+  quadrature_points      = p;
+  weights                = w;
+  is_tensor_product_flag = dim == 1;
 }