Changeset 9236060 for src/prelude

Timestamp:

Aug 14, 2017, 2:03:39 PM (8 years ago)

Author:

Rob Schluntz <rschlunt@…>

Branches:

ADT, aaron-thesis, arm-eh, ast-experimental, cleanup-dtors, deferred_resn, demangler, enum, forall-pointer-decay, jacob/cs343-translation, jenkins-sandbox, master, new-ast, new-ast-unique-expr, new-env, no_list, persistent-indexer, pthread-emulation, qualifiedEnum, resolv-new, with_gc

Children:

74b007ba

Parents:

fd344aa (diff), 54cd58b0 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

Message:

Merge branch 'master' into references

File:

• src/prelude/builtins.c (modified) (1 diff)

src/prelude/builtins.c

@@ +1 @@
+//
+// Cforall Version 1.0.0 Copyright (C) 2016 University of Waterloo
+//
+// The contents of this file are covered under the licence agreement in the
+// file "LICENCE" distributed with Cforall.
+//
+// builtins.c --
+//
+// Author           : Peter A. Buhr
+// Created On       : Fri Jul 21 16:21:03 2017
+// Last Modified By : Andrew Beach
+// Last Modified On : Tus Jul 25 15:33:00 2017
+// Update Count     : 14
+//
+
+// exception implementation
+
 typedef unsigned long long __cfaabi_exception_type_t;
 
+#include "../libcfa/virtual.h"
 #include "../libcfa/exception.h"
+
+// exponentiation operator implementation
+
+extern "C" {
+        float powf( float x, float y );
+        double pow( double x, double y );
+        long double powl( long double x, long double y );
+        float _Complex cpowf( float _Complex x, _Complex float z );
+        double _Complex cpow( double _Complex x, _Complex double z );
+        long double _Complex cpowl( long double _Complex x, _Complex long double z );
+} // extern "C"
+
+static inline float ?\?( float x, float y ) { return powf( x, y ); }
+static inline double ?\?( double x, double y ) { return pow( x, y ); }
+static inline long double ?\?( long double x, long double y ) { return powl( x, y ); }
+static inline float _Complex ?\?( float _Complex x, _Complex float y ) { return cpowf(x, y ); }
+static inline double _Complex ?\?( double _Complex x, _Complex double y ) { return cpow( x, y ); }
+static inline long double _Complex ?\?( long double _Complex x, _Complex long double y ) { return cpowl( x, y ); }
+
+static inline long int ?\?( long int ep, unsigned long int y ) { // disallow negative exponent
+        if ( y == 0 ) return 1;                                                         // base case
+    if ( ep == 2 ) return ep << (y - 1);                                // special case, positive shifting only
+    typeof( ep ) op = 1;                                                                // accumulate odd product
+    for ( ; y > 1; y >>= 1 ) {                                                  // squaring exponentiation, O(log2 y)
+                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
+                ep *= ep;
+        } // for
+    return ep * op;
+} // ?\?
+
+// FIX ME, cannot resolve the "T op = 1".
+
+// static inline forall( otype T | { void ?{}( T * this, one_t ); T ?*?( T, T ); } )
+// T ?\?( T ep, unsigned long int y ) {
+//     if ( y == 0 ) return 1;
+//     T op = 1;
+//     for ( ; y > 1; y >>= 1 ) {                                                       // squaring exponentiation, O(log2 y)
+//              if ( (y & 1) == 1 ) op = op * ep;                               // odd ?
+//              ep = ep * ep;
+//     } // for
+//     return ep * op;
+// } // ?\?
+
+// unsigned computation may be faster and larger
+static inline unsigned long int ?\?( unsigned long int ep, unsigned long int y ) { // disallow negative exponent
+        if ( y == 0 ) return 1;                                                         // base case
+    if ( ep == 2 ) return ep << (y - 1);                                // special case, positive shifting only
+    typeof( ep ) op = 1;                                                                // accumulate odd product
+    for ( ; y > 1; y >>= 1 ) {                                                  // squaring exponentiation, O(log2 y)
+                if ( (y & 1) == 1 ) op *= ep;                                   // odd ?
+                ep *= ep;
+        } // for
+    return ep * op;
+} // ?\?
+
+static inline double ?\?( long int x, signed long int y ) {     // allow negative exponent
+    if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
+    else return 1.0 / x \ (unsigned int)(-y);
+} // ?\?
+
+static inline forall( otype T | { void ?{}( T & this, one_t ); T ?*?( T, T ); double ?/?( double, T ); } )
+double ?\?( T x, signed long int y ) {
+    if ( y >=  0 ) return (double)(x \ (unsigned long int)y);
+    else return 1.0 / x \ (unsigned long int)(-y);
+} // ?\?
+
+static inline long int ?\=?( long int * x, unsigned long int y ) { *x = *x \ y; return *x; }
+static inline unsigned long int ?\=?( unsigned long int * x, unsigned long int y ) { *x = *x \ y; return *x; }
+static inline int ?\=?( int * x, unsigned long int y ) { *x = *x \ y; return *x; }
+static inline unsigned int ?\=?( unsigned int * x, unsigned long int y ) { *x = *x \ y; return *x; }
+
+// Local Variables: //
+// mode: c //
+// tab-width: 4 //
+// End: //