Added a few comments

lib/nBodySim.cpp

@@ -6,7 +6,6 @@
 #include <iomanip>
 
 nBodySim::nBodySim(std::string& datafile) {
-    // reads in data from datafile
     std::string dummy;
     std::ifstream file(datafile);
     file >> dummy >> dummy >> N_;
@@ -41,6 +40,7 @@ nBodySim::nBodySim(std::string& datafile) {
     }
 
     forces_ = new double[3*N_];
+    #pragma omp parallel for
     for (unsigned i=0; i < 3*N_; ++i) {
         forces_[i] = 0.0;
     }

lib/nBodySim.hpp

@@ -5,14 +5,16 @@
 
 class nBodySim {
 public:
-    // constructor & destructor
+    // constructor reads in data from datafile
     nBodySim(std::string& datafile);
+    // destructor deletes arrays 
     ~nBodySim();
-    // methods
+    // loops over all pairs of particles and calculates forces
     void calculateForces();
     double* getMasses() { return masses_; }
     double* getPositions() { return positions_; }
     double* getForces() { return forces_; }
+    // writes data to file, every row is a particle
     void write2file(const double* array, std::string filename, unsigned dim) const;
     unsigned getN() const { return N_; }
 private: