#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