Move epsilon

Author: Sparks29032

src/diffpy/morph/morphs/morphrgrid.py

@@ -55,6 +55,7 @@ def morph(self, x_morph, y_morph, x_target, y_target):
         """Resample arrays onto specified grid."""
         Morph.morph(self, x_morph, y_morph, x_target, y_target)
         rmininc = max(self.x_target_in[0], self.x_morph_in[0])
+        print(self.x_morph_out)
         r_step_target = (self.x_target_in[-1] - self.x_target_in[0]) / (
             len(self.x_target_in) - 1
         )
@@ -66,6 +67,7 @@ def morph(self, x_morph, y_morph, x_target, y_target):
             self.x_target_in[-1] + r_step_target,
             self.x_morph_in[-1] + r_step_morph,
         )
+        print(rmininc, rmaxinc, rstepinc)
         if self.rmin is None or self.rmin < rmininc:
             self.rmin = rmininc
         if self.rmax is None or self.rmax > rmaxinc:
@@ -78,6 +80,7 @@ def morph(self, x_morph, y_morph, x_target, y_target):
         self.x_morph_out = numpy.arange(
             self.rmin, self.rmax - epsilon, self.rstep
         )
+        print((self.rmax - self.rmin) / self.rstep, len(self.x_morph_out))
         self.y_morph_out = numpy.interp(
             self.x_morph_out, self.x_morph_in, self.y_morph_in
         )

src/diffpy/morph/refine.py

@@ -100,6 +100,8 @@ def _residual(self, pvals):
                 diff_length = len(rvec) - self.res_length
                 rvec = list(rvec)
                 for _ in range(diff_length):
+                    # Remove the largest residuals since the larger overlapping
+                    # grid means we are likely closer to the target
                     rvec.pop(max(rvec))
                 rvec = array(rvec)