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

Cardinal B-spline derivatives [CI SKIP]

Browse Source
This commit is contained in:
NAThompson
2019-08-11 14:52:47 -04:00
parent 30f95e4220
commit e327be5887
3 changed files with 148 additions and 51 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
test/cardinal_b_spline_test.cpp
View File

@@ -19,6 +19,7 @@ using boost::multiprecision::float128;
using std::abs;
using boost::math::cardinal_b_spline;
using boost::math::cardinal_b_spline_prime;
using boost::math::forward_cardinal_b_spline;
template<class Real>
@@ -27,6 +28,7 @@ void test_box()
Real t = cardinal_b_spline<0>(Real(1.1));
Real expected = 0;
CHECK_ULP_CLOSE(expected, t, 0);
CHECK_ULP_CLOSE(expected, cardinal_b_spline_prime<0>(Real(1.1)), 0);
t = cardinal_b_spline<0>(Real(-1.1));
expected = 0;
@@ -37,6 +39,8 @@ void test_box()
{
expected = 1;
CHECK_ULP_CLOSE(expected, cardinal_b_spline<0>(t), 0);
expected = 0;
CHECK_ULP_CLOSE(expected, cardinal_b_spline_prime<0>(Real(1.1)), 0);
}
for (t = h; t < 1; t += h)
@@ -65,6 +69,15 @@ void test_hat()
{
std::cerr << " Problem at t = " << t << "\n";
}
if (t < 0) {
CHECK_ULP_CLOSE(Real(1), cardinal_b_spline_prime<1>(t), 0);
}
else if (t == 0) {
CHECK_ULP_CLOSE(Real(0), cardinal_b_spline_prime<1>(t), 0);
}
else if (t > 0) {
CHECK_ULP_CLOSE(Real(-1), cardinal_b_spline_prime<1>(t), 0);
}
}
for (t = 0; t < 2; t += h)