https://github.com/nbodyx/Nbody6
Tip revision: dff3b5c68673d4d4c4c9f54c8ba110b8416098f2 authored by nitadori on 31 May 2020, 13:04:00 UTC
Avoided touching uninitialized SEMIX
Avoided touching uninitialized SEMIX
Tip revision: dff3b5c
jpred2.f
SUBROUTINE JPRED2(I,IPRED,CMX,CMV)
*
*
* Neighbour prediction of single particle.
* ----------------------------------------
*
INCLUDE 'common6.h'
REAL*8 CMX(3),CMV(3)
*
*
* Choose between copy and standard prediction.
IF (IPRED.EQ.0) THEN
DO 10 K = 1,3
CMX(K) = X(K,I)
CMV(K) = XDOT(K,I)
10 CONTINUE
ELSE
S = TIME - T0(I)
S1 = 1.5*S
S2 = 2.0*S
CMX(1) = ((FDOT(1,I)*S + F(1,I))*S + X0DOT(1,I))*S + X0(1,I)
CMX(2) = ((FDOT(2,I)*S + F(2,I))*S + X0DOT(2,I))*S + X0(2,I)
CMX(3) = ((FDOT(3,I)*S + F(3,I))*S + X0DOT(3,I))*S + X0(3,I)
CMV(1) = (FDOT(1,I)*S1 + F(1,I))*S2 + X0DOT(1,I)
CMV(2) = (FDOT(2,I)*S1 + F(2,I))*S2 + X0DOT(2,I)
CMV(3) = (FDOT(3,I)*S1 + F(3,I))*S2 + X0DOT(3,I)
END IF
*
RETURN
*
END