Neko 1.99.1
A portable framework for high-order spectral element flow simulations
Loading...
Searching...
No Matches
prgm_lib.F90
Go to the documentation of this file.
1
2module opencl_prgm_lib
3 use opencl_intf
4 use utils, only : neko_error
5 use, intrinsic :: iso_c_binding, only : c_ptr, c_null_ptr
6 implicit none
7 private
8
9#ifdef HAVE_OPENCL
10
12 type(c_ptr), public, bind(c) :: math_program = c_null_ptr
13
15 type(c_ptr), public, bind(c) :: mathops_program = c_null_ptr
16
18 type(c_ptr), public, bind(c) :: dirichlet_program = c_null_ptr
19
21 type(c_ptr), public, bind(c) :: inflow_program = c_null_ptr
22
24 type(c_ptr), public, bind(c) :: zero_dirichlet_program = c_null_ptr
25
27 type(c_ptr), public, bind(c) :: symmetry_program = c_null_ptr
28
30 type(c_ptr), public, bind(c) :: facet_normal_program = c_null_ptr
31
33 type(c_ptr), public, bind(c) :: neumann_program = c_null_ptr
34
36 type(c_ptr), public, bind(c) :: inhom_dirichlet_program = c_null_ptr
37
39 type(c_ptr), public, bind(c) :: dudxyz_program = c_null_ptr
40
42 type(c_ptr), public, bind(c) :: cdtp_program = c_null_ptr
43
45 type(c_ptr), public, bind(c) :: conv1_program = c_null_ptr
46
48 type(c_ptr), public, bind(c) :: convect_scalar_program = c_null_ptr
49
51 type(c_ptr), public, bind(c) :: set_convect_rst_program = c_null_ptr
52
54 type(c_ptr), public, bind(c) :: cfl_program = c_null_ptr
55
57 type(c_ptr), public, bind(c) :: opgrad_program = c_null_ptr
58
60 type(c_ptr), public, bind(c) :: gs_program = c_null_ptr
61
63 type(c_ptr), public, bind(c) :: ax_helm_program = c_null_ptr
64
66 type(c_ptr), public, bind(c) :: ax_helm_full_program = c_null_ptr
67
69 type(c_ptr), public, bind(c) :: jacobi_program = c_null_ptr
70
72 type(c_ptr), public, bind(c) :: rhs_maker_program = c_null_ptr
73
75 type(c_ptr), public, bind(c) :: pnpn_res_program = c_null_ptr
76
78 type(c_ptr), public, bind(c) :: pnpn_stress_res_program = c_null_ptr
79
81 type(c_ptr), public, bind(c) :: euler_res_program = c_null_ptr
82
84 type(c_ptr), public, bind(c) :: compressible_ops_compute_max_wave_speed_program = c_null_ptr
85 type(c_ptr), public, bind(c) :: compressible_ops_compute_entropy_program = c_null_ptr
86
88 type(c_ptr), public, bind(c) :: fdm_program = c_null_ptr
89
91 type(c_ptr), public, bind(c) :: tensor_program = c_null_ptr
92
94 type(c_ptr), public, bind(c) :: schwarz_program = c_null_ptr
95
97 type(c_ptr), public, bind(c) :: dong_program = c_null_ptr
98
100 type(c_ptr), public, bind(c) :: coef_program = c_null_ptr
101
103 type(c_ptr), public, bind(c) :: scalar_residual_program = c_null_ptr
104
106 type(c_ptr), public, bind(c) :: lambda2_program = c_null_ptr
107
109 type(c_ptr), public, bind(c) :: compute_max_wave_speed_program = c_null_ptr
110
112 type(c_ptr), public, bind(c) :: mapping_program = c_null_ptr
113
115 type(c_ptr), public, bind(c) :: find_rst_legendre_program = c_null_ptr
116
117 public :: opencl_prgm_lib_release
118
119contains
120
121 subroutine opencl_prgm_lib_release
122
123 if (c_associated(math_program)) then
124 if(clreleaseprogram(math_program) .ne. cl_success) then
125 call neko_error('Failed to release program')
126 end if
127 math_program = c_null_ptr
128 end if
129
130 if (c_associated(mathops_program)) then
131 if(clreleaseprogram(mathops_program) .ne. cl_success) then
132 call neko_error('Failed to release program')
133 end if
134 mathops_program = c_null_ptr
135 end if
136
137 if (c_associated(dirichlet_program)) then
138 if(clreleaseprogram(dirichlet_program) .ne. cl_success) then
139 call neko_error('Failed to release program')
140 end if
141 dirichlet_program = c_null_ptr
142 end if
143
144 if (c_associated(inflow_program)) then
145 if(clreleaseprogram(inflow_program) .ne. cl_success) then
146 call neko_error('Failed to release program')
147 end if
148 inflow_program = c_null_ptr
149 end if
150
151 if (c_associated(zero_dirichlet_program)) then
152 if (clreleaseprogram(zero_dirichlet_program) .ne. cl_success) then
153 call neko_error('Failed to release program')
154 end if
155 zero_dirichlet_program = c_null_ptr
156 end if
157
158 if (c_associated(symmetry_program)) then
159 if(clreleaseprogram(symmetry_program) .ne. cl_success) then
160 call neko_error('Failed to release program')
161 end if
162 symmetry_program = c_null_ptr
163 end if
164
165 if (c_associated(facet_normal_program)) then
166 if(clreleaseprogram(facet_normal_program) .ne. cl_success) then
167 call neko_error('Failed to release program')
168 end if
169 facet_normal_program = c_null_ptr
170 end if
171
172 if (c_associated(inhom_dirichlet_program)) then
173 if(clreleaseprogram(inhom_dirichlet_program) .ne. cl_success) then
174 call neko_error('Failed to release program')
175 end if
176 inhom_dirichlet_program = c_null_ptr
177 end if
178
179 if (c_associated(dudxyz_program)) then
180 if(clreleaseprogram(dudxyz_program) .ne. cl_success) then
181 call neko_error('Failed to release program')
182 end if
183 dudxyz_program = c_null_ptr
184 end if
185
186 if (c_associated(cdtp_program)) then
187 if(clreleaseprogram(cdtp_program) .ne. cl_success) then
188 call neko_error('Failed to release program')
189 end if
190 cdtp_program = c_null_ptr
191 end if
192
193 if (c_associated(conv1_program)) then
194 if(clreleaseprogram(conv1_program) .ne. cl_success) then
195 call neko_error('Failed to release program')
196 end if
197 conv1_program = c_null_ptr
198 end if
199
200 if (c_associated(cfl_program)) then
201 if(clreleaseprogram(cfl_program) .ne. cl_success) then
202 call neko_error('Failed to release program')
203 end if
204 cfl_program = c_null_ptr
205 end if
206
207 if (c_associated(opgrad_program)) then
208 if(clreleaseprogram(opgrad_program) .ne. cl_success) then
209 call neko_error('Failed to release program')
210 end if
211 opgrad_program = c_null_ptr
212 end if
213
214 if (c_associated(gs_program)) then
215 if(clreleaseprogram(gs_program) .ne. cl_success) then
216 call neko_error('Failed to release program')
217 end if
218 gs_program = c_null_ptr
219 end if
220
221 if (c_associated(ax_helm_program)) then
222 if(clreleaseprogram(ax_helm_program) .ne. cl_success) then
223 call neko_error('Failed to release program')
224 end if
225 ax_helm_program = c_null_ptr
226 end if
227
228 if (c_associated(ax_helm_full_program)) then
229 if(clreleaseprogram(ax_helm_full_program) .ne. cl_success) then
230 call neko_error('Failed to release program')
231 end if
232 ax_helm_full_program = c_null_ptr
233 end if
234
235 if (c_associated(jacobi_program)) then
236 if(clreleaseprogram(jacobi_program) .ne. cl_success) then
237 call neko_error('Failed to release program')
238 end if
239 jacobi_program = c_null_ptr
240 end if
241
242 if (c_associated(rhs_maker_program)) then
243 if(clreleaseprogram(rhs_maker_program) .ne. cl_success) then
244 call neko_error('Failed to release program')
245 end if
246 rhs_maker_program = c_null_ptr
247 end if
248
249 if (c_associated(pnpn_res_program)) then
250 if(clreleaseprogram(pnpn_res_program) .ne. cl_success) then
251 call neko_error('Failed to release program')
252 end if
253 pnpn_res_program = c_null_ptr
254 end if
255
256 if (c_associated(pnpn_stress_res_program)) then
257 if(clreleaseprogram(pnpn_stress_res_program) .ne. cl_success) then
258 call neko_error('Failed to release program')
259 end if
260 pnpn_stress_res_program = c_null_ptr
261 end if
262
263 if (c_associated(euler_res_program)) then
264 if(clreleaseprogram(euler_res_program) .ne. cl_success) then
265 call neko_error('Failed to release program')
266 end if
267 euler_res_program = c_null_ptr
268 end if
269
270 if (c_associated(compressible_ops_compute_max_wave_speed_program)) then
271 if(clreleaseprogram(compressible_ops_compute_max_wave_speed_program) .ne. cl_success) then
272 call neko_error('Failed to release program')
273 end if
274 compressible_ops_compute_max_wave_speed_program = c_null_ptr
275 end if
276
277 if (c_associated(compressible_ops_compute_entropy_program)) then
278 if(clreleaseprogram(compressible_ops_compute_entropy_program) .ne. cl_success) then
279 call neko_error('Failed to release program')
280 end if
281 compressible_ops_compute_entropy_program = c_null_ptr
282 end if
283
284 if (c_associated(fdm_program)) then
285 if(clreleaseprogram(fdm_program) .ne. cl_success) then
286 call neko_error('Failed to release program')
287 end if
288 fdm_program = c_null_ptr
289 end if
290
291 if (c_associated(tensor_program)) then
292 if(clreleaseprogram(tensor_program) .ne. cl_success) then
293 call neko_error('Failed to release program')
294 end if
295 tensor_program = c_null_ptr
296 end if
297
298 if (c_associated(schwarz_program)) then
299 if(clreleaseprogram(schwarz_program) .ne. cl_success) then
300 call neko_error('Failed to release program')
301 end if
302 schwarz_program = c_null_ptr
303 end if
304
305 if (c_associated(dong_program)) then
306 if(clreleaseprogram(dong_program) .ne. cl_success) then
307 call neko_error('Failed to release program')
308 end if
309 dong_program = c_null_ptr
310 end if
311
312 if (c_associated(coef_program)) then
313 if(clreleaseprogram(coef_program) .ne. cl_success) then
314 call neko_error('Failed to release program')
315 end if
316 coef_program = c_null_ptr
317 end if
318
319 if (c_associated(scalar_residual_program)) then
320 if(clreleaseprogram(scalar_residual_program) .ne. cl_success) then
321 call neko_error('Failed to release program')
322 end if
323 scalar_residual_program = c_null_ptr
324 end if
325
326 if (c_associated(lambda2_program)) then
327 if(clreleaseprogram(lambda2_program) .ne. cl_success) then
328 call neko_error('Failed to release program')
329 end if
330 lambda2_program = c_null_ptr
331 end if
332
333 if (c_associated(compute_max_wave_speed_program)) then
334 if(clreleaseprogram(compute_max_wave_speed_program) .ne. cl_success) then
335 call neko_error('Failed to release program')
336 end if
337 compute_max_wave_speed_program = c_null_ptr
338 end if
339
340 if (c_associated(mapping_program)) then
341 if(clreleaseprogram(mapping_program) .ne. cl_success) then
342 call neko_error('Failed to release program')
343 end if
344 mapping_program = c_null_ptr
345 end if
346
347 if (c_associated(find_rst_legendre_program)) then
348 if(clreleaseprogram(find_rst_legendre_program) .ne. cl_success) then
349 call neko_error('Failed to release program')
350 end if
351 find_rst_legendre_program = c_null_ptr
352 end if
353
354 end subroutine opencl_prgm_lib_release
355
356#endif
357
358end module opencl_prgm_lib
opencl_intf::clreleaseprogram
Definition opencl_intf.F90:362
utils::neko_error
Definition utils.f90:42
opencl_intf
Fortran OpenCL interface.
Definition opencl_intf.F90:34
opencl_intf::cl_success
@ cl_success
Definition opencl_intf.F90:50
opencl_prgm_lib
OpenCL JIT program library.
Definition prgm_lib.F90:2
opencl_prgm_lib::jacobi_program
type(c_ptr), bind(C), public jacobi_program
Device jacobi kernels.
Definition prgm_lib.F90:69
opencl_prgm_lib::opgrad_program
type(c_ptr), bind(C), public opgrad_program
Device Velocity gradient kernels.
Definition prgm_lib.F90:57
opencl_prgm_lib::math_program
type(c_ptr), bind(C), public math_program
Device math kernels.
Definition prgm_lib.F90:12
opencl_prgm_lib::convect_scalar_program
type(c_ptr), bind(C), public convect_scalar_program
Device convective kernels for oifs.
Definition prgm_lib.F90:48
opencl_prgm_lib::symmetry_program
type(c_ptr), bind(C), public symmetry_program
Device Symmetry kernels.
Definition prgm_lib.F90:27
opencl_prgm_lib::inhom_dirichlet_program
type(c_ptr), bind(C), public inhom_dirichlet_program
Device Blasius profile kernel.
Definition prgm_lib.F90:36
opencl_prgm_lib::set_convect_rst_program
type(c_ptr), bind(C), public set_convect_rst_program
Device convect_rst kernels.
Definition prgm_lib.F90:51
opencl_prgm_lib::mathops_program
type(c_ptr), bind(C), public mathops_program
Device mathops kernels.
Definition prgm_lib.F90:15
opencl_prgm_lib::dong_program
type(c_ptr), bind(C), public dong_program
Device dong kernels.
Definition prgm_lib.F90:97
opencl_prgm_lib::conv1_program
type(c_ptr), bind(C), public conv1_program
Device convective kernels.
Definition prgm_lib.F90:45
opencl_prgm_lib::inflow_program
type(c_ptr), bind(C), public inflow_program
Device Inflow kernels.
Definition prgm_lib.F90:21
opencl_prgm_lib::ax_helm_full_program
type(c_ptr), bind(C), public ax_helm_full_program
Device Ax helm full kernels.
Definition prgm_lib.F90:66
opencl_prgm_lib::dudxyz_program
type(c_ptr), bind(C), public dudxyz_program
Device Derivative kernels.
Definition prgm_lib.F90:39
opencl_prgm_lib::coef_program
type(c_ptr), bind(C), public coef_program
Device coef kernels.
Definition prgm_lib.F90:100
opencl_prgm_lib::lambda2_program
type(c_ptr), bind(C), public lambda2_program
Device lambda2 kernels.
Definition prgm_lib.F90:106
opencl_prgm_lib::tensor_program
type(c_ptr), bind(C), public tensor_program
Device tensor kernels.
Definition prgm_lib.F90:91
opencl_prgm_lib::ax_helm_program
type(c_ptr), bind(C), public ax_helm_program
Device Ax helm kernels.
Definition prgm_lib.F90:63
opencl_prgm_lib::cdtp_program
type(c_ptr), bind(C), public cdtp_program
Device kernels.
Definition prgm_lib.F90:42
opencl_prgm_lib::opencl_prgm_lib_release
subroutine, public opencl_prgm_lib_release
Definition prgm_lib.F90:122
opencl_prgm_lib::mapping_program
type(c_ptr), bind(C), public mapping_program
Device filter kernels.
Definition prgm_lib.F90:112
opencl_prgm_lib::schwarz_program
type(c_ptr), bind(C), public schwarz_program
Device schwarz kernels.
Definition prgm_lib.F90:94
opencl_prgm_lib::neumann_program
type(c_ptr), bind(C), public neumann_program
Device Neumann kernels.
Definition prgm_lib.F90:33
opencl_prgm_lib::gs_program
type(c_ptr), bind(C), public gs_program
Device Gather-Scatter kernels.
Definition prgm_lib.F90:60
opencl_prgm_lib::zero_dirichlet_program
type(c_ptr), bind(C), public zero_dirichlet_program
Device zero dirichlet kernels.
Definition prgm_lib.F90:24
opencl_prgm_lib::cfl_program
type(c_ptr), bind(C), public cfl_program
Device CFL kernels.
Definition prgm_lib.F90:54
opencl_prgm_lib::dirichlet_program
type(c_ptr), bind(C), public dirichlet_program
Device Dirichlet kernels.
Definition prgm_lib.F90:18
opencl_prgm_lib::compressible_ops_compute_entropy_program
type(c_ptr), bind(C), public compressible_ops_compute_entropy_program
Definition prgm_lib.F90:85
opencl_prgm_lib::compute_max_wave_speed_program
type(c_ptr), bind(C), public compute_max_wave_speed_program
Device compute_max_wave_speed kernels.
Definition prgm_lib.F90:109
opencl_prgm_lib::find_rst_legendre_program
type(c_ptr), bind(C), public find_rst_legendre_program
Device find rest kernels.
Definition prgm_lib.F90:115
opencl_prgm_lib::pnpn_res_program
type(c_ptr), bind(C), public pnpn_res_program
Device pnpn residual kernels.
Definition prgm_lib.F90:75
opencl_prgm_lib::rhs_maker_program
type(c_ptr), bind(C), public rhs_maker_program
Device rhs_maker kernels.
Definition prgm_lib.F90:72
opencl_prgm_lib::scalar_residual_program
type(c_ptr), bind(C), public scalar_residual_program
Device scalar residual kernels.
Definition prgm_lib.F90:103
opencl_prgm_lib::euler_res_program
type(c_ptr), bind(C), public euler_res_program
Device euler residual kernels.
Definition prgm_lib.F90:81
opencl_prgm_lib::pnpn_stress_res_program
type(c_ptr), bind(C), public pnpn_stress_res_program
Device pnpn residual kernels (stress formulation)
Definition prgm_lib.F90:78
opencl_prgm_lib::facet_normal_program
type(c_ptr), bind(C), public facet_normal_program
Device Facet normal kernels.
Definition prgm_lib.F90:30
opencl_prgm_lib::fdm_program
type(c_ptr), bind(C), public fdm_program
Device fdm kernels.
Definition prgm_lib.F90:88
opencl_prgm_lib::compressible_ops_compute_max_wave_speed_program
type(c_ptr), bind(C), public compressible_ops_compute_max_wave_speed_program
Device compressible ops kernels.
Definition prgm_lib.F90:84
utils
Utilities.
Definition utils.f90:35