Rev 2071 | Go to most recent revision | Only display areas with differences | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 2071 | Rev 2131 | ||
---|---|---|---|
1 | /* |
1 | /* |
2 | * Copyright (c) 2005 Jakub Vana |
2 | * Copyright (c) 2005 Jakub Vana |
3 | * Copyright (c) 2005 Jakub Jermar |
3 | * Copyright (c) 2005 Jakub Jermar |
4 | * All rights reserved. |
4 | * All rights reserved. |
5 | * |
5 | * |
6 | * Redistribution and use in source and binary forms, with or without |
6 | * Redistribution and use in source and binary forms, with or without |
7 | * modification, are permitted provided that the following conditions |
7 | * modification, are permitted provided that the following conditions |
8 | * are met: |
8 | * are met: |
9 | * |
9 | * |
10 | * - Redistributions of source code must retain the above copyright |
10 | * - Redistributions of source code must retain the above copyright |
11 | * notice, this list of conditions and the following disclaimer. |
11 | * notice, this list of conditions and the following disclaimer. |
12 | * - Redistributions in binary form must reproduce the above copyright |
12 | * - Redistributions in binary form must reproduce the above copyright |
13 | * notice, this list of conditions and the following disclaimer in the |
13 | * notice, this list of conditions and the following disclaimer in the |
14 | * documentation and/or other materials provided with the distribution. |
14 | * documentation and/or other materials provided with the distribution. |
15 | * - The name of the author may not be used to endorse or promote products |
15 | * - The name of the author may not be used to endorse or promote products |
16 | * derived from this software without specific prior written permission. |
16 | * derived from this software without specific prior written permission. |
17 | * |
17 | * |
18 | * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR |
18 | * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR |
19 | * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES |
19 | * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES |
20 | * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. |
20 | * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. |
21 | * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT, |
21 | * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT, |
22 | * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT |
22 | * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT |
23 | * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, |
23 | * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, |
24 | * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY |
24 | * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY |
25 | * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
25 | * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
26 | * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF |
26 | * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF |
27 | * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
27 | * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
28 | */ |
28 | */ |
29 | 29 | ||
30 | #if (defined(ia32) || defined(amd64) || defined(ia64) || defined(ia32xen)) |
30 | #if (defined(ia32) || defined(amd64) || defined(ia64) || defined(ia32xen)) |
31 | 31 | ||
32 | #include <print.h> |
32 | #include <print.h> |
33 | #include <debug.h> |
33 | #include <debug.h> |
34 | 34 | ||
35 | #include <test.h> |
35 | #include <test.h> |
36 | #include <atomic.h> |
36 | #include <atomic.h> |
37 | #include <proc/thread.h> |
37 | #include <proc/thread.h> |
38 | 38 | ||
39 | #include <arch.h> |
39 | #include <arch.h> |
40 | #include <arch/arch.h> |
40 | #include <arch/arch.h> |
41 | 41 | ||
42 | 42 | ||
43 | #define THREADS 150 |
43 | #define THREADS 150 |
44 | #define ATTEMPTS 100 |
44 | #define ATTEMPTS 100 |
45 | 45 | ||
46 | #define E_10e8 271828182 |
46 | #define E_10e8 271828182 |
47 | #define PI_10e8 314159265 |
47 | #define PI_10e8 314159265 |
48 | 48 | ||
49 | 49 | ||
50 | #ifdef KERN_ia32_ARCH_H_ |
50 | #ifdef KERN_ia32_ARCH_H_ |
51 | static inline double sqrt(double x) |
51 | static inline double sqrt(double x) |
52 | { |
52 | { |
53 | double v; |
53 | double v; |
54 | 54 | ||
55 | asm ( |
55 | asm ( |
56 | "fsqrt\n" |
56 | "fsqrt\n" |
57 | : "=t" (v) |
57 | : "=t" (v) |
58 | : "0" (x) |
58 | : "0" (x) |
59 | ); |
59 | ); |
60 | 60 | ||
61 | return v; |
61 | return v; |
62 | } |
62 | } |
63 | #endif |
63 | #endif |
64 | 64 | ||
65 | #ifdef KERN_amd64_ARCH_H_ |
65 | #ifdef KERN_amd64_ARCH_H_ |
66 | static inline double sqrt(double x) |
66 | static inline double sqrt(double x) |
67 | { |
67 | { |
68 | double v; |
68 | double v; |
69 | 69 | ||
70 | asm ( |
70 | asm ( |
71 | "fsqrt\n" |
71 | "fsqrt\n" |
72 | : "=t" (v) |
72 | : "=t" (v) |
73 | : "0" (x) |
73 | : "0" (x) |
74 | ); |
74 | ); |
75 | 75 | ||
76 | return v; |
76 | return v; |
77 | } |
77 | } |
78 | #endif |
78 | #endif |
79 | 79 | ||
80 | #ifdef KERN_ia64_ARCH_H_ |
80 | #ifdef KERN_ia64_ARCH_H_ |
81 | 81 | ||
82 | #undef PI_10e8 |
82 | #undef PI_10e8 |
83 | #define PI_10e8 3141592 |
83 | #define PI_10e8 3141592 |
84 | 84 | ||
85 | static inline long double sqrt(long double a) |
85 | static inline long double sqrt(long double a) |
86 | { |
86 | { |
87 | long double x = 1; |
87 | long double x = 1; |
88 | long double lx = 0; |
88 | long double lx = 0; |
89 | 89 | ||
90 | if (a < 0.00000000000000001) |
90 | if (a < 0.00000000000000001) |
91 | return 0; |
91 | return 0; |
92 | 92 | ||
93 | while(x != lx) { |
93 | while(x != lx) { |
94 | lx = x; |
94 | lx = x; |
95 | x = (x + (a / x)) / 2; |
95 | x = (x + (a / x)) / 2; |
96 | } |
96 | } |
97 | 97 | ||
98 | return x; |
98 | return x; |
99 | } |
99 | } |
100 | #endif |
100 | #endif |
101 | 101 | ||
102 | 102 | ||
103 | static atomic_t threads_ok; |
103 | static atomic_t threads_ok; |
104 | static atomic_t threads_fault; |
104 | static atomic_t threads_fault; |
105 | static waitq_t can_start; |
105 | static waitq_t can_start; |
106 | static bool sh_quiet; |
106 | static bool sh_quiet; |
107 | 107 | ||
108 | static void e(void *data) |
108 | static void e(void *data) |
109 | { |
109 | { |
110 | int i; |
110 | int i; |
111 | double e, d, le, f; |
111 | double e, d, le, f; |
112 | 112 | ||
113 | thread_detach(THREAD); |
113 | thread_detach(THREAD); |
114 | 114 | ||
115 | waitq_sleep(&can_start); |
115 | waitq_sleep(&can_start); |
116 | 116 | ||
117 | for (i = 0; i<ATTEMPTS; i++) { |
117 | for (i = 0; i<ATTEMPTS; i++) { |
118 | le = -1; |
118 | le = -1; |
119 | e = 0; |
119 | e = 0; |
120 | f = 1; |
120 | f = 1; |
121 | 121 | ||
122 | for (d = 1; e != le; d *= f, f += 1) { |
122 | for (d = 1; e != le; d *= f, f += 1) { |
123 | le = e; |
123 | le = e; |
124 | e = e + 1 / d; |
124 | e = e + 1 / d; |
125 | } |
125 | } |
126 | 126 | ||
127 | if ((int) (100000000 * e) != E_10e8) { |
127 | if ((int) (100000000 * e) != E_10e8) { |
128 | if (!sh_quiet) |
128 | if (!sh_quiet) |
129 | printf("tid%d: e*10e8=%zd should be %zd\n", THREAD->tid, (unative_t) (100000000 * e), (unative_t) E_10e8); |
129 | printf("tid%d: e*10e8=%zd should be %zd\n", THREAD->tid, (unative_t) (100000000 * e), (unative_t) E_10e8); |
130 | atomic_inc(&threads_fault); |
130 | atomic_inc(&threads_fault); |
131 | break; |
131 | break; |
132 | } |
132 | } |
133 | } |
133 | } |
134 | atomic_inc(&threads_ok); |
134 | atomic_inc(&threads_ok); |
135 | } |
135 | } |
136 | 136 | ||
137 | static void pi(void *data) |
137 | static void pi(void *data) |
138 | { |
138 | { |
139 | int i; |
139 | int i; |
140 | double lpi, pi; |
140 | double lpi, pi; |
141 | double n, ab, ad; |
141 | double n, ab, ad; |
142 | 142 | ||
143 | thread_detach(THREAD); |
143 | thread_detach(THREAD); |
144 | 144 | ||
145 | waitq_sleep(&can_start); |
145 | waitq_sleep(&can_start); |
146 | 146 | ||
147 | for (i = 0; i < ATTEMPTS; i++) { |
147 | for (i = 0; i < ATTEMPTS; i++) { |
148 | lpi = -1; |
148 | lpi = -1; |
149 | pi = 0; |
149 | pi = 0; |
150 | 150 | ||
151 | for (n = 2, ab = sqrt(2); lpi != pi; n *= 2, ab = ad) { |
151 | for (n = 2, ab = sqrt(2); lpi != pi; n *= 2, ab = ad) { |
152 | double sc, cd; |
152 | double sc, cd; |
153 | 153 | ||
154 | sc = sqrt(1 - (ab * ab / 4)); |
154 | sc = sqrt(1 - (ab * ab / 4)); |
155 | cd = 1 - sc; |
155 | cd = 1 - sc; |
156 | ad = sqrt(ab * ab / 4 + cd * cd); |
156 | ad = sqrt(ab * ab / 4 + cd * cd); |
157 | lpi = pi; |
157 | lpi = pi; |
158 | pi = 2 * n * ad; |
158 | pi = 2 * n * ad; |
159 | } |
159 | } |
160 | 160 | ||
161 | #ifdef KERN_ia64_ARCH_H_ |
161 | #ifdef KERN_ia64_ARCH_H_ |
162 | if ((int) (1000000 * pi) != PI_10e8) { |
162 | if ((int) (1000000 * pi) != PI_10e8) { |
163 | if (!sh_quiet) |
163 | if (!sh_quiet) |
164 | printf("tid%d: pi*10e8=%zd should be %zd\n", THREAD->tid, (unative_t) (1000000 * pi), (unative_t) (PI_10e8 / 100)); |
164 | printf("tid%d: pi*10e8=%zd should be %zd\n", THREAD->tid, (unative_t) (1000000 * pi), (unative_t) (PI_10e8 / 100)); |
165 | atomic_inc(&threads_fault); |
165 | atomic_inc(&threads_fault); |
166 | break; |
166 | break; |
167 | } |
167 | } |
168 | #else |
168 | #else |
169 | if ((int) (100000000 * pi) != PI_10e8) { |
169 | if ((int) (100000000 * pi) != PI_10e8) { |
170 | if (!sh_quiet) |
170 | if (!sh_quiet) |
171 | printf("tid%d: pi*10e8=%zd should be %zd\n", THREAD->tid, (unative_t) (100000000 * pi), (unative_t) PI_10e8); |
171 | printf("tid%d: pi*10e8=%zd should be %zd\n", THREAD->tid, (unative_t) (100000000 * pi), (unative_t) PI_10e8); |
172 | atomic_inc(&threads_fault); |
172 | atomic_inc(&threads_fault); |
173 | break; |
173 | break; |
174 | } |
174 | } |
175 | #endif |
175 | #endif |
176 | } |
176 | } |
177 | atomic_inc(&threads_ok); |
177 | atomic_inc(&threads_ok); |
178 | } |
178 | } |
179 | 179 | ||
180 | char * test_fpu1(bool quiet) |
180 | char * test_fpu1(bool quiet) |
181 | { |
181 | { |
182 | unsigned int i, total = 0; |
182 | unsigned int i, total = 0; |
183 | sh_quiet = quiet; |
183 | sh_quiet = quiet; |
184 | 184 | ||
185 | waitq_initialize(&can_start); |
185 | waitq_initialize(&can_start); |
186 | atomic_set(&threads_ok, 0); |
186 | atomic_set(&threads_ok, 0); |
187 | atomic_set(&threads_fault, 0); |
187 | atomic_set(&threads_fault, 0); |
188 | 188 | ||
189 | if (!quiet) |
189 | if (!quiet) |
190 | printf("Creating %d threads... ", 2 * THREADS); |
190 | printf("Creating %d threads... ", 2 * THREADS); |
191 | 191 | ||
192 | for (i = 0; i < THREADS; i++) { |
192 | for (i = 0; i < THREADS; i++) { |
193 | thread_t *t; |
193 | thread_t *t; |
194 | 194 | ||
195 | if (!(t = thread_create(e, NULL, TASK, 0, "e", false))) { |
195 | if (!(t = thread_create(e, NULL, TASK, 0, "e", false))) { |
196 | if (!quiet) |
196 | if (!quiet) |
197 | printf("could not create thread %d\n", 2 * i); |
197 | printf("could not create thread %d\n", 2 * i); |
198 | break; |
198 | break; |
199 | } |
199 | } |
200 | thread_ready(t); |
200 | thread_ready(t); |
201 | total++; |
201 | total++; |
202 | 202 | ||
203 | if (!(t = thread_create(pi, NULL, TASK, 0, "pi", false))) { |
203 | if (!(t = thread_create(pi, NULL, TASK, 0, "pi", false))) { |
204 | if (!quiet) |
204 | if (!quiet) |
205 | printf("could not create thread %d\n", 2 * i + 1); |
205 | printf("could not create thread %d\n", 2 * i + 1); |
206 | break; |
206 | break; |
207 | } |
207 | } |
208 | thread_ready(t); |
208 | thread_ready(t); |
209 | total++; |
209 | total++; |
210 | } |
210 | } |
211 | 211 | ||
212 | if (!quiet) |
212 | if (!quiet) |
213 | printf("ok\n"); |
213 | printf("ok\n"); |
214 | 214 | ||
215 | thread_sleep(1); |
215 | thread_sleep(1); |
216 | waitq_wakeup(&can_start, WAKEUP_ALL); |
216 | waitq_wakeup(&can_start, WAKEUP_ALL); |
217 | 217 | ||
218 | while (atomic_get(&threads_ok) != total) { |
218 | while (atomic_get(&threads_ok) != total) { |
219 | if (!quiet) |
219 | if (!quiet) |
220 | printf("Threads left: %d\n", total - atomic_get(&threads_ok)); |
220 | printf("Threads left: %d\n", total - atomic_get(&threads_ok)); |
221 | thread_sleep(1); |
221 | thread_sleep(1); |
222 | } |
222 | } |
223 | 223 | ||
224 | if (atomic_get(&threads_fault) == 0) |
224 | if (atomic_get(&threads_fault) == 0) |
225 | return NULL; |
225 | return NULL; |
226 | 226 | ||
227 | return "Test failed"; |
227 | return "Test failed"; |
228 | } |
228 | } |
229 | 229 | ||
230 | #endif |
230 | #endif |
231 | 231 |