From 7b10b74bf252dfb1b5a7c3e5b60fa47a5cd4d9a9 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@uliege.be>
Date: Tue, 8 Feb 2022 21:29:15 +0100
Subject: [PATCH] not clear which rule to use for order quad 2: with the
current choice Gauss2 has 7 points and GaussComposite2 has only 4!
---
src/numeric/GaussQuadratureQuad.cpp | 2 ++
1 file changed, 2 insertions(+)
diff --git a/src/numeric/GaussQuadratureQuad.cpp b/src/numeric/GaussQuadratureQuad.cpp
index 2a5d9994db..661bfb128e 100644
--- a/src/numeric/GaussQuadratureQuad.cpp
+++ b/src/numeric/GaussQuadratureQuad.cpp
@@ -90,6 +90,8 @@ IntPt GQQ16[16] = {
{{xq16[12], yq16[12], 0}, pq16[12]}, {{xq16[13], yq16[13], 0}, pq16[13]},
{{xq16[14], yq16[14], 0}, pq16[14]}, {{xq16[15], yq16[15], 0}, pq16[15]}};
+// Not clear which rule to use for order 2 - with this choice Gauss2 has 7
+// points, but only 4 for the tensor rule...
static IntPt *GQQ[3] = {GQQ1, GQQ3, GQQ7};
static int GQQnPt[3] = {1, 3, 7};
static std::vector<IntPt *> GQQGL(40, nullptr);
--
GitLab