diff --git a/src/kick/kick.f90 b/src/kick/kick.f90 index efad4702a..56e74dfcd 100644 --- a/src/kick/kick.f90 +++ b/src/kick/kick.f90 @@ -56,9 +56,6 @@ module pure subroutine kick_getacch_int_tp(self, GMpl, xhp, npl) acc(:) = 0.0_DP do j = 1, npl dx(:) = tp%xh(:,i) - xhp(:, j) - !rji2 = dot_product(dx(:), dx(:)) - !irij3 = 1.0_DP / (rji2 * sqrt(rji2)) - !fac = GMpl(j) * irij3 r2 = dot_product(dx(:), dx(:)) fac = GMpl(j) / (r2 * sqrt(r2)) acc(:) = acc(:) - fac * dx(:)