source: src/libcfa/math @ 6d54c3a

ADTaaron-thesisarm-ehast-experimentalcleanup-dtorsdeferred_resndemanglerenumforall-pointer-decayjacob/cs343-translationjenkins-sandboxnew-astnew-ast-unique-exprnew-envno_listpersistent-indexerpthread-emulationqualifiedEnumresolv-newwith_gc
Last change on this file since 6d54c3a was dab7ac7, checked in by Peter A. Buhr <pabuhr@…>, 7 years ago

inline math routines, remove math.c, add exponential operator to math tests
  • Property mode set to 100644
File size: 21.0 KB
Line 
1//
2// Cforall Version 1.0.0 Copyright (C) 2016 University of Waterloo
3//
4// The contents of this file are covered under the licence agreement in the
5// file "LICENCE" distributed with Cforall.
6//
7// math --
8//
9// Author           : Peter A. Buhr
10// Created On       : Mon Apr 18 23:37:04 2016
11// Last Modified By : Peter A. Buhr
12// Last Modified On : Thu Jul 20 21:45:07 2017
13// Update Count     : 77
14//
15
16#pragma once
17
18#include <math.h>
19#include <complex.h>
20
21static inline float ?%?( float x, float y ) { return fmodf( x, y ); }
22static inline float fmod( float x, float y ) { return fmodf( x, y ); }
23static inline double ?%?( double x, double y ) { return fmod( x, y ); }
24// extern "C" { double fmod( double, double ); }
25static inline long double ?%?( long double x, long double y ) { return fmodl( x, y ); }
26static inline long double fmod( long double x, long double y ) { return fmodl( x, y ); }
27
28static inline float remainder( float x, float y ) { return remainderf( x, y ); }
29// extern "C" { double remainder( double, double ); }
30static inline long double remainder( long double x, long double y ) { return remainderl( x, y ); }
31
32static inline float remquo( float x, float y, int * quo ) { return remquof( x, y, quo ); }
33// extern "C" { double remquo( double x, double y, int * quo ); }
34static inline long double remquo( long double x, long double y, int * quo ) { return remquol( x, y, quo ); }
35static inline [ int, float ] remquo( float x, float y ) { int quo; x = remquof( x, y, &quo ); return [ quo, x ]; }
36static inline [ int, double ] remquo( double x, double y ) { int quo; x = remquo( x, y, &quo ); return [ quo, x ]; }
37static inline [ int, long double ] remquo( long double x, long double y ) { int quo; x = remquol( x, y, &quo ); return [ quo, x ]; }
38
39// alternative name for remquo
40static inline float div( float x, float y, int * quo ) { return remquof( x, y, quo ); }
41static inline double div( double x, double y, int * quo ) { return remquo( x, y, quo ); }
42static inline long double div( long double x, long double y, int * quo ) { return remquol( x, y, quo ); }
43static inline [ int, float ] div( float x, float y ) { int quo; x = remquof( x, y, &quo ); return [ quo, x ]; }
44static inline [ int, double ] div( double x, double y ) { int quo; x = remquo( x, y, &quo ); return [ quo, x ]; }
45static inline [ int, long double ] div( long double x, long double y ) { int quo; x = remquol( x, y, &quo ); return [ quo, x ]; }
46
47static inline float fma( float x, float y, float z ) { return fmaf( x, y, z ); }
48// extern "C" { double fma( double, double, double ); }
49static inline long double fma( long double x, long double y, long double z ) { return fmal( x, y, z ); }
50
51static inline float fdim( float x, float y ) { return fdimf( x, y ); }
52// extern "C" { double fdim( double, double ); }
53static inline long double fdim( long double x, long double y ) { return fdiml( x, y ); }
54
55static inline float nan( const char * tag ) { return nanf( tag ); }
56// extern "C" { double nan( const char * ); }
57static inline long double nan( const char * tag ) { return nanl( tag ); }
58
59//---------------------- Exponential ----------------------
60
61static inline float exp( float x ) { return expf( x ); }
62// extern "C" { double exp( double ); }
63static inline long double exp( long double x ) { return expl( x ); }
64static inline float _Complex exp( float _Complex x ) { return cexpf( x ); }
65static inline double _Complex exp( double _Complex x ) { return cexp( x ); }
66static inline long double _Complex exp( long double _Complex x ) { return cexpl( x ); }
67
68static inline float exp2( float x ) { return exp2f( x ); }
69// extern "C" { double exp2( double ); }
70static inline long double exp2( long double x ) { return exp2l( x ); }
71//static inline float _Complex exp2( float _Complex x ) { return cexp2f( x ); }
72//static inline double _Complex exp2( double _Complex x ) { return cexp2( x ); }
73//static inline long double _Complex exp2( long double _Complex x ) { return cexp2l( x ); }
74
75static inline float expm1( float x ) { return expm1f( x ); }
76// extern "C" { double expm1( double ); }
77static inline long double expm1( long double x ) { return expm1l( x ); }
78
79static inline float pow( float x, float y ) { return powf( x, y ); }
80// extern "C" { double pow( double, double ); }
81static inline long double pow( long double x, long double y ) { return powl( x, y ); }
82static inline float _Complex pow( float _Complex x, float _Complex y ) { return cpowf( x, y ); }
83static inline double _Complex pow( double _Complex x, double _Complex y ) { return cpow( x, y ); }
84static inline long double _Complex pow( long double _Complex x, long double _Complex y ) { return cpowl( x, y ); }
85
86static inline float ?\?( float x, float y ) { return powf( x, y ); }
87static inline double ?\?( double x, double y ) { return pow( x, y ); }
88static inline long double ?\?( long double x, long double y ) { return powl( x, y ); }
89static inline float _Complex ?\?( float _Complex x, _Complex float y ) { return cpowf(x, y ); }
90static inline double _Complex ?\?( double _Complex x, _Complex double y ) { return cpow( x, y ); }
91static inline long double _Complex ?\?( long double _Complex x, _Complex long double y ) { return cpowl( x, y ); }
92
93static inline float ?\=?( float * x, float y ) { *x = *x \ y; return *x; }
94static inline double ?\=?( double * x, double y ) { *x = *x \ y; return *x; }
95static inline long double ?\=?( long double * x, long double y ) { *x = *x \ y; return *x; }
96static inline float _Complex ?\=?( float _Complex * x, _Complex float y ) { *x = *x \ y; return *x; }
97static inline double _Complex ?\=?( double _Complex * x, _Complex double y ) { *x = *x \ y; return *x; }
98static inline long double _Complex ?\=?( long double _Complex * x, _Complex long double y ) { *x = *x \ y; return *x; }
99
100static inline long int ?\?( long int x, unsigned long y ) {     // disallow negative exponent
101    if ( y == 0 ) return 1;
102    if ( x == 2 ) return x << (y - 1);
103    long int prod = 1;
104    for ( unsigned int i = 0; i < y; i += 1 ) {
105                prod = prod * x;
106    } // for
107    return prod;
108}
109static inline double ?\?( long int x, signed long y ) { // allow negative exponent
110    if ( y >=  0 ) return (double)(x \ (unsigned int)y);
111    else return 1.0 / x \ (unsigned int)(-y);
112}
113static inline forall( otype T | { void ?{}( T * this, one_t ); T ?*?( T, T ); } )
114T ?\?( T x, unsigned long y ) {
115    T prod = 1;
116    for ( unsigned int i = 1; i < y; i += 1 ) {
117                prod = prod * x;
118    } // for
119    return prod;
120}
121static inline long int ?\=?( long int * x, unsigned long y ) { *x = *x \ y; return *x; }
122static inline int ?\=?( int * x, unsigned long y ) { *x = *x \ y; return *x; }
123
124//---------------------- Logarithm ----------------------
125
126static inline float log( float x ) { return logf( x ); }
127// extern "C" { double log( double ); }
128static inline long double log( long double x ) { return logl( x ); }
129static inline float _Complex log( float _Complex x ) { return clogf( x ); }
130static inline double _Complex log( double _Complex x ) { return clog( x ); }
131static inline long double _Complex log( long double _Complex x ) { return clogl( x ); }
132
133static inline float log2( float x ) { return log2f( x ); }
134// extern "C" { double log2( double ); }
135static inline long double log2( long double x ) { return log2l( x ); }
136// static inline float _Complex log2( float _Complex x ) { return clog2f( x ); }
137// static inline double _Complex log2( double _Complex x ) { return clog2( x ); }
138// static inline long double _Complex log2( long double _Complex x ) { return clog2l( x ); }
139
140static inline float log10( float x ) { return log10f( x ); }
141// extern "C" { double log10( double ); }
142static inline long double log10( long double x ) { return log10l( x ); }
143// static inline float _Complex log10( float _Complex x ) { return clog10f( x ); }
144// static inline double _Complex log10( double _Complex x ) { return clog10( x ); }
145// static inline long double _Complex log10( long double _Complex x ) { return clog10l( x ); }
146
147static inline float log1p( float x ) { return log1pf( x ); }
148// extern "C" { double log1p( double ); }
149static inline long double log1p( long double x ) { return log1pl( x ); }
150
151static inline int ilogb( float x ) { return ilogbf( x ); }
152// extern "C" { int ilogb( double ); }
153static inline int ilogb( long double x ) { return ilogbl( x ); }
154
155static inline float logb( float x ) { return logbf( x ); }
156// extern "C" { double logb( double ); }
157static inline long double logb( long double x ) { return logbl( x ); }
158
159static inline float sqrt( float x ) { return sqrtf( x ); }
160// extern "C" { double sqrt( double ); }
161static inline long double sqrt( long double x ) { return sqrtl( x ); }
162static inline float _Complex sqrt( float _Complex x ) { return csqrtf( x ); }
163static inline double _Complex sqrt( double _Complex x ) { return csqrt( x ); }
164static inline long double _Complex sqrt( long double _Complex x ) { return csqrtl( x ); }
165
166static inline float cbrt( float x ) { return cbrtf( x ); }
167// extern "C" { double cbrt( double ); }
168static inline long double cbrt( long double x ) { return cbrtl( x ); }
169
170static inline float hypot( float x, float y ) { return hypotf( x, y ); }
171// extern "C" { double hypot( double, double ); }
172static inline long double hypot( long double x, long double y ) { return hypotl( x, y ); }
173
174//---------------------- Trigonometric ----------------------
175
176static inline float sin( float x ) { return sinf( x ); }
177// extern "C" { double sin( double ); }
178static inline long double sin( long double x ) { return sinl( x ); }
179static inline float _Complex sin( float _Complex x ) { return csinf( x ); }
180static inline double _Complex sin( double _Complex x ) { return csin( x ); }
181static inline long double _Complex sin( long double _Complex x ) { return csinl( x ); }
182
183static inline float cos( float x ) { return cosf( x ); }
184// extern "C" { double cos( double ); }
185static inline long double cos( long double x ) { return cosl( x ); }
186static inline float _Complex cos( float _Complex x ) { return ccosf( x ); }
187static inline double _Complex cos( double _Complex x ) { return ccos( x ); }
188static inline long double _Complex cos( long double _Complex x ) { return ccosl( x ); }
189
190static inline float tan( float x ) { return tanf( x ); }
191// extern "C" { double tan( double ); }
192static inline long double tan( long double x ) { return tanl( x ); }
193static inline float _Complex tan( float _Complex x ) { return ctanf( x ); }
194static inline double _Complex tan( double _Complex x ) { return ctan( x ); }
195static inline long double _Complex tan( long double _Complex x ) { return ctanl( x ); }
196
197static inline float asin( float x ) { return asinf( x ); }
198// extern "C" { double asin( double ); }
199static inline long double asin( long double x ) { return asinl( x ); }
200static inline float _Complex asin( float _Complex x ) { return casinf( x ); }
201static inline double _Complex asin( double _Complex x ) { return casin( x ); }
202static inline long double _Complex asin( long double _Complex x ) { return casinl( x ); }
203
204static inline float acos( float x ) { return acosf( x ); }
205// extern "C" { double acos( double ); }
206static inline long double acos( long double x ) { return acosl( x ); }
207static inline float _Complex acos( float _Complex x ) { return cacosf( x ); }
208static inline double _Complex acos( double _Complex x ) { return cacos( x ); }
209static inline long double _Complex acos( long double _Complex x ) { return cacosl( x ); }
210
211static inline float atan( float x ) { return atanf( x ); }
212// extern "C" { double atan( double ); }
213static inline long double atan( long double x ) { return atanl( x ); }
214static inline float _Complex atan( float _Complex x ) { return catanf( x ); }
215static inline double _Complex atan( double _Complex x ) { return catan( x ); }
216static inline long double _Complex atan( long double _Complex x ) { return catanl( x ); }
217
218static inline float atan2( float x, float y ) { return atan2f( x, y ); }
219// extern "C" { double atan2( double, double ); }
220static inline long double atan2( long double x, long double y ) { return atan2l( x, y ); }
221
222// alternative name for atan2
223static inline float atan( float x, float y ) { return atan2f( x, y ); }
224static inline double atan( double x, double y ) { return atan2( x, y ); }
225static inline long double atan( long double x, long double y ) { return atan2l( x, y ); }
226
227//---------------------- Hyperbolic ----------------------
228
229static inline float sinh( float x ) { return sinhf( x ); }
230// extern "C" { double sinh( double ); }
231static inline long double sinh( long double x ) { return sinhl( x ); }
232static inline float _Complex sinh( float _Complex x ) { return csinhf( x ); }
233static inline double _Complex sinh( double _Complex x ) { return csinh( x ); }
234static inline long double _Complex sinh( long double _Complex x ) { return csinhl( x ); }
235
236static inline float cosh( float x ) { return coshf( x ); }
237// extern "C" { double cosh( double ); }
238static inline long double cosh( long double x ) { return coshl( x ); }
239static inline float _Complex cosh( float _Complex x ) { return ccoshf( x ); }
240static inline double _Complex cosh( double _Complex x ) { return ccosh( x ); }
241static inline long double _Complex cosh( long double _Complex x ) { return ccoshl( x ); }
242
243static inline float tanh( float x ) { return tanhf( x ); }
244// extern "C" { double tanh( double ); }
245static inline long double tanh( long double x ) { return tanhl( x ); }
246static inline float _Complex tanh( float _Complex x ) { return ctanhf( x ); }
247static inline double _Complex tanh( double _Complex x ) { return ctanh( x ); }
248static inline long double _Complex tanh( long double _Complex x ) { return ctanhl( x ); }
249
250static inline float asinh( float x ) { return asinhf( x ); }
251// extern "C" { double asinh( double ); }
252static inline long double asinh( long double x ) { return asinhl( x ); }
253static inline float _Complex asinh( float _Complex x ) { return casinhf( x ); }
254static inline double _Complex asinh( double _Complex x ) { return casinh( x ); }
255static inline long double _Complex asinh( long double _Complex x ) { return casinhl( x ); }
256
257static inline float acosh( float x ) { return acoshf( x ); }
258// extern "C" { double acosh( double ); }
259static inline long double acosh( long double x ) { return acoshl( x ); }
260static inline float _Complex acosh( float _Complex x ) { return cacoshf( x ); }
261static inline double _Complex acosh( double _Complex x ) { return cacosh( x ); }
262static inline long double _Complex acosh( long double _Complex x ) { return cacoshl( x ); }
263
264static inline float atanh( float x ) { return atanhf( x ); }
265// extern "C" { double atanh( double ); }
266static inline long double atanh( long double x ) { return atanhl( x ); }
267static inline float _Complex atanh( float _Complex x ) { return catanhf( x ); }
268static inline double _Complex atanh( double _Complex x ) { return catanh( x ); }
269static inline long double _Complex atanh( long double _Complex x ) { return catanhl( x ); }
270
271//---------------------- Error / Gamma ----------------------
272
273static inline float erf( float x ) { return erff( x ); }
274// extern "C" { double erf( double ); }
275static inline long double erf( long double x ) { return erfl( x ); }
276// float _Complex erf( float _Complex );
277// double _Complex erf( double _Complex );
278// long double _Complex erf( long double _Complex );
279
280static inline float erfc( float x ) { return erfcf( x ); }
281// extern "C" { double erfc( double ); }
282static inline long double erfc( long double x ) { return erfcl( x ); }
283// float _Complex erfc( float _Complex );
284// double _Complex erfc( double _Complex );
285// long double _Complex erfc( long double _Complex );
286
287static inline float lgamma( float x ) { return lgammaf( x ); }
288// extern "C" { double lgamma( double ); }
289static inline long double lgamma( long double x ) { return lgammal( x ); }
290static inline float lgamma( float x, int * sign ) { return lgammaf_r( x, sign ); }
291static inline double lgamma( double x, int * sign ) { return lgamma_r( x, sign ); }
292static inline long double lgamma( long double x, int * sign ) { return lgammal_r( x, sign ); }
293
294static inline float tgamma( float x ) { return tgammaf( x ); }
295// extern "C" { double tgamma( double ); }
296static inline long double tgamma( long double x ) { return tgammal( x ); }
297
298//---------------------- Nearest Integer ----------------------
299
300static inline float floor( float x ) { return floorf( x ); }
301// extern "C" { double floor( double ); }
302static inline long double floor( long double x ) { return floorl( x ); }
303
304static inline float ceil( float x ) { return ceilf( x ); }
305// extern "C" { double ceil( double ); }
306static inline long double ceil( long double x ) { return ceill( x ); }
307
308static inline float trunc( float x ) { return truncf( x ); }
309// extern "C" { double trunc( double ); }
310static inline long double trunc( long double x ) { return truncl( x ); }
311
312static inline float rint( float x ) { return rintf( x ); }
313// extern "C" { double rint( double x ); }
314static inline long double rint( long double x ) { return rintl( x ); }
315static inline long int rint( float x ) { return lrintf( x ); }
316static inline long int rint( double x ) { return lrint( x ); }
317static inline long int rint( long double x ) { return lrintl( x ); }
318static inline long long int rint( float x ) { return llrintf( x ); }
319static inline long long int rint( double x ) { return llrint( x ); }
320static inline long long int rint( long double x ) { return llrintl( x ); }
321
322static inline long int lrint( float x ) { return lrintf( x ); }
323// extern "C" { long int lrint( double ); }
324static inline long int lrint( long double x ) { return lrintl( x ); }
325static inline long long int llrint( float x ) { return llrintf( x ); }
326// extern "C" { long long int llrint( double ); }
327static inline long long int llrint( long double x ) { return llrintl( x ); }
328
329static inline float nearbyint( float x ) { return nearbyintf( x ); }
330// extern "C" { double nearbyint( double ); }
331static inline long double nearbyint( long double x ) { return nearbyintl( x ); }
332
333static inline float round( float x ) { return roundf( x ); }
334// extern "C" { double round( double x ); }
335static inline long double round( long double x ) { return roundl( x ); }
336static inline long int round( float x ) { return lroundf( x ); }
337static inline long int round( double x ) { return lround( x ); }
338static inline long int round( long double x ) { return lroundl( x ); }
339static inline long long int round( float x ) { return llroundf( x ); }
340static inline long long int round( double x ) { return llround( x ); }
341static inline long long int round( long double x ) { return llroundl( x ); }
342
343static inline long int lround( float x ) { return lroundf( x ); }
344// extern "C" { long int lround( double ); }
345static inline long int lround( long double x ) { return lroundl( x ); }
346static inline long long int llround( float x ) { return llroundf( x ); }
347// extern "C" { long long int llround( double ); }
348static inline long long int llround( long double x ) { return llroundl( x ); }
349
350//---------------------- Manipulation ----------------------
351
352static inline float copysign( float x, float y ) { return copysignf( x, y ); }
353// extern "C" { double copysign( double, double ); }
354static inline long double copysign( long double x, long double y ) { return copysignl( x, y ); }
355
356static inline float frexp( float x, int * ip ) { return frexpf( x, ip ); }
357// extern "C" { double frexp( double, int * ); }
358static inline long double frexp( long double x, int * ip ) { return frexpl( x, ip ); }
359
360static inline float ldexp( float x, int exp2 ) { return ldexpf( x, exp2 ); }
361// extern "C" { double ldexp( double, int ); }
362static inline long double ldexp( long double x, int exp2 ) { return ldexpl( x, exp2 ); }
363
364static inline [ float, float ] modf( float x ) { float i; x = modff( x, &i ); return [ i, x ]; }
365static inline float modf( float x, float * i ) { return modff( x, i ); }
366static inline [ double, double ] modf( double x ) { double i; x = modf( x, &i ); return [ i, x ]; }
367// extern "C" { double modf( double, double * ); }
368static inline [ long double, long double ] modf( long double x ) { long double i; x = modfl( x, &i ); return [ i, x ]; }
369static inline long double modf( long double x, long double * i ) { return modfl( x, i ); }
370
371static inline float nextafter( float x, float y ) { return nextafterf( x, y ); }
372// extern "C" { double nextafter( double, double ); }
373static inline long double nextafter( long double x, long double y ) { return nextafterl( x, y ); }
374
375static inline float nexttoward( float x, long double y ) { return nexttowardf( x, y ); }
376// extern "C" { double nexttoward( double, long double ); }
377static inline long double nexttoward( long double x, long double y ) { return nexttowardl( x, y ); }
378
379static inline float scalbn( float x, int exp ) { return scalbnf( x, exp ); }
380// extern "C" { double scalbn( double, int ); }
381static inline long double scalbn( long double x, int exp ) { return scalbnl( x, exp ); }
382static inline float scalbn( float x, long int exp ) { return scalblnf( x, exp ); }
383static inline double scalbn( double x, long int exp ) { return scalbln( x, exp ); }
384static inline long double scalbn( long double x, long int exp ) { return scalblnl( x, exp ); }
385
386static inline float scalbln( float x, long int exp ) { return scalblnf( x, exp ); }
387// extern "C" { double scalbln( double, long int ); }
388static inline long double scalbln( long double x, long int exp ) { return scalblnl( x, exp ); }
389
390// Local Variables: //
391// mode: c //
392// tab-width: 4 //
393// End: //
Note: See TracBrowser for help on using the repository browser.