Remove: analyze diffusion_profile.

Closes: espressomd#271
RudolfWeeber · Nov 21, 2017 · 079c95a · 079c95a
1 parent 1d03321
commit 079c95a
Show file tree

Hide file tree

Showing 2 changed files with 0 additions and 66 deletions.
diff --git a/src/core/statistics.cpp b/src/core/statistics.cpp
@@ -912,68 +912,6 @@ void density_profile_av(PartCfg &partCfg, int n_conf, int n_bin, double density,
     rho_ave[i] /= n_conf;
 }
 
-void calc_diffusion_profile(PartCfg &partCfg, int dir, double xmin, double xmax,
-                            int nbins, int n_part, int n_conf, int time,
-                            int type, double *bins) {
-  int i, t, count, index;
-  double tcount = 0;
-  double xpos;
-  double tpos[3];
-  int img_box[3] = {0, 0, 0};
-  // double delta_x = (box_l[0])/((double) nbins);
-
-  /* create and initialize the array of bins */
-
-  // double *bins;
-
-  int *label;
-  label = (int *)Utils::malloc(n_part * sizeof(int));
-
-  /* calculation over last n_conf configurations */
-  t = n_configs - n_conf;
-
-  while (t < n_configs - time) {
-    /* check initial condition */
-    count = 0;
-
-    for (i = 0; i < n_part; i++) {
-      if (partCfg[i].p.type == type) {
-        tpos[0] = configs[t][3 * i];
-        tpos[1] = configs[t][3 * i + 1];
-        tpos[2] = configs[t][3 * i + 2];
-        fold_coordinate(tpos, img_box, dir);
-        xpos = tpos[dir];
-        if (xpos > xmin && xpos < xmax) {
-          label[count] = i;
-        } else
-          label[count] = -1;
-        count++;
-      }
-    }
-
-    /* check at time 'time' */
-    for (i = 0; i < n_part; i++) {
-      if (label[i] > 0) {
-        tpos[0] = configs[t + time][3 * label[i]];
-        tpos[1] = configs[t + time][3 * label[i] + 1];
-        tpos[2] = configs[t + time][3 * label[i] + 2];
-        fold_coordinate(tpos, img_box, dir);
-        xpos = tpos[dir];
-
-        index = (int)(xpos / box_l[dir] * nbins);
-        bins[index]++;
-      }
-    }
-    t++;
-    tcount++;
-  }
-
-  /* normalization */
-  for (i = 0; i < nbins; i++) {
-    bins[i] = bins[i] / (tcount);
-  }
-  free(label);
-}
 
 int calc_cylindrical_average(
     PartCfg &partCfg, std::vector<double> center_,

diff --git a/src/core/statistics.hpp b/src/core/statistics.hpp
@@ -355,10 +355,6 @@ int calc_radial_density_map(int xbins, int ybins, int thetabins, double xrange,
                             IntList *beadids, DoubleList *density_map,
                             DoubleList *density_profile);
 
-void calc_diffusion_profile(int dir, double xmin, double xmax, int nbins,
-                            int n_part, int n_conf, int time, int type,
-                            double *bins);
-
 /** returns the minimal squared distance between two positions in the perhaps
  periodic
     simulation box.