1 #ifndef CLASSIC_FIELDMAP1DDYNAMICFAST_HH
2 #define CLASSIC_FIELDMAP1DDYNAMICFAST_HH
13 double &rBegin,
double &rEnd)
const;
15 double &yIni,
double &yFinal,
16 double &zIni,
double &zFinal)
const;
21 virtual void getOnaxisEz(std::vector<std::pair<double, double>> &eZ);
31 bool checkFileData(std::ifstream &fieldFile,
bool parsingPassed);
33 double onAxisFieldP[],
double onAxisFieldPP[],
34 double onAxisFieldPPP[]);
36 std::vector<double> fieldComponents)
const;
40 double onAxisFieldPP[],
41 double onAxisFieldPPP[]);
43 void normalizeField(
double maxEz, std::vector<double> &fourierCoefs);
44 double readFileData(std::ifstream &fieldFile,
double fieldData[]);
46 std::vector<std::pair<double, double>> &eZ);
48 void scaleField(
double maxEz, std::vector<std::pair<double, double>> &eZ);
gsl_interp_accel * onAxisFieldPPPAccel_m
virtual bool getFieldstrength(const Vector_t &R, Vector_t &E, Vector_t &B) const
virtual void getInfo(Inform *)
virtual void setFrequency(double freq)
void stripFileHeader(std::ifstream &fieldFile)
double zBegin_m
Maximum radius of field.
double readFileData(std::ifstream &fieldFile, double fieldData[])
void prepareForMapCheck(std::vector< double > &fourierCoefs)
FM1DDynamic_fast(std::string aFilename)
virtual void getOnaxisEz(std::vector< std::pair< double, double >> &eZ)
double deltaZ_m
Number of grid points in field input file.
virtual bool isInside(const Vector_t &r) const
void scaleField(double maxEz, std::vector< std::pair< double, double >> &eZ)
virtual bool getFieldDerivative(const Vector_t &R, Vector_t &E, Vector_t &B, const DiffDirection &dir) const
bool checkFileData(std::ifstream &fieldFile, bool parsingPassed)
unsigned int accuracy_m
Field grid point spacing.
gsl_interp_accel * onAxisFieldAccel_m
On axis field third derivative interpolation structure.
unsigned int numberOfGridPoints_m
Field length.
double twoPiOverLambdaSq_m
Field angular frequency (Hz).
void computeInterpolationVectors(double onAxisFieldP[], double onAxisFieldPP[], double onAxisFieldPPP[])
double rBegin_m
2 Pi divided by the field RF wavelength squared.
double length_m
Longitudinal end of field.
void computeFieldOffAxis(const Vector_t &R, Vector_t &E, Vector_t &B, std::vector< double > fieldComponents) const
gsl_spline * onAxisFieldInterpolants_m
On axis field data.
double zEnd_m
Longitudinal start of field.
gsl_spline * onAxisFieldPPInterpolants_m
On axis field first derivative interpolation structure.
gsl_interp_accel * onAxisFieldPPAccel_m
gsl_spline * onAxisFieldPPPInterpolants_m
On axis field second derivative interpolation structure.
double rEnd_m
Minimum radius of field.
void computeFieldOnAxis(double z, std::vector< double > &fieldComponents) const
gsl_interp_accel * onAxisFieldPAccel_m
virtual double getFrequency() const
gsl_spline * onAxisFieldPInterpolants_m
On axis field interpolation structure.
bool readFileHeader(std::ifstream &fieldFile)
std::vector< double > computeFourierCoefficients(double fieldData[])
void normalizeField(double maxEz, std::vector< double > &fourierCoefs)
void computeFieldDerivatives(std::vector< double > fourierCoefs, double onAxisFieldP[], double onAxisFieldPP[], double onAxisFieldPPP[])
virtual void getFieldDimensions(double &zBegin, double &zEnd, double &rBegin, double &rEnd) const