Minor bug in scattered interploant. Results of comparison between con…

…ts energy cuts -- Kun vs real vs Martin
abuts · Mar 2, 2020 · f38c3ce · f38c3ce
1 parent b80df4b
commit f38c3ce
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 3 deletions.
diff --git a/2020April/Working2018.pptx → 2020April/Working2020.pptx b/2020April/Working2018.pptx → 2020April/Working2020.pptx
diff --git a/Data/DFT/build_ScattInt.m b/Data/DFT/build_ScattInt.m
@@ -1,16 +1,20 @@
 function  Interp_array = build_ScattInt(en_pts,qx,qy,qz,en,ses)
 
-nIp = numel(en_pts);
+if all(size(en_pts) == size(qx))
+    en_pts = unique(sort(en_pts));
+    nIp = numel(en_pts);
+else
+    nIp = numel(en_pts);
+end
 if size(ses,2) > 1
     binned = true;
 else
     binned = false;
 end
 Interp_array = cell(nIp,1);
 bin_size = 8;
-e_min = min(en_pts);
 
-e_ind = floor((en-e_min)/bin_size)+1;
+e_ind = floor(en/bin_size);
 n_empties = 0;
 n_binv = max(e_ind);
 for i=1:n_binv