summaryrefslogtreecommitdiff
path: root/lpic/src/diagnostic_phasespace.h
blob: ef50cc5402dd047f0639120583b77574b05307d9 (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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
/*
   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 DIAGNOSTIC_PHASESPACE_H
#define DIAGNOSTIC_PHASESPACE_H

#include <diagnostic_stepper.h>
#include <common.h>
#include <error.h>
#include <parameter.h>
#include <domain.h>
#include <matrix.h>
#include <math.h>
#include <readfile.h>


class input_phasespace {
private:
  char          errname[filename_size];
  readfile      rf;
  void          save( parameter &p );
  char          species_name[filename_size];

public:
  stepper_param stepper;
  int           n_domains;
  int           cells, cells_per_wl;
  int           Q_restart;
  char          restart_file[filename_size];

  input_phasespace( parameter &p, char *species_name_input );
};


//////////////////////////////////////////////////////////////////////////////////////////


class phasespace {
private:
  char             errname[filename_size];
  readfile         rf;
  input_phasespace input;
  int              species;
  char             species_name[filename_size];

public:
  diagnostic_stepper stepper;

  int    dim;
  int    box_cells;
  double box_length;
  double Beta, Gamma;
  double vcut;
  unsigned char **x, **y, **z;
  double vx, vy, vz;
  int bx, bvx, bvy, bvz;
  char *name;
  FILE *file_x, *file_y, *file_z;

  phasespace           ( parameter &p,
			 int species_input, char *species_name_input );
  void write_phasespace( double time, parameter &p, domain *grid );
};


class el_phasespace : public phasespace {
public:

  el_phasespace( parameter &p,
		 int species_input=0, char *species_name_input="el" );
};


class ion_phasespace : public phasespace {
public:

  ion_phasespace( parameter &p,
		  int species_input=1, char *species_name_input="ion" );
};

//////////////////////////////////////////////////////////////////////////////////////////

#endif