summaryrefslogtreecommitdiff
path: root/lib
diff options
context:
space:
mode:
Diffstat (limited to 'lib')
-rw-r--r--lib/Data/Packed/Internal/Matrix.hs4
-rw-r--r--lib/Data/Packed/Internal/Vector.hs4
2 files changed, 4 insertions, 4 deletions
diff --git a/lib/Data/Packed/Internal/Matrix.hs b/lib/Data/Packed/Internal/Matrix.hs
index 861c72a..100b8ea 100644
--- a/lib/Data/Packed/Internal/Matrix.hs
+++ b/lib/Data/Packed/Internal/Matrix.hs
@@ -391,11 +391,11 @@ conjV :: (Storable a, RealFloat a) => Vector (Complex a) -> Vector (Complex a)
391conjV = mapVector conjugate 391conjV = mapVector conjugate
392 392
393-- | creates a complex vector from vectors with real and imaginary parts 393-- | creates a complex vector from vectors with real and imaginary parts
394toComplexV :: Element a => (Vector a, Vector a) -> Vector (Complex a) 394toComplexV :: (RealFloat a, Element a) => (Vector a, Vector a) -> Vector (Complex a)
395toComplexV (r,i) = asComplex $ flatten $ fromColumns [r,i] 395toComplexV (r,i) = asComplex $ flatten $ fromColumns [r,i]
396 396
397-- | the inverse of 'toComplex' 397-- | the inverse of 'toComplex'
398fromComplexV :: Element a => Vector (Complex a) -> (Vector a, Vector a) 398fromComplexV :: (RealFloat a, Element a) => Vector (Complex a) -> (Vector a, Vector a)
399fromComplexV z = (r,i) where 399fromComplexV z = (r,i) where
400 [r,i] = toColumns $ reshape 2 $ asReal z 400 [r,i] = toColumns $ reshape 2 $ asReal z
401 401
diff --git a/lib/Data/Packed/Internal/Vector.hs b/lib/Data/Packed/Internal/Vector.hs
index 6d39c6e..dfa2e86 100644
--- a/lib/Data/Packed/Internal/Vector.hs
+++ b/lib/Data/Packed/Internal/Vector.hs
@@ -266,13 +266,13 @@ takesV ms w | sum ms > dim w = error $ "takesV " ++ show ms ++ " on dim = " ++ (
266--------------------------------------------------------------- 266---------------------------------------------------------------
267 267
268-- | transforms a complex vector into a real vector with alternating real and imaginary parts 268-- | transforms a complex vector into a real vector with alternating real and imaginary parts
269asReal :: Vector (Complex a) -> Vector a 269asReal :: (RealFloat a, Storable a) => Vector (Complex a) -> Vector a
270--asReal v = V { ioff = 2*ioff v, idim = 2*dim v, fptr = castForeignPtr (fptr v) } 270--asReal v = V { ioff = 2*ioff v, idim = 2*dim v, fptr = castForeignPtr (fptr v) }
271asReal v = unsafeFromForeignPtr (castForeignPtr fp) (2*i) (2*n) 271asReal v = unsafeFromForeignPtr (castForeignPtr fp) (2*i) (2*n)
272 where (fp,i,n) = unsafeToForeignPtr v 272 where (fp,i,n) = unsafeToForeignPtr v
273 273
274-- | transforms a real vector into a complex vector with alternating real and imaginary parts 274-- | transforms a real vector into a complex vector with alternating real and imaginary parts
275asComplex :: Vector a -> Vector (Complex a) 275asComplex :: (RealFloat a, Storable a) => Vector a -> Vector (Complex a)
276--asComplex v = V { ioff = ioff v `div` 2, idim = dim v `div` 2, fptr = castForeignPtr (fptr v) } 276--asComplex v = V { ioff = ioff v `div` 2, idim = dim v `div` 2, fptr = castForeignPtr (fptr v) }
277asComplex v = unsafeFromForeignPtr (castForeignPtr fp) (i `div` 2) (n `div` 2) 277asComplex v = unsafeFromForeignPtr (castForeignPtr fp) (i `div` 2) (n `div` 2)
278 where (fp,i,n) = unsafeToForeignPtr v 278 where (fp,i,n) = unsafeToForeignPtr v