24 */
25 /*
26 * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
27 * Use is subject to license terms.
28 */
29
30 #if defined(ELFOBJ)
31 #pragma weak fmin = __fmin
32 #endif
33
34 /*
35 * fmin(x,y) returns the smaller of x and y. If just one of the
36 * arguments is NaN, fmin returns the other argument. If both
37 * arguments are NaN, fmin returns NaN.
38 *
39 * See fmaxf.c for a discussion of implementation trade-offs.
40 */
41
42 #include "libm.h" /* for islessequal macro */
43
44 #include <sys/isa_defs.h>
45
46 double
47 __fmin(double x, double y) {
48 union {
49 unsigned i[2];
50 double d;
51 } xx, yy;
52 unsigned s;
53
54 /* if y is nan, replace it by x */
55 if (y != y)
56 y = x;
57
58 /* if x is nan, replace it by y */
59 if (x != x)
60 x = y;
61
62 /* if x is greater than y or x and y are unordered, replace x by y */
63 #if defined(COMPARISON_MACRO_BUG)
64 if (x != x || x > y)
65 #else
66 if (!islessequal(x, y))
67 #endif
68 x = y;
69
70 /*
71 * now x and y are either both NaN or both numeric; set the
72 * sign of the result if either x or y has its sign set
73 */
74 xx.d = x;
75 yy.d = y;
76 #if defined(_BIG_ENDIAN)
77 s = (xx.i[0] | yy.i[0]) & 0x80000000;
78 xx.i[0] |= s;
79 #else
80 s = (xx.i[1] | yy.i[1]) & 0x80000000;
81 xx.i[1] |= s;
82 #endif
83
84 return (xx.d);
85 }
|
24 */
25 /*
26 * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
27 * Use is subject to license terms.
28 */
29
30 #if defined(ELFOBJ)
31 #pragma weak fmin = __fmin
32 #endif
33
34 /*
35 * fmin(x,y) returns the smaller of x and y. If just one of the
36 * arguments is NaN, fmin returns the other argument. If both
37 * arguments are NaN, fmin returns NaN.
38 *
39 * See fmaxf.c for a discussion of implementation trade-offs.
40 */
41
42 #include "libm.h" /* for islessequal macro */
43
44 #include "fenv_inlines.h"
45 #include <stdio.h>
46 #include <sys/isa_defs.h>
47
48 double
49 __fmin(double x, double y) {
50 union {
51 unsigned i[2];
52 double d;
53 } xx, yy;
54 unsigned s;
55
56 /* if y is nan, replace it by x */
57 if (y != y)
58 y = x;
59
60 /* if x is nan, replace it by y */
61 if (x != x)
62 x = y;
63
64 /* At this point, x and y are either both numeric, or both NaN */
65 if (!isnan(x) && !islessequal(x, y))
66 x = y;
67
68 /*
69 * set the sign of the result if either x or y has its sign set
70 */
71 xx.d = x;
72 yy.d = y;
73 #if defined(_BIG_ENDIAN)
74 s = (xx.i[0] | yy.i[0]) & 0x80000000;
75 xx.i[0] |= s;
76 #else
77 s = (xx.i[1] | yy.i[1]) & 0x80000000;
78 xx.i[1] |= s;
79 #endif
80
81 return (xx.d);
82 }
|