https://github.com/florentrenaud/nbody6tt
Tip revision: 8a4716382ead3ece116c48a4ae5c65f8c9534437 authored by Florent on 29 January 2015, 12:19:28 UTC
Nbody6 - 29 January 2015 (added GPU2/Build/.gitkeep)
Nbody6 - 29 January 2015 (added GPU2/Build/.gitkeep)
Tip revision: 8a47163
poti.f
***
SUBROUTINE POTI(I,POTJ)
*
*
* Potential of one particle.
* --------------------------
*
INCLUDE 'common6.h'
*
*
* Obtain the potential of body #I on host.
POTJ = 0.D0
KDUM = 0
DO 30 JDUM = IFIRST,NTOT
IF (JDUM.EQ.I) GO TO 30
J = JDUM
IF (J.GT.N) THEN
JPAIR = J - N
* Use c.m. approximation for unperturbed binary.
IF (LIST(1,2*JPAIR-1).GT.0) THEN
KDUM = 2*JPAIR - 1
J = KDUM
END IF
END IF
20 RIJ2 = 0.D0
DO 25 K = 1,3
RIJ2 = RIJ2 + (X(K,I) - X(K,J))**2
25 CONTINUE
POTJ = POTJ - BODY(J)/SQRT(RIJ2)
IF (J.EQ.KDUM) THEN
J = J + 1
GO TO 20
END IF
30 CONTINUE
*
RETURN
END
***