Fixed bug in read_fastfood code when filename is supplied

Author: manodeep

Corrfunc/utils.py

@@ -129,7 +129,8 @@ def read_fastfood(filename, return_dtype=None, need_header=None):
                 input_dtype = np.float32 if skip1 // ngal == 4 else np.float
                 array = np.fromfile(f, input_dtype, ngal)
                 skip2 = struct.unpack(bytes_to_native_str(b'@i'), f.read(4))[0]
-                pos[field] = array if dtype is None else dtype(array)
+                pos[field] = array if return_dtype == input_dtype \
+                             else return_dtype(array)
 
         x = pos['x']
         y = pos['y']
@@ -181,7 +182,7 @@ def read_fastfood(filename, return_dtype=None, need_header=None):
             f = read_fastfood if '.ff' in extension else read_ascii
 
             # default return is double
-            x, y, z = f(filebase, np.float)
+            x, y, z = f(filebase, np.float32)
             return x, y, z
 
         raise IOError("Could not locate file {0}".format(filebase))