summaryrefslogtreecommitdiff
path: root/post/src/ft2d.h
blob: e66504143147342fa7f2f8b422754e4124aecd8b (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
/*
   This file is part of LPIC++, a particle-in-cell code for
   simulating the interaction of laser light with plasma.

   Copyright (C) 1994-1997 Roland Lichters

   LPIC++ is free software; you can redistribute it and/or
   modify it under the terms of the GNU General Public License
   as published by the Free Software Foundation; either version 2
   of the License, or (at your option) any later version.

   This program is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.

   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the Free Software
   Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
*/

#ifndef FT2_H
#define FT2_H

#include <common.h>
#include <utilities.h>
#include <math.h>
#include <fstream>
#include <iomanip>
#include <string.h>

class FFT2D {
 public:

// input
  int   periods_input_1, periods_input_2;
  int   steps_pp_input_1, steps_pp_input_2;
  int   screen_input_1, screen_input_2;

// internal
  int   steps_1, steps_2, steps_half_1, steps_half_2, steps_input_1, steps_input_2;
  float dt_1, dt_input_1, dt_2, dt_input_2;
  float df_1, df_2;
  float **local;
  float *data;
  float *frequency_1, *frequency_2;
  float **si;
  float **co;
  float **power;
  int   *nn;

  char errname[filename_size];

      FFT2D            ( int Q_kw, int periods_1, int steps_pp_1, int screen_1,
			           int periods_2, int steps_pp_2, int screen_2 );
  int   steps_ft       ( int steps_input );
  float window         ( float t, int screen_input, int periods_input );
  void  fftn           ( float* data, int *nn, int ndim, int isign );
  void  RealFt         ( float** input );
};

#endif