1 #include <cuda.h>
2 #include <cuda_runtime.h>
3 #include <c10/cuda/CUDAException.h>
4
5 #include <ATen/ATen.h>
6
sigmoid_add_kernel(const float * __restrict__ x,const float * __restrict__ y,float * __restrict__ output,const int size)7 __global__ void sigmoid_add_kernel(
8 const float* __restrict__ x,
9 const float* __restrict__ y,
10 float* __restrict__ output,
11 const int size) {
12 const int index = blockIdx.x * blockDim.x + threadIdx.x;
13 if (index < size) {
14 const float sigmoid_x = 1.0f / (1.0f + __expf(-x[index]));
15 const float sigmoid_y = 1.0f / (1.0f + __expf(-y[index]));
16 output[index] = sigmoid_x + sigmoid_y;
17 }
18 }
19
sigmoid_add_cuda(const float * x,const float * y,float * output,int size)20 void sigmoid_add_cuda(const float* x, const float* y, float* output, int size) {
21 const int threads = 1024;
22 const int blocks = (size + threads - 1) / threads;
23 sigmoid_add_kernel<<<blocks, threads>>>(x, y, output, size);
24 C10_CUDA_KERNEL_LAUNCH_CHECK();
25 }
26