Merge pull request #640 from mperrin/fix_float_comparison

fix float comparison bug for non-square subarrays

poppy/poppy_core.py

@@ -1184,7 +1184,7 @@ def _propagate_mft(self, det):
 
         if not np.isscalar(self.pixelscale.value):
             # check for rectangular arrays
-            if self.pixelscale[0] == self.pixelscale[1]:
+            if np.isclose(self.pixelscale[0], self.pixelscale[1]):
                 self.pixelscale = self.pixelscale[0]
                 # we're in a rectangular array with same scale in both directions, so treat pixelscale as a scalar
             else: