Neko
1.99.1
A portable framework for high-order spectral element flow simulations
Loading...
Searching...
No Matches
inhom_dirichlet.c
Go to the documentation of this file.
1
/*
2
Copyright (c) 2021-2025, The Neko Authors
3
All rights reserved.
4
5
Redistribution and use in source and binary forms, with or without
6
modification, are permitted provided that the following conditions
7
are met:
8
9
* Redistributions of source code must retain the above copyright
10
notice, this list of conditions and the following disclaimer.
11
12
* Redistributions in binary form must reproduce the above
13
copyright notice, this list of conditions and the following
14
disclaimer in the documentation and/or other materials provided
15
with the distribution.
16
17
* Neither the name of the authors nor the names of its
18
contributors may be used to endorse or promote products derived
19
from this software without specific prior written permission.
20
21
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
22
"AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
23
LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
24
FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
25
COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
26
INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
27
BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
28
LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
29
CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
30
LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
31
ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
32
POSSIBILITY OF SUCH DAMAGE.
33
*/
34
35
#ifdef __APPLE__
36
#include <OpenCL/cl.h>
37
#else
38
#include <CL/cl.h>
39
#endif
40
41
#include <stdio.h>
42
#include <
device/device_config.h
>
43
#include <
device/opencl/jit.h
>
44
#include <
device/opencl/prgm_lib.h
>
45
#include <
device/opencl/check.h
>
46
47
#include "inhom_dirichlet_kernel.cl.h"
48
52
void
opencl_inhom_dirichlet_apply_vector
(
void
*msk,
53
void
*
x
,
void
*y,
void
*z,
54
void
*
bla_x
,
void
*
bla_y
,
void
*
bla_z
,
55
int
*m,
cl_command_queue
cmd_queue
) {
56
cl_int
err
;
57
58
if
(
inhom_dirichlet_program
==
NULL
)
59
opencl_kernel_jit
(
inhom_dirichlet_kernel
,
60
(
cl_program
*) &
inhom_dirichlet_program
);
61
62
cl_kernel
kernel
=
clCreateKernel
(
inhom_dirichlet_program
,
63
"inhom_dirichlet_apply_vector_kernel"
, &
err
);
64
CL_CHECK
(
err
);
65
66
CL_CHECK
(
clSetKernelArg
(
kernel
, 0,
sizeof
(
cl_mem
), (
void
*) &msk));
67
CL_CHECK
(
clSetKernelArg
(
kernel
, 1,
sizeof
(
cl_mem
), (
void
*) &
x
));
68
CL_CHECK
(
clSetKernelArg
(
kernel
, 2,
sizeof
(
cl_mem
), (
void
*) &y));
69
CL_CHECK
(
clSetKernelArg
(
kernel
, 3,
sizeof
(
cl_mem
), (
void
*) &z));
70
CL_CHECK
(
clSetKernelArg
(
kernel
, 4,
sizeof
(
cl_mem
), (
void
*) &
bla_x
));
71
CL_CHECK
(
clSetKernelArg
(
kernel
, 5,
sizeof
(
cl_mem
), (
void
*) &
bla_y
));
72
CL_CHECK
(
clSetKernelArg
(
kernel
, 6,
sizeof
(
cl_mem
), (
void
*) &
bla_z
));
73
CL_CHECK
(
clSetKernelArg
(
kernel
, 7,
sizeof
(
int
), m));
74
75
const
int
nb
= ((*m) + 256 - 1) / 256;
76
const
size_t
global_item_size
= 256 *
nb
;
77
const
size_t
local_item_size
= 256;
78
79
CL_CHECK
(
clEnqueueNDRangeKernel
(
cmd_queue
,
kernel
, 1,
NULL
,
80
&
global_item_size
, &
local_item_size
,
81
0,
NULL
,
NULL
));
82
}
83
84
88
void
opencl_inhom_dirichlet_apply_scalar
(
void
*msk,
89
void
*
x
,
90
void
*
bla_x
,
91
int
*m,
92
cl_command_queue
cmd_queue
) {
93
cl_int
err
;
94
95
if
(
inhom_dirichlet_program
==
NULL
)
96
opencl_kernel_jit
(
inhom_dirichlet_kernel
,
97
(
cl_program
*) &
inhom_dirichlet_program
);
98
99
cl_kernel
kernel
=
clCreateKernel
(
inhom_dirichlet_program
,
100
"inhom_dirichlet_apply_scalar_kernel"
, &
err
);
101
CL_CHECK
(
err
);
102
103
CL_CHECK
(
clSetKernelArg
(
kernel
, 0,
sizeof
(
cl_mem
), (
void
*) &msk));
104
CL_CHECK
(
clSetKernelArg
(
kernel
, 1,
sizeof
(
cl_mem
), (
void
*) &
x
));
105
CL_CHECK
(
clSetKernelArg
(
kernel
, 2,
sizeof
(
cl_mem
), (
void
*) &
bla_x
));
106
CL_CHECK
(
clSetKernelArg
(
kernel
, 3,
sizeof
(
int
), m));
107
108
const
int
nb
= ((*m) + 256 - 1) / 256;
109
const
size_t
global_item_size
= 256 *
nb
;
110
const
size_t
local_item_size
= 256;
111
112
CL_CHECK
(
clEnqueueNDRangeKernel
(
cmd_queue
,
kernel
, 1,
NULL
,
113
&
global_item_size
, &
local_item_size
,
114
0,
NULL
,
NULL
));
115
}
x
__global__ void const T *__restrict__ x
Definition
cdtp_kernel.h:106
dirichlet_apply_scalar_kernel
__global__ void dirichlet_apply_scalar_kernel(const int *__restrict__ msk, T *__restrict__ x, const T g, const int m)
Definition
dirichlet_kernel.h:42
device_config.h
opencl_inhom_dirichlet_apply_scalar
void opencl_inhom_dirichlet_apply_scalar(void *msk, void *x, void *bla_x, int *m, cl_command_queue cmd_queue)
Definition
inhom_dirichlet.c:88
opencl_inhom_dirichlet_apply_vector
void opencl_inhom_dirichlet_apply_vector(void *msk, void *x, void *y, void *z, void *bla_x, void *bla_y, void *bla_z, int *m, cl_command_queue cmd_queue)
Definition
inhom_dirichlet.c:52
opencl_kernel_jit
void opencl_kernel_jit(const char *kernel, cl_program *program)
Definition
jit.c:50
jit.h
check.h
CL_CHECK
#define CL_CHECK(err)
Definition
check.h:12
prgm_lib.h
inhom_dirichlet_program
void * inhom_dirichlet_program
src
bc
bcknd
device
opencl
inhom_dirichlet.c
Generated on Mon Sep 1 2025 03:49:14 for Neko by
1.9.8