Expose spline BC to C/C++ (#31)

Author: perazz

include/fitpack_curves_c.h

@@ -59,6 +59,8 @@ FP_FLAG  fitpack_curve_c_all_derivatives(fitpack_curve_c *self, FP_REAL x, FP_RE
 FP_REAL fitpack_curve_c_smoothing      (fitpack_curve_c *self);
 FP_REAL fitpack_curve_c_mse            (fitpack_curve_c *self);
 FP_SIZE  fitpack_curve_c_degree         (fitpack_curve_c *self);
+void    fitpack_curve_c_set_bc         (fitpack_curve_c *self, FP_FLAG bc);
+FP_FLAG fitpack_curve_c_get_bc         (fitpack_curve_c *self);
 
 #ifdef __cplusplus
 }

include/fpCurve.hpp

@@ -133,6 +133,18 @@ class fpCurve
            return fitpack_curve_c_integral(&cptr, from, to);
         }
 
+        // Set spline behavior outside the support
+        void set_bc(FP_FLAG bc)
+        {
+           fitpack_curve_c_set_bc(&cptr, bc);
+        }
+
+        // Get spline behavior outside the support
+        FP_FLAG get_bc()
+        {
+           return fitpack_curve_c_get_bc(&cptr);
+        }
+
         // Get fourier coefficients
         FP_FLAG fourier(const vector<FP_REAL> &alpha,
                               vector<FP_REAL> &A, vector<FP_REAL> &B)

src/fitpack_curves_c.f90

@@ -24,7 +24,8 @@
 ! **************************************************************************************************
 module fitpack_curves_c
      use fitpack_curves, only: fitpack_curve
-     use fitpack_core, only: FP_SIZE, FP_REAL, FP_FLAG, FP_BOOL
+     use fitpack_core, only: FP_SIZE, FP_REAL, FP_FLAG, FP_BOOL, &
+                             OUTSIDE_EXTRAPOLATE, OUTSIDE_ZERO, OUTSIDE_NOT_ALLOWED, OUTSIDE_NEAREST_BND
      use, intrinsic :: iso_c_binding
      implicit none
      private
@@ -48,6 +49,8 @@ module fitpack_curves_c
      public :: fitpack_curve_c_smoothing
      public :: fitpack_curve_c_mse
      public :: fitpack_curve_c_degree
+     public :: fitpack_curve_c_set_bc
+     public :: fitpack_curve_c_get_bc
 
      !> Opaque-pointer C derived type
      type, public, bind(C) :: fitpack_curve_c
@@ -358,5 +361,36 @@ integer(FP_SIZE) function fitpack_curve_c_degree(this) bind(c,name='fitpack_curv
 
      end function fitpack_curve_c_degree
 
+     !> Set spline behavior outside the support
+     subroutine fitpack_curve_c_set_bc(this, bc) bind(c,name='fitpack_curve_c_set_bc')
+        type(fitpack_curve_c), intent(inout) :: this
+        integer(FP_FLAG), intent(in), value :: bc
+
+        type(fitpack_curve), pointer :: fcurve
+
+        !> Get object; allocate it in case
+        call fitpack_curve_c_pointer(this,fcurve)
+
+        !> Validate extrapolation behavior value
+        if (bc /= OUTSIDE_EXTRAPOLATE .and. bc /= OUTSIDE_ZERO .and. &
+            bc /= OUTSIDE_NOT_ALLOWED .and. bc /= OUTSIDE_NEAREST_BND) return
+
+        fcurve%bc = bc
+
+     end subroutine fitpack_curve_c_set_bc
+
+     !> Get spline behavior outside the support
+     integer(FP_FLAG) function fitpack_curve_c_get_bc(this) bind(c,name='fitpack_curve_c_get_bc')
+        type(fitpack_curve_c), intent(inout) :: this
+
+        type(fitpack_curve), pointer :: fcurve
+
+        !> Get object; allocate it in case
+        call fitpack_curve_c_pointer(this,fcurve)
+
+        fitpack_curve_c_get_bc = fcurve%bc
+
+     end function fitpack_curve_c_get_bc
+
 
 end module fitpack_curves_c