Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed

  • libcfa/prelude/builtins.c

    r8a30423 r7579ac0  
    1010// Created On       : Fri Jul 21 16:21:03 2017
    1111// Last Modified By : Peter A. Buhr
    12 // Last Modified On : Tue Mar 26 23:10:36 2019
    13 // Update Count     : 95
     12// Last Modified On : Sun Mar 10 10:52:50 2019
     13// Update Count     : 31
    1414//
    1515
    … …  
    1818typedef unsigned long long __cfaabi_abi_exception_type_t;
    1919
    20 #include <limits.h>                                                                             // CHAR_BIT
    2120#include "../src/virtual.h"
    2221#include "../src/exception.h"
    … …  
    2726// increment/decrement unification
    2827
    29 static inline {
    30         forall( dtype DT | { DT & ?+=?( DT &, one_t ); } )
    31         DT & ++?( DT & x ) { return x += 1; }
     28static inline forall( dtype T | { T & ?+=?( T &, one_t ); } )
     29T & ++? ( T & x ) { return x += 1; }
    3230
    33         forall( dtype DT | sized(DT) | { void ?{}( DT &, DT ); void ^?{}( DT & ); DT & ?+=?( DT &, one_t ); } )
    34         DT & ?++( DT & x ) { DT tmp = x; x += 1; return tmp; }
     31static inline forall( dtype T | sized(T) | { void ?{}( T &, T ); void ^?{}( T & ); T & ?+=?( T &, one_t ); } )
     32T & ?++ ( T & x ) { T tmp = x; x += 1; return tmp; }
    3533
    36         forall( dtype DT | { DT & ?-=?( DT &, one_t ); } )
    37         DT & --?( DT & x ) { return x -= 1; }
     34static inline forall( dtype T | { T & ?-=?( T &, one_t ); } )
     35T & --? ( T & x ) { return x -= 1; }
    3836
    39         forall( dtype DT | sized(DT) | { void ?{}( DT &, DT ); void ^?{}( DT & ); DT & ?-=?( DT &, one_t ); } )
    40         DT & ?--( DT & x ) { DT tmp = x; x -= 1; return tmp; }
    41 } // distribution
     37static inline forall( dtype T | sized(T) | { void ?{}( T &, T ); void ^?{}( T & ); T & ?-=?( T &, one_t ); } )
     38T & ?-- ( T & x ) { T tmp = x; x -= 1; return tmp; }
    4239
    4340// universal typed pointer constant
    44 // Compiler issue: there is a problem with anonymous types that do not have a size.
    45 static inline forall( dtype DT | sized(DT) ) DT * intptr( uintptr_t addr ) { return (DT *)addr; }
     41
     42static inline forall( dtype T ) T * intptr( uintptr_t addr ) { return (T *)addr; }
    4643
    4744// exponentiation operator implementation
    … …  
    5653} // extern "C"
    5754
    58 static inline {
    59         float ?\?( float x, float y ) { return powf( x, y ); }
    60         double ?\?( double x, double y ) { return pow( x, y ); }
    61         long double ?\?( long double x, long double y ) { return powl( x, y ); }
    62         float _Complex ?\?( float _Complex x, _Complex float y ) { return cpowf(x, y ); }
    63         double _Complex ?\?( double _Complex x, _Complex double y ) { return cpow( x, y ); }
    64         long double _Complex ?\?( long double _Complex x, _Complex long double y ) { return cpowl( x, y ); }
    65 } // distribution
     55static inline float ?\?( float x, float y ) { return powf( x, y ); }
     56static inline double ?\?( double x, double y ) { return pow( x, y ); }
     57static inline long double ?\?( long double x, long double y ) { return powl( x, y ); }
     58static inline float _Complex ?\?( float _Complex x, _Complex float y ) { return cpowf(x, y ); }
     59static inline double _Complex ?\?( double _Complex x, _Complex double y ) { return cpow( x, y ); }
     60static inline long double _Complex ?\?( long double _Complex x, _Complex long double y ) { return cpowl( x, y ); }
    6661
    67 #define __CFA_BASE_COMP_1__() if ( ep == 1 ) return 1
    68 #define __CFA_BASE_COMP_2__() if ( ep == 2 ) return ep << (y - 1)
    69 #define __CFA_EXP_OVERFLOW__() if ( y >= sizeof(y) * CHAR_BIT ) return 0
     62static inline long int ?\?( long int ep, unsigned long int y ) { // disallow negative exponent
     63        if ( y == 0 ) return 1;                                                         // base case
     64        if ( ep == 2 ) return ep << (y - 1);                            // special case, positive shifting only
     65        typeof( ep ) op = 1;                                                            // accumulate odd product
     66        for ( ; y > 1; y >>= 1 ) {                                                      // squaring exponentiation, O(log2 y)
     67                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
     68                ep *= ep;
     69        } // for
     70        return ep * op;
     71} // ?\?
    7072
    71 #define __CFA_EXP__() \
    72         if ( y == 0 ) return 1;                                                         /* convention */ \
    73         __CFA_BASE_COMP_1__();                                                          /* base case */ \
    74         __CFA_BASE_COMP_2__();                                                          /* special case, positive shifting for integral types */ \
    75         __CFA_EXP_OVERFLOW__();                                                         /* immediate overflow, negative exponent > 2^size-1 */ \
    76         typeof(ep) op = 1;                                                                      /* accumulate odd product */ \
    77         for ( ; y > 1; y >>= 1 ) {                                                      /* squaring exponentiation, O(log2 y) */ \
    78                 if ( (y & 1) == 1 ) op = op * ep;                               /* odd ? */ \
    79                 ep = ep * ep; \
    80         } \
    81         return ep * op
     73static inline forall( otype T | { void ?{}( T & this, one_t ); T ?*?( T, T ); } )
     74T ?\?( T ep, unsigned long int y ) {
     75        if ( y == 0 ) return 1;
     76        T op = 1;
     77        for ( ; y > 1; y >>= 1 ) {                                                      // squaring exponentiation, O(log2 y)
     78                if ( (y & 1) == 1 ) op = op * ep;                               // odd ?
     79                ep = ep * ep;
     80        } // for
     81        return ep * op;
     82} // ?\?
    8283
    83 static inline {
    84         long int ?\?( int ep, unsigned int y ) { __CFA_EXP__(); }
    85         long int ?\?( long int ep, unsigned long int y ) { __CFA_EXP__(); }
    86         // unsigned computation may be faster and larger
    87         unsigned long int ?\?( unsigned int ep, unsigned int y ) { __CFA_EXP__(); }
    88         unsigned long int ?\?( unsigned long int ep, unsigned long int y ) { __CFA_EXP__(); }
    89 } // distribution
     84// unsigned computation may be faster and larger
     85static inline unsigned long int ?\?( unsigned long int ep, unsigned long int y ) { // disallow negative exponent
     86        if ( y == 0 ) return 1;                                                         // base case
     87        if ( ep == 2 ) return ep << (y - 1);                            // special case, positive shifting only
     88        typeof( ep ) op = 1;                                                            // accumulate odd product
     89        for ( ; y > 1; y >>= 1 ) {                                                      // squaring exponentiation, O(log2 y)
     90                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
     91                ep *= ep;
     92        } // for
     93        return ep * op;
     94} // ?\?
    9095
    91 #undef __CFA_BASE_COMP_1__
    92 #undef __CFA_BASE_COMP_2__
    93 #undef __CFA_EXP_OVERFLOW__
    94 #define __CFA_BASE_COMP_1__()
    95 #define __CFA_BASE_COMP_2__()
    96 #define __CFA_EXP_OVERFLOW__()
     96static inline double ?\?( long int x, signed long int y ) {     // allow negative exponent
     97        if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
     98        else return 1.0 / x \ (unsigned int)(-y);
     99} // ?\?
    97100
    98 static inline forall( otype OT | { void ?{}( OT & this, one_t ); OT ?*?( OT, OT ); } ) {
    99         OT ?\?( OT ep, unsigned int y ) { __CFA_EXP__(); }
    100         OT ?\?( OT ep, unsigned long int y ) { __CFA_EXP__(); }
    101 } // distribution
     101// FIXME (x \ (unsigned long int)y) relies on X ?\?(T, unsigned long) a function that is neither
     102// defined, nor passed as an assertion parameter. Without user-defined conversions, cannot specify
     103// X as a type that casts to double, yet it doesn't make sense to write functions with that type
     104// signature where X is double.
    102105
    103 #undef __CFA_BASE_COMP_1__
    104 #undef __CFA_BASE_COMP_2__
    105 #undef __CFA_EXP_OVERFLOW__
     106// static inline forall( otype T | { void ?{}( T & this, one_t ); T ?*?( T, T ); double ?/?( double, T ); } )
     107// double ?\?( T x, signed long int y ) {
     108//     if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
     109//     else return 1.0 / x \ (unsigned long int)(-y);
     110// } // ?\?
    106111
    107 static inline {
    108         long int ?\=?( long int & x, unsigned long int y ) { x = x \ y; return x; }
    109         unsigned long int ?\=?( unsigned long int & x, unsigned long int y ) { x = x \ y; return x; }
    110         int ?\=?( int & x, unsigned long int y ) { x = x \ y; return x; }
    111         unsigned int ?\=?( unsigned int & x, unsigned long int y ) { x = x \ y; return x; }
    112 } // distribution
     112static inline long int ?\=?( long int & x, unsigned long int y ) { x = x \ y; return x; }
     113static inline unsigned long int ?\=?( unsigned long int & x, unsigned long int y ) { x = x \ y; return x; }
     114static inline int ?\=?( int & x, unsigned long int y ) { x = x \ y; return x; }
     115static inline unsigned int ?\=?( unsigned int & x, unsigned long int y ) { x = x \ y; return x; }
    113116
    114117// Local Variables: //
Note: See TracChangeset for help on using the changeset viewer.