Skip to content
Snippets Groups Projects
Commit 16216cd3 authored by Christian Engwer's avatar Christian Engwer
Browse files

better epsilon

[[Imported from SVN: r4376]]
parent 4d8d4c37
Branches
Tags
No related merge requests found
......@@ -11,7 +11,7 @@
double success = true;
template<class ctype, int dim>
bool checkQuadrature(Dune::GeometryType t, int p)
void checkQuadrature(Dune::GeometryType t, int p)
{
double volume = 0;
// Quadratures
......@@ -27,7 +27,7 @@ bool checkQuadrature(Dune::GeometryType t, int p)
}
if (std::abs(volume -
Dune::ReferenceElements<ctype, dim>::general(t).volume())
> 10*std::numeric_limits<double>::epsilon())
> 20*std::numeric_limits<double>::epsilon())
{
std::cerr << "Error: Quadrature for " << t << " and order=" << p
<< " does not sum to volume of RefElem" << std::endl;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment