Added optimizations (/ to sqrt)

TimGoTheCreator · web-flow · commit b210e7db16fe · 2026-02-02T20:10:09.000+03:00
diff --git a/src/gravity/octree.h b/src/gravity/octree.h
@@ -115,8 +115,8 @@ inline void bhAccel(Octree* node, const Particle& p, real theta,
     // Correct Barnes–Hut acceptance criterion
     if (node->leaf || (node->size / dist) < theta)
     {
-        real invDist  = real(1) / dist;
-        real invDist3 = invDist * invDist * invDist;
+        real invDist = std::sqrt(real(1.0) / distSq);
+real invDist3 = invDist * invDist * invDist;
 
         real fac = G * node->m * invDist3;
 

Original file line number	Diff line number	Diff line change
`@@ -115,8 +115,8 @@ inline void bhAccel(Octree* node, const Particle& p, real theta,`
`115`	`115`	`// Correct Barnes–Hut acceptance criterion`
`116`	`116`	`if (node->leaf \|\| (node->size / dist) < theta)`
`117`	`117`	`{`
`118`		`- real invDist = real(1) / dist;`
`119`		`- real invDist3 = invDist * invDist * invDist;`
	`118`	`+ real invDist = std::sqrt(real(1.0) / distSq);`
	`119`	`+real invDist3 = invDist * invDist * invDist;`
`120`	`120`
`121`	`121`	`real fac = G * node->m * invDist3;`
`122`	`122`