35 "An error occured when reading the fieldmap '" +
Filename_m +
"'");
82 const double RR2 =
R(0) *
R(0) +
R(1) *
R(1);
90 const double BfieldR = -bzp / 2. + bzppp / 16. * RR2;
92 B(0) += BfieldR *
R(0);
93 B(1) += BfieldR *
R(1);
94 B(2) += bz - bzpp * RR2 / 4.;
124 std::string tmpString;
129 passed = interpretLine<std::string, int>(file, tmpString, tmpInt);
131 passed = interpretLine<std::string, int, std::string>(file, tmpString, tmpInt, tmpString);
134 if (tmpString !=
"TRUE" &&
135 tmpString !=
"FALSE")
137 "The third string on the first line of 1D field "
138 "maps has to be either TRUE or FALSE");
147 std::string tmpString;
151 interpretLine<std::string, int>(file, tmpString, accuracy);
153 interpretLine<std::string, int, std::string>(file, tmpString, accuracy, tmpString);
virtual void setFrequency(double freq)
virtual ~_Astra1DMagnetoStatic_fast()
virtual double getFrequency() const
_Astra1DMagnetoStatic_fast(const std::string &filename)
static std::string typeset_msg(const std::string &msg, const std::string &title)
std::shared_ptr< _Astra1DMagnetoStatic_fast > Astra1DMagnetoStatic_fast
gsl_spline * onAxisInterpolants_m[4]
std::string toUpper(const std::string &str)
Inform & endl(Inform &inf)
std::vector< double > interpolateFieldData(std::vector< double > &samplingPoints)
bool readFileHeader(std::ifstream &file)
virtual void getFieldDimensions(double &zBegin, double &zEnd) const
int stripFileHeader(std::ifstream &file)
void normalizeFieldData(double maxEz)
double readFieldData(std::ifstream &file)
static Astra1DMagnetoStatic_fast create(const std::string &filename)
std::vector< double > getEvenlyDistributedSamplingPoints()
void disableFieldmapWarning()
virtual void getInfo(Inform *)
virtual bool getFieldDerivative(const Vector_t &R, Vector_t &E, Vector_t &B, const DiffDirection &dir) const
gsl_interp_accel * onAxisAccel_m[4]
Inform & level3(Inform &inf)
constexpr double e
The value of .
void computeFieldDerivatives(std::vector< double > &fourierComponents, int accuracy)
std::vector< double > computeFourierCoefficients(int accuracy, std::vector< double > &evenSampling)
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)
bool determineNumSamplingPoints(std::ifstream &file)