diff options
Diffstat (limited to 'lib/Numeric/GSL/Polynomials.hs')
-rw-r--r-- | lib/Numeric/GSL/Polynomials.hs | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/lib/Numeric/GSL/Polynomials.hs b/lib/Numeric/GSL/Polynomials.hs index 42694f0..e663711 100644 --- a/lib/Numeric/GSL/Polynomials.hs +++ b/lib/Numeric/GSL/Polynomials.hs | |||
@@ -47,7 +47,8 @@ polySolve = toList . polySolve' . fromList | |||
47 | polySolve' :: Vector Double -> Vector (Complex Double) | 47 | polySolve' :: Vector Double -> Vector (Complex Double) |
48 | polySolve' v | dim v > 1 = unsafePerformIO $ do | 48 | polySolve' v | dim v > 1 = unsafePerformIO $ do |
49 | r <- createVector (dim v-1) | 49 | r <- createVector (dim v-1) |
50 | c_polySolve // vec v // vec r // check "polySolve" [v] | 50 | ww2 withVector v withVector r $ \ v r -> |
51 | c_polySolve // v // r // check "polySolve" | ||
51 | return r | 52 | return r |
52 | | otherwise = error "polySolve on a polynomial of degree zero" | 53 | | otherwise = error "polySolve on a polynomial of degree zero" |
53 | 54 | ||