1 files changed, 31 insertions, 0 deletions

diff --git a/packages/sundials/src/helpers.c b/packages/sundials/src/helpers.c
index d51310c..eab5ac9 100644
--- a/packages/sundials/src/helpers.c
+++ b/packages/sundials/src/helpers.c
@@ -8,6 +8,12 @@
 #include <sundials/sundials_types.h>       /* definition of type realtype          */
 #include <sundials/sundials_math.h>
 
+#include "farkode.h"
+
+#include <HsFFI.h>
+#include "Main_stub.h"
+
+
 /* Check function return value...
     opt == 0 means SUNDIALS function allocates memory so check if
              returned NULL pointer
@@ -56,6 +62,31 @@ int f(realtype t, N_Vector y, N_Vector ydot, void *user_data)
   return 0;                                   /* return with success */
 }
 
+int FARK_IMP_FUN(realtype *T, realtype *Y, realtype *YDOT,
+                 long int *IPAR, realtype *RPAR, int *IER) {
+  realtype t = *T;
+  realtype u = Y[0];
+  realtype lamda = -100.0;
+  YDOT[0] = singleEq(t, u);
+  return 0;
+}
+
+/* C interface to user-supplied FORTRAN function FARKIFUN; see
+   farkode.h for further details */
+int FARKfi(realtype t, N_Vector y, N_Vector ydot, void *user_data) {
+
+  int ier;
+  realtype *ydata, *dydata;
+  FARKUserData ARK_userdata;
+  ydata  = N_VGetArrayPointer(y);
+  dydata = N_VGetArrayPointer(ydot);
+  ARK_userdata = (FARKUserData) user_data;
+
+  FARK_IMP_FUN(&t, ydata, dydata, ARK_userdata->ipar,
+               ARK_userdata->rpar, &ier);
+  return(ier);
+}
+
 /* Jacobian routine to compute J(t,y) = df/dy. */
 int Jac(realtype t, N_Vector y, N_Vector fy, SUNMatrix J,
         void *user_data, N_Vector tmp1, N_Vector tmp2, N_Vector tmp3)