Disable debug prints

atcollab · Sep 11, 2023 · cbecf0f · cbecf0f
1 parent 542cc44
commit cbecf0f
Show file tree

Hide file tree

Showing 2 changed files with 16 additions and 16 deletions.
diff --git a/atintegrators/driftkickrad.c b/atintegrators/driftkickrad.c
@@ -1,4 +1,4 @@
-#define ZUT
+#undef DPRINT
 #undef NOP
 /***********************************************************************
  Note: in the US convention the transverse multipole field is written as:
@@ -143,18 +143,18 @@ the polynomial terms in PolynomB.
    B2P = B2perp(ImSum, ReSum +irho, irho, x , xpr, y ,ypr);
 
    dp_0 = r[4];
-   #ifdef ZUT
+   #ifdef DPRINT
    double dr4 = SQR(1+r[4])*B2P*(1 + x*irho + (SQR(xpr)+SQR(ypr))/2 );
    atPrintf("B2P: % 15.12g\t% 15.12g\t% 15.12g\t% 15.12g\n", x, xpr, B2P, dr4);
-   #endif
+   #endif /*DPRINT*/
    #ifndef NOP
    r[4] = r[4] - CRAD*SQR(1+r[4])*B2P*(1 + x*irho + (SQR(xpr)+SQR(ypr))/2 )*L;
 
    /* recalculate momentums from angles after losing energy for radiation 	*/
    p_norm = 1/(1+r[4]);
    r[1] = xpr/p_norm;
    r[3] = ypr/p_norm;
-   #endif /*ZUT*/
+   #endif /*NOP*/
 
    r[1] -=  L*(ReSum-(dp_0-r[0]*irho)*irho);
    r[3] +=  L*ImSum;
@@ -198,18 +198,18 @@ static void strthinkickrad(double* r, const double* A, const double* B, double L
    B2P = StrB2perp(ImSum, ReSum , x , xpr, y ,ypr);
 
    dp_0 = r[4];
-   #ifdef ZUT
+   #ifdef DPRINT
    double dr4 = SQR(1+r[4])*B2P*(1 + x*irho + (SQR(xpr)+SQR(ypr))/2 );
    atPrintf("B2P: % 15.12g\t% 15.12g\t% 15.12g\t% 15.12g\n", x, xpr, B2P, dr4);
-   #endif
+   #endif /*DPRINT */
    #ifndef NOP
    r[4] = r[4] - CRAD*SQR(1+r[4])*B2P*(1 + x*irho + (SQR(xpr)+SQR(ypr))/2 )*L;
 
    /* recalculate momentums from angles after losing energy for radiation 	*/
    p_norm = 1/(1+r[4]);
    r[1] = xpr/p_norm;
    r[3] = ypr/p_norm;
-   #endif /*ZUT*/
+   #endif /*NOP*/
 
    r[1] -=  L*(ReSum-(dp_0-r[0]*irho)*irho);
    r[3] +=  L*ImSum;

diff --git a/atintegrators/exactkickrad.c b/atintegrators/exactkickrad.c
@@ -1,4 +1,4 @@
-#define ZUT
+#undef DPRINT
 #undef NOP
 /***********************************************************************
  Note: in the US convention the transverse multipole field is written as:
@@ -133,10 +133,10 @@ the polynomial terms in PolynomB.
 
    B2P = B2perp(ImSum, ReSum+irho, irho, x , xpr, y ,ypr);
 
-   #ifdef ZUT
+   #ifdef DPRINT
    double dr4 = SQR(1+r[4]) * B2P * (1.0+x*irho);
    atPrintf("B2P: % 15.12g\t% 15.12g\t% 15.12g\t% 15.12g\n", x, xpr, B2P, dr4);
-   #endif
+   #endif /*DPRINT */
    #ifndef NOP
    r[4] = r[4] - CRAD * SQR(1+r[4]) * B2P * (1.0+x*irho) * L;
 //   r[4] = r[4] - CRAD*SQR(1+r[4])*B2P*(1 + x*irho + (SQR(xpr)+SQR(ypr))/2 )*L;
@@ -183,18 +183,18 @@ static void ex_strthinkickrad(double* r, const double* A, const double* B, doubl
 
    B2P = StrB2perp(ImSum, ReSum+B0 , x , xpr, y ,ypr);
 
-   #ifdef ZUT
+   #ifdef DPRINT
    double dr4 = SQR(1+r[4]) * B2P;
    atPrintf("B2P: % 15.12g\t% 15.12g\t% 15.12g\t% 15.12g\n", x, xpr, B2P, dr4);
-   #endif
+   #endif /*DPRINT */
    #ifndef NOP
    r[4] = r[4] - CRAD * SQR(1+r[4]) * B2P * L;
 
    /* recalculate momentums from angles after losing energy for radiation 	*/
    p_norm = 1/(1+r[4]);
    r[1] = xpr/p_norm;
    r[3] = ypr/p_norm;
-   #endif /*ZUT*/
+   #endif /*NOP*/
 
    r[1] -=  L*ReSum;
    r[3] +=  L*ImSum;
@@ -227,18 +227,18 @@ static void ex2_bndthinkickrad(double* r, double* A, double* B, double K, double
 
    B2P = B2perp(ImSum, ReSum+irho, irho, x , xpr, y ,ypr);
 
-   #ifdef ZUT
+   #ifdef DPRINT
    double dr4 = SQR(1+r[4]) * B2P * (1.0+x*irho);
    atPrintf("B2P: % 15.12g\t% 15.12g\t% 15.12g\t% 15.12g\n", x, xpr, B2P, dr4);
-   #endif
+   #endif /*DPRINT */
    #ifndef NOP
    r[4] = r[4] - CRAD * SQR(1+r[4]) * B2P * (1.0+x*irho) * L;
 //   r[4] = r[4] - CRAD*SQR(1+r[4])*B2P*(1 + x*irho + (SQR(xpr)+SQR(ypr))/2 )*L;
    /* recalculate momentums from angles after losing energy for radiation 	*/
    p_norm = 1/(1+r[4]);
    r[1] = xpr/p_norm;
    r[3] = ypr/p_norm;
-   #endif /*ZUT*/
+   #endif /*NOP*/
 
    r[1] -=  0.5*K*ReSum;
    r[3] +=  0.5*K*ImSum;