Commit d27204ac authored by ph's avatar ph

add transport solver

parent d54bcd0b
...@@ -14,8 +14,8 @@ wbord = 0. ...@@ -14,8 +14,8 @@ wbord = 0.
# grid size # grid size
nx = 50 nx = 512
ny = 50 ny = 512
Lx = 1. Lx = 1.
Ly = 1. Ly = 1.
...@@ -148,16 +148,12 @@ while t < Tmax: ...@@ -148,16 +148,12 @@ while t < Tmax:
plt.colorbar() plt.colorbar()
plt.pause(0.01) plt.pause(0.01)
# in case of an odd number of iterations, perform a last exchange
#if iter % 2 != 0:
# wn_gpu, wnp1_gpu = wnp1_gpu, wn_gpu
# prg.K2(queue, (taille, ), None, y_gpu).wait() # prg.K2(queue, (taille, ), None, y_gpu).wait()
# prg.K3(queue, (taille, ), None , x_gpu, y_gpu, z_gpu, wait_for = [event]).wait() # prg.K3(queue, (taille, ), None , x_gpu, y_gpu, z_gpu, wait_for = [event]).wait()
#
cl.enqueue_copy(queue, wn_cpu, wn_gpu).wait() cl.enqueue_copy(queue, wn_cpu, wn_gpu).wait()
#
wplot = np.reshape(wn_cpu,(nx, ny)) wplot = np.reshape(wn_cpu,(nx, ny))
plt.clf() plt.clf()
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment