@@ -42,7 +42,7 @@ function relax!(F, UX, UY, UZ, nx, ny, nz, deltaU, t1D, t2D, t3D, sSQU, chunkid,
42
42
for l = 1 : size (F,4 )
43
43
density = density + F[i,j,k,l]
44
44
end
45
- fs = Array ( Float64, 6 )
45
+ fs = Vector { Float64} ( 6 )
46
46
for l = 1 : 6
47
47
fs[l] = 0.0
48
48
for m = 1 : 5
@@ -135,7 +135,7 @@ function perf_lattice_boltzmann(n)
135
135
136
136
CI = [0 : matsize: matsize* 19 ;]'
137
137
138
- BOUND = Array ( Float64, nx,ny,nz)
138
+ BOUND = Array { Float64} ( nx, ny, nz)
139
139
140
140
for i= 1 : nx, j= 1 : ny, k= 1 : nz
141
141
BOUND[i,j,k] = ((i- 5 )^ 2 + (j- 6 )^ 2 + (k- 7 )^ 2 ) < 6
@@ -149,14 +149,14 @@ function perf_lattice_boltzmann(n)
149
149
TO_REFLECT = [ON+ CI[2 ] ON+ CI[3 ] ON+ CI[4 ] ON+ CI[5 ] ON+ CI[6 ] ON+ CI[7 ] ON+ CI[8 ] ON+ CI[9 ] ON+ CI[10 ] ON+ CI[11 ] ON+ CI[12 ] ON+ CI[13 ] ON+ CI[14 ] ON+ CI[15 ] ON+ CI[16 ] ON+ CI[17 ] ON+ CI[18 ] ON+ CI[19 ]]
150
150
REFLECTED = [ON+ CI[3 ] ON+ CI[2 ] ON+ CI[5 ] ON+ CI[4 ] ON+ CI[7 ] ON+ CI[6 ] ON+ CI[11 ] ON+ CI[10 ] ON+ CI[9 ] ON+ CI[8 ] ON+ CI[15 ] ON+ CI[14 ] ON+ CI[13 ] ON+ CI[12 ] ON+ CI[19 ] ON+ CI[18 ] ON+ CI[17 ] ON+ CI[16 ]]
151
151
152
- UX = Array ( Float64, nx,ny,nz)
153
- UY = Array ( Float64, nx,ny,nz)
154
- UZ = Array ( Float64, nx,ny,nz)
155
- U = Array ( Float64, 12 ,nchunk)
156
- t1D = Array ( Float64, nx,ny,nz)
157
- t2D = Array ( Float64, nx,ny,nz)
158
- t3D = Array ( Float64, nx,ny,nz)
159
- sSQU = Array ( Float64, nx,ny,nz)
152
+ UX = Array { Float64} ( nx, ny, nz)
153
+ UY = Array { Float64} ( nx, ny, nz)
154
+ UZ = Array { Float64} ( nx, ny, nz)
155
+ U = Matrix { Float64} ( 12 , nchunk)
156
+ t1D = Array { Float64} ( nx, ny, nz)
157
+ t2D = Array { Float64} ( nx, ny, nz)
158
+ t3D = Array { Float64} ( nx, ny, nz)
159
+ sSQU = Array { Float64} ( nx, ny, nz)
160
160
161
161
avu = 1
162
162
prevavu = 1
0 commit comments