#include <stdexcept>
#include <catch2/catch.hpp>
#include "CubicInterpolation.hpp"
namespace datasketches {
TEST_CASE("hll tables: interpolation exception", "[hll_tables]") {
REQUIRE_THROWS_AS(CubicInterpolation<>::usingXAndYTables(-1.0), std::invalid_argument);
REQUIRE_THROWS_AS(CubicInterpolation<>::usingXAndYTables(1e12), std::invalid_argument);
}
TEST_CASE("hll tables: check corner case", "[hll_tables]") {
int len = 10;
double xArr[] = {1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0};
double yArr[] = {2.0, 4.0, 8.0, 16.0, 32.0, 64.0, 128.0, 256.0, 512.0, 1024.0};
double x = xArr[len - 1];
double y = CubicInterpolation<>::usingXAndYTables(xArr, yArr, len, x);
double yExp = yArr[len - 1];
REQUIRE(y == yExp);
}
}