source: src/prelude/builtins.c @ d2887f7

aaron-thesisarm-ehcleanup-dtorsdeferred_resndemanglerjacob/cs343-translationjenkins-sandboxnew-astnew-ast-unique-exprnew-envno_listpersistent-indexerresolv-newwith_gc
Last change on this file since d2887f7 was d2887f7, checked in by Aaron Moss <a3moss@…>, 4 years ago

Default increment/decrement operators defined in terms of one_t

  • Property mode set to 100644
File size: 4.8 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 : Peter A. Buhr
12// Last Modified On : Thu Feb  8 12:47:59 2018
13// Update Count     : 19
14//
15
16// exception implementation
17
18typedef unsigned long long __cfaabi_abi_exception_type_t;
19
20#include "../libcfa/virtual.h"
21#include "../libcfa/exception.h"
22
23void exit( int status, const char fmt[], ... ) __attribute__ (( format(printf, 2, 3), __nothrow__, __leaf__, __noreturn__ ));
24void abort ( const char fmt[], ... ) __attribute__ (( format(printf, 1, 2), __nothrow__, __leaf__, __noreturn__ ));
25
26// increment/decrement unification
27
28static inline forall( dtype T | { T& ?+=?( T&, one_t ); } )
29T& ++? ( T& x ) { return x += 1; }
30
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; }
33
34static inline forall( dtype T | { T& ?-=?( T&, one_t ); } )
35T& --? ( T& x ) { return x -= 1; }
36
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; }
39
40// exponentiation operator implementation
41
42extern "C" {
43        float powf( float x, float y );
44        double pow( double x, double y );
45        long double powl( long double x, long double y );
46        float _Complex cpowf( float _Complex x, _Complex float z );
47        double _Complex cpow( double _Complex x, _Complex double z );
48        long double _Complex cpowl( long double _Complex x, _Complex long double z );
49} // extern "C"
50
51static inline float ?\?( float x, float y ) { return powf( x, y ); }
52static inline double ?\?( double x, double y ) { return pow( x, y ); }
53static inline long double ?\?( long double x, long double y ) { return powl( x, y ); }
54static inline float _Complex ?\?( float _Complex x, _Complex float y ) { return cpowf(x, y ); }
55static inline double _Complex ?\?( double _Complex x, _Complex double y ) { return cpow( x, y ); }
56static inline long double _Complex ?\?( long double _Complex x, _Complex long double y ) { return cpowl( x, y ); }
57
58static inline long int ?\?( long int ep, unsigned long int y ) { // disallow negative exponent
59        if ( y == 0 ) return 1;                                                         // base case
60    if ( ep == 2 ) return ep << (y - 1);                                // special case, positive shifting only
61    typeof( ep ) op = 1;                                                                // accumulate odd product
62    for ( ; y > 1; y >>= 1 ) {                                                  // squaring exponentiation, O(log2 y)
63                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
64                ep *= ep;
65        } // for
66    return ep * op;
67} // ?\?
68
69// FIX ME, cannot resolve the "T op = 1".
70
71// static inline forall( otype T | { void ?{}( T * this, one_t ); T ?*?( T, T ); } )
72// T ?\?( T ep, unsigned long int y ) {
73//     if ( y == 0 ) return 1;
74//     T op = 1;
75//     for ( ; y > 1; y >>= 1 ) {                                                       // squaring exponentiation, O(log2 y)
76//              if ( (y & 1) == 1 ) op = op * ep;                               // odd ?
77//              ep = ep * ep;
78//     } // for
79//     return ep * op;
80// } // ?\?
81
82// unsigned computation may be faster and larger
83static inline unsigned long int ?\?( unsigned long int ep, unsigned long int y ) { // disallow negative exponent
84        if ( y == 0 ) return 1;                                                         // base case
85    if ( ep == 2 ) return ep << (y - 1);                                // special case, positive shifting only
86    typeof( ep ) op = 1;                                                                // accumulate odd product
87    for ( ; y > 1; y >>= 1 ) {                                                  // squaring exponentiation, O(log2 y)
88                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
89                ep *= ep;
90        } // for
91    return ep * op;
92} // ?\?
93
94static inline double ?\?( long int x, signed long int y ) {     // allow negative exponent
95    if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
96    else return 1.0 / x \ (unsigned int)(-y);
97} // ?\?
98
99// FIXME (x \ (unsigned long int)y) relies on X ?\?(T, unsigned long) a function that is neither
100// defined, nor passed as an assertion parameter. Without user-defined conversions, cannot specify
101// X as a type that casts to double, yet it doesn't make sense to write functions with that type
102// signature where X is double.
103
104// static inline forall( otype T | { void ?{}( T & this, one_t ); T ?*?( T, T ); double ?/?( double, T ); } )
105// double ?\?( T x, signed long int y ) {
106//     if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
107//     else return 1.0 / x \ (unsigned long int)(-y);
108// } // ?\?
109
110static inline long int ?\=?( long int & x, unsigned long int y ) { x = x \ y; return x; }
111static inline unsigned long int ?\=?( unsigned long int & x, unsigned long int y ) { x = x \ y; return x; }
112static inline int ?\=?( int & x, unsigned long int y ) { x = x \ y; return x; }
113static inline unsigned int ?\=?( unsigned int & x, unsigned long int y ) { x = x \ y; return x; }
114
115// Local Variables: //
116// mode: c //
117// tab-width: 4 //
118// End: //
Note: See TracBrowser for help on using the repository browser.