provisional fix for interpolation

1 files changed, 12 insertions, 0 deletions

diff --git a/packages/gsl/src/Numeric/GSL/Interpolation.hs b/packages/gsl/src/Numeric/GSL/Interpolation.hs
index d060468..6f02405 100644
--- a/packages/gsl/src/Numeric/GSL/Interpolation.hs
+++ b/packages/gsl/src/Numeric/GSL/Interpolation.hs
@@ -1,3 +1,5 @@
+{-# LANGUAGE MagicHash, UnboxedTuples #-}
+
 {- |
 Module      :  Numeric.GSL.Interpolation
 Copyright   :  (c) Matthew Peddie 2015
@@ -40,6 +42,10 @@ import Foreign.Storable(peek)
 import Numeric.GSL.Internal
 import System.IO.Unsafe(unsafePerformIO)
 
+-- FIXME
+import qualified Data.Vector.Storable as S
+import GHC.Base (IO(..), realWorld#)
+
 data InterpolationMethod = Linear
                          | Polynomial
                          | CSpline
@@ -59,6 +65,12 @@ methodToInt AkimaPeriodic = 5
 dim :: Numeric t => Vector t -> Int
 dim = size
 
+-- FIXME
+appVector f x = unsafeInlinePerformIO (S.unsafeWith x (return . f))
+
+unsafeInlinePerformIO (IO f) = case f realWorld# of
+    (# _, x #) -> x
+
 applyCFun hsname cname fun mth xs ys x
   | dim xs /= dim ys = error $
                          "Error: Vectors of unequal sizes " ++