This MR fixes several shortcomings in the quadrature test. It changes all
test-quadrature and use ADL for the lookup of
pow to allow other than standard floating point types to be used in quadrature rules. The test is extended to
In the previous implementation of the
test-quadrature.hh and in the quadrature rules were some problems:
- In several places explicit use of
- qualified names for
- Specialization for 1d rules to be constructed from string if
numeric_limits<ctype>::is_specialized. This is is probably to correct, since, e.g.,
Dune::Float128is specialized but should be constructed from string instead of from
Dune::GMPFielddoes not work in the current implementation (I have not found the reason why)
- The test fails if LAPACK is not enabled
These problems are cleaned up:
ctypeeverywhere and provide replacement for
numeric_limits<ctype>::epsilonif not available.
using std::abs; using std::powover qualified names
- Specialization for 1d rules is done with
std::is_floating_point<ctype>::value. This should be correct, since in this specialization the numbers are constructed from floating-point numbers in
GMPFieldis not fixed, but
Float128is tested instead.
GaussJacobi_n_0test is disabled in case LAPACK is not found.
Closes #24 (closed)