Changeset 32 in TRACY3 for trunk/tracy/tracy/src/t2ring.cc

Timestamp:

Apr 9, 2014, 3:50:11 PM (10 years ago)

Author:

zhangj

Message:

active the transport of the twiss functions and orbits of the transfer line.

File:

• trunk/tracy/tracy/src/t2ring.cc (modified) (7 diffs)

trunk/tracy/tracy/src/t2ring.cc

@@ -196 +196 @@
 
    Purpose:
-     Computes Twiss parameters alpha and beta from the matrix Ascr
+     Assign the Twiss parameters alpha and beta from the matrix Ascr
+       to alpha & beta.
        
        M oneturn matrix    (A and M are symplectic)
@@ -218 +219 @@
 
    Return:
-       none
+       alpha, beta.
 
    Global variables:
@@ -233 +234 @@
 void getprm(Matrix &Ascr, Vector2 &alpha, Vector2 &beta)
 {
-  int  i, j;
+  int  i=0, j=0;
 
   for (i = 1; i <= 2; i++) {
@@ -253 +254 @@
        i0     first element
        i1     last element
-       ring   true if a ring
+       ring   true if a ring, false is a transfer line
        chroma true if compute chromaticities
        dP     energy offset
@@ -278 +279 @@
                   bool ring, double dP)
 {
-  long int  i;
-  int       j, k;
+  long int  i=0L;
+  int       j=0, k=0;
   Vector2   nu1, dnu;
-  Vector    xref;
-  Matrix    Ascr0, Ascr1;
+  Vector    xref;                   /*6*1 vector*/
+  Matrix    Ascr0, Ascr1;           /*6*6 matrix*/
   CellType  *cellp;
 
   if (dP != globval.dPparticle) Cell_SetdP(dP);
 
-  /* Init */  
+  /* Initize the phase */  
   for (j = 0; j <= 1; j++)
-    nu1[j] = 0.0;
-
-  /* get alpha beta for i0 */  
+    nu1[j] = 0.0;  /*phase advance*/
+
+  /* get alpha & beta for i0 */  
   cellp = &Cell[i0];
-  getprm(Ascr, cellp->Alpha, cellp->Beta);
-  memcpy(cellp->Nu, nu1, sizeof(Vector2));
+  getprm(Ascr, cellp->Alpha, cellp->Beta);  /*initialize the alpha & beta functions */
+  memcpy(cellp->Nu, nu1, sizeof(Vector2));  /*initialize the phase advance*/
   CopyMat(n+1, Ascr, Ascr0); CopyVec(n+2L, globval.CODvect, xref);
 
@@ -301 +302 @@
     /* Ascr1=Elem_M*Ascr0 */
     /* xref =Elem(xref)   */
-    Elem_Pass_M(i, xref, Ascr1); 
+    Elem_Pass_M(i, xref, Ascr1);   /* transport orbit & optic matrix */
 
     cellp = &Cell[i];
     /* get alpha and beta for element i */
     getprm(Ascr1, cellp->Alpha, cellp->Beta);
+    /*save the orbit at the end of the lattice elment*/
+    CopyVec(ss_dim, xref,cellp->CODvect);
 
     for (j = 0; j <= 1; j++) {
@@ -861 +864 @@
                 Vector2 &etap, Vector &codvect)
 {
-  long i, j, lastpos;
-  double sb;
-  Matrix Ascr;
+  long i=0L, j=0L, lastpos=0L;
+  double sb=0.0;
+  Matrix Ascr; /* 6*6 transfer matrix*/
 
   UnitMat(6L, Ascr);