@@ -13,9 +13,9 @@ using Plots, Printf
13
13
macro qx (ix, iy) esc (:(- D_dx * (C[$ ix+ 1 , $ iy+ 1 ] - C[$ ix, $ iy+ 1 ]))) end
14
14
macro qy (ix, iy) esc (:(- D_dy * (C[$ ix+ 1 , $ iy+ 1 ] - C[$ ix+ 1 , $ iy]))) end
15
15
16
- @parallel_indices (ix, iy) function compute! (C2, C, D_dx, D_dy, dt, _dx, _dy, size_C1_2, size_C2_2)
16
+ @parallel_indices (ix, iy) function compute! (C2, C, D_dx, D_dy, dt, size_C1_2, size_C2_2)
17
17
if (ix <= size_C1_2 && iy <= size_C2_2)
18
- C2[ix+ 1 , iy+ 1 ] = C[ix+ 1 , iy+ 1 ] - dt * ((@qx (ix + 1 , iy) - @qx (ix, iy)) * _dx + (@qy (ix, iy + 1 ) - @qy (ix, iy)) * _dy )
18
+ C2[ix+ 1 , iy+ 1 ] = C[ix+ 1 , iy+ 1 ] - dt * ((@qx (ix + 1 , iy) - @qx (ix, iy)) * D_dx + (@qy (ix, iy + 1 ) - @qy (ix, iy)) * D_dy )
19
19
end
20
20
return
21
21
end
35
35
xc, yc = LinRange (dx / 2 , Lx - dx / 2 , nx), LinRange (dy / 2 , Ly - dy / 2 , ny)
36
36
D_dx = D / dx
37
37
D_dy = D / dy
38
- _dx, _dy = 1.0 / dx, 1.0 / dy
39
38
# Array initialisation
40
39
C = Data. Array (exp .(.- (xc .- Lx / 2 ) .^ 2 .- (yc' .- Ly / 2 ) .^ 2 ))
41
40
C2 = copy (C)
49
48
# Time loop
50
49
for it = 1 : nt
51
50
if (it == 11 ) t_tic = Base. time (); niter = 0 end
52
- @parallel compute! (C2, C, D_dx, D_dy, dt, _dx, _dy, size_C1_2, size_C2_2)
51
+ @parallel compute! (C2, C, D_dx, D_dy, dt, size_C1_2, size_C2_2)
53
52
C, C2 = C2, C # pointer swap
54
53
niter += 1
55
54
if do_visu && (it % nout == 0 )
0 commit comments