diff options
Diffstat (limited to 'lib/Numeric/GSL/Fourier.hs')
-rw-r--r-- | lib/Numeric/GSL/Fourier.hs | 47 |
1 files changed, 0 insertions, 47 deletions
diff --git a/lib/Numeric/GSL/Fourier.hs b/lib/Numeric/GSL/Fourier.hs deleted file mode 100644 index 86aedd6..0000000 --- a/lib/Numeric/GSL/Fourier.hs +++ /dev/null | |||
@@ -1,47 +0,0 @@ | |||
1 | {-# LANGUAGE ForeignFunctionInterface #-} | ||
2 | ----------------------------------------------------------------------------- | ||
3 | {- | | ||
4 | Module : Numeric.GSL.Fourier | ||
5 | Copyright : (c) Alberto Ruiz 2006 | ||
6 | License : GPL-style | ||
7 | |||
8 | Maintainer : Alberto Ruiz (aruiz at um dot es) | ||
9 | Stability : provisional | ||
10 | Portability : uses ffi | ||
11 | |||
12 | Fourier Transform. | ||
13 | |||
14 | <http://www.gnu.org/software/gsl/manual/html_node/Fast-Fourier-Transforms.html#Fast-Fourier-Transforms> | ||
15 | |||
16 | -} | ||
17 | ----------------------------------------------------------------------------- | ||
18 | module Numeric.GSL.Fourier ( | ||
19 | fft, | ||
20 | ifft | ||
21 | ) where | ||
22 | |||
23 | import Data.Packed.Internal | ||
24 | import Data.Complex | ||
25 | import Foreign.C.Types | ||
26 | import System.IO.Unsafe (unsafePerformIO) | ||
27 | |||
28 | genfft code v = unsafePerformIO $ do | ||
29 | r <- createVector (dim v) | ||
30 | app2 (c_fft code) vec v vec r "fft" | ||
31 | return r | ||
32 | |||
33 | foreign import ccall unsafe "gsl-aux.h fft" c_fft :: CInt -> TCVCV | ||
34 | |||
35 | |||
36 | {- | Fast 1D Fourier transform of a 'Vector' @(@'Complex' 'Double'@)@ using /gsl_fft_complex_forward/. It uses the same scaling conventions as GNU Octave. | ||
37 | |||
38 | >>> fft (fromList [1,2,3,4]) | ||
39 | fromList [10.0 :+ 0.0,(-2.0) :+ 2.0,(-2.0) :+ 0.0,(-2.0) :+ (-2.0)] | ||
40 | |||
41 | -} | ||
42 | fft :: Vector (Complex Double) -> Vector (Complex Double) | ||
43 | fft = genfft 0 | ||
44 | |||
45 | -- | The inverse of 'fft', using /gsl_fft_complex_inverse/. | ||
46 | ifft :: Vector (Complex Double) -> Vector (Complex Double) | ||
47 | ifft = genfft 1 | ||