1 | // -*- Mode: C -*-
|
---|
2 | //
|
---|
3 | // Cforall Version 1.0.0 Copyright (C) 2016 University of Waterloo
|
---|
4 | //
|
---|
5 | // The contents of this file are covered under the licence agreement in the
|
---|
6 | // file "LICENCE" distributed with Cforall.
|
---|
7 | //
|
---|
8 | // math.c --
|
---|
9 | //
|
---|
10 | // Author : Peter A. Buhr
|
---|
11 | // Created On : Tue Apr 19 22:23:08 2016
|
---|
12 | // Last Modified By : Peter A. Buhr
|
---|
13 | // Last Modified On : Thu Apr 21 07:56:49 2016
|
---|
14 | // Update Count : 51
|
---|
15 | //
|
---|
16 |
|
---|
17 | #include "math"
|
---|
18 |
|
---|
19 | //---------------------------------------
|
---|
20 |
|
---|
21 | extern "C" {
|
---|
22 | #include <math.h>
|
---|
23 | #include <complex.h>
|
---|
24 | } // extern "C"
|
---|
25 |
|
---|
26 | float cos( float x ) { return cosf( x ); }
|
---|
27 | long double cos( long double x ) { return cosl( x ); }
|
---|
28 | float _Complex cos( float _Complex x ) { return ccosf( x ); }
|
---|
29 | double _Complex cos( double _Complex x ) { return ccos( x ); }
|
---|
30 | long double _Complex cos( long double _Complex x ) { return ccosl( x ); }
|
---|
31 |
|
---|
32 | float cosh( float x ) { return coshf( x ); }
|
---|
33 | long double cosh( long double x ) { return coshl( x ); }
|
---|
34 | float _Complex cosh( float _Complex x ) { return ccoshf( x ); }
|
---|
35 | double _Complex cosh( double _Complex x ) { return ccosh( x ); }
|
---|
36 | long double _Complex cosh( long double _Complex x ) { return ccoshl( x ); }
|
---|
37 |
|
---|
38 | float acos( float x ) { return acosf( x ); }
|
---|
39 | long double acos( long double x ) { return acosl( x ); }
|
---|
40 | float _Complex acos( float _Complex x ) { return cacosf( x ); }
|
---|
41 | double _Complex acos( double _Complex x ) { return cacos( x ); }
|
---|
42 | long double _Complex acos( long double _Complex x ) { return cacosl( x ); }
|
---|
43 |
|
---|
44 | float acosh( float x ) { return acoshf( x ); }
|
---|
45 | long double acosh( long double x ) { return acoshl( x ); }
|
---|
46 | float _Complex acosh( float _Complex x ) { return cacoshf( x ); }
|
---|
47 | double _Complex acosh( double _Complex x ) { return cacosh( x ); }
|
---|
48 | long double _Complex acosh( long double _Complex x ) { return cacoshl( x ); }
|
---|
49 |
|
---|
50 | float sin( float x ) { return sinf( x ); }
|
---|
51 | long double sin( long double x ) { return sinl( x ); }
|
---|
52 | float _Complex sin( float _Complex x ) { return csinf( x ); }
|
---|
53 | double _Complex sin( double _Complex x ) { return csin( x ); }
|
---|
54 | long double _Complex sin( long double _Complex x ) { return csinl( x ); }
|
---|
55 |
|
---|
56 | float sinh( float x ) { return sinhf( x ); }
|
---|
57 | long double sinh( long double x ) { return sinhl( x ); }
|
---|
58 | float _Complex sinh( float _Complex x ) { return csinhf( x ); }
|
---|
59 | double _Complex sinh( double _Complex x ) { return csinh( x ); }
|
---|
60 | long double _Complex sinh( long double _Complex x ) { return csinhl( x ); }
|
---|
61 |
|
---|
62 | float asin( float x ) { return asinf( x ); }
|
---|
63 | long double asin( long double x ) { return asinl( x ); }
|
---|
64 | float _Complex asin( float _Complex x ) { return casinf( x ); }
|
---|
65 | double _Complex asin( double _Complex x ) { return casin( x ); }
|
---|
66 | long double _Complex asin( long double _Complex x ) { return casinl( x ); }
|
---|
67 |
|
---|
68 | float asinh( float x ) { return asinhf( x ); }
|
---|
69 | long double asinh( long double x ) { return asinhl( x ); }
|
---|
70 | float _Complex asinh( float _Complex x ) { return casinhf( x ); }
|
---|
71 | double _Complex asinh( double _Complex x ) { return casinh( x ); }
|
---|
72 | long double _Complex asinh( long double _Complex x ) { return casinhl( x ); }
|
---|
73 |
|
---|
74 | float tan( float x ) { return tanf( x ); }
|
---|
75 | long double tan( long double x ) { return tanl( x ); }
|
---|
76 | float _Complex tan( float _Complex x ) { return ctanf( x ); }
|
---|
77 | double _Complex tan( double _Complex x ) { return ctan( x ); }
|
---|
78 | long double _Complex tan( long double _Complex x ) { return ctanl( x ); }
|
---|
79 |
|
---|
80 | float tanh( float x ) { return tanhf( x ); }
|
---|
81 | long double tanh( long double x ) { return tanhl( x ); }
|
---|
82 | float _Complex tanh( float _Complex x ) { return ctanhf( x ); }
|
---|
83 | double _Complex tanh( double _Complex x ) { return ctanh( x ); }
|
---|
84 | long double _Complex tanh( long double _Complex x ) { return ctanhl( x ); }
|
---|
85 |
|
---|
86 | float atan( float x ) { return atanf( x ); }
|
---|
87 | long double atan( long double x ) { return atanl( x ); }
|
---|
88 | float _Complex atan( float _Complex x ) { return catanf( x ); }
|
---|
89 | double _Complex atan( double _Complex x ) { return catan( x ); }
|
---|
90 | long double _Complex atan( long double _Complex x ) { return catanl( x ); }
|
---|
91 |
|
---|
92 | float atan( float x, float y ) { return atan2f( x, y ); }
|
---|
93 | double atan( double x, double y ) { return atan2( x, y ); }
|
---|
94 | long double atan( long double x, long double y ) { return atan2l( x, y ); }
|
---|
95 |
|
---|
96 | float atanh( float x ) { return atanhf( x ); }
|
---|
97 | long double atanh( long double x ) { return atanhl( x ); }
|
---|
98 | float _Complex atanh( float _Complex x ) { return catanhf( x ); }
|
---|
99 | double _Complex atanh( double _Complex x ) { return catanh( x ); }
|
---|
100 | long double _Complex atanh( long double _Complex x ) { return catanhl( x ); }
|
---|
101 |
|
---|
102 | float exp( float x ) { return expf( x ); }
|
---|
103 | long double exp( long double x ) { return expl( x ); }
|
---|
104 | float _Complex exp( float _Complex x ) { return cexpf( x ); }
|
---|
105 | double _Complex exp( double _Complex x ) { return cexp( x ); }
|
---|
106 | long double _Complex exp( long double _Complex x ) { return cexpl( x ); }
|
---|
107 |
|
---|
108 | float log( float x ) { return logf( x ); }
|
---|
109 | long double log( long double x ) { return logl( x ); }
|
---|
110 | float _Complex log( float _Complex x ) { return clogf( x ); }
|
---|
111 | double _Complex log( double _Complex x ) { return clog( x ); }
|
---|
112 | long double _Complex log( long double _Complex x ) { return clogl( x ); }
|
---|
113 |
|
---|
114 | float pow( float x, float y ) { return powf( x, y ); }
|
---|
115 | long double pow( long double x, long double y ) { return powl( x, y ); }
|
---|
116 | float _Complex pow( float _Complex x, float _Complex y ) { return cpowf( x, y ); }
|
---|
117 | double _Complex pow( double _Complex x, double _Complex y ) { return cpow( x, y ); }
|
---|
118 | long double _Complex pow( long double _Complex x, long double _Complex y ) { return cpowl( x, y ); }
|
---|
119 |
|
---|
120 | float sqrt( float x ) { return sqrtf( x ); }
|
---|
121 | long double sqrt( long double x ) { return sqrtl( x ); }
|
---|
122 | float _Complex sqrt( float _Complex x ) { return csqrtf( x ); }
|
---|
123 | double _Complex sqrt( double _Complex x ) { return csqrt( x ); }
|
---|
124 | long double _Complex sqrt( long double _Complex x ) { return csqrtl( x ); }
|
---|
125 |
|
---|
126 | float fabs( float x ) { return fabsf( x ); }
|
---|
127 | long double fabs( long double x ) { return fabsl( x ); }
|
---|
128 | float cabs( float _Complex x ) { return cabsf( x ); }
|
---|
129 | long double cabs( long double _Complex x ) { return cabsl( x ); }
|
---|
130 |
|
---|
131 | float floor( float x ) { return floorf( x ); }
|
---|
132 | long double floor( long double x ) { return floorl( x ); }
|
---|
133 |
|
---|
134 | float ceil( float x ) { return ceilf( x ); }
|
---|
135 | long double ceil( long double x ) { return ceill( x ); }
|
---|
136 |
|
---|
137 | //---------------------------------------
|
---|
138 |
|
---|
139 | float cbrt( float x ) { return cbrtf( x ); }
|
---|
140 | long double cbrt( long double x ) { return cbrtl( x ); }
|
---|
141 |
|
---|
142 | float copysign( float x, float y ) { return copysignf( x, y ); }
|
---|
143 | long double copysign( long double x, long double y ) { return copysignl( x, y ); }
|
---|
144 |
|
---|
145 | float erf( float x ) { return erff( x ); }
|
---|
146 | long double erf( long double x ) { return erfl( x ); }
|
---|
147 |
|
---|
148 | float erfc( float x ) { return erfcf( x ); }
|
---|
149 | long double erfc( long double x ) { return erfcl( x ); }
|
---|
150 |
|
---|
151 | float exp2( float x ) { return exp2f( x ); }
|
---|
152 | long double exp2( long double x ) { return exp2l( x ); }
|
---|
153 |
|
---|
154 | float expm1( float x ) { return expm1f( x ); }
|
---|
155 | long double expm1( long double x ) { return expm1l( x ); }
|
---|
156 |
|
---|
157 | float fdim( float x, float y ) { return fdimf( x, y ); }
|
---|
158 | long double fdim( long double x, long double y ) { return fdiml( x, y ); }
|
---|
159 |
|
---|
160 | float fma( float x, float y, float z ) { return fmaf( x, y, z ); }
|
---|
161 | long double fma( long double x, long double y, long double z ) { return fmal( x, y, z ); }
|
---|
162 |
|
---|
163 | float ?%?( float x, float y ) { return fmodf( x, y ); }
|
---|
164 | float fmod( float x, float y ) { return fmodf( x, y ); }
|
---|
165 | double ?%?( double x, double y ) { return fmod( x, y ); }
|
---|
166 | long double ?%?( long double x, long double y ) { return fmodl( x, y ); }
|
---|
167 | long double fmod( long double x, long double y ) { return fmodl( x, y ); }
|
---|
168 |
|
---|
169 | //[ int, float ] div( float x ) { int i; x = frexpf( x, &i ); return [ i, x ]; }
|
---|
170 | float frexp( float x, int *ip ) { return frexpf( x, ip ); }
|
---|
171 | //[ int, double ] div( double x ) { int i; x = frexp( x, &i ); return [ i, x ]; }
|
---|
172 | //[ int, long double ] div( long double x ) { int i; x = frexpl( x, &i ); return [ i, x ]; }
|
---|
173 | long double frexp( long double x, int *ip ) { return frexpl( x, ip ); }
|
---|
174 |
|
---|
175 | float hypot( float x, float y ) { return hypotf( x, y ); }
|
---|
176 | long double hypot( long double x, long double y ) { return hypotl( x, y ); }
|
---|
177 |
|
---|
178 | int ilogb( float x ) { return ilogbf( x ); }
|
---|
179 | int ilogb( long double x ) { return ilogbl( x ); }
|
---|
180 |
|
---|
181 | float ldexp( float x, int exp2 ) { return ldexpf( x, exp2 ); }
|
---|
182 | long double ldexp( long double x, int exp2 ) { return ldexpl( x, exp2 ); }
|
---|
183 |
|
---|
184 | float lgamma( float x ) { return lgammaf( x ); }
|
---|
185 | long double lgamma( long double x ) { return lgammal( x ); }
|
---|
186 |
|
---|
187 | float rint( float x ) { return rintf( x ); }
|
---|
188 | long double rint( long double x ) { return rintl( x ); }
|
---|
189 | long int rint( float x ) { return lrintf( x ); }
|
---|
190 | long int rint( double x ) { return lrint( x ); }
|
---|
191 | long int rint( long double x ) { return lrintl( x ); }
|
---|
192 | long long int rint( float x ) { return llrintf( x ); }
|
---|
193 | long long int rint( double x ) { return llrint( x ); }
|
---|
194 | long long int rint( long double x ) { return llrintl( x ); }
|
---|
195 |
|
---|
196 | long int lrint( float x ) { return lrintf( x ); }
|
---|
197 | long int lrint( long double x ) { return lrintl( x ); }
|
---|
198 | long long int llrint( float x ) { return llrintf( x ); }
|
---|
199 | long long int llrint( long double x ) { return llrintl( x ); }
|
---|
200 |
|
---|
201 | float round( float x ) { return roundf( x ); }
|
---|
202 | long double round( long double x ) { return roundl( x ); }
|
---|
203 | long int round( float x ) { return lroundf( x ); }
|
---|
204 | long int round( double x ) { return lround( x ); }
|
---|
205 | long int round( long double x ) { return lroundl( x ); }
|
---|
206 | long long int round( float x ) { return llroundf( x ); }
|
---|
207 | long long int round( double x ) { return llround( x ); }
|
---|
208 | long long int round( long double x ) { return llroundl( x ); }
|
---|
209 |
|
---|
210 | long int lround( float x ) { return lroundf( x ); }
|
---|
211 | long int lround( long double x ) { return lroundl( x ); }
|
---|
212 | long long int llround( float x ) { return llroundf( x ); }
|
---|
213 | long long int llround( long double x ) { return llroundl( x ); }
|
---|
214 |
|
---|
215 | float log10( float x ) { return log10f( x ); }
|
---|
216 | long double log10( long double x ) { return log10l( x ); }
|
---|
217 |
|
---|
218 | float log1p( float x ) { return log1pf( x ); }
|
---|
219 | long double log1p( long double x ) { return log1pl( x ); }
|
---|
220 |
|
---|
221 | float log2( float x ) { return log2f( x ); }
|
---|
222 | long double log2( long double x ) { return log2l( x ); }
|
---|
223 |
|
---|
224 | float logb( float x ) { return logbf( x ); }
|
---|
225 | long double logb( long double x ) { return logbl( x ); }
|
---|
226 |
|
---|
227 | float nearbyint( float x ) { return nearbyintf( x ); }
|
---|
228 | long double nearbyint( long double x ) { return nearbyintl( x ); }
|
---|
229 |
|
---|
230 | float nextafter( float x, float y ) { return nextafterf( x, y ); }
|
---|
231 | long double nextafter( long double x, long double y ) { return nextafterl( x, y ); }
|
---|
232 |
|
---|
233 | float nexttoward( float x, long double y ) { return nexttowardf( x, y ); }
|
---|
234 | long double nexttoward( long double x, long double y ) { return nexttowardl( x, y ); }
|
---|
235 |
|
---|
236 | float remainder( float x, float y ) { return remainderf( x, y ); }
|
---|
237 | long double remainder( long double x, long double y ) { return remainderl( x, y ); }
|
---|
238 |
|
---|
239 | //[ int, float ] remquo( float x, float y ) { int quo; x = remquof( x, y, &quo ); return [ quo, x ]; }
|
---|
240 | float remquo( float x, float y, int *quo ) { return remquof( x, y, quo ); }
|
---|
241 | //[ int, double ] remquo( double x, double y ) { int quo; x = remquo( x, y, &quo ); return [ quo, x ]; }
|
---|
242 | //[ int, long double ] remquo( long double x, long double y ) { int quo; x = remquol( x, y, &quo ); return [ quo, x ]; }
|
---|
243 | long double remquo( long double x, long double y, int *quo ) { return remquol( x, y, quo ); }
|
---|
244 |
|
---|
245 | float scalbn( float x, int exp ) { return scalbnf( x, exp ); }
|
---|
246 | long double scalbn( long double x, int exp ) { return scalbnl( x, exp ); }
|
---|
247 | float scalbn( float x, long int exp ) { return scalblnf( x, exp ); }
|
---|
248 | double scalbn( double x, long int exp ) { return scalbln( x, exp ); }
|
---|
249 | long double scalbn( long double x, long int exp ) { return scalblnl( x, exp ); }
|
---|
250 |
|
---|
251 | float scalbln( float x, long int exp ) { return scalblnf( x, exp ); }
|
---|
252 | long double scalbln( long double x, long int exp ) { return scalblnl( x, exp ); }
|
---|
253 |
|
---|
254 | float tgamma( float x ) { return tgammaf( x ); }
|
---|
255 | long double tgamma( long double x ) { return tgammal( x ); }
|
---|
256 |
|
---|
257 | float trunc( float x ) { return truncf( x ); }
|
---|
258 | long double trunc( long double x ) { return truncl( x ); }
|
---|
259 |
|
---|
260 | // Local Variables: //
|
---|
261 | // mode: c //
|
---|
262 | // tab-width: 4 //
|
---|
263 | // End: //
|
---|