diff --git a/corsika/detail/framework/geometry/LeapFrogTrajectory.inl b/corsika/detail/framework/geometry/LeapFrogTrajectory.inl index 7293a7f5969aa999e4aa14e6ef1422850b7ba3b9..af31f8209eedad47bd9de9cbf1344c17c9ef758e 100644 --- a/corsika/detail/framework/geometry/LeapFrogTrajectory.inl +++ b/corsika/detail/framework/geometry/LeapFrogTrajectory.inl @@ -48,6 +48,7 @@ namespace corsika { inline TimeType LeapFrogTrajectory::getDuration(double const u) const { TimeType const step = timeStep_ * u; double const correction = 1; + // the eventual (delta-L to L) correction factor is: // (initialDirection_ + initialDirection_.cross(magneticfield_) * step * // k_).getNorm(); return step / 2 * (correction + 1); @@ -66,7 +67,7 @@ namespace corsika { /* initial attempt to calculate delta-L from assumed full-leap-frog-length L: - Note: often return 0. Not good enough yet. + Note: often return 0. Not good enough yet. LengthType const L = initialVelocity_.getNorm() * limit; // distance double const a = (initialVelocity_.cross(magneticfield_) * k_).getSquaredNorm() / 4 /