fa3f3bf9792730ac2c3096f9e0f56e3c6e7322c0
[kopensolaris-gnu/glibc.git] / sysdeps / m68k / fpu / __math.h
1 /* Copyright (C) 1991, 1992 Free Software Foundation, Inc.
2 This file is part of the GNU C Library.
3
4 The GNU C Library is free software; you can redistribute it and/or
5 modify it under the terms of the GNU Library General Public License as
6 published by the Free Software Foundation; either version 2 of the
7 License, or (at your option) any later version.
8
9 The GNU C Library is distributed in the hope that it will be useful,
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12 Library General Public License for more details.
13
14 You should have received a copy of the GNU Library General Public
15 License along with the GNU C Library; see the file COPYING.LIB.  If
16 not, write to the Free Software Foundation, Inc., 675 Mass Ave,
17 Cambridge, MA 02139, USA.  */
18
19 #ifdef  __GNUC__
20
21 #include <sys/cdefs.h>
22
23 #ifdef  __NO_MATH_INLINES
24 #define __m81_u(x)      __CONCAT(__,x)
25 #else
26 #define __m81_u(x)      x
27 #define __MATH_INLINES  1
28 #endif
29
30 #define __inline_mathop2(func, op)                                            \
31   extern __inline __const double                                              \
32   __m81_u(func)(double __mathop_x)                                            \
33   {                                                                           \
34     double __result;                                                          \
35     __asm("f" __STRING(op) "%.x %1, %0" : "=f" (__result) : "f" (__mathop_x));\
36     return __result;                                                          \
37   }
38 #define __inline_mathop(op)             __inline_mathop2(op, op)
39
40 __inline_mathop(acos)
41 __inline_mathop(asin)
42 __inline_mathop(atan)
43 __inline_mathop(cos)
44 __inline_mathop(sin)
45 __inline_mathop(tan)
46 __inline_mathop(cosh)
47 __inline_mathop(sinh)
48 __inline_mathop(tanh)
49 __inline_mathop2(exp, etox)
50 __inline_mathop2(fabs, abs)
51 __inline_mathop(log10)
52 __inline_mathop2(log, logn)
53 __inline_mathop2(floor, intrz)
54 __inline_mathop(sqrt)
55
56 #ifdef  __USE_MISC
57 __inline_mathop2(rint, int)
58 __inline_mathop2(expm1, etoxm1)
59 __inline_mathop2(log1p, lognp1)
60 __inline_mathop(atanh)
61 #endif
62
63 extern __inline __const double
64 __m81_u(__drem)(double __x, double __y)
65 {
66   double __result;
67   __asm("frem%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x));
68   return __result;
69 }
70
71 extern __inline __const double
72 __m81_u(ldexp)(double __x, int __e)
73 {
74   double __result;
75   double __double_e = (double) __e;
76   __asm("fscale%.x %1, %0" : "=f" (__result) : "f" (__double_e), "0" (__x));
77   return __result;
78 }
79
80 extern __inline __const double
81 __m81_u(fmod)(double __x, double __y)
82 {
83   double __result;
84   __asm("fmod%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x));
85   return __result;
86 }
87
88 extern __inline double
89 __m81_u(frexp)(double __value, int *__expptr)
90 {
91   double __mantissa, __exponent;
92   __asm("fgetexp%.x %1, %0" : "=f" (__exponent) : "f" (__value));
93   __asm("fgetman%.x %1, %0" : "=f" (__mantissa) : "f" (__value));
94   *__expptr = (int) __exponent;
95   return __mantissa;
96 }
97
98 extern __inline __const double
99 __m81_u(pow)(double __x, double __y)
100 {
101   double __result;
102   if (__y == 0.0 || __x == 1.0)
103     __result = 1.0;
104   else if (__y == 1.0)
105     __result = __x;
106   else if (__y == 2.0)
107     __result = __x * __x;
108   else if (__x == 10.0)
109     __asm("ftentox%.x %1, %0" : "=f" (__result) : "f" (__y));
110   else if (__x == 2.0)
111     __asm("ftwotox%.x %1, %0" : "=f" (__result) : "f" (__y));
112   else
113     __result = __m81_u(exp)(__y * __m81_u(log)(__x));
114   return __result;
115 }
116
117 extern __inline __const double
118 __m81_u(ceil)(double __x)
119 {
120   double __result;
121   unsigned long int __ctrl_reg;
122   __asm("fmove%.l fpcr, %0" : "=g" (__ctrl_reg));
123   /* Set rounding towards positive infinity.  */
124   __asm("fmove%.l %0, fpcr" : /* No outputs.  */ : "g" (__ctrl_reg | 0x30));
125   /* Convert X to an integer, using +Inf rounding.  */
126   __asm("fint%.x %1, %0" : "=f" (__result) : "f" (__x));
127   /* Restore the previous rounding mode.  */
128   __asm("fmove%.l %0, fpcr" : /* No outputs.  */ : "g" (__ctrl_reg));
129   return __result;
130 }
131
132 extern __inline double
133 __m81_u(modf)(double __value, double *__iptr)
134 {
135   double __modf_int = __m81_u(floor)(__value);
136   *__iptr = __modf_int;
137   return __value - __modf_int;
138 }
139
140 extern __inline int
141 __m81_u(__isinf)(double __value)
142 {
143   /* There is no branch-condition for infinity,
144      so we must extract and examine the condition codes manually.  */
145   unsigned long int __fpsr;
146   __asm("ftst%.x %1\n"
147         "fmove%.l fpsr, %0" : "=g" (__fpsr) : "f" (__value));
148   return (__fpsr & (2 << (3 * 8))) ? (__value < 0 ? -1 : 1) : 0;
149 }
150
151 extern __inline int
152 __m81_u(__isnan)(double __value)
153 {
154   char __result;
155   __asm("ftst%.x %1\n"
156         "fsun %0" : "=g" (__result) : "f" (__value));
157   return __result;
158 }
159
160 #endif  /* GCC.  */