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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
| /*
|
| Copyright (c) 2009, 2010, 2011 STMicroelectronics
| Written by Christophe Lyon
|
| Permission is hereby granted, free of charge, to any person obtaining a copy
| of this software and associated documentation files (the "Software"), to deal
| in the Software without restriction, including without limitation the rights
| to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
| copies of the Software, and to permit persons to whom the Software is
| furnished to do so, subject to the following conditions:
|
| The above copyright notice and this permission notice shall be included in
| all copies or substantial portions of the Software.
|
| THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
| IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
| FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
| AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
| LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
| OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
| THE SOFTWARE.
|
| */
|
| #if defined(__arm__) || defined(__aarch64__)
| #include <arm_neon.h>
| #else
| #include "stm-arm-neon.h"
| #endif
|
| #include "stm-arm-neon-ref.h"
| #include <math.h>
|
| #define TEST_MSG "VRECPS/VRECPSQ"
| void exec_vrecps(void)
| {
| int i;
|
| /* Basic test: y=vrecps(x), then store the result. */
| #define TEST_VRECPS(Q, T1, T2, W, N) \
| VECT_VAR(vector_res, T1, W, N) = \
| vrecps##Q##_##T2##W(VECT_VAR(vector, T1, W, N), \
| VECT_VAR(vector2, T1, W, N)); \
| vst1##Q##_##T2##W(VECT_VAR(result, T1, W, N), \
| VECT_VAR(vector_res, T1, W, N))
|
| /* With ARM RVCT, we need to declare variables before any executable
| statement */
|
| /* No need for integer variants */
| DECL_VARIABLE(vector, float, 32, 2);
| DECL_VARIABLE(vector, float, 32, 4);
|
| DECL_VARIABLE(vector2, float, 32, 2);
| DECL_VARIABLE(vector2, float, 32, 4);
|
| DECL_VARIABLE(vector_res, float, 32, 2);
| DECL_VARIABLE(vector_res, float, 32, 4);
|
| clean_results ();
|
| /* Choose init value arbitrarily */
| VDUP(vector, , float, f, 32, 2, 12.9f);
| VDUP(vector, q, float, f, 32, 4, 9.2f);
|
| VDUP(vector2, , float, f, 32, 2, 8.9f);
| VDUP(vector2, q, float, f, 32, 4, 3.2f);
|
| /* Apply the operator */
| TEST_VRECPS(, float, f, 32, 2);
| TEST_VRECPS(q, float, f, 32, 4);
|
| fprintf (ref_file, "\n%s output:\n", TEST_MSG);
| fprintf (gcc_tests_file, "\n%s output:\n", TEST_MSG);
| DUMP_FP(TEST_MSG, float, 32, 2, PRIx32);
| DUMP_FP(TEST_MSG, float, 32, 4, PRIx32);
|
|
| /* Test FP variants with special input values (NaN) */
| VDUP(vector, , float, f, 32, 2, NAN);
| VDUP(vector2, q, float, f, 32, 4, NAN);
|
| /* Apply the operator */
| TEST_VRECPS(, float, f, 32, 2);
| TEST_VRECPS(q, float, f, 32, 4);
|
| fprintf (ref_file, "\n%s output:\n", TEST_MSG " FP special (NAN) and normal values");
| DUMP_FP(TEST_MSG, float, 32, 2, PRIx32);
| DUMP_FP(TEST_MSG, float, 32, 4, PRIx32);
|
|
| /* Test FP variants with special input values (infinity, 0) */
| VDUP(vector, , float, f, 32, 2, HUGE_VALF);
| VDUP(vector, q, float, f, 32, 4, 0.0f);
| VDUP(vector2, q, float, f, 32, 4, 3.2f); /* Restore a normal value */
|
| /* Apply the operator */
| TEST_VRECPS(, float, f, 32, 2);
| TEST_VRECPS(q, float, f, 32, 4);
|
| fprintf (ref_file, "\n%s output:\n", TEST_MSG " FP special (infinity, 0) and normal values");
| DUMP_FP(TEST_MSG, float, 32, 2, PRIx32);
| DUMP_FP(TEST_MSG, float, 32, 4, PRIx32);
|
|
| /* Test FP variants with only special input values (infinity, 0) */
| VDUP(vector, , float, f, 32, 2, HUGE_VALF);
| VDUP(vector, q, float, f, 32, 4, 0.0f);
| VDUP(vector2, , float, f, 32, 2, 0.0f);
| VDUP(vector2, q, float, f, 32, 4, HUGE_VALF);
|
| /* Apply the operator */
| TEST_VRECPS(, float, f, 32, 2);
| TEST_VRECPS(q, float, f, 32, 4);
|
| fprintf (ref_file, "\n%s output:\n", TEST_MSG " FP special (infinity, 0)");
| DUMP_FP(TEST_MSG, float, 32, 2, PRIx32);
| DUMP_FP(TEST_MSG, float, 32, 4, PRIx32);
| }
|
|
