diff --git a/tracy/tracy/src/t2elem.cc b/tracy/tracy/src/t2elem.cc index 75fdd491be35cf752fce3ed63d87f39f509705cf..41310c50710fecb0d567b7c749b895a3b258b905 100644 --- a/tracy/tracy/src/t2elem.cc +++ b/tracy/tracy/src/t2elem.cc @@ -1760,8 +1760,9 @@ void Insertion_Pass(CellType &Cell, ss_vect<T> &X) { switch (elemp->ID->Pmethod) { case Meth_Second: - for (i = 1; i <= Nslice; i++) { - Drift(LN/2.0, X); + Drift(0.5*LN, X); + + for (i = 1; i <= Nslice; i++) { // fprintf(stdout, "Second order L= %f, N= %d LN=%f S1=%f S2 = %f\n", //elemp->PL, Nslice, LN, elemp->ID->scaling1, elemp->ID->scaling1); @@ -1806,7 +1807,7 @@ void Insertion_Pass(CellType &Cell, ss_vect<T> &X) { if (i != Nslice) Drift(LN, X); else - Drift(LN/2.0, X); + Drift(0.5*LN, X); } // for break; case Meth_Fourth: /* 4-th order integrator */