src/libnptm/Configuration.cpp

@@ -297,9 +297,10 @@ ScattererConfiguration::ScattererConfiguration(
   exdc = ex;
   wp = w;
   xip = x;
-  if (variable_name == "XIV") scale_vec = scale_vector;
-  else {
-    scale_vec = new double[number_of_scales]();
+  scale_vec = new double[number_of_scales]();
+  if (variable_name == "XIV") {
+    for (int xi = 0; xi < nxi; xi++) scale_vec[xi] = scale_vector[xi];
+  } else {
     const double pi2 = 2.0 * acos(-1.0);
     const double evc = 6.5821188e-16;
     for (int si = 0; si < number_of_scales; si++) {
@@ -312,6 +313,7 @@ ScattererConfiguration::ScattererConfiguration(
     }
   }
 }
+
 ScattererConfiguration::ScattererConfiguration(const  ScattererConfiguration& rhs)
 {
   number_of_spheres = rhs.number_of_spheres;
@@ -642,6 +644,7 @@ ScattererConfiguration* ScattererConfiguration::from_dedfb(string dedfb_file_nam
 							      _xip
   );
   delete[] file_lines;
+  delete[] variable_vector;
   return config;
 }
 

src/sphere/sphere.cpp

@@ -225,7 +225,7 @@ void sphere(string config_file, string data_file, string output_path) {
 	  c1->rc[i116][j115] = sconf->get_rcf(i116, j115) * c1->ros[i116];
 	}
       }
-      gcs += c1->gcsv[iogi];
+      gcs += c1->gcsv[iogi - 1];
     }
     double ****zpv = new double***[l_max]; //[l_max][3][2][2]; // Matrix: dim[LM x 3 x 2 x 2]
     for (int zi = 0; zi < l_max; zi++) {