We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent bb0a461 commit 3670d6fCopy full SHA for 3670d6f
1 file changed
src/struct/particle.h
@@ -25,13 +25,14 @@ struct alignas(32) Particle {
25
inline void Gravity(Particle& a, Particle& b, real dt)
26
{
27
constexpr real G = real(1.0);
28
- real eps2 = pairSoftening(a.m, b.m); // must be epsilon^2
+
29
+ real eps = pairSoftening(a.m, b.m);
30
31
real dx = b.x - a.x;
32
real dy = b.y - a.y;
33
real dz = b.z - a.z;
34
- real r2 = dx*dx + dy*dy + dz*dz + eps2;
35
+ real r2 = dx*dx + dy*dy + dz*dz + eps*eps;
36
37
real invR2 = real(1.0) / r2;
38
real invR = std::sqrt(invR2);
0 commit comments