Free cookie consent management tool by TermsFeed Policy Generator

source: branches/HeuristicLab.Problems.GaussianProcessTuning/ILNumerics.2.14.4735.573/Native/FFT/ILFFTW3FFT.cs @ 11826

Last change on this file since 11826 was 9102, checked in by gkronber, 12 years ago

#1967: ILNumerics source for experimentation
File size: 46.1 KB
Line 
1///
2///    This file is part of ILNumerics Community Edition.
3///
4///    ILNumerics Community Edition - high performance computing for applications.
5///    Copyright (C) 2006 - 2012 Haymo Kutschbach, http://ilnumerics.net
6///
7///    ILNumerics Community Edition is free software: you can redistribute it and/or modify
8///    it under the terms of the GNU General Public License version 3 as published by
9///    the Free Software Foundation.
10///
11///    ILNumerics Community Edition is distributed in the hope that it will be useful,
12///    but WITHOUT ANY WARRANTY; without even the implied warranty of
13///    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14///    GNU General Public License for more details.
15///
16///    You should have received a copy of the GNU General Public License
17///    along with ILNumerics Community Edition. See the file License.txt in the root
18///    of your distribution package. If not, see <http://www.gnu.org/licenses/>.
19///
20///    In addition this software uses the following components and/or licenses:
21///
22///    =================================================================================
23///    The Open Toolkit Library License
24///   
25///    Copyright (c) 2006 - 2009 the Open Toolkit library.
26///   
27///    Permission is hereby granted, free of charge, to any person obtaining a copy
28///    of this software and associated documentation files (the "Software"), to deal
29///    in the Software without restriction, including without limitation the rights to
30///    use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
31///    the Software, and to permit persons to whom the Software is furnished to do
32///    so, subject to the following conditions:
33///
34///    The above copyright notice and this permission notice shall be included in all
35///    copies or substantial portions of the Software.
36///
37///    =================================================================================
38///   
39
40#pragma warning disable 1591
41
42using System;
43using System.Collections.Generic;
44using System.Text;
45using System.Security;
46using System.Runtime.InteropServices;
47using ILNumerics.Exceptions; 
48using ILNumerics.Misc;
49
50namespace ILNumerics.Native {
51   
52    /// <summary>
53    /// Wrapper for FFT interface using FFTW3 libs
54    /// </summary>
55    public unsafe class ILFFTW3FFT : IILFFT, IDisposable {
56
57        #region pinvoke definitions
58       
59        /// <summary>
60        /// This struct is used to define (n-dimensional) transform sizes
61        /// </summary>
62        /// <remarks>This struct is only user in C-API. It is NOT used for the Fortran interface (i.e. it's not used in here)!!! </remarks>
63        [StructLayout(LayoutKind.Sequential)]
64        private struct fftw_iodim {
65            /// <summary>
66            /// length of dimension
67            /// </summary>
68            public int n;
69            /// <summary>
70            /// input stride
71            /// </summary>
72            public int iS;
73            /// <summary>
74            /// output stride
75            /// </summary>
76            public int oS;
77        }
78
79        #region DOUBLE PRECISION INTERFACE ********* COMPLEX
80        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
81    private static extern void dfftw_plan_dft_1d_(ref IntPtr plan, ref int n, [In,Out] complex* input, [In,Out] complex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
82        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
83    private static extern void dfftw_plan_dft_2d_(ref IntPtr plan, ref int nx, ref int ny, [In,Out] complex* input, [In,Out] complex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
84        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
85    private static extern void dfftw_plan_dft_3d_(ref IntPtr plan, ref int nx, ref int ny, ref int nz, [In,Out] complex* input, [In,Out] complex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
86        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
87    private static extern void dfftw_plan_dft_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] complex* input, [In,Out] complex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
88       
89        /****************************************************************** REAL-COMPLEX ***/
90        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
91    private static extern void dfftw_plan_dft_r2c_1d_(ref IntPtr plan, ref int n, [In,Out] double* input, [In,Out] complex* output, ref FFTW_PLANCREATION flags);
92        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
93    private static extern void dfftw_plan_dft_r2c_2d_(ref IntPtr plan, ref int nx, ref int ny, [In,Out] double* input, [In,Out] complex* output, ref FFTW_PLANCREATION flags);
94        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
95    private static extern void dfftw_plan_dft_r2c_3d_(ref IntPtr plan, ref int nx, ref int ny, ref int nz, [In,Out] double* input, [In,Out] complex* output, ref FFTW_PLANCREATION flags);
96        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
97    private static extern void dfftw_plan_dft_r2c_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] double* input, [In,Out] complex* output, ref FFTW_PLANCREATION flags);
98
99        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
100    private static extern void dfftw_plan_dft_c2r_1d_(ref IntPtr plan, ref int n, [In,Out] complex* input, [In,Out] double* output, ref FFTW_PLANCREATION flags);
101        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
102    private static extern void dfftw_plan_dft_c2r_2d_(ref IntPtr plan, ref int nx, ref int ny, [In,Out] complex* input, [In,Out] double* output, ref FFTW_PLANCREATION flags);
103        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
104    private static extern void dfftw_plan_dft_c2r_3d_(ref IntPtr plan, ref int nx, ref int ny, ref int nz, [In,Out] complex* input, [In,Out] double* output, ref FFTW_PLANCREATION flags);
105        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
106    private static extern void dfftw_plan_dft_c2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] complex* input, [In,Out] double* output, ref FFTW_PLANCREATION flags);
107
108        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
109    private static extern void dfftw_plan_r2r_1d_(ref IntPtr plan, ref int n, [In,Out] double* input, [In,Out] double* output,  ref FFTW_PLANCREATION flags);
110        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
111    private static extern void dfftw_plan_r2r_2d_(ref IntPtr plan, ref int nx, ref int ny, [In,Out] double* input, [In,Out] double* output, ref FFTW_PLANCREATION flags);
112        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
113    private static extern void dfftw_plan_r2r_3d_(ref IntPtr plan, ref int nx, ref int ny, ref int nz, [In,Out] double* input, [In,Out] double* output, ref FFTW_PLANCREATION flags);
114        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
115    private static extern void dfftw_plan_r2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] double* input, [In,Out] double* output, ref FFTW_PLANCREATION flags);
116
117        /************************************ ADVANCED INTERFACE *************************************/
118        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
119        private static extern void dfftw_plan_many_dft_(ref IntPtr plan, ref int rank, [In,Out] int[] n, ref int howmany, [In,Out] complex* input, [In,Out] int[] inembed, ref int istride, ref int idist, [In,Out] complex* output, [In,Out] int[] onembed, ref int ostride, ref int odist, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);   
120        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
121        private static extern void dfftw_plan_many_dft_r2c_(ref IntPtr plan, ref int rank, [In,Out] int [] n, ref int howmany, [In,Out] double* input, [In,Out] int[] inembed, ref int istride, ref int idist, [In,Out] complex* output, [In,Out] int[] onembed, ref int ostride, ref int odist, ref FFTW_PLANCREATION flags);       
122        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
123        private static extern void dfftw_plan_many_dft_c2r_(ref IntPtr plan, ref int rank, [In,Out] int [] n, ref int howmany, [In,Out] complex* input, [In,Out] int[] inembed, ref int istride, ref int idist, [In,Out] double* output, [In,Out] int[] onembed, ref int ostride, ref int odist, ref FFTW_PLANCREATION flags);       
124        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
125        private static extern void dfftw_plan_many_r2r_(    ref IntPtr plan, ref int rank, [In,Out] int [] n, ref int howmany, [In,Out] double [] input, [In,Out] int[] inembed, ref int istride, ref int idist, [In,Out] double* output, [In,Out] int[] onembed, ref int ostride, ref int odist, ref FFTW_KIND kind, ref FFTW_PLANCREATION flags);       
126       
127        /************************************ GURU INTERFACE *************************************/
128        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
129    private static extern void dfftw_plan_guru_dft_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] complex* input, [In,Out] complex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
130        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
131        private static extern void dfftw_plan_guru_split_dft_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] double* ri, [In,Out] double* ii, [In,Out] double* ro, [In,Out] double* io, ref FFTW_PLANCREATION flags);
132        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
133    private static extern void dfftw_plan_guru_dft_r2c_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] double* input, [In,Out] complex* output, ref FFTW_PLANCREATION flags);
134        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
135    private static extern void dfftw_plan_guru_dft_c2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] complex* input, [In,Out] double* output, ref FFTW_PLANCREATION flags);
136        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
137        private static extern void dfftw_plan_guru_split_dft_r2c_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] double* ri, [In,Out] double* ro, [In,Out] double* io, ref FFTW_PLANCREATION flags);
138        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
139        private static extern void dfftw_plan_guru_split_dft_c2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] double* ri, [In,Out] double* ii, [In,Out] double* ro, ref FFTW_PLANCREATION flags);
140        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
141    private static extern void dfftw_plan_guru_r2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] double* input, [In,Out] double* output, ref FFTW_KIND kind, ref FFTW_PLANCREATION flags);
142        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
143    private static extern void dfftw_execute_dft_(ref IntPtr plan, [In,Out] complex* input, [In,Out] complex* output);
144        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
145    private static extern void dfftw_execute_split_dft_(ref IntPtr plan, [In,Out] double* ri, [In,Out] double* ii, [In,Out] double* ro, [In,Out] double* io);
146        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
147    private static extern void dfftw_execute_dft_r2c_(ref IntPtr plan, [In,Out] double* ri, [In,Out] complex* ii);
148        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
149    private static extern void dfftw_execute_split_dft_r2c_(ref IntPtr plan, [In,Out] double* input, [In,Out] double* ro, [In,Out] double* ri);
150        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
151    private static extern void dfftw_execute_dft_c2r_(ref IntPtr plan, [In,Out] complex* input, [In,Out] double* output);
152        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
153    private static extern void dfftw_execute_split_dft_c2r_(ref IntPtr plan, [In,Out] double* ri, [In,Out] double* ii, [In,Out] double* output);
154        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
155    private static extern void dfftw_execute_r2r_(ref IntPtr plan, [In,Out] double* input, [In,Out] double* output);
156
157        /************************************ BASIC PLAN EXECUTION & MISC ******************************/
158        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
159    private static extern void dfftw_execute_(ref IntPtr plan);
160        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
161    private static extern void dfftw_destroy_plan_(ref IntPtr plan);
162        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
163    private static extern void dfftw_cleanup_();
164        [DllImport("libfftw3-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
165    private static extern void fftw_set_timelimit(ref double seconds);
166        #endregion DOUBLE PRECISION INTERFACE
167
168        #region SINGLE PRECISION INTERFACE
169        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
170    private static extern void sfftw_plan_dft_1d_(ref IntPtr plan, ref int n, [In,Out] fcomplex* input, [In,Out] fcomplex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
171        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
172    private static extern void sfftw_plan_dft_2d_(ref IntPtr plan, ref int nx, ref int ny, [In,Out] fcomplex* input, [In,Out] fcomplex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
173        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
174    private static extern void sfftw_plan_dft_3d_(ref IntPtr plan, ref int nx, ref int ny, ref int nz, [In,Out] fcomplex* input, [In,Out] fcomplex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
175        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
176    private static extern void sfftw_plan_dft_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] fcomplex* input, [In,Out] fcomplex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
177       
178        /****************************************************************** REAL-COMPLEX ***/
179        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
180    private static extern void sfftw_plan_dft_r2c_1d_(ref IntPtr plan, ref int n, [In,Out] float* input, [In,Out] fcomplex* output, ref FFTW_PLANCREATION flags);
181        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
182    private static extern void sfftw_plan_dft_r2c_2d_(ref IntPtr plan, ref int nx, ref int ny, [In,Out] float* input, [In,Out] fcomplex* output, ref FFTW_PLANCREATION flags);
183        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
184    private static extern void sfftw_plan_dft_r2c_3d_(ref IntPtr plan, ref int nx, ref int ny, ref int nz, [In,Out] float* input, [In,Out] fcomplex* output, ref FFTW_PLANCREATION flags);
185        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
186    private static extern void sfftw_plan_dft_r2c_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] float* input, [In,Out] fcomplex* output, ref FFTW_PLANCREATION flags);
187
188        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
189    private static extern void sfftw_plan_dft_c2r_1d_(ref IntPtr plan, ref int n, [In,Out] fcomplex* input, [In,Out] float* output, ref FFTW_PLANCREATION flags);
190        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
191    private static extern void sfftw_plan_dft_c2r_2d_(ref IntPtr plan, ref int nx, ref int ny, [In,Out] fcomplex* input, [In,Out] float* output, ref FFTW_PLANCREATION flags);
192        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
193    private static extern void sfftw_plan_dft_c2r_3d_(ref IntPtr plan, ref int nx, ref int ny, ref int nz, [In,Out] fcomplex* input, [In,Out] float* output, ref FFTW_PLANCREATION flags);
194        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
195    private static extern void sfftw_plan_dft_c2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] fcomplex* input, [In,Out] float* output, ref FFTW_PLANCREATION flags);
196
197        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
198    private static extern void sfftw_plan_r2r_1d_(ref IntPtr plan, ref int n, [In,Out] float* input, [In,Out] float* output,  ref FFTW_PLANCREATION flags);
199        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
200    private static extern void sfftw_plan_r2r_2d_(ref IntPtr plan, ref int nx, ref int ny, [In,Out] float* input, [In,Out] float* output, ref FFTW_PLANCREATION flags);
201        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
202    private static extern void sfftw_plan_r2r_3d_(ref IntPtr plan, ref int nx, ref int ny, ref int nz, [In,Out] float* input, [In,Out] float* output, ref FFTW_PLANCREATION flags);
203        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
204    private static extern void sfftw_plan_r2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] float* input, [In,Out] float* output, ref FFTW_PLANCREATION flags);
205
206        /************************************ ADVANCED INTERFACE *************************************/
207        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
208        private static extern void sfftw_plan_many_dft_(ref IntPtr plan, ref int rank, [In,Out] int[] n, ref int howmany, [In,Out] fcomplex* input, [In,Out] int[] inembed, ref int istride, ref int idist, [In,Out] fcomplex* output, [In,Out] int[] onembed, ref int ostride, ref int odist, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);   
209        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
210        private static extern void sfftw_plan_many_dft_r2c_(ref IntPtr plan, ref int rank, [In,Out] int [] n, ref int howmany, [In,Out] float* input, [In,Out] int[] inembed, ref int istride, ref int idist, [In,Out] fcomplex* output, [In,Out] int[] onembed, ref int ostride, ref int odist, ref FFTW_PLANCREATION flags);       
211        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
212        private static extern void sfftw_plan_many_dft_c2r_(ref IntPtr plan, ref int rank, [In,Out] int [] n, ref int howmany, [In,Out] fcomplex* input, [In,Out] int[] inembed, ref int istride, ref int idist, [In,Out] float* output, [In,Out] int[] onembed, ref int ostride, ref int odist, ref FFTW_PLANCREATION flags);       
213        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
214        private static extern void sfftw_plan_many_r2r_(    ref IntPtr plan, ref int rank, [In,Out] int [] n, ref int howmany, [In,Out] float [] input, [In,Out] int[] inembed, ref int istride, ref int idist, [In,Out] float* output, [In,Out] int[] onembed, ref int ostride, ref int odist, ref FFTW_KIND kind, ref FFTW_PLANCREATION flags);       
215       
216        /************************************ GURU INTERFACE *************************************/
217        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
218    private static extern void sfftw_plan_guru_dft_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] fcomplex* input, [In,Out] fcomplex* output, ref FFTW_DIRECTION sign, ref FFTW_PLANCREATION flags);
219        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
220        private static extern void sfftw_plan_guru_split_dft_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] float* ri, [In,Out] float* ii, [In,Out] float* ro, [In,Out] float* io, ref FFTW_PLANCREATION flags);
221        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
222    private static extern void sfftw_plan_guru_dft_r2c_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] float* input, [In,Out] fcomplex* output, ref FFTW_PLANCREATION flags);
223        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
224    private static extern void sfftw_plan_guru_dft_c2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] fcomplex* input, [In,Out] float* output, ref FFTW_PLANCREATION flags);
225        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
226        private static extern void sfftw_plan_guru_split_dft_r2c_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] float* ri, [In,Out] float* ro, [In,Out] float* io, ref FFTW_PLANCREATION flags);
227        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
228        private static extern void sfftw_plan_guru_split_dft_c2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] float* ri, [In,Out] float* ii, [In,Out] float* ro, ref FFTW_PLANCREATION flags);
229        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
230    private static extern void sfftw_plan_guru_r2r_(ref IntPtr plan, ref int rank, [In,Out] int[] n, [In,Out] int[] iS,  [In,Out] int[] iO, ref int howmany_rank, [In,Out] int[] hn, [In,Out] int[] hiS,  [In,Out] int[] hiO, [In,Out] float* input, [In,Out] float* output, ref FFTW_KIND kind, ref FFTW_PLANCREATION flags);
231        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
232    private static extern void sfftw_execute_dft_(ref IntPtr plan, [In,Out] fcomplex* input, [In,Out] fcomplex* output);
233        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
234    private static extern void sfftw_execute_split_dft_(ref IntPtr plan, [In,Out] float* ri, [In,Out] float* ii, [In,Out] float* ro, [In,Out] float* io);
235        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
236    private static extern void sfftw_execute_dft_r2c_(ref IntPtr plan, [In,Out] float* ri, [In,Out] fcomplex* ii);
237        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
238    private static extern void sfftw_execute_split_dft_r2c_(ref IntPtr plan, [In,Out] float* input, [In,Out] float* ro, [In,Out] float* ri);
239        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
240    private static extern void sfftw_execute_dft_c2r_(ref IntPtr plan, [In,Out] fcomplex* input, [In,Out] float* output);
241        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
242    private static extern void sfftw_execute_split_dft_c2r_(ref IntPtr plan, [In,Out] float* ri, [In,Out] float* ii, [In,Out] float* output);
243        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
244    private static extern void sfftw_execute_r2r_(ref IntPtr plan, [In,Out] float* input, [In,Out] float* output);
245
246        /************************************ BASIC PLAN EXECUTION & MISC ******************************/
247        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
248    private static extern void sfftw_execute_(ref IntPtr plan);
249        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
250    private static extern void sfftw_destroy_plan_(ref IntPtr plan);
251        [DllImport("libfftw3f-3",CallingConvention =CallingConvention.Cdecl),SuppressUnmanagedCodeSecurity]
252    private static extern void sfftw_cleanup_();
253        #endregion SINGLE PRECISION INTERFACE
254        #endregion
255
256        #region FFTW CONSTANTS
257
258        [FlagsAttribute]
259        public enum FFTW_PLANCREATION : uint {
260            MEASURE = 0u,
261            DESTROY_INPUT = (1U << 0),
262            UNALIGNED = (1U << 1),
263            CONSERVE_MEMORY = (1U << 2),
264            EXHAUSTIVE = (1U << 3), /* NO_EXHAUSTIVE is default */
265            PRESERVE_INPUT = (1U << 4), /* cancels FFTW_DESTROY_INPUT */
266            PATIENT = (1U << 5), /* IMPATIENT is default */
267            ESTIMATE = (1U << 6)
268        }
269       
270        public enum FFTW_DIRECTION {
271            FORWARD = -1,
272            BACKWARDS = 1
273        }
274
275        public enum FFTW_KIND {
276            R2HC=0,
277            HC2R=1,
278            DHT=2,
279            REDFT00=3,
280            REDFT01=4,
281            REDFT10=5,
282            REDFT11=6,
283            RODFT00=7,
284            RODFT01=8,
285            RODFT10=9,
286            RODFT11=10
287        }
288        #endregion
289
290        #region attributes
291        object _lockobject = new object();
292        int[] m_n,m_is,m_os, m_hn, m_his, m_hos;
293        //List<int> m_tmpDims = new List<int>(10);
294        #endregion
295
296        #region constructor
297        public ILFFTW3FFT () {
298            // fftw internally caches the plans, so ....
299            m_n = new int[10];
300            m_is = new int[10];
301            m_os = new int[10];
302            m_hn = new int[2];
303            m_his = new int[2];
304            m_hos = new int[2];
305        }
306        #endregion
307
308        #region IILFFT Member - 1-D
309
310
311
312        public ILRetArray< complex >  FFTForward1D (ILInArray< double > A, int dim) {
313            return  FFTForward (A,dim,1);
314        }
315
316#region HYCALPER AUTO GENERATED CODE
317
318
319        public ILRetArray< fcomplex >  FFTBackward1D (ILInArray< fcomplex > A, int dim) {
320            return  FFTBackward (A,dim,1);
321        }
322
323        public ILRetArray< complex >  FFTBackward1D (ILInArray< complex > A, int dim) {
324            return  FFTBackward (A,dim,1);
325        }
326
327        public ILRetArray< fcomplex >  FFTForward1D (ILInArray< fcomplex > A, int dim) {
328            return  FFTForward (A,dim,1);
329        }
330
331        public ILRetArray< fcomplex >  FFTForward1D (ILInArray< float > A, int dim) {
332            return  FFTForward (A,dim,1);
333        }
334
335        public ILRetArray< complex >  FFTForward1D (ILInArray< complex > A, int dim) {
336            return  FFTForward (A,dim,1);
337        }
338
339#endregion HYCALPER AUTO GENERATED CODE
340
341
342
343        public ILRetArray< double > FFTBackwSym1D(ILInArray< complex > A, int dim) {
344            return ILMath.real(FFTBackward1D(A,dim));
345        }
346
347#region HYCALPER AUTO GENERATED CODE
348
349
350        public ILRetArray< float > FFTBackwSym1D(ILInArray< fcomplex > A, int dim) {
351            return ILMath.real(FFTBackward1D(A,dim));
352        }
353
354#endregion HYCALPER AUTO GENERATED CODE
355
356        #endregion
357
358        #region IILFFT Member - n-D
359
360
361
362        public ILRetArray< complex >  FFTForward (ILInArray< double > A, int nDims) {
363            return  FFTForward (A,0,nDims);
364        }
365
366        public ILRetArray< complex>  FFTForward(ILInArray< double> A, int firstDim, int nDims) {
367            using (ILScope.Enter(A)) {
368                if (object.Equals(A, null) || nDims < 0 || firstDim < 0)
369                    throw new ILArgumentException("invalid parameter!");
370                if (A.IsEmpty) return ILArray<  complex>.empty(A.Size);
371                if (A.IsScalar || (A.Size[firstDim] == 1 && nDims == 1))
372                    return  ILMath.tocomplex(A);
373                if (nDims > A.Size.NumberOfDimensions)
374                    nDims = A.Size.NumberOfDimensions;
375                // prepare output array
376                ILArray< complex> ret =  ILMath.tocomplex(A);
377                IntPtr plan = IntPtr.Zero;
378                int hrank = 2;
379                FFTW_DIRECTION dir =  FFTW_DIRECTION.FORWARD;
380                lock (_lockobject) {
381                    FFTW_PLANCREATION flags = FFTW_PLANCREATION.ESTIMATE;
382                    if (nDims > m_n.Length) resizeCache(nDims);
383                    for (int i = 0; i < nDims; i++) {
384                        m_n[i] = A.Size[i + firstDim];
385                        m_is[i] = A.Size.SequentialIndexDistance(i + firstDim);
386                        m_os[i] = m_is[i];
387                    }
388                    m_hn[0] = A.Size.SequentialIndexDistance(firstDim);
389                    m_his[0] = 1;
390                    m_hos[0] = 1;
391                    m_his[1] = A.Size.SequentialIndexDistance(nDims + firstDim);
392                    m_hn[1] = A.Size.NumberOfElements / m_his[1];
393                    m_hos[1] = m_his[1];
394                    fixed ( complex* retArr = ret.GetArrayForWrite()) {
395                       
396                        dfftw_plan_guru_dft_(ref plan, ref nDims, m_n, m_is, m_os,
397                                      ref hrank, m_hn, m_his, m_hos,
398                                      retArr, retArr,
399                                      ref dir, ref flags);
400                    }
401                }
402                if (plan == IntPtr.Zero)
403                    throw new ILInvalidOperationException("error creating plan for fftw3 (guru interface)");
404                dfftw_execute_(ref plan);
405                dfftw_destroy_plan_(ref plan);
406               
407
408                return ret;
409            }
410        }
411
412#region HYCALPER AUTO GENERATED CODE
413
414
415        public ILRetArray< fcomplex >  FFTBackward (ILInArray< fcomplex > A, int nDims) {
416            return  FFTBackward (A,0,nDims);
417        }
418
419        public ILRetArray< fcomplex>  FFTBackward(ILInArray< fcomplex> A, int firstDim, int nDims) {
420            using (ILScope.Enter(A)) {
421                if (object.Equals(A, null) || nDims < 0 || firstDim < 0)
422                    throw new ILArgumentException("invalid parameter!");
423                if (A.IsEmpty) return ILArray<  fcomplex>.empty(A.Size);
424                if (A.IsScalar || (A.Size[firstDim] == 1 && nDims == 1))
425                    return  A.C;
426                if (nDims > A.Size.NumberOfDimensions)
427                    nDims = A.Size.NumberOfDimensions;
428                // prepare output array
429                ILArray< fcomplex> ret =  A.C;
430                IntPtr plan = IntPtr.Zero;
431                int hrank = 2;
432                FFTW_DIRECTION dir =  FFTW_DIRECTION.BACKWARDS;
433                lock (_lockobject) {
434                    FFTW_PLANCREATION flags = FFTW_PLANCREATION.ESTIMATE;
435                    if (nDims > m_n.Length) resizeCache(nDims);
436                    for (int i = 0; i < nDims; i++) {
437                        m_n[i] = A.Size[i + firstDim];
438                        m_is[i] = A.Size.SequentialIndexDistance(i + firstDim);
439                        m_os[i] = m_is[i];
440                    }
441                    m_hn[0] = A.Size.SequentialIndexDistance(firstDim);
442                    m_his[0] = 1;
443                    m_hos[0] = 1;
444                    m_his[1] = A.Size.SequentialIndexDistance(nDims + firstDim);
445                    m_hn[1] = A.Size.NumberOfElements / m_his[1];
446                    m_hos[1] = m_his[1];
447                    fixed ( fcomplex* retArr = ret.GetArrayForWrite()) {
448                       
449                        sfftw_plan_guru_dft_(ref plan, ref nDims, m_n, m_is, m_os,
450                                      ref hrank, m_hn, m_his, m_hos,
451                                      retArr, retArr,
452                                      ref dir, ref flags);
453                    }
454                }
455                if (plan == IntPtr.Zero)
456                    throw new ILInvalidOperationException("error creating plan for fftw3 (guru interface)");
457                dfftw_execute_(ref plan);
458                dfftw_destroy_plan_(ref plan);
459                ret /= new fcomplex(m_his[1]/m_hn[0],0.0f);
460                return ret;
461            }
462        }
463
464        public ILRetArray< complex >  FFTBackward (ILInArray< complex > A, int nDims) {
465            return  FFTBackward (A,0,nDims);
466        }
467
468        public ILRetArray< complex>  FFTBackward(ILInArray< complex> A, int firstDim, int nDims) {
469            using (ILScope.Enter(A)) {
470                if (object.Equals(A, null) || nDims < 0 || firstDim < 0)
471                    throw new ILArgumentException("invalid parameter!");
472                if (A.IsEmpty) return ILArray<  complex>.empty(A.Size);
473                if (A.IsScalar || (A.Size[firstDim] == 1 && nDims == 1))
474                    return  A.C;
475                if (nDims > A.Size.NumberOfDimensions)
476                    nDims = A.Size.NumberOfDimensions;
477                // prepare output array
478                ILArray< complex> ret =  A.C;
479                IntPtr plan = IntPtr.Zero;
480                int hrank = 2;
481                FFTW_DIRECTION dir =  FFTW_DIRECTION.BACKWARDS;
482                lock (_lockobject) {
483                    FFTW_PLANCREATION flags = FFTW_PLANCREATION.ESTIMATE;
484                    if (nDims > m_n.Length) resizeCache(nDims);
485                    for (int i = 0; i < nDims; i++) {
486                        m_n[i] = A.Size[i + firstDim];
487                        m_is[i] = A.Size.SequentialIndexDistance(i + firstDim);
488                        m_os[i] = m_is[i];
489                    }
490                    m_hn[0] = A.Size.SequentialIndexDistance(firstDim);
491                    m_his[0] = 1;
492                    m_hos[0] = 1;
493                    m_his[1] = A.Size.SequentialIndexDistance(nDims + firstDim);
494                    m_hn[1] = A.Size.NumberOfElements / m_his[1];
495                    m_hos[1] = m_his[1];
496                    fixed ( complex* retArr = ret.GetArrayForWrite()) {
497                       
498                        dfftw_plan_guru_dft_(ref plan, ref nDims, m_n, m_is, m_os,
499                                      ref hrank, m_hn, m_his, m_hos,
500                                      retArr, retArr,
501                                      ref dir, ref flags);
502                    }
503                }
504                if (plan == IntPtr.Zero)
505                    throw new ILInvalidOperationException("error creating plan for fftw3 (guru interface)");
506                dfftw_execute_(ref plan);
507                dfftw_destroy_plan_(ref plan);
508                ret /= new complex(m_his[1]/m_hn[0],0.0f);
509                return ret;
510            }
511        }
512
513        public ILRetArray< fcomplex >  FFTForward (ILInArray< fcomplex > A, int nDims) {
514            return  FFTForward (A,0,nDims);
515        }
516
517        public ILRetArray< fcomplex>  FFTForward(ILInArray< fcomplex> A, int firstDim, int nDims) {
518            using (ILScope.Enter(A)) {
519                if (object.Equals(A, null) || nDims < 0 || firstDim < 0)
520                    throw new ILArgumentException("invalid parameter!");
521                if (A.IsEmpty) return ILArray<  fcomplex>.empty(A.Size);
522                if (A.IsScalar || (A.Size[firstDim] == 1 && nDims == 1))
523                    return  A.C;
524                if (nDims > A.Size.NumberOfDimensions)
525                    nDims = A.Size.NumberOfDimensions;
526                // prepare output array
527                ILArray< fcomplex> ret =  A.C;
528                IntPtr plan = IntPtr.Zero;
529                int hrank = 2;
530                FFTW_DIRECTION dir =  FFTW_DIRECTION.FORWARD;
531                lock (_lockobject) {
532                    FFTW_PLANCREATION flags = FFTW_PLANCREATION.ESTIMATE;
533                    if (nDims > m_n.Length) resizeCache(nDims);
534                    for (int i = 0; i < nDims; i++) {
535                        m_n[i] = A.Size[i + firstDim];
536                        m_is[i] = A.Size.SequentialIndexDistance(i + firstDim);
537                        m_os[i] = m_is[i];
538                    }
539                    m_hn[0] = A.Size.SequentialIndexDistance(firstDim);
540                    m_his[0] = 1;
541                    m_hos[0] = 1;
542                    m_his[1] = A.Size.SequentialIndexDistance(nDims + firstDim);
543                    m_hn[1] = A.Size.NumberOfElements / m_his[1];
544                    m_hos[1] = m_his[1];
545                    fixed ( fcomplex* retArr = ret.GetArrayForWrite()) {
546                       
547                        sfftw_plan_guru_dft_(ref plan, ref nDims, m_n, m_is, m_os,
548                                      ref hrank, m_hn, m_his, m_hos,
549                                      retArr, retArr,
550                                      ref dir, ref flags);
551                    }
552                }
553                if (plan == IntPtr.Zero)
554                    throw new ILInvalidOperationException("error creating plan for fftw3 (guru interface)");
555                dfftw_execute_(ref plan);
556                dfftw_destroy_plan_(ref plan);
557               
558                return ret;
559            }
560        }
561
562        public ILRetArray< fcomplex >  FFTForward (ILInArray< float > A, int nDims) {
563            return  FFTForward (A,0,nDims);
564        }
565
566        public ILRetArray< fcomplex>  FFTForward(ILInArray< float> A, int firstDim, int nDims) {
567            using (ILScope.Enter(A)) {
568                if (object.Equals(A, null) || nDims < 0 || firstDim < 0)
569                    throw new ILArgumentException("invalid parameter!");
570                if (A.IsEmpty) return ILArray<  fcomplex>.empty(A.Size);
571                if (A.IsScalar || (A.Size[firstDim] == 1 && nDims == 1))
572                    return  ILMath.tofcomplex(A);
573                if (nDims > A.Size.NumberOfDimensions)
574                    nDims = A.Size.NumberOfDimensions;
575                // prepare output array
576                ILArray< fcomplex> ret =  ILMath.tofcomplex(A);
577                IntPtr plan = IntPtr.Zero;
578                int hrank = 2;
579                FFTW_DIRECTION dir =  FFTW_DIRECTION.FORWARD;
580                lock (_lockobject) {
581                    FFTW_PLANCREATION flags = FFTW_PLANCREATION.ESTIMATE;
582                    if (nDims > m_n.Length) resizeCache(nDims);
583                    for (int i = 0; i < nDims; i++) {
584                        m_n[i] = A.Size[i + firstDim];
585                        m_is[i] = A.Size.SequentialIndexDistance(i + firstDim);
586                        m_os[i] = m_is[i];
587                    }
588                    m_hn[0] = A.Size.SequentialIndexDistance(firstDim);
589                    m_his[0] = 1;
590                    m_hos[0] = 1;
591                    m_his[1] = A.Size.SequentialIndexDistance(nDims + firstDim);
592                    m_hn[1] = A.Size.NumberOfElements / m_his[1];
593                    m_hos[1] = m_his[1];
594                    fixed ( fcomplex* retArr = ret.GetArrayForWrite()) {
595                       
596                        sfftw_plan_guru_dft_(ref plan, ref nDims, m_n, m_is, m_os,
597                                      ref hrank, m_hn, m_his, m_hos,
598                                      retArr, retArr,
599                                      ref dir, ref flags);
600                    }
601                }
602                if (plan == IntPtr.Zero)
603                    throw new ILInvalidOperationException("error creating plan for fftw3 (guru interface)");
604                dfftw_execute_(ref plan);
605                dfftw_destroy_plan_(ref plan);
606               
607                return ret;
608            }
609        }
610
611        public ILRetArray< complex >  FFTForward (ILInArray< complex > A, int nDims) {
612            return  FFTForward (A,0,nDims);
613        }
614
615        public ILRetArray< complex>  FFTForward(ILInArray< complex> A, int firstDim, int nDims) {
616            using (ILScope.Enter(A)) {
617                if (object.Equals(A, null) || nDims < 0 || firstDim < 0)
618                    throw new ILArgumentException("invalid parameter!");
619                if (A.IsEmpty) return ILArray<  complex>.empty(A.Size);
620                if (A.IsScalar || (A.Size[firstDim] == 1 && nDims == 1))
621                    return  A.C;
622                if (nDims > A.Size.NumberOfDimensions)
623                    nDims = A.Size.NumberOfDimensions;
624                // prepare output array
625                ILArray< complex> ret =  A.C;
626                IntPtr plan = IntPtr.Zero;
627                int hrank = 2;
628                FFTW_DIRECTION dir =  FFTW_DIRECTION.FORWARD;
629                lock (_lockobject) {
630                    FFTW_PLANCREATION flags = FFTW_PLANCREATION.ESTIMATE;
631                    if (nDims > m_n.Length) resizeCache(nDims);
632                    for (int i = 0; i < nDims; i++) {
633                        m_n[i] = A.Size[i + firstDim];
634                        m_is[i] = A.Size.SequentialIndexDistance(i + firstDim);
635                        m_os[i] = m_is[i];
636                    }
637                    m_hn[0] = A.Size.SequentialIndexDistance(firstDim);
638                    m_his[0] = 1;
639                    m_hos[0] = 1;
640                    m_his[1] = A.Size.SequentialIndexDistance(nDims + firstDim);
641                    m_hn[1] = A.Size.NumberOfElements / m_his[1];
642                    m_hos[1] = m_his[1];
643                    fixed ( complex* retArr = ret.GetArrayForWrite()) {
644                       
645                        dfftw_plan_guru_dft_(ref plan, ref nDims, m_n, m_is, m_os,
646                                      ref hrank, m_hn, m_his, m_hos,
647                                      retArr, retArr,
648                                      ref dir, ref flags);
649                    }
650                }
651                if (plan == IntPtr.Zero)
652                    throw new ILInvalidOperationException("error creating plan for fftw3 (guru interface)");
653                dfftw_execute_(ref plan);
654                dfftw_destroy_plan_(ref plan);
655               
656                return ret;
657            }
658        }
659
660#endregion HYCALPER AUTO GENERATED CODE
661       
662
663       
664        public ILRetArray< double > FFTBackwSym(ILInArray< complex > A, int nDims) {
665            return ILMath.real(FFTBackward(A,nDims)); 
666        }
667
668#region HYCALPER AUTO GENERATED CODE
669
670       
671        public ILRetArray< float > FFTBackwSym(ILInArray< fcomplex > A, int nDims) {
672            return ILMath.real(FFTBackward(A,nDims)); 
673        }
674
675#endregion HYCALPER AUTO GENERATED CODE
676
677        #endregion
678       
679        #region IILFFT Member - Misc
680       
681        public bool CachePlans {
682            get {
683                return true;
684            }
685        }
686
687        public void FreePlans() {
688            FreeAllDescriptors();
689        }
690
691        public bool SpeedyHermitian {
692            get { return false; }
693        }
694
695        #endregion
696
697        #region private helper
698        private void FreeAllDescriptors() {
699            try {
700                dfftw_cleanup_();
701                sfftw_cleanup_();
702            } catch (Exception) { }
703        }
704        private void resizeCache(int size) {
705            if (m_n.Length >= size) return;
706            m_n = new int[size];
707            m_is = new int[size];
708            m_os = new int[size];
709        }
710        #endregion
711
712        #region IDisposable Member
713
714        public void Dispose() {
715            Dispose(true);
716            GC.SuppressFinalize(this);
717        }
718
719        public virtual void Dispose (bool manual) {
720            if (manual) {
721                FreeAllDescriptors();             
722            }
723        }
724
725        ~ILFFTW3FFT () {
726            Dispose(true);     
727        }
728        #endregion
729
730    }
731}
Note: See TracBrowser for help on using the repository browser.