source: TRACY3/trunk/tracy/tracy/inc/t2elem.h @ 11

/* Tracy-3

   J. Bengtsson, BNL 2007

*/

extern double  cl_rad, q_fluct, I2, I4, I5;


void getelem(long i, CellType *cellrec);

void putelem(long i, CellType *cellrec);


int GetnKid(const int Fnum1);

long Elem_GetPos(const int Fnum1, const int Knum1);


template<typename T>
void GtoL(ss_vect<T> &X, Vector2 &S, Vector2 &R,
          const double c0, const double c1, const double s1);

template<typename T>
void LtoG(ss_vect<T> &X, Vector2 &S, Vector2 &R,
          double c0, double c1, double s1);

template<typename T>
void p_rot(double phi, ss_vect<T> &x);


template<typename T>
void get_B2(const double h_ref, const T B[], const ss_vect<T> &xp,
            T &B2_perp, T &B2_par);

template<typename T>
void radiate(ss_vect<T> &x, const double L, const double h_ref, const T B[]);

template<typename T>
void Drift(double L,double h_bend, ss_vect<T> &x);

template<typename T>
void Drift(double L, ss_vect<T> &x);

template<typename T>
void bend_fringe(double hb, ss_vect<T> &x);

template<typename T>
void bend_fringe(double hb, double gap, ss_vect<T> &x);

template<typename T>
static void bendCurvature(double irho, double H, ss_vect<T> &x); 

template<typename T>
static void EdgeFocus(double irho, double phi, double gap, ss_vect<T> &x, bool Entrance);

template<typename T>
void quad_fringe(double b2, ss_vect<T> &x);


template<typename T>
void Drift_Pass(CellType &Cell, ss_vect<T> &x);

template<typename T>
void dipole_kick(double L, double h_ref, double h_bend, ss_vect<T> &x);
template<typename T>
void multipole_kick(int Order, double MB[], double L, double h_bend, ss_vect<T> &x);

template<typename T>
void thin_kick(int Order, double MB[], double L, double h_bend, double h_ref,
               ss_vect<T> &x);

template<typename T>
void Mpole_Pass(CellType &Cell, ss_vect<T> &x);

template<typename T>
void Marker_Pass(CellType &Cell, ss_vect<T> &X);

template<typename T>
void Cav_Pass(CellType &Cell, ss_vect<T> &X);

template<typename T>
void Wiggler_pass_EF(const elemtype &elem, ss_vect<T> &x);

template<typename T>
void Wiggler_pass_EF2(int nstep, double L, double kxV, double kxH, double kz, 
                      double BoBrhoV, double BoBrhoH, double phi,
                      ss_vect<T> &x);

template<typename T>
void Wiggler_pass_EF3(const elemtype &elem, ss_vect<T> &x);

template<typename T>
void Wiggler_Pass(CellType &Cell, ss_vect<T> &X);

template<typename T>
void FieldMap_Pass(CellType &Cell, ss_vect<T> &X);

template<typename T>
void Insertion_Pass(CellType &Cell, ss_vect<T> &x);

template<typename T>
void sol_pass(const elemtype &elem, ss_vect<T> &x);

template<typename T>
void Solenoid_Pass(CellType &Cell, ss_vect<T> &x);


void Drift_SetMatrix(int Fnum1, int Knum1);

void Insertion_SetMatrix(int Fnum1, int Knum1);


void Mpole_SetPB(int Fnum1, int Knum1, int Order);

void Wiggler_SetPB(int Fnum1, int Knum1, int Order);


void GtoL_M(Matrix &X, Vector2 &T);

void LtoG_M(Matrix &X, Vector2 &T);

void Drift_Pass_M(CellType &Cell, Vector &xref, Matrix &X);

void thin_kick_M(int Order, double MB[], double L, double irho,
                 Vector &xref, Matrix &x);

void Mpole_Pass_M(CellType &Cell, Vector &xref, Matrix &x);

void Wiggler_Pass_M(CellType &Cell, Vector &xref, Matrix &x);

void Insertion_Pass_M(CellType &Cell, Vector &xref, Matrix &M);


void MulLsMat(Matrix &A, Matrix &B);

void LinsTrans(Matrix &A, Vector &b);


void Drift_Alloc(elemtype *Elem);

void Mpole_Alloc(elemtype *Elem);

void Cav_Alloc(elemtype *Elem);

void Wiggler_Alloc(elemtype *Elem);

void FieldMap_Alloc(elemtype *Elem, const bool alloc_fm);

void Insertion_Alloc(elemtype *Elem);

void Spreader_Alloc(elemtype *Elem);

void Recombiner_Alloc(elemtype *Elem);

void Solenoid_Alloc(elemtype *Elem);


void SI_init(void);

void Drift_Init(int Fnum1);

void Mpole_Init(int Fnum1);

void Wiggler_Init(int Fnum1);

void FieldMap_Init(int Fnum1);

void Cav_Init(int Fnum1);

void Marker_Init(int Fnum1);

void Insertion_Init(int Fnum1);

void Spreader_Init(int Fnum1);

void Recombiner_Init(int Fnum1);

void Solenoid_Init(int Fnum1);


void get_B(const char *file_name, FieldMapType *FM);

double Elem_GetKval(int Fnum1, int Knum1, int Order);

void Mpole_SetdS(int Fnum1, int Knum1);

void Mpole_SetdT(int Fnum1, int Knum1);