8 #include "gsl/gsl_interp.h"
9 #include "gsl/gsl_spline.h"
10 #include "gsl/gsl_fft_real.h"
18 FourCoefs_m(nullptr) {
21 int skippedValues = 0;
22 std::string tmpString;
30 bool parsing_passed =
true;
32 parsing_passed = interpretLine<std::string, int>(file,
36 parsing_passed = interpretLine<std::string, int, std::string>(file,
42 if (tmpString !=
"TRUE" &&
45 "The third string on the first line of 1D field "
46 "maps has to be either TRUE or FALSE");
51 parsing_passed = parsing_passed &&
53 parsing_passed = parsing_passed &&
54 interpretLine<double, double>(file,
zbegin_m, tmpDouble);
57 while(!file.eof() && parsing_passed) {
58 parsing_passed = interpretLine<double, double>(file,
zend_m, tmpDouble,
false);
59 if (zend_m - tmpDouble2 > 1e-10) {
61 }
else if (parsing_passed) {
69 if (!parsing_passed && !file.eof()) {
73 "An error occured when reading the fieldmap '" +
Filename_m +
"'");
102 bool parsing_passed =
true;
104 std::string tmpString;
113 gsl_spline *Ez_interpolant = gsl_spline_alloc(gsl_interp_cspline,
num_gridpz_m);
114 gsl_interp_accel *Ez_accel = gsl_interp_accel_alloc();
116 gsl_fft_real_wavetable *
real = gsl_fft_real_wavetable_alloc(2 *
num_gridpz_m);
117 gsl_fft_real_workspace *
work = gsl_fft_real_workspace_alloc(2 *
num_gridpz_m);
128 parsing_passed = interpretLine<double, double>(in, zvals[i], RealValues[i]);
131 if (zvals[i] - tmpDouble > 1
e-10) {
132 if (
std::abs(RealValues[i]) > Ez_max) {
135 tmpDouble = zvals[i];
141 gsl_spline_init(Ez_interpolant, zvals, RealValues,
num_gridpz_m);
148 RealValues[ii] = gsl_spline_eval(Ez_interpolant, z, Ez_accel);
150 RealValues[ii ++] = RealValues[num_gridpz_m - 1];
154 RealValues[i] = RealValues[ii];
157 gsl_fft_real_transform(RealValues, 1, 2 * num_gridpz_m, real, work);
164 for (
int i = 1; i < 2 *
accuracy_m - 1; ++ i) {
168 gsl_spline_free(Ez_interpolant);
169 gsl_interp_accel_free(Ez_accel);
171 gsl_fft_real_workspace_free(work);
172 gsl_fft_real_wavetable_free(real);
190 const double RR2 =
R(0) *
R(0) +
R(1) *
R(1);
200 for (
int l = 1; l <
accuracy_m ; ++ l, n += 2) {
202 double somefactor = 1.0;
203 double coskzl =
cos(kz * l);
204 double sinkzl =
sin(kz * l);
206 somefactor *= somefactor_base;
208 somefactor *= somefactor_base;
210 somefactor *= somefactor_base;
217 const double EfieldR = -(ezp / 2. + fp * RR2);
220 E(0) += EfieldR *
R(0);
221 E(1) += EfieldR *
R(1);
222 E(2) += ez + 4. * f * RR2;
223 B(0) += -BfieldT *
R(1);
224 B(1) += BfieldT *
R(0);
268 std::string tmpString;
272 interpretLine<std::string, int>(in, tmpString, tmpInt);
273 interpretLine<double>(in, tmpDouble);
276 interpretLine<double, double>(in, F[i].first, F[i].second);
277 if (
std::abs(F[i].second) > Ez_max) {
284 F[i].second /= Ez_max;
constexpr double c
The velocity of light in m/s.
constexpr double two_pi
The value of .
std::shared_ptr< _Astra1DDynamic > Astra1DDynamic
PETE_TUTree< FnAbs, typename T::PETE_Expr_t > abs(const PETE_Expr< T > &l)
_Astra1DDynamic(const std::string &filename)
double *__restrict__ FourCoefs_m
virtual void getOnaxisEz(std::vector< std::pair< double, double > > &F)
static std::string typeset_msg(const std::string &msg, const std::string &title)
constexpr double pi
The value of .
std::string toUpper(const std::string &str)
constexpr double Vpm2MVpm
Inform & endl(Inform &inf)
virtual double getFrequency() const
virtual void setFrequency(double freq)
static Astra1DDynamic create(const std::string &filename)
virtual void getFieldDimensions(double &zBegin, double &zEnd) const
virtual void getInfo(Inform *)
void disableFieldmapWarning()
Tps< T > cos(const Tps< T > &x)
Cosine.
and that you know you can do these things To protect your we need to make restrictions that forbid anyone to deny you these rights or to ask you to surrender the rights These restrictions translate to certain responsibilities for you if you distribute copies of the or if you modify it For if you distribute copies of such a whether gratis or for a you must give the recipients all the rights that you have You must make sure that receive or can get the source code And you must show them these terms so they know their rights We protect your rights with two distribute and or modify the software for each author s protection and we want to make certain that everyone understands that there is no warranty for this free software If the software is modified by someone else and passed we want its recipients to know that what they have is not the so that any problems introduced by others will not reflect on the original authors reputations any free program is threatened constantly by software patents We wish to avoid the danger that redistributors of a free program will individually obtain patent in effect making the program proprietary To prevent we have made it clear that any patent must be licensed for everyone s free use or not licensed at all The precise terms and conditions for distribution and modification follow GNU GENERAL PUBLIC LICENSE TERMS AND CONDITIONS FOR DISTRIBUTION AND MODIFICATION This License applies to any program or other work which contains a notice placed by the copyright holder saying it may be distributed under the terms of this General Public License The refers to any such program or work
void getLine(std::ifstream &in, std::string &buffer)
virtual ~_Astra1DDynamic()
Inform & level3(Inform &inf)
constexpr double e
The value of .
virtual bool getFieldDerivative(const Vector_t &R, Vector_t &E, Vector_t &B, const DiffDirection &dir) const
Tps< T > sin(const Tps< T > &x)
Sine.
FLieGenerator< T, N > real(const FLieGenerator< std::complex< T >, N > &)
Take real part of a complex generator.
virtual bool getFieldstrength(const Vector_t &R, Vector_t &E, Vector_t &B) const