On 05/14/2015 12:51 AM, briand@aracnet.com wrote:
Well what I was_really_ trying to do was to get to a Vector of (Complex Double).
Interestingly, the following, where ys is type Vector.Storable (Complex Double),
ys_array = listArray (0, V.length ys -1) (V.toList ys) y = V.fromList $ map (\z -> realPart z :+ imagPart z) (elems (dft ys_array))
worked just fine, but it's not obvious to me that it should have. Seems like the realPart/imagPart should have failed to work because the argument is not Complex Double it's Complex FFTWReal, so I should have ended up with (Complex FFTWReal) when I really wanted (Complex Double) and fail.
BTW, it seems to me that it would be awfully nice if FFT delivered results as Vectors instead of the Array thingy. It looked to me like V.Storable is something which is supposed to ease the marshalling to and from the C world.
Would it be a whole lot of work to change the interface to Vector.Storable ?
Among other advantages, it makes interfacing to HMatrix easier. Hi Brian,
I think the problem is in the way of seeing the types and type classes. First, the type of `ys` is not `Vector.Storable (Complex Double)`, it is Vector (Complex Double). `FFTWReal` is not a type, it is a type class. It is used as a constraint for types, for example, the type of `dft`: dft :: (Ix i, Shapable i, FFTWReal r) => CArray i (Complex r) -> CArray i (Complex r) this means that the type `r` must be instance o FFTWReal type class. So, the type of `dft ys_array` is `CArray i (Complex Double)`, since the type of `ys_array` is `CArray i (Complex Double)`, that means that the `r` in the type of `dft` is `Double`. So, the `map` is unnecessary. You can try: ``` vectorToArray v = listArray (0, V.length v - 1) (V.toList v) arrayToVector = V.fromList . elems fft :: FFTWReal r => V.Vector (Complex r) -> V.Vector (Complex r) fft = arrayToVector . dft . vectorToArray ``` Hope it can help! -- Leza Morais Lutonda, Lemol-C Electronics & Telecommunications Engineer Software Development and Architecture Enthusiast http://lemol.github.io 50 Aniversario de la Cujae. Inaugurada por Fidel el 2 de diciembre de 1964 http://cujae.edu.cu