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