def MC_int(c, s, dist): x = np.empty([s], dtype=np.double, dist=dist) y = np.empty([s], dtype=np.double, dist=dist) sum=0.0 np.timer_reset() np.evalflush() start=time.time() for i in range(c): np.ufunc_random(x,x) np.ufunc_random(y,y) np.square(x,x) np.square(y,y) np.add(x,y,x) z = np.less_equal(x, 1) sum += np.add.reduce(z)*4.0/s sum = sum / c np.evalflush() stop=time.time() print 'Pi: ', sum, ' with ', s,' samples in sec: ', stop-start, if dist: print "(Dist) notes: %s"%sys.argv[4] else: print "(Non-Dist) notes: %s"%sys.argv[4]
def MC_int(c, s, dist): x = np.empty([s], dtype=np.double, dist=dist) y = np.empty([s], dtype=np.double, dist=dist) sum = 0.0 np.timer_reset() np.evalflush() start = time.time() for i in range(c): np.ufunc_random(x, x) np.ufunc_random(y, y) np.square(x, x) np.square(y, y) np.add(x, y, x) z = np.less_equal(x, 1) sum += np.add.reduce(z) * 4.0 / s sum = sum / c np.evalflush() stop = time.time() print 'Pi: ', sum, ' with ', s, ' samples in sec: ', stop - start, if dist: print "(Dist) notes: %s" % sys.argv[4] else: print "(Non-Dist) notes: %s" % sys.argv[4]
BOUND = np.zeros((nx,ny,nz), dtype=float, dist=DIST) BOUNDi = np.zeros((nx,ny,nz), dtype=float, dist=DIST) BOUNDi += 1#ones if not NO_OBST: for i in xrange(nx): for j in xrange(ny): for k in xrange(nz): if ((i-4)**2+(j-5)**2+(k-6)**2) < 6: BOUND[i,j,k] += 1.0 BOUNDi[i,j,k] += 0.0 BOUND[:,0,:] += 1.0 BOUNDi[:,0,:] *= 0.0 np.timer_reset() np.evalflush() t1 = time.time() while ts<ITER: ##Propagate / Streaming step T[:] = F #nearest-neighbours F[1,:,:,0] = T[1,:,:,-1] F[1,:,:,1:] = T[1,:,:,:-1] F[2,:,:,:-1] = T[2,:,:,1:] F[2,:,:,-1] = T[2,:,:,0] F[3,:,0,:] = T[3,:,-1,:] F[3,:,1:,:] = T[3,:,:-1,:] F[4,:,:-1,:] = T[4,:,1:,:] F[4,:,-1,:] = T[4,:,0,:] F[5,0,:,:] = T[5,-1,:,:] F[5,1:,:,:] = T[5,:-1,:,:]
seed = int(a) elif o in ("-f", "--file"): script_list.append(a) elif o in ("-e", "--exclude"): exclude_list.append(a) else: assert False, "unhandled option" if len(script_list) == 0: script_list = os.listdir(\ os.path.dirname(os.path.abspath(__file__))) if np.RANK == 0: print "*"*100 print "*"*31, "Testing Distributed Numerical Python", "*"*31 np.evalflush(barrier=True) for i in xrange(len(script_list)): f = script_list[i] if f.startswith("test_") and f.endswith("py")\ and f not in exclude_list: m = f[:-3]#Remove ".py" m = __import__(m) np.datalayout(None) np.evalflush(barrier=True) if np.RANK == 0: print "*"*100 print "Testing %s"%f np.evalflush(barrier=True) err = False msg = "" r1 = 0; r2 = 0