boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-01-19 04:22:09 +00:00
Code Issues Projects Releases Wiki Activity

Additional test for quintic_hermite

Browse Source
This commit is contained in:
NAThompson
2020-01-26 14:19:33 +08:00
parent 8213c19b30
commit a035ffca87
1 changed files with 33 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

33
test/quintic_hermite_test.cpp
View File

@@ -112,6 +112,35 @@ void test_cubic()
}
}
template<typename Real>
void test_quartic()
{
std::vector<Real> x{0,1,2,3, 4,5,6,7,8,9, 10, 11};
std::vector<Real> y(x.size());
for (size_t i = 0; i < y.size(); ++i)
{
y[i] = x[i]*x[i]*x[i]*x[i]/24;
}
std::vector<Real> dydx(x.size());
for (size_t i = 0; i < y.size(); ++i) {
dydx[i] = x[i]*x[i]*x[i]/6;
}
std::vector<Real> d2ydx2(x.size());
for (size_t i = 0; i < y.size(); ++i) {
d2ydx2[i] = x[i]*x[i]/2;
}
auto qh = quintic_hermite(std::move(x), std::move(y), std::move(dydx), std::move(d2ydx2));
for (Real t = 0; t <= 9; t += 0.125) {
CHECK_ULP_CLOSE(Real(t*t*t*t)/24, qh(t), 43);
}
}
template<typename Real>
void test_interpolation_condition()
{
@@ -152,18 +181,21 @@ int main()
test_linear<float>();
test_quadratic<float>();
test_cubic<float>();
test_quartic<float>();
test_interpolation_condition<float>();
test_constant<double>();
test_linear<double>();
test_quadratic<double>();
test_cubic<double>();
test_quartic<double>();
test_interpolation_condition<double>();
test_constant<long double>();
test_linear<long double>();
test_quadratic<long double>();
test_cubic<long double>();
test_quartic<long double>();
test_interpolation_condition<long double>();
#ifdef BOOST_HAS_FLOAT128
@@ -171,6 +203,7 @@ int main()
test_linear<float128>();
test_quadratic<float128>();
test_cubic<float128>();
test_quartic<float128>();
#endif
return boost::math::test::report_errors();