##
## p c h i p . R Test suite
##
pchip <- pracma::pchip
x <- c(1, 2, 3, 4, 5, 6)
y <- c(16, 18, 21, 17, 15, 12)
xs <- c(1.5, 2.5, 3.5, 4.5, 5.5)
ys <- pchip(x, y, xs)
# ys <- interp1(x, y, xs, method="cubic") # the same
# 16.88750 19.80000 19.33333 15.96667 13.63750
yml <- c(16.887499999999999, 19.800000000000001, 19.333333333333332,
15.966666666666667, 13.637499999999999)
all.equal(ys, yml, tolerance = 1e-7)