8#include "gsl/gsl_interp.h"
9#include "gsl/gsl_spline.h"
10#include "gsl/gsl_fft_real.h"
17 FourCoefs_m(nullptr) {
20 int skippedValues = 0;
21 std::string tmpString;
29 bool parsing_passed =
true;
31 parsing_passed = interpretLine<std::string, int>(file, tmpString,
accuracy_m);
33 parsing_passed = interpretLine<std::string, int, std::string>(file,
39 if (tmpString !=
"TRUE" &&
42 "The third string on the first line of 1D field "
43 "maps has to be either TRUE or FALSE");
48 parsing_passed = parsing_passed &&
49 interpretLine<double, double>(file,
zbegin_m, tmpDouble);
52 while(!file.eof() && parsing_passed) {
53 parsing_passed = interpretLine<double, double>(file,
zend_m, tmpDouble,
false);
54 if (
zend_m - tmpDouble2 > 1
e-10) {
56 }
else if (parsing_passed) {
64 if (!parsing_passed && !file.eof()) {
68 "An error occured when reading the fieldmap '" +
Filename_m +
"'");
89 bool parsing_passed =
true;
91 std::string tmpString;
100 gsl_spline *Ez_interpolant = gsl_spline_alloc(gsl_interp_cspline,
num_gridpz_m);
101 gsl_interp_accel *Ez_accel = gsl_interp_accel_alloc();
103 gsl_fft_real_wavetable *
real = gsl_fft_real_wavetable_alloc(2 *
num_gridpz_m);
104 gsl_fft_real_workspace *work = gsl_fft_real_workspace_alloc(2 *
num_gridpz_m);
113 parsing_passed = interpretLine<double, double>(in, zvals[i], RealValues[i]);
116 if (zvals[i] - tmpDouble > 1
e-10) {
117 if (
std::abs(RealValues[i]) > Ez_max) {
120 tmpDouble = zvals[i];
126 gsl_spline_init(Ez_interpolant, zvals, RealValues,
num_gridpz_m);
133 RealValues[ii] = gsl_spline_eval(Ez_interpolant, z, Ez_accel);
139 RealValues[i] = RealValues[ii];
152 for (
int i = 1; i < 2 *
accuracy_m - 1; i++) {
156 gsl_fft_real_workspace_free(work);
157 gsl_fft_real_wavetable_free(
real);
159 gsl_spline_free(Ez_interpolant);
160 gsl_interp_accel_free(Ez_accel);
182 const double RR2 =
R(0) *
R(0) +
R(1) *
R(1);
194 double somefactor = 1.0;
195 double coskzl =
cos(kz * l);
196 double sinkzl =
sin(kz * l);
198 somefactor *= somefactor_base;
200 somefactor *= somefactor_base;
202 somefactor *= somefactor_base;
207 const double EfieldR = -ezp / 2. + ezppp / 16. * RR2;
209 E(0) += EfieldR *
R(0);
210 E(1) += EfieldR *
R(1);
211 E(2) += ez - ezpp * RR2 / 4.;
FLieGenerator< T, N > real(const FLieGenerator< std::complex< T >, N > &)
Take real part of a complex generator.
Tps< T > cos(const Tps< T > &x)
Cosine.
Tps< T > sin(const Tps< T > &x)
Sine.
PETE_TUTree< FnAbs, typename T::PETE_Expr_t > abs(const PETE_Expr< T > &l)
Inform & endl(Inform &inf)
Inform & level3(Inform &inf)
constexpr double two_pi
The value of.
constexpr double e
The value of.
constexpr double pi
The value of.
constexpr double MVpm2Vpm
std::string toUpper(const std::string &str)
virtual bool getFieldstrength(const Vector_t &R, Vector_t &E, Vector_t &B) const
virtual double getFrequency() const
Astra1DElectroStatic(std::string aFilename)
virtual void setFrequency(double freq)
virtual void getFieldDimensions(double &zBegin, double &zEnd) const
virtual bool getFieldDerivative(const Vector_t &R, Vector_t &E, Vector_t &B, const DiffDirection &dir) const
virtual void getInfo(Inform *)
void disableFieldmapWarning()
static std::string typeset_msg(const std::string &msg, const std::string &title)
void getLine(std::ifstream &in, std::string &buffer)