refactor: Correct spline naming/types in python module

Signed-off-by: Sietze van Buuren <s.van.buuren@gmail.com>

Author: swvanbuuren

include/cubic_splines_2d.hpp

@@ -50,7 +50,7 @@ class AkimaSpline2D : public CubicInterpND<T, N>
 
 
 
-template <typename T, cip::BoundaryConditionType BC=cip::BoundaryConditionType::NotAKnot, std::size_t N=2>
+template <typename T, cip::BoundaryConditionType BC=cip::BoundaryConditionType::Natural, std::size_t N=2>
 class NaturalSpline2D : public CubicInterpND<T, N>
 {
     using Vector = std::vector<T>;
@@ -71,4 +71,25 @@ class NaturalSpline2D : public CubicInterpND<T, N>
 };
 
 
+template <typename T, cip::BoundaryConditionType BC=cip::BoundaryConditionType::Periodic, std::size_t N=2>
+class NaturalPeriodicSpline2D : public CubicInterpND<T, N>
+{
+    using Vector = std::vector<T>;
+    using Mdspan1D = std::mdspan<T, std::dextents<std::size_t, 1>, std::layout_stride>;
+    using VectorN = cip::VectorN<T, N>;
+public:
+    NaturalPeriodicSpline2D(const Vector &_x, const Vector &_y, const VectorN &_f)
+    : CubicInterpND<T, N>(_f, _x, _y)
+    {
+        this->build(_f, _x, _y);
+    }
+    ~NaturalPeriodicSpline2D() {}
+
+    Vector calc_slopes(const Vector &x, const Mdspan1D &f) const override
+    {
+        return natural_spline_slopes<T, BC>(x, f);
+    }
+};
+
+
 } // namespace cip

include/cubic_splines_3d.hpp

@@ -7,7 +7,7 @@
 namespace cip {
 
 
-template <typename T, cip::BoundaryConditionType BC=cip::BoundaryConditionType::NotAKnot, std::size_t N=3>
+template <typename T, cip::BoundaryConditionType BC=cip::BoundaryConditionType::Natural, std::size_t N=3>
 class NaturalSpline3D : public CubicInterpND<T, N>
 {
     using Vector = std::vector<T>;

src/cubinterpp_py_module.cpp

@@ -60,6 +60,10 @@ PYBIND11_MODULE(cubinterpp_py, m) {
         .def(py::init<DoubleVector, DoubleVector, DoubleVector2>())
         .def("eval", &cip::NaturalSpline2D<double>::eval<double, double>, py::return_value_policy::reference_internal);
 
+    py::class_<cip::NaturalPeriodicSpline2D<double>>(m, "NaturalPeriodicSpline2D")
+        .def(py::init<DoubleVector, DoubleVector, DoubleVector2>())
+        .def("eval", &cip::NaturalPeriodicSpline2D<double>::eval<double, double>, py::return_value_policy::reference_internal);
+
     py::class_<cip::NaturalSpline3D<double>>(m, "NaturalSpline3D")
         .def(py::init<DoubleVector, DoubleVector, DoubleVector, DoubleVector3>())
         .def("eval", &cip::NaturalSpline3D<double>::eval<double, double, double>, py::return_value_policy::reference_internal);