-
-int fourier_calc(t_fourier *fourier) {
-
- dprintf(fourier->outfd,"[fourier] %d dimensional %c-%cft calculation ...\n",
- fourier->dim,
- ((fourier->type&FWD)?'f':'b'),
- ((fourier->type&DFT)?'d':'f'));
-
- if(fourier->type&DFT) {
- switch(fourier->dim) {
- case 1:
- fourier_dft_1d(fourier);
- break;
- case 2:
- fourier_dft_2d(fourier);
- break;
- case 3:
- fourier_dft_3d(fourier);
- break;
- default:
- dprintf(fourier->outfd,"[fourier] dimension failure\n");
- return F_DIM_FAILURE;
- }
- return F_SUCCESS;
- }
- else {
- dprintf(fourier->outfd,"[fourier] fft not supported by now\n");
- return F_NOT_SUPPORTED;
- }
-}