source: src/prelude/builtins.c @ df7a162

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

Allow builtin objects to be constructed

  • Property mode set to 100644
File size: 4.2 KB
RevLine 
[a5f0529]1//
[a1edafa]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.
[a5f0529]6//
7// builtins.c --
8//
[a1edafa]9// Author           : Peter A. Buhr
10// Created On       : Fri Jul 21 16:21:03 2017
[169d944]11// Last Modified By : Peter A. Buhr
12// Last Modified On : Thu Feb  8 12:47:59 2018
13// Update Count     : 19
[a5f0529]14//
[a1edafa]15
16// exception implementation
17
[36982fc]18typedef unsigned long long __cfaabi_abi_exception_type_t;
[fa4805f]19
[a5f0529]20#include "../libcfa/virtual.h"
[fa4805f]21#include "../libcfa/exception.h"
[a1edafa]22
[169d944]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
[a1edafa]26// exponentiation operator implementation
27
28extern "C" {
29        float powf( float x, float y );
30        double pow( double x, double y );
31        long double powl( long double x, long double y );
32        float _Complex cpowf( float _Complex x, _Complex float z );
33        double _Complex cpow( double _Complex x, _Complex double z );
34        long double _Complex cpowl( long double _Complex x, _Complex long double z );
35} // extern "C"
36
37static inline float ?\?( float x, float y ) { return powf( x, y ); }
38static inline double ?\?( double x, double y ) { return pow( x, y ); }
39static inline long double ?\?( long double x, long double y ) { return powl( x, y ); }
40static inline float _Complex ?\?( float _Complex x, _Complex float y ) { return cpowf(x, y ); }
41static inline double _Complex ?\?( double _Complex x, _Complex double y ) { return cpow( x, y ); }
42static inline long double _Complex ?\?( long double _Complex x, _Complex long double y ) { return cpowl( x, y ); }
43
[fe97a7d8]44static inline long int ?\?( long int ep, unsigned long int y ) { // disallow negative exponent
[a1edafa]45        if ( y == 0 ) return 1;                                                         // base case
[df7a162]46        if ( ep == 2 ) return ep << (y - 1);                            // special case, positive shifting only
47        typeof( ep ) op = 1;                                                            // accumulate odd product
48        for ( ; y > 1; y >>= 1 ) {                                                      // squaring exponentiation, O(log2 y)
[fe97a7d8]49                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
50                ep *= ep;
51        } // for
[df7a162]52        return ep * op;
[a1edafa]53} // ?\?
54
[df7a162]55static inline forall( otype T | { void ?{}( T & this, one_t ); T ?*?( T, T ); } )
56T ?\?( T ep, unsigned long int y ) {
57        if ( y == 0 ) return 1;
58        T op = 1;
59        for ( ; y > 1; y >>= 1 ) {                                                      // squaring exponentiation, O(log2 y)
60                if ( (y & 1) == 1 ) op = op * ep;                               // odd ?
61                ep = ep * ep;
62        } // for
63        return ep * op;
64} // ?\?
[a1edafa]65
66// unsigned computation may be faster and larger
[fe97a7d8]67static inline unsigned long int ?\?( unsigned long int ep, unsigned long int y ) { // disallow negative exponent
[a1edafa]68        if ( y == 0 ) return 1;                                                         // base case
[df7a162]69        if ( ep == 2 ) return ep << (y - 1);                            // special case, positive shifting only
70        typeof( ep ) op = 1;                                                            // accumulate odd product
71        for ( ; y > 1; y >>= 1 ) {                                                      // squaring exponentiation, O(log2 y)
[fe97a7d8]72                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
73                ep *= ep;
74        } // for
[df7a162]75        return ep * op;
[a1edafa]76} // ?\?
77
[fe97a7d8]78static inline double ?\?( long int x, signed long int y ) {     // allow negative exponent
[df7a162]79        if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
80        else return 1.0 / x \ (unsigned int)(-y);
[a1edafa]81} // ?\?
82
[36982fc]83// FIXME (x \ (unsigned long int)y) relies on X ?\?(T, unsigned long) a function that is neither
84// defined, nor passed as an assertion parameter. Without user-defined conversions, cannot specify
85// X as a type that casts to double, yet it doesn't make sense to write functions with that type
[e472d54]86// signature where X is double.
87
88// static inline forall( otype T | { void ?{}( T & this, one_t ); T ?*?( T, T ); double ?/?( double, T ); } )
89// double ?\?( T x, signed long int y ) {
90//     if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
91//     else return 1.0 / x \ (unsigned long int)(-y);
92// } // ?\?
[a1edafa]93
[1ae06fa]94static inline long int ?\=?( long int & x, unsigned long int y ) { x = x \ y; return x; }
95static inline unsigned long int ?\=?( unsigned long int & x, unsigned long int y ) { x = x \ y; return x; }
96static inline int ?\=?( int & x, unsigned long int y ) { x = x \ y; return x; }
97static inline unsigned int ?\=?( unsigned int & x, unsigned long int y ) { x = x \ y; return x; }
[a1edafa]98
99// Local Variables: //
100// mode: c //
101// tab-width: 4 //
102// End: //
Note: See TracBrowser for help on using the repository browser.