xref: /freebsd/contrib/netbsd-tests/lib/libm/t_log.c (revision 63d1fd5970ec814904aa0f4580b10a0d302d08b2)
1 /* $NetBSD: t_log.c,v 1.13 2015/02/09 19:39:48 martin Exp $ */
2 
3 /*-
4  * Copyright (c) 2011 The NetBSD Foundation, Inc.
5  * All rights reserved.
6  *
7  * This code is derived from software contributed to The NetBSD Foundation
8  * by Jukka Ruohonen.
9  *
10  * Redistribution and use in source and binary forms, with or without
11  * modification, are permitted provided that the following conditions
12  * are met:
13  * 1. Redistributions of source code must retain the above copyright
14  *    notice, this list of conditions and the following disclaimer.
15  * 2. Redistributions in binary form must reproduce the above copyright
16  *    notice, this list of conditions and the following disclaimer in the
17  *    documentation and/or other materials provided with the distribution.
18  *
19  * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
20  * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
21  * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
22  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
23  * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
24  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
25  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
26  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
27  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
28  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
29  * POSSIBILITY OF SUCH DAMAGE.
30  */
31 #include <sys/cdefs.h>
32 __RCSID("$NetBSD: t_log.c,v 1.13 2015/02/09 19:39:48 martin Exp $");
33 
34 #include <atf-c.h>
35 
36 #include <math.h>
37 #include <stdio.h>
38 #include <string.h>
39 
40 /*
41  * log10(3)
42  */
43 ATF_TC(log10_base);
44 ATF_TC_HEAD(log10_base, tc)
45 {
46 	atf_tc_set_md_var(tc, "descr", "Test log10(10) == 1");
47 }
48 
49 ATF_TC_BODY(log10_base, tc)
50 {
51 	ATF_CHECK(log10(10.0) == 1.0);
52 }
53 
54 ATF_TC(log10_nan);
55 ATF_TC_HEAD(log10_nan, tc)
56 {
57 	atf_tc_set_md_var(tc, "descr", "Test log10(NaN) == NaN");
58 }
59 
60 ATF_TC_BODY(log10_nan, tc)
61 {
62 	const double x = 0.0L / 0.0L;
63 
64 	ATF_CHECK(isnan(x) != 0);
65 	ATF_CHECK(isnan(log10(x)) != 0);
66 }
67 
68 ATF_TC(log10_inf_neg);
69 ATF_TC_HEAD(log10_inf_neg, tc)
70 {
71 	atf_tc_set_md_var(tc, "descr", "Test log10(-Inf) == NaN");
72 }
73 
74 ATF_TC_BODY(log10_inf_neg, tc)
75 {
76 	const double x = -1.0L / 0.0L;
77 	const double y = log10(x);
78 
79 	ATF_CHECK(isnan(y) != 0);
80 }
81 
82 ATF_TC(log10_inf_pos);
83 ATF_TC_HEAD(log10_inf_pos, tc)
84 {
85 	atf_tc_set_md_var(tc, "descr", "Test log10(+Inf) == +Inf");
86 }
87 
88 ATF_TC_BODY(log10_inf_pos, tc)
89 {
90 	const double x = 1.0L / 0.0L;
91 
92 	ATF_CHECK(log10(x) == x);
93 }
94 
95 ATF_TC(log10_one_pos);
96 ATF_TC_HEAD(log10_one_pos, tc)
97 {
98 	atf_tc_set_md_var(tc, "descr", "Test log10(1.0) == +0.0");
99 }
100 
101 ATF_TC_BODY(log10_one_pos, tc)
102 {
103 	const double x = log10(1.0);
104 	const double y = 0.0L;
105 
106 	ATF_CHECK(x == y);
107 	ATF_CHECK(signbit(x) == 0);
108 	ATF_CHECK(signbit(y) == 0);
109 }
110 
111 ATF_TC(log10_zero_neg);
112 ATF_TC_HEAD(log10_zero_neg, tc)
113 {
114 	atf_tc_set_md_var(tc, "descr", "Test log10(-0.0) == -HUGE_VAL");
115 }
116 
117 ATF_TC_BODY(log10_zero_neg, tc)
118 {
119 	const double x = -0.0L;
120 
121 	ATF_CHECK(log10(x) == -HUGE_VAL);
122 }
123 
124 ATF_TC(log10_zero_pos);
125 ATF_TC_HEAD(log10_zero_pos, tc)
126 {
127 	atf_tc_set_md_var(tc, "descr", "Test log10(+0.0) == -HUGE_VAL");
128 }
129 
130 ATF_TC_BODY(log10_zero_pos, tc)
131 {
132 	const double x = 0.0L;
133 
134 	ATF_CHECK(log10(x) == -HUGE_VAL);
135 }
136 
137 /*
138  * log10f(3)
139  */
140 ATF_TC(log10f_base);
141 ATF_TC_HEAD(log10f_base, tc)
142 {
143 	atf_tc_set_md_var(tc, "descr", "Test log10f(10) == 1");
144 }
145 
146 ATF_TC_BODY(log10f_base, tc)
147 {
148 	ATF_CHECK(log10f(10.0) == 1.0);
149 }
150 
151 ATF_TC(log10f_nan);
152 ATF_TC_HEAD(log10f_nan, tc)
153 {
154 	atf_tc_set_md_var(tc, "descr", "Test log10f(NaN) == NaN");
155 }
156 
157 ATF_TC_BODY(log10f_nan, tc)
158 {
159 	const float x = 0.0L / 0.0L;
160 
161 	ATF_CHECK(isnan(x) != 0);
162 	ATF_CHECK(isnan(log10f(x)) != 0);
163 }
164 
165 ATF_TC(log10f_inf_neg);
166 ATF_TC_HEAD(log10f_inf_neg, tc)
167 {
168 	atf_tc_set_md_var(tc, "descr", "Test log10f(-Inf) == NaN");
169 }
170 
171 ATF_TC_BODY(log10f_inf_neg, tc)
172 {
173 	const float x = -1.0L / 0.0L;
174 	const float y = log10f(x);
175 
176 	ATF_CHECK(isnan(y) != 0);
177 }
178 
179 ATF_TC(log10f_inf_pos);
180 ATF_TC_HEAD(log10f_inf_pos, tc)
181 {
182 	atf_tc_set_md_var(tc, "descr", "Test log10f(+Inf) == +Inf");
183 }
184 
185 ATF_TC_BODY(log10f_inf_pos, tc)
186 {
187 	const float x = 1.0L / 0.0L;
188 
189 	ATF_CHECK(log10f(x) == x);
190 }
191 
192 ATF_TC(log10f_one_pos);
193 ATF_TC_HEAD(log10f_one_pos, tc)
194 {
195 	atf_tc_set_md_var(tc, "descr", "Test log10f(1.0) == +0.0");
196 }
197 
198 ATF_TC_BODY(log10f_one_pos, tc)
199 {
200 	const float x = log10f(1.0);
201 	const float y = 0.0L;
202 
203 	ATF_CHECK(x == y);
204 	ATF_CHECK(signbit(x) == 0);
205 	ATF_CHECK(signbit(y) == 0);
206 }
207 
208 ATF_TC(log10f_zero_neg);
209 ATF_TC_HEAD(log10f_zero_neg, tc)
210 {
211 	atf_tc_set_md_var(tc, "descr", "Test log10f(-0.0) == -HUGE_VALF");
212 }
213 
214 ATF_TC_BODY(log10f_zero_neg, tc)
215 {
216 	const float x = -0.0L;
217 
218 	ATF_CHECK(log10f(x) == -HUGE_VALF);
219 }
220 
221 ATF_TC(log10f_zero_pos);
222 ATF_TC_HEAD(log10f_zero_pos, tc)
223 {
224 	atf_tc_set_md_var(tc, "descr", "Test log10f(+0.0) == -HUGE_VALF");
225 }
226 
227 ATF_TC_BODY(log10f_zero_pos, tc)
228 {
229 	const float x = 0.0L;
230 
231 	ATF_CHECK(log10f(x) == -HUGE_VALF);
232 }
233 
234 /*
235  * log1p(3)
236  */
237 ATF_TC(log1p_nan);
238 ATF_TC_HEAD(log1p_nan, tc)
239 {
240 	atf_tc_set_md_var(tc, "descr", "Test log1p(NaN) == NaN");
241 }
242 
243 ATF_TC_BODY(log1p_nan, tc)
244 {
245 	const double x = 0.0L / 0.0L;
246 
247 	ATF_CHECK(isnan(x) != 0);
248 	ATF_CHECK(isnan(log1p(x)) != 0);
249 }
250 
251 ATF_TC(log1p_inf_neg);
252 ATF_TC_HEAD(log1p_inf_neg, tc)
253 {
254 	atf_tc_set_md_var(tc, "descr", "Test log1p(-Inf) == NaN");
255 }
256 
257 ATF_TC_BODY(log1p_inf_neg, tc)
258 {
259 	const double x = -1.0L / 0.0L;
260 	const double y = log1p(x);
261 
262 	if (isnan(y) == 0) {
263 		atf_tc_expect_fail("PR lib/45362");
264 		atf_tc_fail("log1p(-Inf) != NaN");
265 	}
266 }
267 
268 ATF_TC(log1p_inf_pos);
269 ATF_TC_HEAD(log1p_inf_pos, tc)
270 {
271 	atf_tc_set_md_var(tc, "descr", "Test log1p(+Inf) == +Inf");
272 }
273 
274 ATF_TC_BODY(log1p_inf_pos, tc)
275 {
276 	const double x = 1.0L / 0.0L;
277 
278 	ATF_CHECK(log1p(x) == x);
279 }
280 
281 ATF_TC(log1p_one_neg);
282 ATF_TC_HEAD(log1p_one_neg, tc)
283 {
284 	atf_tc_set_md_var(tc, "descr", "Test log1p(-1.0) == -HUGE_VAL");
285 }
286 
287 ATF_TC_BODY(log1p_one_neg, tc)
288 {
289 	const double x = log1p(-1.0);
290 
291 	if (x != -HUGE_VAL) {
292 		atf_tc_expect_fail("PR lib/45362");
293 		atf_tc_fail("log1p(-1.0) != -HUGE_VAL");
294 	}
295 }
296 
297 ATF_TC(log1p_zero_neg);
298 ATF_TC_HEAD(log1p_zero_neg, tc)
299 {
300 	atf_tc_set_md_var(tc, "descr", "Test log1p(-0.0) == -0.0");
301 }
302 
303 ATF_TC_BODY(log1p_zero_neg, tc)
304 {
305 	const double x = -0.0L;
306 
307 	ATF_CHECK(log1p(x) == x);
308 }
309 
310 ATF_TC(log1p_zero_pos);
311 ATF_TC_HEAD(log1p_zero_pos, tc)
312 {
313 	atf_tc_set_md_var(tc, "descr", "Test log1p(+0.0) == +0.0");
314 }
315 
316 ATF_TC_BODY(log1p_zero_pos, tc)
317 {
318 	const double x = 0.0L;
319 
320 	ATF_CHECK(log1p(x) == x);
321 }
322 
323 /*
324  * log1pf(3)
325  */
326 ATF_TC(log1pf_nan);
327 ATF_TC_HEAD(log1pf_nan, tc)
328 {
329 	atf_tc_set_md_var(tc, "descr", "Test log1pf(NaN) == NaN");
330 }
331 
332 ATF_TC_BODY(log1pf_nan, tc)
333 {
334 	const float x = 0.0L / 0.0L;
335 
336 	ATF_CHECK(isnan(x) != 0);
337 	ATF_CHECK(isnan(log1pf(x)) != 0);
338 }
339 
340 ATF_TC(log1pf_inf_neg);
341 ATF_TC_HEAD(log1pf_inf_neg, tc)
342 {
343 	atf_tc_set_md_var(tc, "descr", "Test log1pf(-Inf) == NaN");
344 }
345 
346 ATF_TC_BODY(log1pf_inf_neg, tc)
347 {
348 	const float x = -1.0L / 0.0L;
349 	const float y = log1pf(x);
350 
351 	if (isnan(y) == 0) {
352 		atf_tc_expect_fail("PR lib/45362");
353 		atf_tc_fail("log1pf(-Inf) != NaN");
354 	}
355 }
356 
357 ATF_TC(log1pf_inf_pos);
358 ATF_TC_HEAD(log1pf_inf_pos, tc)
359 {
360 	atf_tc_set_md_var(tc, "descr", "Test log1pf(+Inf) == +Inf");
361 }
362 
363 ATF_TC_BODY(log1pf_inf_pos, tc)
364 {
365 	const float x = 1.0L / 0.0L;
366 
367 	ATF_CHECK(log1pf(x) == x);
368 }
369 
370 ATF_TC(log1pf_one_neg);
371 ATF_TC_HEAD(log1pf_one_neg, tc)
372 {
373 	atf_tc_set_md_var(tc, "descr", "Test log1pf(-1.0) == -HUGE_VALF");
374 }
375 
376 ATF_TC_BODY(log1pf_one_neg, tc)
377 {
378 	const float x = log1pf(-1.0);
379 
380 	if (x != -HUGE_VALF) {
381 		atf_tc_expect_fail("PR lib/45362");
382 		atf_tc_fail("log1pf(-1.0) != -HUGE_VALF");
383 	}
384 }
385 
386 ATF_TC(log1pf_zero_neg);
387 ATF_TC_HEAD(log1pf_zero_neg, tc)
388 {
389 	atf_tc_set_md_var(tc, "descr", "Test log1pf(-0.0) == -0.0");
390 }
391 
392 ATF_TC_BODY(log1pf_zero_neg, tc)
393 {
394 	const float x = -0.0L;
395 
396 	ATF_CHECK(log1pf(x) == x);
397 }
398 
399 ATF_TC(log1pf_zero_pos);
400 ATF_TC_HEAD(log1pf_zero_pos, tc)
401 {
402 	atf_tc_set_md_var(tc, "descr", "Test log1pf(+0.0) == +0.0");
403 }
404 
405 ATF_TC_BODY(log1pf_zero_pos, tc)
406 {
407 	const float x = 0.0L;
408 
409 	ATF_CHECK(log1pf(x) == x);
410 }
411 
412 /*
413  * log2(3)
414  */
415 ATF_TC(log2_base);
416 ATF_TC_HEAD(log2_base, tc)
417 {
418 	atf_tc_set_md_var(tc, "descr", "Test log2(2) == 1");
419 }
420 
421 ATF_TC_BODY(log2_base, tc)
422 {
423 	ATF_CHECK(log2(2.0) == 1.0);
424 }
425 
426 ATF_TC(log2_nan);
427 ATF_TC_HEAD(log2_nan, tc)
428 {
429 	atf_tc_set_md_var(tc, "descr", "Test log2(NaN) == NaN");
430 }
431 
432 ATF_TC_BODY(log2_nan, tc)
433 {
434 	const double x = 0.0L / 0.0L;
435 
436 	ATF_CHECK(isnan(x) != 0);
437 	ATF_CHECK(isnan(log2(x)) != 0);
438 }
439 
440 ATF_TC(log2_inf_neg);
441 ATF_TC_HEAD(log2_inf_neg, tc)
442 {
443 	atf_tc_set_md_var(tc, "descr", "Test log2(-Inf) == NaN");
444 }
445 
446 ATF_TC_BODY(log2_inf_neg, tc)
447 {
448 	const double x = -1.0L / 0.0L;
449 	const double y = log2(x);
450 
451 	ATF_CHECK(isnan(y) != 0);
452 }
453 
454 ATF_TC(log2_inf_pos);
455 ATF_TC_HEAD(log2_inf_pos, tc)
456 {
457 	atf_tc_set_md_var(tc, "descr", "Test log2(+Inf) == +Inf");
458 }
459 
460 ATF_TC_BODY(log2_inf_pos, tc)
461 {
462 	const double x = 1.0L / 0.0L;
463 
464 	ATF_CHECK(log2(x) == x);
465 }
466 
467 ATF_TC(log2_one_pos);
468 ATF_TC_HEAD(log2_one_pos, tc)
469 {
470 	atf_tc_set_md_var(tc, "descr", "Test log2(1.0) == +0.0");
471 }
472 
473 ATF_TC_BODY(log2_one_pos, tc)
474 {
475 	const double x = log2(1.0);
476 	const double y = 0.0L;
477 
478 	ATF_CHECK(x == y);
479 	ATF_CHECK(signbit(x) == 0);
480 	ATF_CHECK(signbit(y) == 0);
481 }
482 
483 ATF_TC(log2_zero_neg);
484 ATF_TC_HEAD(log2_zero_neg, tc)
485 {
486 	atf_tc_set_md_var(tc, "descr", "Test log2(-0.0) == -HUGE_VAL");
487 }
488 
489 ATF_TC_BODY(log2_zero_neg, tc)
490 {
491 	const double x = -0.0L;
492 
493 	ATF_CHECK(log2(x) == -HUGE_VAL);
494 }
495 
496 ATF_TC(log2_zero_pos);
497 ATF_TC_HEAD(log2_zero_pos, tc)
498 {
499 	atf_tc_set_md_var(tc, "descr", "Test log2(+0.0) == -HUGE_VAL");
500 }
501 
502 ATF_TC_BODY(log2_zero_pos, tc)
503 {
504 	const double x = 0.0L;
505 
506 	ATF_CHECK(log2(x) == -HUGE_VAL);
507 }
508 
509 /*
510  * log2f(3)
511  */
512 ATF_TC(log2f_base);
513 ATF_TC_HEAD(log2f_base, tc)
514 {
515 	atf_tc_set_md_var(tc, "descr", "Test log2f(2) == 1");
516 }
517 
518 ATF_TC_BODY(log2f_base, tc)
519 {
520 	ATF_CHECK(log2f(2.0) == 1.0);
521 }
522 
523 ATF_TC(log2f_nan);
524 ATF_TC_HEAD(log2f_nan, tc)
525 {
526 	atf_tc_set_md_var(tc, "descr", "Test log2f(NaN) == NaN");
527 }
528 
529 ATF_TC_BODY(log2f_nan, tc)
530 {
531 	const float x = 0.0L / 0.0L;
532 
533 	ATF_CHECK(isnan(x) != 0);
534 	ATF_CHECK(isnan(log2f(x)) != 0);
535 }
536 
537 ATF_TC(log2f_inf_neg);
538 ATF_TC_HEAD(log2f_inf_neg, tc)
539 {
540 	atf_tc_set_md_var(tc, "descr", "Test log2f(-Inf) == NaN");
541 }
542 
543 ATF_TC_BODY(log2f_inf_neg, tc)
544 {
545 	const float x = -1.0L / 0.0L;
546 	const float y = log2f(x);
547 
548 	ATF_CHECK(isnan(y) != 0);
549 }
550 
551 ATF_TC(log2f_inf_pos);
552 ATF_TC_HEAD(log2f_inf_pos, tc)
553 {
554 	atf_tc_set_md_var(tc, "descr", "Test log2f(+Inf) == +Inf");
555 }
556 
557 ATF_TC_BODY(log2f_inf_pos, tc)
558 {
559 	const float x = 1.0L / 0.0L;
560 
561 	ATF_CHECK(log2f(x) == x);
562 }
563 
564 ATF_TC(log2f_one_pos);
565 ATF_TC_HEAD(log2f_one_pos, tc)
566 {
567 	atf_tc_set_md_var(tc, "descr", "Test log2f(1.0) == +0.0");
568 }
569 
570 ATF_TC_BODY(log2f_one_pos, tc)
571 {
572 	const float x = log2f(1.0);
573 	const float y = 0.0L;
574 
575 	ATF_CHECK(x == y);
576 	ATF_CHECK(signbit(x) == 0);
577 	ATF_CHECK(signbit(y) == 0);
578 }
579 
580 ATF_TC(log2f_zero_neg);
581 ATF_TC_HEAD(log2f_zero_neg, tc)
582 {
583 	atf_tc_set_md_var(tc, "descr", "Test log2f(-0.0) == -HUGE_VALF");
584 }
585 
586 ATF_TC_BODY(log2f_zero_neg, tc)
587 {
588 	const float x = -0.0L;
589 
590 	ATF_CHECK(log2f(x) == -HUGE_VALF);
591 }
592 
593 ATF_TC(log2f_zero_pos);
594 ATF_TC_HEAD(log2f_zero_pos, tc)
595 {
596 	atf_tc_set_md_var(tc, "descr", "Test log2f(+0.0) == -HUGE_VALF");
597 }
598 
599 ATF_TC_BODY(log2f_zero_pos, tc)
600 {
601 	const float x = 0.0L;
602 
603 	ATF_CHECK(log2f(x) == -HUGE_VALF);
604 }
605 
606 /*
607  * log(3)
608  */
609 ATF_TC(log_base);
610 ATF_TC_HEAD(log_base, tc)
611 {
612 	atf_tc_set_md_var(tc, "descr", "Test log(e) == 1");
613 }
614 
615 ATF_TC_BODY(log_base, tc)
616 {
617 	const double eps = 1.0e-38;
618 
619 	if (fabs(log(M_E) - 1.0) > eps)
620 		atf_tc_fail_nonfatal("log(e) != 1");
621 }
622 
623 ATF_TC(log_nan);
624 ATF_TC_HEAD(log_nan, tc)
625 {
626 	atf_tc_set_md_var(tc, "descr", "Test log(NaN) == NaN");
627 }
628 
629 ATF_TC_BODY(log_nan, tc)
630 {
631 	const double x = 0.0L / 0.0L;
632 
633 	ATF_CHECK(isnan(x) != 0);
634 	ATF_CHECK(isnan(log(x)) != 0);
635 }
636 
637 ATF_TC(log_inf_neg);
638 ATF_TC_HEAD(log_inf_neg, tc)
639 {
640 	atf_tc_set_md_var(tc, "descr", "Test log(-Inf) == NaN");
641 }
642 
643 ATF_TC_BODY(log_inf_neg, tc)
644 {
645 	const double x = -1.0L / 0.0L;
646 	const double y = log(x);
647 
648 	ATF_CHECK(isnan(y) != 0);
649 }
650 
651 ATF_TC(log_inf_pos);
652 ATF_TC_HEAD(log_inf_pos, tc)
653 {
654 	atf_tc_set_md_var(tc, "descr", "Test log(+Inf) == +Inf");
655 }
656 
657 ATF_TC_BODY(log_inf_pos, tc)
658 {
659 	const double x = 1.0L / 0.0L;
660 
661 	ATF_CHECK(log(x) == x);
662 }
663 
664 ATF_TC(log_one_pos);
665 ATF_TC_HEAD(log_one_pos, tc)
666 {
667 	atf_tc_set_md_var(tc, "descr", "Test log(1.0) == +0.0");
668 }
669 
670 ATF_TC_BODY(log_one_pos, tc)
671 {
672 	const double x = log(1.0);
673 	const double y = 0.0L;
674 
675 	ATF_CHECK(x == y);
676 	ATF_CHECK(signbit(x) == 0);
677 	ATF_CHECK(signbit(y) == 0);
678 }
679 
680 ATF_TC(log_zero_neg);
681 ATF_TC_HEAD(log_zero_neg, tc)
682 {
683 	atf_tc_set_md_var(tc, "descr", "Test log(-0.0) == -HUGE_VAL");
684 }
685 
686 ATF_TC_BODY(log_zero_neg, tc)
687 {
688 	const double x = -0.0L;
689 
690 	ATF_CHECK(log(x) == -HUGE_VAL);
691 }
692 
693 ATF_TC(log_zero_pos);
694 ATF_TC_HEAD(log_zero_pos, tc)
695 {
696 	atf_tc_set_md_var(tc, "descr", "Test log(+0.0) == -HUGE_VAL");
697 }
698 
699 ATF_TC_BODY(log_zero_pos, tc)
700 {
701 	const double x = 0.0L;
702 
703 	ATF_CHECK(log(x) == -HUGE_VAL);
704 }
705 
706 /*
707  * logf(3)
708  */
709 ATF_TC(logf_base);
710 ATF_TC_HEAD(logf_base, tc)
711 {
712 	atf_tc_set_md_var(tc, "descr", "Test logf(e) == 1");
713 }
714 
715 ATF_TC_BODY(logf_base, tc)
716 {
717 	const float eps = 1.0e-7;
718 
719 	if (fabsf(logf(M_E) - 1.0f) > eps)
720 		atf_tc_fail_nonfatal("logf(e) != 1");
721 }
722 
723 ATF_TC(logf_nan);
724 ATF_TC_HEAD(logf_nan, tc)
725 {
726 	atf_tc_set_md_var(tc, "descr", "Test logf(NaN) == NaN");
727 }
728 
729 ATF_TC_BODY(logf_nan, tc)
730 {
731 	const float x = 0.0L / 0.0L;
732 
733 	ATF_CHECK(isnan(x) != 0);
734 	ATF_CHECK(isnan(logf(x)) != 0);
735 }
736 
737 ATF_TC(logf_inf_neg);
738 ATF_TC_HEAD(logf_inf_neg, tc)
739 {
740 	atf_tc_set_md_var(tc, "descr", "Test logf(-Inf) == NaN");
741 }
742 
743 ATF_TC_BODY(logf_inf_neg, tc)
744 {
745 	const float x = -1.0L / 0.0L;
746 	const float y = logf(x);
747 
748 	ATF_CHECK(isnan(y) != 0);
749 }
750 
751 ATF_TC(logf_inf_pos);
752 ATF_TC_HEAD(logf_inf_pos, tc)
753 {
754 	atf_tc_set_md_var(tc, "descr", "Test logf(+Inf) == +Inf");
755 }
756 
757 ATF_TC_BODY(logf_inf_pos, tc)
758 {
759 	const float x = 1.0L / 0.0L;
760 
761 	ATF_CHECK(logf(x) == x);
762 }
763 
764 ATF_TC(logf_one_pos);
765 ATF_TC_HEAD(logf_one_pos, tc)
766 {
767 	atf_tc_set_md_var(tc, "descr", "Test logf(1.0) == +0.0");
768 }
769 
770 ATF_TC_BODY(logf_one_pos, tc)
771 {
772 	const float x = logf(1.0);
773 	const float y = 0.0L;
774 
775 	ATF_CHECK(x == y);
776 	ATF_CHECK(signbit(x) == 0);
777 	ATF_CHECK(signbit(y) == 0);
778 }
779 
780 ATF_TC(logf_zero_neg);
781 ATF_TC_HEAD(logf_zero_neg, tc)
782 {
783 	atf_tc_set_md_var(tc, "descr", "Test logf(-0.0) == -HUGE_VALF");
784 }
785 
786 ATF_TC_BODY(logf_zero_neg, tc)
787 {
788 	const float x = -0.0L;
789 
790 	ATF_CHECK(logf(x) == -HUGE_VALF);
791 }
792 
793 ATF_TC(logf_zero_pos);
794 ATF_TC_HEAD(logf_zero_pos, tc)
795 {
796 	atf_tc_set_md_var(tc, "descr", "Test logf(+0.0) == -HUGE_VALF");
797 }
798 
799 ATF_TC_BODY(logf_zero_pos, tc)
800 {
801 	const float x = 0.0L;
802 
803 	ATF_CHECK(logf(x) == -HUGE_VALF);
804 }
805 
806 ATF_TP_ADD_TCS(tp)
807 {
808 
809 	ATF_TP_ADD_TC(tp, log10_base);
810 	ATF_TP_ADD_TC(tp, log10_nan);
811 	ATF_TP_ADD_TC(tp, log10_inf_neg);
812 	ATF_TP_ADD_TC(tp, log10_inf_pos);
813 	ATF_TP_ADD_TC(tp, log10_one_pos);
814 	ATF_TP_ADD_TC(tp, log10_zero_neg);
815 	ATF_TP_ADD_TC(tp, log10_zero_pos);
816 
817 	ATF_TP_ADD_TC(tp, log10f_base);
818 	ATF_TP_ADD_TC(tp, log10f_nan);
819 	ATF_TP_ADD_TC(tp, log10f_inf_neg);
820 	ATF_TP_ADD_TC(tp, log10f_inf_pos);
821 	ATF_TP_ADD_TC(tp, log10f_one_pos);
822 	ATF_TP_ADD_TC(tp, log10f_zero_neg);
823 	ATF_TP_ADD_TC(tp, log10f_zero_pos);
824 
825 	ATF_TP_ADD_TC(tp, log1p_nan);
826 	ATF_TP_ADD_TC(tp, log1p_inf_neg);
827 	ATF_TP_ADD_TC(tp, log1p_inf_pos);
828 	ATF_TP_ADD_TC(tp, log1p_one_neg);
829 	ATF_TP_ADD_TC(tp, log1p_zero_neg);
830 	ATF_TP_ADD_TC(tp, log1p_zero_pos);
831 
832 	ATF_TP_ADD_TC(tp, log1pf_nan);
833 	ATF_TP_ADD_TC(tp, log1pf_inf_neg);
834 	ATF_TP_ADD_TC(tp, log1pf_inf_pos);
835 	ATF_TP_ADD_TC(tp, log1pf_one_neg);
836 	ATF_TP_ADD_TC(tp, log1pf_zero_neg);
837 	ATF_TP_ADD_TC(tp, log1pf_zero_pos);
838 
839 	ATF_TP_ADD_TC(tp, log2_base);
840 	ATF_TP_ADD_TC(tp, log2_nan);
841 	ATF_TP_ADD_TC(tp, log2_inf_neg);
842 	ATF_TP_ADD_TC(tp, log2_inf_pos);
843 	ATF_TP_ADD_TC(tp, log2_one_pos);
844 	ATF_TP_ADD_TC(tp, log2_zero_neg);
845 	ATF_TP_ADD_TC(tp, log2_zero_pos);
846 
847 	ATF_TP_ADD_TC(tp, log2f_base);
848 	ATF_TP_ADD_TC(tp, log2f_nan);
849 	ATF_TP_ADD_TC(tp, log2f_inf_neg);
850 	ATF_TP_ADD_TC(tp, log2f_inf_pos);
851 	ATF_TP_ADD_TC(tp, log2f_one_pos);
852 	ATF_TP_ADD_TC(tp, log2f_zero_neg);
853 	ATF_TP_ADD_TC(tp, log2f_zero_pos);
854 
855 	ATF_TP_ADD_TC(tp, log_base);
856 	ATF_TP_ADD_TC(tp, log_nan);
857 	ATF_TP_ADD_TC(tp, log_inf_neg);
858 	ATF_TP_ADD_TC(tp, log_inf_pos);
859 	ATF_TP_ADD_TC(tp, log_one_pos);
860 	ATF_TP_ADD_TC(tp, log_zero_neg);
861 	ATF_TP_ADD_TC(tp, log_zero_pos);
862 
863 	ATF_TP_ADD_TC(tp, logf_base);
864 	ATF_TP_ADD_TC(tp, logf_nan);
865 	ATF_TP_ADD_TC(tp, logf_inf_neg);
866 	ATF_TP_ADD_TC(tp, logf_inf_pos);
867 	ATF_TP_ADD_TC(tp, logf_one_pos);
868 	ATF_TP_ADD_TC(tp, logf_zero_neg);
869 	ATF_TP_ADD_TC(tp, logf_zero_pos);
870 
871 	return atf_no_error();
872 }
873