From b8396cc443180c216fb79974a642bab754d0a7e6 Mon Sep 17 00:00:00 2001
From: Klaus <klaus.boehnlein@tu-dresden.de>
Date: Fri, 17 May 2024 16:19:22 +0200
Subject: [PATCH] test different QuadRule
---
.../discretekirchhoffbendingenergyprestrained.hh | 13 +++++++++++--
1 file changed, 11 insertions(+), 2 deletions(-)
diff --git a/dune/microstructure/energies/discretekirchhoffbendingenergyprestrained.hh b/dune/microstructure/energies/discretekirchhoffbendingenergyprestrained.hh
index a512bbd2..def8d5b9 100644
--- a/dune/microstructure/energies/discretekirchhoffbendingenergyprestrained.hh
+++ b/dune/microstructure/energies/discretekirchhoffbendingenergyprestrained.hh
@@ -592,19 +592,28 @@ namespace Dune::GFE
*/
RT harmonicEnergy = 0;
- // Gauss-Quadrature:
+ // Gauss-Legendre-Quadrature (default):
const auto &quadRule = QuadratureRules<double, gridDim>::rule(lagrangeLFE_.type(), quadOrder);
+ //TEST Gauß-Jacobi-Quadrature formula:
+ // const auto &quadRule = QuadratureRules<double, gridDim>::rule(lagrangeLFE_.type(), quadOrder, QuadratureType::GaussJacobi_1_0);
+
+
// Trapezoidal-rule:
// std::vector<Dune::QuadraturePoint<double,2>> quadRule = { {{0.0,0.0}, 1.0/6.0}, {{1.0,0.0}, 1.0/6.0}, {{0.0,1.0}, 1.0/6.0} };
// Trapezoidal-rule + evaluation on edge centers (only for testing purposes):
// std::vector<Dune::QuadraturePoint<double,2>> quadRule = { {{0.0,0.0}, 1.0/6.0}, {{0.5,0.0}, 1.0/6.0}, {{1.0,0.0}, 1.0/6.0}, {{0.0,0.5}, 1.0/6.0}, {{0.5,0.5}, 1.0/6.0}, {{0.0,1.0}, 1.0/6.0} };
// int count = 0;
- // std::cout << "Number of quadrature points: " << quadRule.size() << std::endl;
+
+
+ std::cout << "Number of quadrature points: " << quadRule.size() << std::endl;
+
+
for (auto&& quadPoint : quadRule)
{
const auto& quadPos = quadPoint.position();
const auto integrationElement = geometry.integrationElement(quadPos);
+ std::cout << "quadPos:" << quadPos << std::endl;
// exit(0);
--
GitLab