/* Tracy-2

   J. Bengtsson, CBP, LBL      1990 - 1994   Pascal version
                 SLS, PSI      1995 - 1997
   M. Boege      SLS, PSI      1998          C translation
   L. Nadolski   SOLEIL        2002          Link to NAFF, Radia field maps
   J. Bengtsson  NSLS-II, BNL  2004 -        

*/

// Intel compiler
#ifndef __APPLE_CC__
#pragma warning (disable:193)
#pragma warning (disable:981)
#pragma warning (disable:1572)
#endif
// LSN

#define DOF    (ss_dim/2)

// defined in math.h
//#define M_PI   3.14159265358979323846  // pi

#ifndef LONG_MAX
# define LONG_MAX   ((long)(((unsigned long) -1) >> 1))
# define LONG_MIN   (~LONG_MAX)
#endif

#define debug           false

#define S_SIZE 200  // max size for file name of a lattice file

#define nv_             6

#define mnp             400

// maximum number of LEGO blocks (Cell_nLoc)
#define Cell_nLocMax    8000

// maximum number of families for Elem_NFam
#define Elem_nFamMax    3000

// maximum number of kids */
#define nKidMax         1000

#define DBNameLen       39

#define NameLength      150  // maximum length of identifiers
#define SymbolLength    15   // maximum length of element name

#define blankname       "               "

#define maxincl         5
#define maxfil          10
#define bigvectmax      4096

#define HOMmax          30

#define maxtransfmat    2000
#define maxkicks        5


const double  c0    = 2.99792458e8;             // speed of light in vacuum
const double  q_e   = 1.60217733e-19;                // electron charge
const double  m_e   = 0.51099906e6;             // electron rest mass [eV]
const double  mu_0  = 4.0*M_PI*1e-7;            // permittivity of free space
const double  eps_0 = 1.0/(sqr(c0)*mu_0);       // permeability of free space
const double  r_e   = q_e/(4.0*M_PI*eps_0*m_e); // classical electron radius
const double  h_bar = 6.58211899e-16;           // reduced Planck constant [eV]

// const double  max_ampl = 10.0; // [m]   nsrl-ii
const double  max_ampl = 1.0; // [m]    soleil


enum pthicktype { thick = 0, thin = 1 };

enum PartsKind { drift = 0, Wigl = 1, Mpole = 2, Cavity = 3, marker = 4,
		 undef = 5, Insertion = 6, FieldMap = 7,
                 Spreader = 8, Recombiner = 9, Solenoid = 10 };

enum { All = 0, Dip = 1, Quad = 2, Sext = 3, Oct = 4, Dec = 5, Dodec = 6 };

enum { Horizontal = 1, Vertical = 2 };

enum { Meth_Linear = 0, Meth_First = 1, Meth_Second = 2, Meth_Fourth = 4,
       Meth_genfun = 5 };

typedef double svdarray[mnp];
typedef double svdmat[mnp][mnp];

typedef char str80[80];

typedef char alfa_[NameLength];

typedef long   iVector2[2];
typedef double Vector2[2];
typedef double Vector3[3];

#define fitvectmax      200
typedef long   fitvect[fitvectmax];

typedef char partsName[NameLength];

typedef char DBNameType[DBNameLen];

typedef double mpolArray[HOMmax+HOMmax+1];


typedef struct statusrec{
  bool tuneflag, chromflag, codflag, mapflag, passflag, overflag, chambre;
  int lossplane; /* lost in: horizontal    1
		                     vertical      2
			                 longitudinal  3 */
} statusrec;


extern bool  stable;
extern bool  ErrFlag;
extern bool  trace, cellconcat, traceID;

extern double  Fdrift1, Fkick1, Fdrift2, Fkick2, crad, cfluc;

extern str80  finame,   /* input  data file  */
              foname,   /* output data file */
              fname;    /* temp file name */

extern FILE  *fi,       /* lattice input  file  */
             *fo,       /* lattice output file */
             *psin[],   /* program input file */
             *psout,              /* program output file*/
             *prr[];   /* prr[1] : input, prr[2] : output */

extern statusrec status;


extern int P_eof(FILE *f);

extern int P_eoln(FILE *f);

extern void GDiag(int n, double C, Matrix &A, Matrix &Ainv, Matrix &R,
		  Matrix &M, double &Omega, double &Yalphac);
extern void NormEigenVec(Matrix &Vr, Matrix &Vi, double *wr, double *wi,
			 Matrix &t6a);

extern void EigenVal(double (*Ai)[mnp], long n, double *wr, double *wi);

extern void t2init(void);

/****** for orbit correction *********/
extern void prt_gcmat(int bpm, int corr, int plane);

extern void gcmat(int bpm, int corr, int plane);

extern void lsoc(int niter, int bpm, int corr, int plane);


/**** same as asctime in C without the \n at the end****/
char *asctime2(const struct tm *timeptr);
struct tm* GetTime(void);
uint32_t stampstop(uint32_t start);
uint32_t stampstart(void);