Fixed Bugs

Fixed eps softening + Newton's 3rd law.

Author: TimGoTheCreator

src/struct/particle.h

@@ -16,23 +16,22 @@ inline void Gravity(Particle& a, Particle& b, real dt)
     real dx = b.x - a.x;
     real dy = b.y - a.y;
     real dz = b.z - a.z;
-    real eps = real(1e-6);
+    real eps = real(1e6);
     real distSq = dx*dx + dy*dy + dz*dz + eps*eps;
     real dist = std::sqrt(distSq);
 
     real invDist = real(1) / dist;
     real invDist3 = invDist * invDist * invDist;
-        
+
     real ax = G * b.m * dx * invDist3;
     real ay = G * b.m * dy * invDist3;
     real az = G * b.m * dz * invDist3;
 
-    real inv_bm = real(1) / b.m;
-    real scale = -a.m * inv_bm;
+    real bx = -G * a.m * dx * invDist3;
+    real by = -G * a.m * dy * invDist3;
+    real bz = -G * a.m * dz * invDist3;
 
-    real bx = ax * scale;
-    real by = ay * scale;
-    real bz = az * scale;
+    
 
     a.vx += ax * dt;
     a.vy += ay * dt;
@@ -42,3 +41,4 @@ inline void Gravity(Particle& a, Particle& b, real dt)
     b.vy += by * dt;
     b.vz += bz * dt;
 }
+