Neko 1.99.1
A portable framework for high-order spectral element flow simulations
Loading...
Searching...
No Matches
sx_conv1.f90
Go to the documentation of this file.
1! Copyright (c) 2021-2024, The Neko Authors
2! All rights reserved.
3!
4! Redistribution and use in source and binary forms, with or without
5! modification, are permitted provided that the following conditions
6! are met:
7!
8! * Redistributions of source code must retain the above copyright
9! notice, this list of conditions and the following disclaimer.
10!
11! * Redistributions in binary form must reproduce the above
12! copyright notice, this list of conditions and the following
13! disclaimer in the documentation and/or other materials provided
14! with the distribution.
15!
16! * Neither the name of the authors nor the names of its
17! contributors may be used to endorse or promote products derived
18! from this software without specific prior written permission.
19!
20! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
21! "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
22! LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
23! FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
24! COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
25! INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
26! BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
27! LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
28! CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
29! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
30! ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
31! POSSIBILITY OF SUCH DAMAGE.
32!
34submodule(opr_sx) sx_conv1
35 implicit none
36
37contains
38
39 module subroutine opr_sx_conv1(du, u, vx, vy, vz, xh, coef, nelv)
40 type(space_t), intent(in) :: Xh
41 type(coef_t), intent(in) :: coef
42 integer, intent(in) :: nelv
43 real(kind=rp), intent(inout) :: du(xh%lxyz, nelv)
44 real(kind=rp), intent(inout) :: u(xh%lx, xh%ly, xh%lz, nelv)
45 real(kind=rp), intent(inout) :: vx(xh%lx, xh%ly, xh%lz, nelv)
46 real(kind=rp), intent(inout) :: vy(xh%lx, xh%ly, xh%lz, nelv)
47 real(kind=rp), intent(inout) :: vz(xh%lx, xh%ly, xh%lz, nelv)
48
49 select case (xh%lx)
50 case (14)
51 call sx_conv1_lx14(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
52 coef%drdx, coef%dsdx, coef%dtdx, &
53 coef%drdy, coef%dsdy, coef%dtdy, &
54 coef%drdz, coef%dsdz, coef%dtdz, &
55 coef%jacinv, nelv)
56 case (13)
57 call sx_conv1_lx13(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
58 coef%drdx, coef%dsdx, coef%dtdx, &
59 coef%drdy, coef%dsdy, coef%dtdy, &
60 coef%drdz, coef%dsdz, coef%dtdz, &
61 coef%jacinv, nelv)
62 case (12)
63 call sx_conv1_lx12(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
64 coef%drdx, coef%dsdx, coef%dtdx, &
65 coef%drdy, coef%dsdy, coef%dtdy, &
66 coef%drdz, coef%dsdz, coef%dtdz, &
67 coef%jacinv, nelv)
68 case (11)
69 call sx_conv1_lx11(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
70 coef%drdx, coef%dsdx, coef%dtdx, &
71 coef%drdy, coef%dsdy, coef%dtdy, &
72 coef%drdz, coef%dsdz, coef%dtdz, &
73 coef%jacinv, nelv)
74 case (10)
75 call sx_conv1_lx10(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
76 coef%drdx, coef%dsdx, coef%dtdx, &
77 coef%drdy, coef%dsdy, coef%dtdy, &
78 coef%drdz, coef%dsdz, coef%dtdz, &
79 coef%jacinv, nelv)
80 case (9)
81 call sx_conv1_lx9(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
82 coef%drdx, coef%dsdx, coef%dtdx, &
83 coef%drdy, coef%dsdy, coef%dtdy, &
84 coef%drdz, coef%dsdz, coef%dtdz, &
85 coef%jacinv, nelv)
86 case (8)
87 call sx_conv1_lx8(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
88 coef%drdx, coef%dsdx, coef%dtdx, &
89 coef%drdy, coef%dsdy, coef%dtdy, &
90 coef%drdz, coef%dsdz, coef%dtdz, &
91 coef%jacinv, nelv)
92 case (7)
93 call sx_conv1_lx7(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
94 coef%drdx, coef%dsdx, coef%dtdx, &
95 coef%drdy, coef%dsdy, coef%dtdy, &
96 coef%drdz, coef%dsdz, coef%dtdz, &
97 coef%jacinv, nelv)
98 case (6)
99 call sx_conv1_lx6(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
100 coef%drdx, coef%dsdx, coef%dtdx, &
101 coef%drdy, coef%dsdy, coef%dtdy, &
102 coef%drdz, coef%dsdz, coef%dtdz, &
103 coef%jacinv, nelv)
104 case (5)
105 call sx_conv1_lx5(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
106 coef%drdx, coef%dsdx, coef%dtdx, &
107 coef%drdy, coef%dsdy, coef%dtdy, &
108 coef%drdz, coef%dsdz, coef%dtdz, &
109 coef%jacinv, nelv)
110 case (4)
111 call sx_conv1_lx4(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
112 coef%drdx, coef%dsdx, coef%dtdx, &
113 coef%drdy, coef%dsdy, coef%dtdy, &
114 coef%drdz, coef%dsdz, coef%dtdz, &
115 coef%jacinv, nelv)
116 case (3)
117 call sx_conv1_lx3(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
118 coef%drdx, coef%dsdx, coef%dtdx, &
119 coef%drdy, coef%dsdy, coef%dtdy, &
120 coef%drdz, coef%dsdz, coef%dtdz, &
121 coef%jacinv, nelv)
122 case (2)
123 call sx_conv1_lx2(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
124 coef%drdx, coef%dsdx, coef%dtdx, &
125 coef%drdy, coef%dsdy, coef%dtdy, &
126 coef%drdz, coef%dsdz, coef%dtdz, &
127 coef%jacinv, nelv)
128 case default
129 call sx_conv1_lx(du, u, vx, vy, vz, xh%dx, xh%dy, xh%dz, &
130 coef%drdx, coef%dsdx, coef%dtdx, &
131 coef%drdy, coef%dsdy, coef%dtdy, &
132 coef%drdz, coef%dsdz, coef%dtdz, &
133 coef%jacinv, nelv, xh%lx)
134 end select
135
136 end subroutine opr_sx_conv1
137
138 subroutine sx_conv1_lx(du, u, vx, vy, vz, dx, dy, dz, &
139 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
140 jacinv, nelv, lx)
141 integer, intent(in) :: nelv, lx
142 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
143 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
144 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
145 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
146 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
147 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
148 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
149 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
150 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
151 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
152 real(kind=rp) :: wr, ws, wt
153 integer :: e, i, j, k, jj, kk
154
155 do i = 1, lx
156 do jj = 1, lx * lx * nelv
157 wr = 0d0
158 do kk = 1, lx
159 wr = wr + dx(i, kk) * u(kk, jj,1,1)
160 end do
161 dudr(i, jj,1,1) = wr
162 end do
163 end do
164
165 do k = 1, lx
166 do i = 1, lx
167 do j = 1, lx
168 do e = 1, nelv
169 ws = 0d0
170 !NEC$ unroll_completely
171 do kk = 1, lx
172 ws = ws + dy(j, kk) * u(i, kk,k,e)
173 end do
174 duds(i,j,k,e) = ws
175 end do
176 end do
177 end do
178 end do
179
180 do j = 1, lx
181 do i = 1, lx
182 do k = 1, lx
183 do e = 1, nelv
184 wt = 0d0
185 !NEC$ unroll_completely
186 do kk = 1, lx
187 wt = wt + dz(k, kk) * u(i,j, kk,e)
188 end do
189 dudt(i,j,k,e) = wt
190 end do
191 end do
192 end do
193 end do
194
195 do i = 1, nelv * lx * lx * lx
196 du(i,1,1,1) = jacinv(i,1,1,1) * &
197 (vx(i,1,1,1) * ( &
198 drdx(i,1,1,1)*dudr(i,1,1,1) &
199 + dsdx(i,1,1,1)*duds(i,1,1,1) &
200 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
201 + vy(i,1,1,1)*( &
202 drdy(i,1,1,1)*dudr(i,1,1,1) &
203 + dsdy(i,1,1,1)*duds(i,1,1,1) &
204 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
205 + vz(i,1,1,1)*( &
206 drdz(i,1,1,1)*dudr(i,1,1,1) &
207 + dsdz(i,1,1,1)*duds(i,1,1,1) &
208 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
209 end do
210
211 end subroutine sx_conv1_lx
212
213 subroutine sx_conv1_lx14(du, u, vx, vy, vz, dx, dy, dz, &
214 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
215 jacinv, nelv)
216 integer, parameter :: lx = 14
217 integer, intent(in) :: nelv
218 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
219 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
220 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
221 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
222 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
223 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
224 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
225 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
226 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
227 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
228 real(kind=rp) :: wr, ws, wt
229 integer :: e, i, j, k, jj, kk
230
231 do i = 1, lx
232 do jj = 1, lx * lx * nelv
233 wr = 0d0
234 do kk = 1, lx
235 wr = wr + dx(i, kk) * u(kk, jj,1,1)
236 end do
237 dudr(i, jj,1,1) = wr
238 end do
239 end do
240
241 do k = 1, lx
242 do i = 1, lx
243 do j = 1, lx
244 do e = 1, nelv
245 ws = 0d0
246 !NEC$ unroll_completely
247 do kk = 1, lx
248 ws = ws + dy(j, kk) * u(i, kk,k,e)
249 end do
250 duds(i,j,k,e) = ws
251 end do
252 end do
253 end do
254 end do
255
256 do j = 1, lx
257 do i = 1, lx
258 do k = 1, lx
259 do e = 1, nelv
260 wt = 0d0
261 !NEC$ unroll_completely
262 do kk = 1, lx
263 wt = wt + dz(k, kk) * u(i,j, kk,e)
264 end do
265 dudt(i,j,k,e) = wt
266 end do
267 end do
268 end do
269 end do
270
271 do i = 1, nelv * lx * lx * lx
272 du(i,1,1,1) = jacinv(i,1,1,1) * &
273 (vx(i,1,1,1) * ( &
274 drdx(i,1,1,1)*dudr(i,1,1,1) &
275 + dsdx(i,1,1,1)*duds(i,1,1,1) &
276 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
277 + vy(i,1,1,1)*( &
278 drdy(i,1,1,1)*dudr(i,1,1,1) &
279 + dsdy(i,1,1,1)*duds(i,1,1,1) &
280 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
281 + vz(i,1,1,1)*( &
282 drdz(i,1,1,1)*dudr(i,1,1,1) &
283 + dsdz(i,1,1,1)*duds(i,1,1,1) &
284 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
285 end do
286
287 end subroutine sx_conv1_lx14
288
289 subroutine sx_conv1_lx13(du, u, vx, vy, vz, dx, dy, dz, &
290 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
291 jacinv, nelv)
292 integer, parameter :: lx = 13
293 integer, intent(in) :: nelv
294 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
295 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
296 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
297 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
298 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
299 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
300 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
301 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
302 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
303 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
304 real(kind=rp) :: wr, ws, wt
305 integer :: e, i, j, k, jj, kk
306
307 do i = 1, lx
308 do jj = 1, lx * lx * nelv
309 wr = 0d0
310 do kk = 1, lx
311 wr = wr + dx(i, kk) * u(kk, jj,1,1)
312 end do
313 dudr(i, jj,1,1) = wr
314 end do
315 end do
316
317 do k = 1, lx
318 do i = 1, lx
319 do j = 1, lx
320 do e = 1, nelv
321 ws = 0d0
322 !NEC$ unroll_completely
323 do kk = 1, lx
324 ws = ws + dy(j, kk) * u(i, kk,k,e)
325 end do
326 duds(i,j,k,e) = ws
327 end do
328 end do
329 end do
330 end do
331
332 do j = 1, lx
333 do i = 1, lx
334 do k = 1, lx
335 do e = 1, nelv
336 wt = 0d0
337 !NEC$ unroll_completely
338 do kk = 1, lx
339 wt = wt + dz(k, kk) * u(i,j, kk,e)
340 end do
341 dudt(i,j,k,e) = wt
342 end do
343 end do
344 end do
345 end do
346
347 do i = 1, nelv * lx * lx * lx
348 du(i,1,1,1) = jacinv(i,1,1,1) * &
349 (vx(i,1,1,1) * ( &
350 drdx(i,1,1,1)*dudr(i,1,1,1) &
351 + dsdx(i,1,1,1)*duds(i,1,1,1) &
352 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
353 + vy(i,1,1,1)*( &
354 drdy(i,1,1,1)*dudr(i,1,1,1) &
355 + dsdy(i,1,1,1)*duds(i,1,1,1) &
356 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
357 + vz(i,1,1,1)*( &
358 drdz(i,1,1,1)*dudr(i,1,1,1) &
359 + dsdz(i,1,1,1)*duds(i,1,1,1) &
360 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
361 end do
362
363 end subroutine sx_conv1_lx13
364
365 subroutine sx_conv1_lx12(du, u, vx, vy, vz, dx, dy, dz, &
366 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
367 jacinv, nelv)
368 integer, parameter :: lx = 12
369 integer, intent(in) :: nelv
370 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
371 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
372 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
373 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
374 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
375 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
376 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
377 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
378 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
379 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
380 real(kind=rp) :: wr, ws, wt
381 integer :: e, i, j, k, jj, kk
382
383 do i = 1, lx
384 do jj = 1, lx * lx * nelv
385 wr = 0d0
386 do kk = 1, lx
387 wr = wr + dx(i, kk) * u(kk, jj,1,1)
388 end do
389 dudr(i, jj,1,1) = wr
390 end do
391 end do
392
393 do k = 1, lx
394 do i = 1, lx
395 do j = 1, lx
396 do e = 1, nelv
397 ws = 0d0
398 !NEC$ unroll_completely
399 do kk = 1, lx
400 ws = ws + dy(j, kk) * u(i, kk,k,e)
401 end do
402 duds(i,j,k,e) = ws
403 end do
404 end do
405 end do
406 end do
407
408 do j = 1, lx
409 do i = 1, lx
410 do k = 1, lx
411 do e = 1, nelv
412 wt = 0d0
413 !NEC$ unroll_completely
414 do kk = 1, lx
415 wt = wt + dz(k, kk) * u(i,j, kk,e)
416 end do
417 dudt(i,j,k,e) = wt
418 end do
419 end do
420 end do
421 end do
422
423 do i = 1, nelv * lx * lx * lx
424 du(i,1,1,1) = jacinv(i,1,1,1) * &
425 (vx(i,1,1,1) * ( &
426 drdx(i,1,1,1)*dudr(i,1,1,1) &
427 + dsdx(i,1,1,1)*duds(i,1,1,1) &
428 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
429 + vy(i,1,1,1)*( &
430 drdy(i,1,1,1)*dudr(i,1,1,1) &
431 + dsdy(i,1,1,1)*duds(i,1,1,1) &
432 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
433 + vz(i,1,1,1)*( &
434 drdz(i,1,1,1)*dudr(i,1,1,1) &
435 + dsdz(i,1,1,1)*duds(i,1,1,1) &
436 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
437 end do
438
439 end subroutine sx_conv1_lx12
440
441 subroutine sx_conv1_lx11(du, u, vx, vy, vz, dx, dy, dz, &
442 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
443 jacinv, nelv)
444 integer, parameter :: lx = 11
445 integer, intent(in) :: nelv
446 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
447 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
448 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
449 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
450 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
451 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
452 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
453 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
454 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
455 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
456 real(kind=rp) :: wr, ws, wt
457 integer :: e, i, j, k, jj, kk
458
459 do i = 1, lx
460 do jj = 1, lx * lx * nelv
461 wr = 0d0
462 do kk = 1, lx
463 wr = wr + dx(i, kk) * u(kk, jj,1,1)
464 end do
465 dudr(i, jj,1,1) = wr
466 end do
467 end do
468
469 do k = 1, lx
470 do i = 1, lx
471 do j = 1, lx
472 do e = 1, nelv
473 ws = 0d0
474 !NEC$ unroll_completely
475 do kk = 1, lx
476 ws = ws + dy(j, kk) * u(i, kk,k,e)
477 end do
478 duds(i,j,k,e) = ws
479 end do
480 end do
481 end do
482 end do
483
484 do j = 1, lx
485 do i = 1, lx
486 do k = 1, lx
487 do e = 1, nelv
488 wt = 0d0
489 !NEC$ unroll_completely
490 do kk = 1, lx
491 wt = wt + dz(k, kk) * u(i,j, kk,e)
492 end do
493 dudt(i,j,k,e) = wt
494 end do
495 end do
496 end do
497 end do
498
499 do i = 1, nelv * lx * lx * lx
500 du(i,1,1,1) = jacinv(i,1,1,1) * &
501 (vx(i,1,1,1) * ( &
502 drdx(i,1,1,1)*dudr(i,1,1,1) &
503 + dsdx(i,1,1,1)*duds(i,1,1,1) &
504 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
505 + vy(i,1,1,1)*( &
506 drdy(i,1,1,1)*dudr(i,1,1,1) &
507 + dsdy(i,1,1,1)*duds(i,1,1,1) &
508 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
509 + vz(i,1,1,1)*( &
510 drdz(i,1,1,1)*dudr(i,1,1,1) &
511 + dsdz(i,1,1,1)*duds(i,1,1,1) &
512 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
513 end do
514
515 end subroutine sx_conv1_lx11
516
517 subroutine sx_conv1_lx10(du, u, vx, vy, vz, dx, dy, dz, &
518 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
519 jacinv, nelv)
520 integer, parameter :: lx = 10
521 integer, intent(in) :: nelv
522 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
523 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
524 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
525 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
526 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
527 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
528 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
529 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
530 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
531 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
532 real(kind=rp) :: wr, ws, wt
533 integer :: e, i, j, k, jj, kk
534
535 do i = 1, lx
536 do jj = 1, lx * lx * nelv
537 wr = 0d0
538 do kk = 1, lx
539 wr = wr + dx(i, kk) * u(kk, jj,1,1)
540 end do
541 dudr(i, jj,1,1) = wr
542 end do
543 end do
544
545 do k = 1, lx
546 do i = 1, lx
547 do j = 1, lx
548 do e = 1, nelv
549 ws = 0d0
550 !NEC$ unroll_completely
551 do kk = 1, lx
552 ws = ws + dy(j, kk) * u(i, kk,k,e)
553 end do
554 duds(i,j,k,e) = ws
555 end do
556 end do
557 end do
558 end do
559
560 do j = 1, lx
561 do i = 1, lx
562 do k = 1, lx
563 do e = 1, nelv
564 wt = 0d0
565 !NEC$ unroll_completely
566 do kk = 1, lx
567 wt = wt + dz(k, kk) * u(i,j, kk,e)
568 end do
569 dudt(i,j,k,e) = wt
570 end do
571 end do
572 end do
573 end do
574
575 do i = 1, nelv * lx * lx * lx
576 du(i,1,1,1) = jacinv(i,1,1,1) * &
577 (vx(i,1,1,1) * ( &
578 drdx(i,1,1,1)*dudr(i,1,1,1) &
579 + dsdx(i,1,1,1)*duds(i,1,1,1) &
580 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
581 + vy(i,1,1,1)*( &
582 drdy(i,1,1,1)*dudr(i,1,1,1) &
583 + dsdy(i,1,1,1)*duds(i,1,1,1) &
584 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
585 + vz(i,1,1,1)*( &
586 drdz(i,1,1,1)*dudr(i,1,1,1) &
587 + dsdz(i,1,1,1)*duds(i,1,1,1) &
588 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
589 end do
590
591 end subroutine sx_conv1_lx10
592
593 subroutine sx_conv1_lx9(du, u, vx, vy, vz, dx, dy, dz, &
594 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
595 jacinv, nelv)
596 integer, parameter :: lx = 9
597 integer, intent(in) :: nelv
598 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
599 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
600 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
601 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
602 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
603 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
604 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
605 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
606 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
607 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
608 real(kind=rp) :: wr, ws, wt
609 integer :: e, i, j, k, jj, kk
610
611 do i = 1, lx
612 do jj = 1, lx * lx * nelv
613 wr = 0d0
614 do kk = 1, lx
615 wr = wr + dx(i, kk) * u(kk, jj,1,1)
616 end do
617 dudr(i, jj,1,1) = wr
618 end do
619 end do
620
621 do k = 1, lx
622 do i = 1, lx
623 do j = 1, lx
624 do e = 1, nelv
625 ws = 0d0
626 !NEC$ unroll_completely
627 do kk = 1, lx
628 ws = ws + dy(j, kk) * u(i, kk,k,e)
629 end do
630 duds(i,j,k,e) = ws
631 end do
632 end do
633 end do
634 end do
635
636 do j = 1, lx
637 do i = 1, lx
638 do k = 1, lx
639 do e = 1, nelv
640 wt = 0d0
641 !NEC$ unroll_completely
642 do kk = 1, lx
643 wt = wt + dz(k, kk) * u(i,j, kk,e)
644 end do
645 dudt(i,j,k,e) = wt
646 end do
647 end do
648 end do
649 end do
650
651 do i = 1, nelv * lx * lx * lx
652 du(i,1,1,1) = jacinv(i,1,1,1) * &
653 (vx(i,1,1,1) * ( &
654 drdx(i,1,1,1)*dudr(i,1,1,1) &
655 + dsdx(i,1,1,1)*duds(i,1,1,1) &
656 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
657 + vy(i,1,1,1)*( &
658 drdy(i,1,1,1)*dudr(i,1,1,1) &
659 + dsdy(i,1,1,1)*duds(i,1,1,1) &
660 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
661 + vz(i,1,1,1)*( &
662 drdz(i,1,1,1)*dudr(i,1,1,1) &
663 + dsdz(i,1,1,1)*duds(i,1,1,1) &
664 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
665 end do
666
667 end subroutine sx_conv1_lx9
668
669 subroutine sx_conv1_lx8(du, u, vx, vy, vz, dx, dy, dz, &
670 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
671 jacinv, nelv)
672 integer, parameter :: lx = 8
673 integer, intent(in) :: nelv
674 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
675 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
676 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
677 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
678 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
679 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
680 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
681 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
682 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
683 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
684 real(kind=rp) :: wr, ws, wt
685 integer :: e, i, j, k, jj, kk
686
687 do i = 1, lx
688 do jj = 1, lx * lx * nelv
689 wr = 0d0
690 do kk = 1, lx
691 wr = wr + dx(i, kk) * u(kk, jj,1,1)
692 end do
693 dudr(i, jj,1,1) = wr
694 end do
695 end do
696
697 do k = 1, lx
698 do i = 1, lx
699 do j = 1, lx
700 do e = 1, nelv
701 ws = 0d0
702 !NEC$ unroll_completely
703 do kk = 1, lx
704 ws = ws + dy(j, kk) * u(i, kk,k,e)
705 end do
706 duds(i,j,k,e) = ws
707 end do
708 end do
709 end do
710 end do
711
712 do j = 1, lx
713 do i = 1, lx
714 do k = 1, lx
715 do e = 1, nelv
716 wt = 0d0
717 !NEC$ unroll_completely
718 do kk = 1, lx
719 wt = wt + dz(k, kk) * u(i,j, kk,e)
720 end do
721 dudt(i,j,k,e) = wt
722 end do
723 end do
724 end do
725 end do
726
727 do i = 1, nelv * lx * lx * lx
728 du(i,1,1,1) = jacinv(i,1,1,1) * &
729 (vx(i,1,1,1) * ( &
730 drdx(i,1,1,1)*dudr(i,1,1,1) &
731 + dsdx(i,1,1,1)*duds(i,1,1,1) &
732 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
733 + vy(i,1,1,1)*( &
734 drdy(i,1,1,1)*dudr(i,1,1,1) &
735 + dsdy(i,1,1,1)*duds(i,1,1,1) &
736 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
737 + vz(i,1,1,1)*( &
738 drdz(i,1,1,1)*dudr(i,1,1,1) &
739 + dsdz(i,1,1,1)*duds(i,1,1,1) &
740 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
741 end do
742
743 end subroutine sx_conv1_lx8
744
745 subroutine sx_conv1_lx7(du, u, vx, vy, vz, dx, dy, dz, &
746 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
747 jacinv, nelv)
748 integer, parameter :: lx = 7
749 integer, intent(in) :: nelv
750 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
751 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
752 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
753 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
754 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
755 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
756 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
757 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
758 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
759 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
760 real(kind=rp) :: wr, ws, wt
761 integer :: e, i, j, k, jj, kk
762
763 do i = 1, lx
764 do jj = 1, lx * lx * nelv
765 wr = 0d0
766 do kk = 1, lx
767 wr = wr + dx(i, kk) * u(kk, jj,1,1)
768 end do
769 dudr(i, jj,1,1) = wr
770 end do
771 end do
772
773 do k = 1, lx
774 do i = 1, lx
775 do j = 1, lx
776 do e = 1, nelv
777 ws = 0d0
778 !NEC$ unroll_completely
779 do kk = 1, lx
780 ws = ws + dy(j, kk) * u(i, kk,k,e)
781 end do
782 duds(i,j,k,e) = ws
783 end do
784 end do
785 end do
786 end do
787
788 do j = 1, lx
789 do i = 1, lx
790 do k = 1, lx
791 do e = 1, nelv
792 wt = 0d0
793 !NEC$ unroll_completely
794 do kk = 1, lx
795 wt = wt + dz(k, kk) * u(i,j, kk,e)
796 end do
797 dudt(i,j,k,e) = wt
798 end do
799 end do
800 end do
801 end do
802
803 do i = 1, nelv * lx * lx * lx
804 du(i,1,1,1) = jacinv(i,1,1,1) * &
805 (vx(i,1,1,1) * ( &
806 drdx(i,1,1,1)*dudr(i,1,1,1) &
807 + dsdx(i,1,1,1)*duds(i,1,1,1) &
808 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
809 + vy(i,1,1,1)*( &
810 drdy(i,1,1,1)*dudr(i,1,1,1) &
811 + dsdy(i,1,1,1)*duds(i,1,1,1) &
812 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
813 + vz(i,1,1,1)*( &
814 drdz(i,1,1,1)*dudr(i,1,1,1) &
815 + dsdz(i,1,1,1)*duds(i,1,1,1) &
816 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
817 end do
818
819 end subroutine sx_conv1_lx7
820
821 subroutine sx_conv1_lx6(du, u, vx, vy, vz, dx, dy, dz, &
822 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
823 jacinv, nelv)
824 integer, parameter :: lx = 6
825 integer, intent(in) :: nelv
826 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
827 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
828 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
829 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
830 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
831 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
832 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
833 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
834 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
835 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
836 real(kind=rp) :: wr, ws, wt
837 integer :: e, i, j, k, jj, kk
838
839 do i = 1, lx
840 do jj = 1, lx * lx * nelv
841 wr = 0d0
842 do kk = 1, lx
843 wr = wr + dx(i, kk) * u(kk, jj,1,1)
844 end do
845 dudr(i, jj,1,1) = wr
846 end do
847 end do
848
849 do k = 1, lx
850 do i = 1, lx
851 do j = 1, lx
852 do e = 1, nelv
853 ws = 0d0
854 !NEC$ unroll_completely
855 do kk = 1, lx
856 ws = ws + dy(j, kk) * u(i, kk,k,e)
857 end do
858 duds(i,j,k,e) = ws
859 end do
860 end do
861 end do
862 end do
863
864 do j = 1, lx
865 do i = 1, lx
866 do k = 1, lx
867 do e = 1, nelv
868 wt = 0d0
869 !NEC$ unroll_completely
870 do kk = 1, lx
871 wt = wt + dz(k, kk) * u(i,j, kk,e)
872 end do
873 dudt(i,j,k,e) = wt
874 end do
875 end do
876 end do
877 end do
878
879 do i = 1, nelv * lx * lx * lx
880 du(i,1,1,1) = jacinv(i,1,1,1) * &
881 (vx(i,1,1,1) * ( &
882 drdx(i,1,1,1)*dudr(i,1,1,1) &
883 + dsdx(i,1,1,1)*duds(i,1,1,1) &
884 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
885 + vy(i,1,1,1)*( &
886 drdy(i,1,1,1)*dudr(i,1,1,1) &
887 + dsdy(i,1,1,1)*duds(i,1,1,1) &
888 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
889 + vz(i,1,1,1)*( &
890 drdz(i,1,1,1)*dudr(i,1,1,1) &
891 + dsdz(i,1,1,1)*duds(i,1,1,1) &
892 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
893 end do
894
895 end subroutine sx_conv1_lx6
896
897 subroutine sx_conv1_lx5(du, u, vx, vy, vz, dx, dy, dz, &
898 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
899 jacinv, nelv)
900 integer, parameter :: lx = 5
901 integer, intent(in) :: nelv
902 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
903 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
904 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
905 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
906 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
907 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
908 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
909 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
910 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
911 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
912 real(kind=rp) :: wr, ws, wt
913 integer :: e, i, j, k, jj, kk
914
915 do i = 1, lx
916 do jj = 1, lx * lx * nelv
917 wr = 0d0
918 do kk = 1, lx
919 wr = wr + dx(i, kk) * u(kk, jj,1,1)
920 end do
921 dudr(i, jj,1,1) = wr
922 end do
923 end do
924
925 do k = 1, lx
926 do i = 1, lx
927 do j = 1, lx
928 do e = 1, nelv
929 ws = 0d0
930 !NEC$ unroll_completely
931 do kk = 1, lx
932 ws = ws + dy(j, kk) * u(i, kk,k,e)
933 end do
934 duds(i,j,k,e) = ws
935 end do
936 end do
937 end do
938 end do
939
940 do j = 1, lx
941 do i = 1, lx
942 do k = 1, lx
943 do e = 1, nelv
944 wt = 0d0
945 !NEC$ unroll_completely
946 do kk = 1, lx
947 wt = wt + dz(k, kk) * u(i,j, kk,e)
948 end do
949 dudt(i,j,k,e) = wt
950 end do
951 end do
952 end do
953 end do
954
955 do i = 1, nelv * lx * lx * lx
956 du(i,1,1,1) = jacinv(i,1,1,1) * &
957 (vx(i,1,1,1) * ( &
958 drdx(i,1,1,1)*dudr(i,1,1,1) &
959 + dsdx(i,1,1,1)*duds(i,1,1,1) &
960 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
961 + vy(i,1,1,1)*( &
962 drdy(i,1,1,1)*dudr(i,1,1,1) &
963 + dsdy(i,1,1,1)*duds(i,1,1,1) &
964 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
965 + vz(i,1,1,1)*( &
966 drdz(i,1,1,1)*dudr(i,1,1,1) &
967 + dsdz(i,1,1,1)*duds(i,1,1,1) &
968 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
969 end do
970
971 end subroutine sx_conv1_lx5
972
973 subroutine sx_conv1_lx4(du, u, vx, vy, vz, dx, dy, dz, &
974 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
975 jacinv, nelv)
976 integer, parameter :: lx = 4
977 integer, intent(in) :: nelv
978 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
979 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
980 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
981 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
982 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
983 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
984 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
985 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
986 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
987 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
988 real(kind=rp) :: wr, ws, wt
989 integer :: e, i, j, k, jj, kk
990
991 do i = 1, lx
992 do jj = 1, lx * lx * nelv
993 wr = 0d0
994 do kk = 1, lx
995 wr = wr + dx(i, kk) * u(kk, jj,1,1)
996 end do
997 dudr(i, jj,1,1) = wr
998 end do
999 end do
1000
1001 do k = 1, lx
1002 do i = 1, lx
1003 do j = 1, lx
1004 do e = 1, nelv
1005 ws = 0d0
1006 !NEC$ unroll_completely
1007 do kk = 1, lx
1008 ws = ws + dy(j, kk) * u(i, kk,k,e)
1009 end do
1010 duds(i,j,k,e) = ws
1011 end do
1012 end do
1013 end do
1014 end do
1015
1016 do j = 1, lx
1017 do i = 1, lx
1018 do k = 1, lx
1019 do e = 1, nelv
1020 wt = 0d0
1021 !NEC$ unroll_completely
1022 do kk = 1, lx
1023 wt = wt + dz(k, kk) * u(i,j, kk,e)
1024 end do
1025 dudt(i,j,k,e) = wt
1026 end do
1027 end do
1028 end do
1029 end do
1030
1031 do i = 1, nelv * lx * lx * lx
1032 du(i,1,1,1) = jacinv(i,1,1,1) * &
1033 (vx(i,1,1,1) * ( &
1034 drdx(i,1,1,1)*dudr(i,1,1,1) &
1035 + dsdx(i,1,1,1)*duds(i,1,1,1) &
1036 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
1037 + vy(i,1,1,1)*( &
1038 drdy(i,1,1,1)*dudr(i,1,1,1) &
1039 + dsdy(i,1,1,1)*duds(i,1,1,1) &
1040 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
1041 + vz(i,1,1,1)*( &
1042 drdz(i,1,1,1)*dudr(i,1,1,1) &
1043 + dsdz(i,1,1,1)*duds(i,1,1,1) &
1044 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
1045 end do
1046
1047 end subroutine sx_conv1_lx4
1048
1049 subroutine sx_conv1_lx3(du, u, vx, vy, vz, dx, dy, dz, &
1050 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
1051 jacinv, nelv)
1052 integer, parameter :: lx = 3
1053 integer, intent(in) :: nelv
1054 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
1055 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
1056 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
1057 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
1058 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
1059 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
1060 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
1061 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
1062 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
1063 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
1064 real(kind=rp) :: wr, ws, wt
1065 integer :: e, i, j, k, jj, kk
1066
1067 do i = 1, lx
1068 do jj = 1, lx * lx * nelv
1069 wr = 0d0
1070 do kk = 1, lx
1071 wr = wr + dx(i, kk) * u(kk, jj,1,1)
1072 end do
1073 dudr(i, jj,1,1) = wr
1074 end do
1075 end do
1076
1077 do k = 1, lx
1078 do i = 1, lx
1079 do j = 1, lx
1080 do e = 1, nelv
1081 ws = 0d0
1082 !NEC$ unroll_completely
1083 do kk = 1, lx
1084 ws = ws + dy(j, kk) * u(i, kk,k,e)
1085 end do
1086 duds(i,j,k,e) = ws
1087 end do
1088 end do
1089 end do
1090 end do
1091
1092 do j = 1, lx
1093 do i = 1, lx
1094 do k = 1, lx
1095 do e = 1, nelv
1096 wt = 0d0
1097 !NEC$ unroll_completely
1098 do kk = 1, lx
1099 wt = wt + dz(k, kk) * u(i,j, kk,e)
1100 end do
1101 dudt(i,j,k,e) = wt
1102 end do
1103 end do
1104 end do
1105 end do
1106
1107 do i = 1, nelv * lx * lx * lx
1108 du(i,1,1,1) = jacinv(i,1,1,1) * &
1109 (vx(i,1,1,1) * ( &
1110 drdx(i,1,1,1)*dudr(i,1,1,1) &
1111 + dsdx(i,1,1,1)*duds(i,1,1,1) &
1112 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
1113 + vy(i,1,1,1)*( &
1114 drdy(i,1,1,1)*dudr(i,1,1,1) &
1115 + dsdy(i,1,1,1)*duds(i,1,1,1) &
1116 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
1117 + vz(i,1,1,1)*( &
1118 drdz(i,1,1,1)*dudr(i,1,1,1) &
1119 + dsdz(i,1,1,1)*duds(i,1,1,1) &
1120 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
1121 end do
1122
1123 end subroutine sx_conv1_lx3
1124
1125 subroutine sx_conv1_lx2(du, u, vx, vy, vz, dx, dy, dz, &
1126 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, &
1127 jacinv, nelv)
1128 integer, parameter :: lx = 2
1129 integer, intent(in) :: nelv
1130 real(kind=rp), dimension(lx, lx, lx, nelv), intent(inout) :: du
1131 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: u, vx, vy, vz
1132 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdx, dsdx, dtdx
1133 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdy, dsdy, dtdy
1134 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: drdz, dsdz, dtdz
1135 real(kind=rp), dimension(lx, lx, lx, nelv), intent(in) :: jacinv
1136 real(kind=rp), dimension(lx, lx), intent(in) :: dx, dy, dz
1137 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudr
1138 real(kind=rp), dimension(lx, lx, lx, nelv) :: duds
1139 real(kind=rp), dimension(lx, lx, lx, nelv) :: dudt
1140 real(kind=rp) :: wr, ws, wt
1141 integer :: e, i, j, k, jj, kk
1142
1143 do i = 1, lx
1144 do jj = 1, lx * lx * nelv
1145 wr = 0d0
1146 do kk = 1, lx
1147 wr = wr + dx(i, kk) * u(kk, jj,1,1)
1148 end do
1149 dudr(i, jj,1,1) = wr
1150 end do
1151 end do
1152
1153 do k = 1, lx
1154 do i = 1, lx
1155 do j = 1, lx
1156 do e = 1, nelv
1157 ws = 0d0
1158 !NEC$ unroll_completely
1159 do kk = 1, lx
1160 ws = ws + dy(j, kk) * u(i, kk,k,e)
1161 end do
1162 duds(i,j,k,e) = ws
1163 end do
1164 end do
1165 end do
1166 end do
1167
1168 do j = 1, lx
1169 do i = 1, lx
1170 do k = 1, lx
1171 do e = 1, nelv
1172 wt = 0d0
1173 !NEC$ unroll_completely
1174 do kk = 1, lx
1175 wt = wt + dz(k, kk) * u(i,j, kk,e)
1176 end do
1177 dudt(i,j,k,e) = wt
1178 end do
1179 end do
1180 end do
1181 end do
1182
1183 do i = 1, nelv * lx * lx * lx
1184 du(i,1,1,1) = jacinv(i,1,1,1) * &
1185 (vx(i,1,1,1) * ( &
1186 drdx(i,1,1,1)*dudr(i,1,1,1) &
1187 + dsdx(i,1,1,1)*duds(i,1,1,1) &
1188 + dtdx(i,1,1,1)*dudt(i,1,1,1)) &
1189 + vy(i,1,1,1)*( &
1190 drdy(i,1,1,1)*dudr(i,1,1,1) &
1191 + dsdy(i,1,1,1)*duds(i,1,1,1) &
1192 + dtdy(i,1,1,1)*dudt(i,1,1,1)) &
1193 + vz(i,1,1,1)*( &
1194 drdz(i,1,1,1)*dudr(i,1,1,1) &
1195 + dsdz(i,1,1,1)*duds(i,1,1,1) &
1196 + dtdz(i,1,1,1)*dudt(i,1,1,1)))
1197 end do
1198
1199 end subroutine sx_conv1_lx2
1200
1201end submodule sx_conv1
opr_sx
Operators SX-Aurora backend.
Definition opr_sx.f90:2