13 #ifndef trapezoid_test_suite_hpp
14 #define trapezoid_test_suite_hpp
26 template<
typename data__>
61 data__ x1(1.0), x2(3.0), F_exact(std::exp(x2)-std::exp(x1)), F_quad;
63 std::vector<data__> x(npts), f(npts);
66 for (n=0; n<npts; ++n)
67 x[n]=(n)*(x2-x1)/(npts-1)+x1;
70 for (n=0; n<npts; ++n)
76 F_quad=quad(x[1]-x[0], f.begin(), f.end());
78 if (
typeid(data__)==
typeid(float))
80 TEST_ASSERT_DELTA(1, F_quad/F_exact, 2e-5);
84 TEST_ASSERT_DELTA(1, F_quad/F_exact, 4e-7);
90 data__ x1(1.0), x2(3.0), F_exact(std::exp(x2)-std::exp(x1)), F_quad;
92 std::vector<data__> x(npts), f(npts);
95 for (n=0; n<npts; ++n)
99 for (n=0; n<npts; ++n)
105 F_quad=quad(x.begin(), f.begin(), f.end());
107 TEST_ASSERT_DELTA(1, F_quad/F_exact, 5e-7);
void nonuniform_points_test()
Definition: trapezoid_test_suite.hpp:88
void AddTests(const long double &)
Definition: trapezoid_test_suite.hpp:42
Definition: trapezoid.hpp:25
trapezoid_test_suite()
Definition: trapezoid_test_suite.hpp:49
void uniform_points_test()
Definition: trapezoid_test_suite.hpp:59
void AddTests(const double &)
Definition: trapezoid_test_suite.hpp:36
~trapezoid_test_suite()
Definition: trapezoid_test_suite.hpp:54
void AddTests(const float &)
Definition: trapezoid_test_suite.hpp:30
Definition: trapezoid_test_suite.hpp:27