sfsqrt_opt.S
1.77 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
//===----------------------Hexagon builtin routine ------------------------===//
//
// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
// See https://llvm.org/LICENSE.txt for license information.
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
//
//===----------------------------------------------------------------------===//
.macro FUNCTION_BEGIN name
.text
.p2align 5
.globl \name
.type \name, @function
\name:
.endm
.macro FUNCTION_END name
.size \name, . - \name
.endm
#define RIN r0
#define S r0
#define H r1
#define D r2
#define E r3
#define HALF r4
#define R r5
#define Q6_ALIAS(TAG) .global __qdsp_##TAG ; .set __qdsp_##TAG, __hexagon_##TAG
#define FAST_ALIAS(TAG) .global __hexagon_fast_##TAG ; .set __hexagon_fast_##TAG, __hexagon_##TAG
#define FAST2_ALIAS(TAG) .global __hexagon_fast2_##TAG ; .set __hexagon_fast2_##TAG, __hexagon_##TAG
FUNCTION_BEGIN __hexagon_sqrtf
{
E,p0 = sfinvsqrta(RIN)
R = sffixupr(RIN)
HALF = ##0x3f000000 // 0.5
r1:0 = combine(#0,#0) // clear S/H
}
{
S += sfmpy(E,R):lib // S0
H += sfmpy(E,HALF):lib // H0
D = HALF
E = R
}
{
D -= sfmpy(S,H):lib // d0
p1 = sfclass(R,#1) // is zero?
//E -= sfmpy(S,S):lib // e0
}
{
S += sfmpy(S,D):lib // S1
H += sfmpy(H,D):lib // H1
D = HALF
E = R
}
{
D -= sfmpy(S,H):lib // d0
E -= sfmpy(S,S):lib // e0
}
{
S += sfmpy(H,E):lib // S2
H += sfmpy(H,D):lib // H2
D = HALF
E = R
}
{
//D -= sfmpy(S,H):lib // d2
E -= sfmpy(S,S):lib // e2
if (p1) r0 = or(r0,R) // sqrt(-0.0) = -0.0
}
{
S += sfmpy(H,E,p0):scale // S3
jumpr r31
}
FUNCTION_END __hexagon_sqrtf
Q6_ALIAS(sqrtf)
FAST_ALIAS(sqrtf)
FAST2_ALIAS(sqrtf)