fixed memory leaks

Author: mlietzow

src/MathFunctions_t.cpp

@@ -67,10 +67,7 @@ TEST_CASE("CMathFunctions::calcWVMie::BH", "[MathFunctions][CMathFunctions]")
     dcomplex ri {1.55, 0.0};
     double x {PIx2 * 0.525 / 0.6328};
     double qext, qabs, qsca, gsca {};
-    double* S11 = new double[nr_theta];
-    double* S12 = new double[nr_theta];
-    double* S33 = new double[nr_theta];
-    double* S34 = new double[nr_theta];
+    double S11[nr_theta], S12[nr_theta], S33[nr_theta], S34[nr_theta];
 
     dlist theta(nr_theta);
     for(uint th = 0; th < nr_theta; th++) {
@@ -121,10 +118,7 @@ TEST_CASE("CMathFunctions::calcWVMie::W1", "[MathFunctions][CMathFunctions]")
     dcomplex ri {1.5, 0.0};
     auto x = GENERATE(10.0, 100.0, 1000.0, 5000.0);
     double qext, qabs, qsca, gsca {};
-    double *S11 = new double[nr_theta];
-    double *S12 = new double[nr_theta];
-    double *S33 = new double[nr_theta];
-    double *S34 = new double[nr_theta];
+    double S11[nr_theta], S12[nr_theta], S33[nr_theta], S34[nr_theta];
 
     dlist theta(nr_theta);
     for(uint th = 0; th < nr_theta; th++) {
@@ -232,10 +226,7 @@ TEST_CASE("CMathFunctions::calcWVMie::W2", "[MathFunctions][CMathFunctions]")
     dcomplex ri {1.5, 0.1};
     auto x = GENERATE(10.0, 100.0, 1000.0, 5000.0);
     double qext, qabs, qsca, gsca {};
-    double *S11 = new double[nr_theta];
-    double *S12 = new double[nr_theta];
-    double *S33 = new double[nr_theta];
-    double *S34 = new double[nr_theta];
+    double S11[nr_theta], S12[nr_theta], S33[nr_theta], S34[nr_theta];
 
     dlist theta(nr_theta);
     for(uint th = 0; th < nr_theta; th++) {

src/Pipeline.cpp

@@ -1610,7 +1610,7 @@ pos_1.normalized();
         cell_1->updateData(data_off + 1, -F1.Y());
         cell_1->updateData(data_off + 2, -F1.Z());
 
-        /*for(long c_2 = c_1 + 1; c_2 < long(max_cells - 1); c_2++)
+        for(long c_2 = c_1 + 1; c_2 < long(max_cells - 1); c_2++)
         {
         cell_oc * cell_2 = (cell_oc*) grid->getCellFromIndex(c_2);
         double dens_2 = grid->getGasNumberDensity(cell_2);