ncnn

1
// Tencent is pleased to support the open source community by making ncnn available.
2
//
3
// Copyright (C) 2019 THL A29 Limited, a Tencent company. All rights reserved.
4
//
5
// Licensed under the BSD 3-Clause License (the "License"); you may not use this file except
6
// in compliance with the License. You may obtain a copy of the License at
7
//
8
// https://opensource.org/licenses/BSD-3-Clause
9
//
10
// Unless required by applicable law or agreed to in writing, software distributed
11
// under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR
12
// CONDITIONS OF ANY KIND, either express or implied. See the License for the
13
// specific language governing permissions and limitations under the License.
14

15
#version 450
16

17
#if NCNN_fp16_storage
18
#extension GL_EXT_shader_16bit_storage: require
19
#endif
20
#if NCNN_fp16_arithmetic
21
#extension GL_EXT_shader_explicit_arithmetic_types_float16: require
22
#endif
23

24
#define shape_constant_id_offset 0
25
layout (constant_id = shape_constant_id_offset + 0) const int dims = 0;
26
layout (constant_id = shape_constant_id_offset + 1) const int w = 0;
27
layout (constant_id = shape_constant_id_offset + 2) const int h = 0;
28
layout (constant_id = shape_constant_id_offset + 3) const int c = 0;
29
layout (constant_id = shape_constant_id_offset + 4) const int cstep = 0;
30

31
#if NCNN_image_shader
32
layout (binding = 0) uniform unfp sampler3D bottom_blob_3d;
33
layout (binding = 1, imfmtc1) writeonly uniform unfp image3D top_blob_3d;
34
#else
35
layout (binding = 0) buffer bottom_top_blob { sfp bottom_top_blob_data[]; };
36
#endif
37

38
layout (push_constant) uniform parameter
39
{
40
    int dims;
41
    int w;
42
    int h;
43
    int c;
44
    int cstep;
45
} p;
46

47
float erf(float x)
48
{
49
    float a1 =  0.254829592f;
50
    float a2 = -0.284496736f;
51
    float a3 =  1.421413741f;
52
    float a4 = -1.453152027f;
53
    float a5 =  1.061405429f;
54
    float p  =  0.3275911f;
55
    float s = sign(x);
56
    float x_abs = abs(x);
57
    float t = 1.0f/(1.0f + p*x_abs);
58
    float y = 1.0f - (((((a5*t + a4)*t) + a3)*t + a2)*t + a1)*t*exp(-x_abs*x_abs);
59
    return s * y;
60
}
61

62
void main()
63
{
64
    int gx = int(gl_GlobalInvocationID.x);
65
    int gy = int(gl_GlobalInvocationID.y);
66
    int gz = int(gl_GlobalInvocationID.z);
67

68
    if (gx >= psc(w) || gy >= psc(h) || gz >= psc(c))
69
        return;
70

71
#if NCNN_image_shader
72
    afp v = image3d_ld1(bottom_blob_3d, ivec3(gx, gy, gz));
73
#else
74
    const int gi = gz * psc(cstep) + gy * psc(w) + gx;
75

76
    afp v = buffer_ld1(bottom_top_blob_data, gi);
77
#endif
78

79
    v = erf(v);
80

81
#if NCNN_image_shader
82
    image3d_st1(top_blob_3d, ivec3(gx, gy, gz), v);
83
#else
84
    buffer_st1(bottom_top_blob_data, gi, v);
85
#endif
86
}
87