OPAL (Object Oriented Parallel Accelerator Library)  2024.1
OPAL
Astra1D_fast.h
Go to the documentation of this file.
1 #ifndef CLASSIC_AstraFIELDMAP1DFAST_HH
2 #define CLASSIC_AstraFIELDMAP1DFAST_HH
3 
4 #include "Fields/Fieldmap.h"
5 
6 class _Astra1D_fast: public _Fieldmap {
7 
8 public:
9  virtual ~_Astra1D_fast();
10 
11  virtual bool getFieldstrength(const Vector_t &R, Vector_t &E, Vector_t &B) const = 0;
12  virtual bool getFieldDerivative(const Vector_t &R, Vector_t &E, Vector_t &B, const DiffDirection &dir) const = 0;
13  virtual void getFieldDimensions(double &zBegin, double &zEnd) const = 0;
14  virtual void getFieldDimensions(double &xIni, double &xFinal, double &yIni, double &yFinal, double &zIni, double &zFinal) const = 0;
15  virtual void swap() = 0;
16  virtual void getInfo(Inform *) = 0;
17  virtual double getFrequency() const = 0;
18  virtual void setFrequency(double freq) = 0;
19  virtual void getOnaxisEz(std::vector<std::pair<double, double> > & F);
20 
21  virtual bool isInside(const Vector_t &r) const;
22 protected:
23  _Astra1D_fast(const std::string& filename);
24 
25  virtual void readMap();
26  virtual void freeMap();
27 
28  bool determineNumSamplingPoints(std::ifstream &file);
29  double readFieldData(std::ifstream &file);
30  void normalizeFieldData(double maxEz);
31  std::vector<double> getEvenlyDistributedSamplingPoints();
32  std::vector<double> interpolateFieldData(std::vector<double> &samplingPoints);
33  void computeFieldDerivatives(std::vector<double> & fourierComponents, int accuracy);
34  std::vector<double> computeFourierCoefficients(int accuracy, std::vector<double> &evenSampling);
35 
36 
37  double* onAxisField_m;
38  double* zvals_m;
39  gsl_spline *onAxisInterpolants_m[4];
40  gsl_interp_accel *onAxisAccel_m[4];
41 
42  double hz_m;
43 
44  double zbegin_m;
45  double zend_m;
46  double length_m;
47 
50 
51  friend class _Fieldmap;
52 };
53 
54 inline bool _Astra1D_fast::isInside(const Vector_t &r) const
55 {
56  return r(2) >= zbegin_m && r(2) < zend_m;
57 }
58 
59 using Astra1D_fast = std::shared_ptr<_Astra1D_fast>;
60 
61 #endif
virtual ~_Astra1D_fast()
int numHeaderLines_m
Definition: Astra1D_fast.h:49
DiffDirection
Definition: Fieldmap.h:55
virtual void getInfo(Inform *)=0
double length_m
Definition: Astra1D_fast.h:46
std::shared_ptr< _Astra1D_fast > Astra1D_fast
Definition: Astra1D_fast.h:59
virtual void getFieldDimensions(double &zBegin, double &zEnd) const =0
virtual double getFrequency() const =0
double zbegin_m
Definition: Astra1D_fast.h:44
virtual bool getFieldDerivative(const Vector_t &R, Vector_t &E, Vector_t &B, const DiffDirection &dir) const =0
gsl_spline * onAxisInterpolants_m[4]
Definition: Astra1D_fast.h:39
std::vector< double > interpolateFieldData(std::vector< double > &samplingPoints)
double * onAxisField_m
Definition: Astra1D_fast.h:37
virtual bool getFieldstrength(const Vector_t &R, Vector_t &E, Vector_t &B) const =0
virtual void swap()=0
virtual void freeMap()
void normalizeFieldData(double maxEz)
virtual void getOnaxisEz(std::vector< std::pair< double, double > > &F)
double readFieldData(std::ifstream &file)
Definition: Inform.h:42
virtual bool isInside(const Vector_t &r) const
Definition: Astra1D_fast.h:54
std::vector< double > getEvenlyDistributedSamplingPoints()
virtual void setFrequency(double freq)=0
_Astra1D_fast(const std::string &filename)
gsl_interp_accel * onAxisAccel_m[4]
Definition: Astra1D_fast.h:40
void computeFieldDerivatives(std::vector< double > &fourierComponents, int accuracy)
std::vector< double > computeFourierCoefficients(int accuracy, std::vector< double > &evenSampling)
virtual void readMap()
double * zvals_m
Definition: Astra1D_fast.h:38
bool determineNumSamplingPoints(std::ifstream &file)