我只学了几个月的python,完全是C语言的新手,我从网上得到了一个C代码,我很想学习它.但我只懂python语言,所以有人可以帮助将以下代码翻译成python会很棒.提前致谢!
for(i=0; i<n; i++) { /* Foreach particle "i" ... */
ax=0.0;
ay=0.0;
az=0.0;
for(j=0; j<n; j++) { /* Loop over all particles "j" */
dx=x[j]-x[i];
dy=y[j]-y[i];
dz=yz[j]-z[i];
invr = 1.0/sqrt(dx*dx + dy*dy + dz*dz + eps);
invr3 = invr*invr*invr;
f=m[j]*invr3;
ax += f*dx; /* accumulate the acceleration from gravitational attraction */
ay += f*dy;
az += f*dx;
}
xnew[i] = x[i] + dt*vx[i] + 0.5*dt*dt*ax; /* update position of particle "i" */
ynew[i] = y[i] + dt*vy[i] + 0.5*dt*dt*ay;
znew[i] = z[i] + dt*vz[i] + 0.5*dt*dt*az;
vx[i] += dt*ax; /* update velocity of particle "i" */
vy[i] += dt*ay;
vz[i] += dt*az;
}
Run Code Online (Sandbox Code Playgroud)
再次感谢!
归档时间: |
|
查看次数: |
297 次 |
最近记录: |