max performance

lib/nBodySim.cpp

 #include "nBodySim.hpp"
 #include "tree.hpp"
+#include "utils.hpp"
 #include <omp.h>
 #include <cmath>
 #include <string>
@@ -74,13 +75,14 @@ void nBodySim::runSimulation(double dt, unsigned nSteps) {
         write2file(positions_, positionsFile, 3, skip);
         energyFile << calculateTotalEnergy() << std::endl;
 
-        Vector COM = getCenterOfMass();
+        // Vector COM = getCenterOfMass();
 
         // Integrate with leapfrog
         auto start = std::chrono::high_resolution_clock::now();
         kick(dt/2.0);
         drift(dt);
         treeUpdateForces();
+        // double dev = verifyForces();
         kick(dt/2.0);
         updateTree();
         auto stop = std::chrono::high_resolution_clock::now();
@@ -88,7 +90,7 @@ void nBodySim::runSimulation(double dt, unsigned nSteps) {
 
         std::cout << std::right << "Step "    << std::setw(5)  << i
                                 // << ": MFM = " << std::setw(12) << calculateMeanForceMagnitude()
-                                // << ", DEV = " << std::setw(12) << verifyForces()
+                                // << ", DEV = " << std::setw(12) << dev
                                 // << ", MCD = " << std::setw(8)  << calculateMeanCenterDistance()
                                 // << ", MID = " << std::setw(8)  << calculateMeanInterparticleDistance()
                                 // << ", COM = " << "(" << COM[0] << ", " << COM[1] << ", " << COM[2] << ")"

main.cpp

@@ -18,7 +18,7 @@ int main(int argc, char** argv) {
             return 1;
     }
 
-    double dt = 0.000005;
+    double dt = 5e-7;
     unsigned nSteps = 240;
 
     std::cout << "dt = " << dt << std::endl;