source: src/prelude/builtins.c @ 9236060

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 9236060 was 9236060, checked in by Rob Schluntz <rschlunt@…>, 7 years ago

Merge branch 'master' into references

  • Property mode set to 100644
File size: 3.7 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// builtins.c --
8//
9// Author           : Peter A. Buhr
10// Created On       : Fri Jul 21 16:21:03 2017
11// Last Modified By : Andrew Beach
12// Last Modified On : Tus Jul 25 15:33:00 2017
13// Update Count     : 14
14//
15
16// exception implementation
17
18typedef unsigned long long __cfaabi_exception_type_t;
19
20#include "../libcfa/virtual.h"
21#include "../libcfa/exception.h"
22
23// exponentiation operator implementation
24
25extern "C" {
26        float powf( float x, float y );
27        double pow( double x, double y );
28        long double powl( long double x, long double y );
29        float _Complex cpowf( float _Complex x, _Complex float z );
30        double _Complex cpow( double _Complex x, _Complex double z );
31        long double _Complex cpowl( long double _Complex x, _Complex long double z );
32} // extern "C"
33
34static inline float ?\?( float x, float y ) { return powf( x, y ); }
35static inline double ?\?( double x, double y ) { return pow( x, y ); }
36static inline long double ?\?( long double x, long double y ) { return powl( x, y ); }
37static inline float _Complex ?\?( float _Complex x, _Complex float y ) { return cpowf(x, y ); }
38static inline double _Complex ?\?( double _Complex x, _Complex double y ) { return cpow( x, y ); }
39static inline long double _Complex ?\?( long double _Complex x, _Complex long double y ) { return cpowl( x, y ); }
40
41static inline long int ?\?( long int ep, unsigned long int y ) { // disallow negative exponent
42        if ( y == 0 ) return 1;                                                         // base case
43    if ( ep == 2 ) return ep << (y - 1);                                // special case, positive shifting only
44    typeof( ep ) op = 1;                                                                // accumulate odd product
45    for ( ; y > 1; y >>= 1 ) {                                                  // squaring exponentiation, O(log2 y)
46                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
47                ep *= ep;
48        } // for
49    return ep * op;
50} // ?\?
51
52// FIX ME, cannot resolve the "T op = 1".
53
54// static inline forall( otype T | { void ?{}( T * this, one_t ); T ?*?( T, T ); } )
55// T ?\?( T ep, unsigned long int y ) {
56//     if ( y == 0 ) return 1;
57//     T op = 1;
58//     for ( ; y > 1; y >>= 1 ) {                                                       // squaring exponentiation, O(log2 y)
59//              if ( (y & 1) == 1 ) op = op * ep;                               // odd ?
60//              ep = ep * ep;
61//     } // for
62//     return ep * op;
63// } // ?\?
64
65// unsigned computation may be faster and larger
66static inline unsigned long int ?\?( unsigned long int ep, unsigned long int y ) { // disallow negative exponent
67        if ( y == 0 ) return 1;                                                         // base case
68    if ( ep == 2 ) return ep << (y - 1);                                // special case, positive shifting only
69    typeof( ep ) op = 1;                                                                // accumulate odd product
70    for ( ; y > 1; y >>= 1 ) {                                                  // squaring exponentiation, O(log2 y)
71                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
72                ep *= ep;
73        } // for
74    return ep * op;
75} // ?\?
76
77static inline double ?\?( long int x, signed long int y ) {     // allow negative exponent
78    if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
79    else return 1.0 / x \ (unsigned int)(-y);
80} // ?\?
81
82static inline forall( otype T | { void ?{}( T & this, one_t ); T ?*?( T, T ); double ?/?( double, T ); } )
83double ?\?( T x, signed long int y ) {
84    if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
85    else return 1.0 / x \ (unsigned long int)(-y);
86} // ?\?
87
88static inline long int ?\=?( long int * x, unsigned long int y ) { *x = *x \ y; return *x; }
89static inline unsigned long int ?\=?( unsigned long int * x, unsigned long int y ) { *x = *x \ y; return *x; }
90static inline int ?\=?( int * x, unsigned long int y ) { *x = *x \ y; return *x; }
91static inline unsigned int ?\=?( unsigned int * x, unsigned long int y ) { *x = *x \ y; return *x; }
92
93// Local Variables: //
94// mode: c //
95// tab-width: 4 //
96// End: //
Note: See TracBrowser for help on using the repository browser.