from time import time import numpy as np import refine size=420 a=np.arange(size**3).astype('f8').reshape((size,size,size)) ti=time() b=refine.refine(a,2,2,2) tf=time() dt=tf-ti print a print b print dt raw_input("")