softening now a single scalar for every particle

lib/nBodySim.cpp

@@ -29,9 +29,8 @@ nBodySim::nBodySim(std::string datafile) {
         }
     }
 
-    softening_ = new double[nParticles_];
     for (unsigned i=0; i < nParticles_; ++i) {
-        file >> softening_[i];
+        file >> softening_;
     }
     
     potential_ = new double[nParticles_];
@@ -52,7 +51,6 @@ nBodySim::~nBodySim() {
     delete[] masses_;
     delete[] positions_;
     delete[] velocities_;
-    delete[] softening_;
     delete[] potential_;
     delete[] forces_;
 }
@@ -223,7 +221,7 @@ void nBodySim::saveState2file(unsigned step, std::string filename) const {
     }
     // save softening
     for (unsigned i=0; i < nParticles_; ++i) {
-        file << softening_[i] << std::endl;
+        file << softening_ << std::endl;
     }
     // save potential
     for (unsigned i=0; i < nParticles_; ++i) {

lib/nBodySim.hpp

@@ -37,7 +37,7 @@ private:
     double* masses_;
     double* positions_;
     double* velocities_;
-    double* softening_;
+    double softening_;
     double* potential_;
     double* forces_;
     double meanInterparticleDistance_;