Commit 7384272f authored by Mulas, Giacomo's avatar Mulas, Giacomo
Browse files

bump up one more allocation from size configurations to nspheres

parent 71a70bde

src/libnptm/Configuration.cpp

+6 −6
Original line number Original line Diff line number Diff line
@@ -289,9 +289,9 @@ ScattererConfiguration::ScattererConfiguration(const ScattererConfiguration& rh 
  wp = rhs.wp;

  wp = rhs.wp;

  xip = rhs.xip;

  xip = rhs.xip;

  iog_vec = new int[number_of_spheres]();

  iog_vec = new int[number_of_spheres]();

  radii_of_spheres = new double[configurations]();

  radii_of_spheres = new double[number_of_spheres]();

  nshl_vec = new int[configurations]();

  nshl_vec = new int[number_of_spheres]();

  rcf = new double*[configurations];

  rcf = new double*[number_of_spheres];

  scale_vec = new double[number_of_scales]();

  scale_vec = new double[number_of_scales]();

  dc0_matrix = new dcomplex**[configurations];

  dc0_matrix = new dcomplex**[configurations];

  for (int si = 0; si < number_of_scales; si++) scale_vec[si] = rhs.scale_vec[si];

  for (int si = 0; si < number_of_scales; si++) scale_vec[si] = rhs.scale_vec[si];
@@ -639,9 +639,9 @@ ScattererConfiguration* ScattererConfiguration::from_hdf5(string file_name) { 
      if (iog[ci] < ci + 1) continue;

      if (iog[ci] < ci + 1) continue;

      configuration_count++;

      configuration_count++;

    }

    }

    nshl_vector = new int[configuration_count]();

    nshl_vector = new int[nsph]();

    ros_vector = new double[configuration_count]();

    ros_vector = new double[nsph]();

    rcf_vector = new double*[configuration_count];

    rcf_vector = new double*[nsph];

    for (int i115 = 1; i115 <= nsph; i115++) {

    for (int i115 = 1; i115 <= nsph; i115++) {

      if (iog[i115 - 1] < i115) continue;

      if (iog[i115 - 1] < i115) continue;

      str_name = "NSHL_" + to_string(iog[i115 - 1]);

      str_name = "NSHL_" + to_string(iog[i115 - 1]);