source: src/prelude/builtins.c@ 300d75b

ADT aaron-thesis arm-eh ast-experimental cleanup-dtors deferred_resn demangler enum forall-pointer-decay jacob/cs343-translation jenkins-sandbox new-ast new-ast-unique-expr new-env no_list persistent-indexer pthread-emulation qualifiedEnum resolv-new with_gc
Last change on this file since 300d75b was 36982fc, checked in by Thierry Delisle <tdelisle@…>, 8 years ago

Renamed internal stuff to cfaabi_...

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