29 if(!parsing_passed && !file.eof()) {
33 "An error occured when reading the fieldmap '" +
Filename_m +
"'");
77 const double RR2 =
R(0) *
R(0) +
R(1) *
R(1);
85 const double EfieldR = -ezp / 2. + ezppp / 16. * RR2;
87 E(0) += EfieldR *
R(0);
88 E(1) += EfieldR *
R(1);
89 E(2) += ez - ezpp * RR2 / 4.;
118 std::string tmpString;
123 passed = interpreteLine<std::string, int>(file, tmpString, tmpInt);
125 passed = interpreteLine<std::string, int, std::string>(file, tmpString, tmpInt, tmpString);
128 if (tmpString !=
"TRUE" &&
129 tmpString !=
"FALSE")
131 "The third string on the first line of 1D field "
132 "maps has to be either TRUE or FALSE");
141 std::string tmpString;
145 interpreteLine<std::string, int>(file, tmpString, accuracy);
147 interpreteLine<std::string, int, std::string>(file, tmpString, accuracy, tmpString);
constexpr double e
The value of .
gsl_spline * onAxisInterpolants_m[4]
virtual void getFieldDimensions(double &zBegin, double &zEnd, double &rBegin, double &rEnd) const
std::string toUpper(const std::string &str)
void disableFieldmapWarning()
virtual void getInfo(Inform *)
std::vector< double > getEvenlyDistributedSamplingPoints()
static std::string typeset_msg(const std::string &msg, const std::string &title)
bool determineNumSamplingPoints(std::ifstream &file)
Astra1DElectroStatic_fast(std::string aFilename)
std::vector< double > interpolateFieldData(std::vector< double > &samplingPoints)
double readFieldData(std::ifstream &file)
bool readFileHeader(std::ifstream &file)
gsl_interp_accel * onAxisAccel_m[4]
void computeFieldDerivatives(std::vector< double > &fourierComponents, int accuracy)
virtual void setFrequency(double freq)
virtual double getFrequency() const
std::vector< double > computeFourierCoefficients(int accuracy, std::vector< double > &evenSampling)
virtual bool getFieldDerivative(const Vector_t &R, Vector_t &E, Vector_t &B, const DiffDirection &dir) const
~Astra1DElectroStatic_fast()
Inform & level3(Inform &inf)
virtual bool getFieldstrength(const Vector_t &R, Vector_t &E, Vector_t &B) const
void checkMap(unsigned int accuracy, std::pair< double, double > fieldDimensions, double deltaZ, const std::vector< double > &fourierCoefficients, gsl_spline *splineCoefficients, gsl_interp_accel *splineAccelerator)
int stripFileHeader(std::ifstream &file)
void normalizeFieldData(double maxEz)
Inform & endl(Inform &inf)