| 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 : Peter A. Buhr
 | 
|---|
| 12 | // Last Modified On : Sat Aug 14 08:45:54 2021
 | 
|---|
| 13 | // Update Count     : 133
 | 
|---|
| 14 | //
 | 
|---|
| 15 | 
 | 
|---|
| 16 | #define __cforall_builtins__
 | 
|---|
| 17 | 
 | 
|---|
| 18 | // type that wraps a pointer and a destructor-like function - used in generating implicit destructor calls for struct members in user-defined functions
 | 
|---|
| 19 | // Note: needs to occur early, because it is used to generate destructor calls during code generation
 | 
|---|
| 20 | forall(T &)
 | 
|---|
| 21 | struct __Destructor {
 | 
|---|
| 22 |         T * object;
 | 
|---|
| 23 |         void (*dtor)(T *);
 | 
|---|
| 24 | };
 | 
|---|
| 25 | 
 | 
|---|
| 26 | // defined destructor in the case that non-generated code wants to use __Destructor
 | 
|---|
| 27 | forall(T &)
 | 
|---|
| 28 | static inline void ^?{}(__Destructor(T) & x) {
 | 
|---|
| 29 |         if (x.object && x.dtor) {
 | 
|---|
| 30 |                 x.dtor(x.object);
 | 
|---|
| 31 |         }
 | 
|---|
| 32 | }
 | 
|---|
| 33 | 
 | 
|---|
| 34 | // easy interface into __Destructor's destructor for easy codegen purposes
 | 
|---|
| 35 | extern "C" {
 | 
|---|
| 36 |         forall(T &)
 | 
|---|
| 37 |         static inline void __destroy_Destructor(__Destructor(T) * dtor) {
 | 
|---|
| 38 |                 ^(*dtor){};
 | 
|---|
| 39 |         }
 | 
|---|
| 40 | }
 | 
|---|
| 41 | 
 | 
|---|
| 42 | // exception implementation
 | 
|---|
| 43 | 
 | 
|---|
| 44 | typedef unsigned long long __cfaabi_abi_exception_type_t;
 | 
|---|
| 45 | 
 | 
|---|
| 46 | #include <limits.h>                                                                             // CHAR_BIT
 | 
|---|
| 47 | #include "../src/virtual.h"
 | 
|---|
| 48 | #include "../src/exception.h"
 | 
|---|
| 49 | 
 | 
|---|
| 50 | void exit( int status, const char fmt[], ... ) __attribute__ (( format(printf, 2, 3), __nothrow__, __leaf__, __noreturn__ ));
 | 
|---|
| 51 | void abort( const char fmt[], ... ) __attribute__ (( format(printf, 1, 2), __nothrow__, __leaf__, __noreturn__ ));
 | 
|---|
| 52 | 
 | 
|---|
| 53 | forall(T &)
 | 
|---|
| 54 | static inline T & identity(T & i) {
 | 
|---|
| 55 |         return i;
 | 
|---|
| 56 | }
 | 
|---|
| 57 | 
 | 
|---|
| 58 | // generator support
 | 
|---|
| 59 | struct generator$ {
 | 
|---|
| 60 |         inline int;
 | 
|---|
| 61 | };
 | 
|---|
| 62 | 
 | 
|---|
| 63 | static inline void  ?{}(generator$ & this) { ((int&)this) = 0; }
 | 
|---|
| 64 | static inline void ^?{}(generator$ &) {}
 | 
|---|
| 65 | 
 | 
|---|
| 66 | trait is_generator(T &) {
 | 
|---|
| 67 |       void main(T & this);
 | 
|---|
| 68 |       generator$ * get_generator(T & this);
 | 
|---|
| 69 | };
 | 
|---|
| 70 | 
 | 
|---|
| 71 | forall(T & | is_generator(T))
 | 
|---|
| 72 | static inline T & resume(T & gen) {
 | 
|---|
| 73 |         main(gen);
 | 
|---|
| 74 |         return gen;
 | 
|---|
| 75 | }
 | 
|---|
| 76 | 
 | 
|---|
| 77 | // implicit increment, decrement if += defined, and implicit not if != defined
 | 
|---|
| 78 | 
 | 
|---|
| 79 | // C11 reference manual Section 6.5.16 (page101): "An assignment expression has the value of the left operand after the
 | 
|---|
| 80 | // assignment, but is not an lvalue." Hence, return a value not a reference.
 | 
|---|
| 81 | static inline {
 | 
|---|
| 82 |         forall( T | { T ?+=?( T &, one_t ); } )
 | 
|---|
| 83 |         T ++?( T & x ) { return x += 1; }
 | 
|---|
| 84 | 
 | 
|---|
| 85 |         forall( T | { T ?+=?( T &, one_t ); } )
 | 
|---|
| 86 |         T ?++( T & x ) { T tmp = x; x += 1; return tmp; }
 | 
|---|
| 87 | 
 | 
|---|
| 88 |         forall( T | { T ?-=?( T &, one_t ); } )
 | 
|---|
| 89 |         T --?( T & x ) { return x -= 1; }
 | 
|---|
| 90 | 
 | 
|---|
| 91 |         forall( T | { T ?-=?( T &, one_t ); } )
 | 
|---|
| 92 |         T ?--( T & x ) { T tmp = x; x -= 1; return tmp; }
 | 
|---|
| 93 | 
 | 
|---|
| 94 |         forall( T | { int ?!=?( T, zero_t ); } )
 | 
|---|
| 95 |         int !?( T & x ) { return !( x != 0 ); }
 | 
|---|
| 96 | } // distribution
 | 
|---|
| 97 | 
 | 
|---|
| 98 | // universal typed pointer constant
 | 
|---|
| 99 | static inline forall( DT & ) DT * intptr( uintptr_t addr ) { return (DT *)addr; }
 | 
|---|
| 100 | static inline forall( ftype FT ) FT * intptr( uintptr_t addr ) { return (FT *)addr; }
 | 
|---|
| 101 | 
 | 
|---|
| 102 | #if defined(__SIZEOF_INT128__)
 | 
|---|
| 103 | // constructor for 128-bit numbers (all constants are unsigned as +/- are operators)
 | 
|---|
| 104 | static inline void ?{}( unsigned int128 & this, unsigned long int h, unsigned long int l ) {
 | 
|---|
| 105 |         this = (unsigned int128)h << 64 | (unsigned int128)l;
 | 
|---|
| 106 | } // ?{}
 | 
|---|
| 107 | #endif // __SIZEOF_INT128__
 | 
|---|
| 108 | 
 | 
|---|
| 109 | // for-control index constraints
 | 
|---|
| 110 | // forall( T | { void ?{}( T &, zero_t ); void ?{}( T &, one_t ); T ?+=?( T &, T ); T ?-=?( T &, T ); int ?<?( T, T ); } )
 | 
|---|
| 111 | // static inline T __for_control_index_constraints__( T t ) { return t; }
 | 
|---|
| 112 | 
 | 
|---|
| 113 | // exponentiation operator implementation
 | 
|---|
| 114 | 
 | 
|---|
| 115 | extern "C" {
 | 
|---|
| 116 |         float powf( float x, float y );
 | 
|---|
| 117 |         double pow( double x, double y );
 | 
|---|
| 118 |         long double powl( long double x, long double y );
 | 
|---|
| 119 |         float _Complex cpowf( float _Complex x, _Complex float z );
 | 
|---|
| 120 |         double _Complex cpow( double _Complex x, _Complex double z );
 | 
|---|
| 121 |         long double _Complex cpowl( long double _Complex x, _Complex long double z );
 | 
|---|
| 122 | } // extern "C"
 | 
|---|
| 123 | 
 | 
|---|
| 124 | static inline {
 | 
|---|
| 125 |         float ?\?( float x, float y ) { return powf( x, y ); }
 | 
|---|
| 126 |         double ?\?( double x, double y ) { return pow( x, y ); }
 | 
|---|
| 127 |         long double ?\?( long double x, long double y ) { return powl( x, y ); }
 | 
|---|
| 128 |         float _Complex ?\?( float _Complex x, _Complex float y ) { return cpowf(x, y ); }
 | 
|---|
| 129 |         double _Complex ?\?( double _Complex x, _Complex double y ) { return cpow( x, y ); }
 | 
|---|
| 130 |         long double _Complex ?\?( long double _Complex x, _Complex long double y ) { return cpowl( x, y ); }
 | 
|---|
| 131 | } // distribution
 | 
|---|
| 132 | 
 | 
|---|
| 133 | #define __CFA_BASE_COMP_1__() if ( x == 1 ) return 1
 | 
|---|
| 134 | #define __CFA_BASE_COMP_2__() if ( x == 2 ) return x << (y - 1)
 | 
|---|
| 135 | #define __CFA_EXP_OVERFLOW__() if ( y >= sizeof(y) * CHAR_BIT ) return 0
 | 
|---|
| 136 | 
 | 
|---|
| 137 | #define __CFA_EXP__() \
 | 
|---|
| 138 |         if ( y == 0 ) return 1;                                                         /* convention */ \
 | 
|---|
| 139 |         __CFA_BASE_COMP_1__();                                                          /* base case */ \
 | 
|---|
| 140 |         __CFA_BASE_COMP_2__();                                                          /* special case, positive shifting for integral types */ \
 | 
|---|
| 141 |         __CFA_EXP_OVERFLOW__();                                                         /* immediate overflow, negative exponent > 2^size-1 */ \
 | 
|---|
| 142 |         typeof(x) op = 1;                                                                       /* accumulate odd product */ \
 | 
|---|
| 143 |         for ( ; y > 1; y >>= 1 ) {                                                      /* squaring exponentiation, O(log2 y) */ \
 | 
|---|
| 144 |                 if ( (y & 1) == 1 ) op = op * x;                                /* odd ? */ \
 | 
|---|
| 145 |                 x = x * x; \
 | 
|---|
| 146 |         } \
 | 
|---|
| 147 |         return x * op
 | 
|---|
| 148 | 
 | 
|---|
| 149 | static inline {
 | 
|---|
| 150 |         long int ?\?( int x, unsigned int y ) { __CFA_EXP__(); }
 | 
|---|
| 151 |         long int ?\?( long int x, unsigned long int y ) { __CFA_EXP__(); }
 | 
|---|
| 152 |         long long int ?\?( long long int x, unsigned long long int y ) { __CFA_EXP__(); }
 | 
|---|
| 153 |         // unsigned computation may be faster and larger
 | 
|---|
| 154 |         unsigned long int ?\?( unsigned int x, unsigned int y ) { __CFA_EXP__(); }
 | 
|---|
| 155 |         unsigned long int ?\?( unsigned long int x, unsigned long int y ) { __CFA_EXP__(); }
 | 
|---|
| 156 |         unsigned long long int ?\?( unsigned long long int x, unsigned long long int y ) { __CFA_EXP__(); }
 | 
|---|
| 157 | } // distribution
 | 
|---|
| 158 | 
 | 
|---|
| 159 | #undef __CFA_BASE_COMP_1__
 | 
|---|
| 160 | #undef __CFA_BASE_COMP_2__
 | 
|---|
| 161 | #undef __CFA_EXP_OVERFLOW__
 | 
|---|
| 162 | #define __CFA_BASE_COMP_1__()
 | 
|---|
| 163 | #define __CFA_BASE_COMP_2__()
 | 
|---|
| 164 | #define __CFA_EXP_OVERFLOW__()
 | 
|---|
| 165 | 
 | 
|---|
| 166 | static inline forall( OT | { void ?{}( OT & this, one_t ); OT ?*?( OT, OT ); } ) {
 | 
|---|
| 167 |         OT ?\?( OT x, unsigned int y ) { __CFA_EXP__(); }
 | 
|---|
| 168 |         OT ?\?( OT x, unsigned long int y ) { __CFA_EXP__(); }
 | 
|---|
| 169 |         OT ?\?( OT x, unsigned long long int y ) { __CFA_EXP__(); }
 | 
|---|
| 170 | } // distribution
 | 
|---|
| 171 | 
 | 
|---|
| 172 | #undef __CFA_BASE_COMP_1__
 | 
|---|
| 173 | #undef __CFA_BASE_COMP_2__
 | 
|---|
| 174 | #undef __CFA_EXP_OVERFLOW__
 | 
|---|
| 175 | 
 | 
|---|
| 176 | static inline {
 | 
|---|
| 177 |         long int ?\=?( int & x, unsigned int y ) { x = x \ y; return x; }
 | 
|---|
| 178 |         long int ?\=?( long int & x, unsigned long int y ) { x = x \ y; return x; }
 | 
|---|
| 179 |         long long int ?\=?( long long int & x, unsigned long long int y ) { x = x \ y; return x; }
 | 
|---|
| 180 |         unsigned long int ?\=?( unsigned int & x, unsigned int y ) { x = x \ y; return x; }
 | 
|---|
| 181 |         unsigned long int ?\=?( unsigned long int & x, unsigned long int y ) { x = x \ y; return x; }
 | 
|---|
| 182 |         unsigned long long int ?\=?( unsigned long long int & x, unsigned long long int y ) { x = x \ y; return x; }
 | 
|---|
| 183 | } // distribution
 | 
|---|
| 184 | 
 | 
|---|
| 185 | // Local Variables: //
 | 
|---|
| 186 | // mode: c //
 | 
|---|
| 187 | // tab-width: 4 //
 | 
|---|
| 188 | // End: //
 | 
|---|