Skip to content

#113 Quad weights don't sum up to element volume

Metadata

Property Value
Reported by Oliver Sander (oliver.sander@tu-dresden.de)
Reported at Mar 20, 2006 13:12
Type Bug Report
Version Git (pre2.4) [autotools]
Operating System Unspecified / All
Closed by Christian Engwer (christi@conan.iwr.uni-heidelberg.de)
Closed at Mar 29, 2006 22:35
Closed in version 1.0
Resolution Fixed
Comment

Description

The quadrature weights for the first-order quadrature rule for tetrahedra don't sum up to the reference element volume (1/6). Instead they sum up to one. Here is a test program

#include <config.h>

#include <dune/quadrature/quadraturerules.hh>

int main () { Dune::GeometryType type(Dune::GeometryType::simplex, 3); const Dune::QuadratureRule<double, 3>& quad = Dune::QuadratureRules<double, 3>::rule(type, 1); std::cout << "size: " << quad.size() << std::endl; std::cout << "only weight: " << quad[0].weight() << std::endl; }

Result: size: 1 only weight: 1