Warning: This file is not a C or C++ file. It does not have highlighting.

1/*===---- velintrin_approx.h - VEL intrinsics helper for VE ----------------===
2 *
3 * Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
4 * See https://llvm.org/LICENSE.txt for license information.
5 * SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
6 *
7 *===-----------------------------------------------------------------------===
8 */
9#ifndef __VEL_INTRIN_APPROX_H__
10#define __VEL_INTRIN_APPROX_H__
11
12static inline __vr _vel_approx_vfdivs_vvvl(__vr v0, __vr v1, int l) {
13 float s0;
14 __vr v2, v3, v4, v5;
15 v5 = _vel_vrcps_vvl(v1, l);
16 s0 = 1.0;
17 v4 = _vel_vfnmsbs_vsvvl(s0, v1, v5, l);
18 v3 = _vel_vfmads_vvvvl(v5, v5, v4, l);
19 v2 = _vel_vfmuls_vvvl(v0, v3, l);
20 v4 = _vel_vfnmsbs_vvvvl(v0, v2, v1, l);
21 v2 = _vel_vfmads_vvvvl(v2, v5, v4, l);
22 v0 = _vel_vfnmsbs_vvvvl(v0, v2, v1, l);
23 v0 = _vel_vfmads_vvvvl(v2, v3, v0, l);
24 return v0;
25}
26
27static inline __vr _vel_approx_pvfdiv_vvvl(__vr v0, __vr v1, int l) {
28 float s0;
29 __vr v2, v3, v4, v5;
30 v5 = _vel_pvrcp_vvl(v1, l);
31 s0 = 1.0;
32 v4 = _vel_pvfnmsb_vsvvl(s0, v1, v5, l);
33 v3 = _vel_pvfmad_vvvvl(v5, v5, v4, l);
34 v2 = _vel_pvfmul_vvvl(v0, v3, l);
35 v4 = _vel_pvfnmsb_vvvvl(v0, v2, v1, l);
36 v2 = _vel_pvfmad_vvvvl(v2, v5, v4, l);
37 v0 = _vel_pvfnmsb_vvvvl(v0, v2, v1, l);
38 v0 = _vel_pvfmad_vvvvl(v2, v3, v0, l);
39 return v0;
40}
41
42static inline __vr _vel_approx_vfdivs_vsvl(float s0, __vr v0, int l) {
43 float s1;
44 __vr v1, v2, v3, v4;
45 v4 = _vel_vrcps_vvl(v0, l);
46 s1 = 1.0;
47 v2 = _vel_vfnmsbs_vsvvl(s1, v0, v4, l);
48 v2 = _vel_vfmads_vvvvl(v4, v4, v2, l);
49 v1 = _vel_vfmuls_vsvl(s0, v2, l);
50 v3 = _vel_vfnmsbs_vsvvl(s0, v1, v0, l);
51 v1 = _vel_vfmads_vvvvl(v1, v4, v3, l);
52 v3 = _vel_vfnmsbs_vsvvl(s0, v1, v0, l);
53 v0 = _vel_vfmads_vvvvl(v1, v2, v3, l);
54 return v0;
55}
56
57static inline __vr _vel_approx_vfdivs_vvsl(__vr v0, float s0, int l) {
58 float s1;
59 __vr v1, v2;
60 s1 = 1.0f / s0;
61 v1 = _vel_vfmuls_vsvl(s1, v0, l);
62 v2 = _vel_vfnmsbs_vvsvl(v0, s0, v1, l);
63 v0 = _vel_vfmads_vvsvl(v1, s1, v2, l);
64 return v0;
65}
66
67static inline __vr _vel_approx_vfdivd_vsvl(double s0, __vr v0, int l) {
68 __vr v1, v2, v3;
69 v2 = _vel_vrcpd_vvl(v0, l);
70 double s1 = 1.0;
71 v3 = _vel_vfnmsbd_vsvvl(s1, v0, v2, l);
72 v2 = _vel_vfmadd_vvvvl(v2, v2, v3, l);
73 v1 = _vel_vfnmsbd_vsvvl(s1, v0, v2, l);
74 v1 = _vel_vfmadd_vvvvl(v2, v2, v1, l);
75 v1 = _vel_vaddul_vsvl(1, v1, l);
76 v3 = _vel_vfnmsbd_vsvvl(s1, v0, v1, l);
77 v3 = _vel_vfmadd_vvvvl(v1, v1, v3, l);
78 v1 = _vel_vfmuld_vsvl(s0, v3, l);
79 v0 = _vel_vfnmsbd_vsvvl(s0, v1, v0, l);
80 v0 = _vel_vfmadd_vvvvl(v1, v3, v0, l);
81 return v0;
82}
83
84static inline __vr _vel_approx_vfsqrtd_vvl(__vr v0, int l) {
85 double s0, s1;
86 __vr v1, v2, v3;
87 v2 = _vel_vrsqrtdnex_vvl(v0, l);
88 v1 = _vel_vfmuld_vvvl(v0, v2, l);
89 s0 = 1.0;
90 s1 = 0.5;
91 v3 = _vel_vfnmsbd_vsvvl(s0, v1, v2, l);
92 v3 = _vel_vfmuld_vsvl(s1, v3, l);
93 v2 = _vel_vfmadd_vvvvl(v2, v2, v3, l);
94 v1 = _vel_vfmuld_vvvl(v0, v2, l);
95 v3 = _vel_vfnmsbd_vsvvl(s0, v1, v2, l);
96 v3 = _vel_vfmuld_vsvl(s1, v3, l);
97 v0 = _vel_vfmadd_vvvvl(v1, v1, v3, l);
98 return v0;
99}
100
101static inline __vr _vel_approx_vfsqrts_vvl(__vr v0, int l) {
102 float s0, s1;
103 __vr v1, v2, v3;
104 v0 = _vel_vcvtds_vvl(v0, l);
105 v2 = _vel_vrsqrtdnex_vvl(v0, l);
106 v1 = _vel_vfmuld_vvvl(v0, v2, l);
107 s0 = 1.0;
108 s1 = 0.5;
109 v3 = _vel_vfnmsbd_vsvvl(s0, v1, v2, l);
110 v3 = _vel_vfmuld_vsvl(s1, v3, l);
111 v2 = _vel_vfmadd_vvvvl(v2, v2, v3, l);
112 v1 = _vel_vfmuld_vvvl(v0, v2, l);
113 v3 = _vel_vfnmsbd_vsvvl(s0, v1, v2, l);
114 v3 = _vel_vfmuld_vsvl(s1, v3, l);
115 v0 = _vel_vfmadd_vvvvl(v1, v1, v3, l);
116 v0 = _vel_vcvtsd_vvl(v0, l);
117 return v0;
118}
119
120#endif
121

Warning: This file is not a C or C++ file. It does not have highlighting.

source code of clang/lib/Headers/velintrin_approx.h