def get_vorticity(particle,cluster):
particle.wi = arange(particle.ni)*1.
for ic in range(cluster.n):
particle,cluster = get_trunc(particle,cluster,ic)
ista = cluster.ista[ic]
iend = cluster.iend[ic]
for ip in range(ista,iend+1):
dx = particle.xi[ip]-cluster.xjt
dy = particle.yi[ip]-cluster.yjt
particle.wi[ip] = dot(cluster.gjt,exp(-(dx**2+dy**2)/\
(2*particle.sigma**2))/(2*pi*particle.sigma**2))
return particle,cluster
t[0] = tic-toc
# solve rbf using c wrapper
if parameter.wrapper == 1:
rbf_solver(particle.xi,particle.yi,particle.gj,particle.ei,
particle.sigma,parameter.vis)
toc = tic
tic = cpu_time()
# generate clusters
particle,cluster = get_cluster(particle,cluster)
# set up reference data
if parameter.wrapper == 1:
particle,cluster = get_trunc(particle,cluster)
particle,cluster = get_vorticity(particle,cluster)
particle.gj = exp(-(particle.xj**2+particle.yj**2)/\
(4*parameter.vis*parameter.t))/(pi*4*parameter.vis*parameter.t)*particle.h**2
particle.ei = particle.wi
particle.wi = particle.gj/particle.h**2
toc = tic
tic = cpu_time()
t[0] += tic-toc
# RBF interpolation
it = -1; iconv = 0; err = []; grid.r = []
while iconv < 5:
it = it+1
