Index: src/examples/math.c
===================================================================
--- src/examples/math.c	(revision add371c15f16e8ece211a386eec66e2e03a21afa)
+++ src/examples/math.c	(revision add371c15f16e8ece211a386eec66e2e03a21afa)
@@ -0,0 +1,144 @@
+// 
+// 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.
+// 
+// math.c -- 
+// 
+// Author           : Peter A. Buhr
+// Created On       : Fri Apr 22 14:59:21 2016
+// Last Modified By : Peter A. Buhr
+// Last Modified On : Sun Apr 24 13:24:20 2016
+// Update Count     : 70
+// 
+
+#include <fstream>
+#include <math>
+
+int main( void ) {
+	float f;
+	double d;
+	long double l;
+
+	sout | "fabs:" | fabs( -1.0F ) | fabs( -1.0D ) | fabs( -1.0L ) | cabs( -1.0F+1.0FI ) | cabs( -1.0D+1.0DI ) | cabs( -1.0DL+1.0LI ) | endl;
+	sout | "fmod:" | 5.0F % -2.0F | fmod( 5.0F, -2.0F ) | 5.0D % -2.0D | fmod( 5.0D, -2.0D ) | 5.0L % -2.0L | fmod( 5.0L, -2.0L ) | endl;
+	sout | "remainder:" | remainder( 2.0F, 3.0F ) | remainder( 2.0D, 3.0D ) | remainder( 2.0L, 3.0L ) | endl;
+	int quot;
+	f = remquo( 3.6F, 0.5F, &quot );
+	sout | "remquo:" | quot | f;
+	d = remquo( 3.6D, 0.5F, &quot );
+	sout | quot | d;
+	l = remquo( 3.6L, 0.5L, &quot );
+	sout | quot | l | endl;
+	f = div( 3.6F, 0.5F, &quot );
+	sout | "div:" | quot | f;
+	d = div( 3.6D, 0.5F, &quot );
+	sout | quot | d;
+	l = div( 3.6L, 0.5L, &quot );
+	sout | quot | l | endl;
+	sout | "fma:" | fma( 3.0F, -1.0F, 1.0F ) | fma( 3.0D, -1.0D, 1.0D ) | fma( 3.0L, -1.0L, , 1.0L ) | endl;
+	sout | "fdim:" | fdim( 1.0F, -1.0F ) | fdim( 1.0D, -1.0D ) | fdim( 1.0L, -1.0L ) | endl;
+	sout | "nan:" | (float)nan( "" ) | (double)nan( "" ) | (long double)nan( "" ) | endl;
+
+	//---------------------- Exponential ----------------------
+
+	sout | "exp:" | exp( 1.0F ) | exp( 1.0D ) | exp( 1.0L ) | exp( 1.0F+1.0FI ) | exp( 1.0D+1.0DI ) | exp( 1.0DL+1.0LI ) | endl;
+	sout | "exp2:" | exp2( 1.0F ) | exp2( 1.0D ) | exp2( 1.0L ) | endl;
+	sout | "expm1:" | expm1( 1.0F ) | expm1( 1.0D ) | expm1( 1.0L ) | endl;
+	sout | "log:" | log( 1.0F ) | log( 1.0D ) | log( 1.0L ) | log( 1.0F+1.0FI ) | log( 1.0D+1.0DI ) | log( 1.0DL+1.0LI ) | endl;
+	sout | "log2:" | log2( 8.0F ) | log2( 8.0D ) | log2( 8.0L ) | endl;
+	sout | "log10:" | log10( 100.0F ) | log10( 100.0D ) | log10( 100.0L ) | endl;
+	sout | "log1p:" | log1p( 1.0F ) | log1p( 1.0D ) | log1p( 1.0L ) | endl;
+	sout | "ilogb:" | ilogb( 1.0F ) | ilogb( 1.0D ) | ilogb( 1.0L ) | endl;
+	sout | "logb:" | logb( 8.0F ) | logb( 8.0D ) | logb( 8.0L ) | endl;
+
+	//---------------------- Power ----------------------
+
+	sout | "sqrt:" | sqrt( 1.0F ) | sqrt( 1.0D ) | sqrt( 1.0L ) | sqrt( 1.0F+1.0FI ) | sqrt( 1.0D+1.0DI ) | sqrt( 1.0DL+1.0LI ) | endl;
+	sout | "cbrt:" | cbrt( 27.0F ) | cbrt( 27.0D ) | cbrt( 27.0L ) | endl;
+	sout | "hypot:" | hypot( 1.0F, -1.0F ) | hypot( 1.0D, -1.0D ) | hypot( 1.0L, -1.0L ) | endl;
+	sout | "pow:" | pow( 1.0F, 1.0F ) | pow( 1.0D, 1.0D ) | pow( 1.0L, 1.0L ) | pow( 1.0F+1.0FI, 1.0F+1.0FI ) | pow( 1.0D+1.0DI, 1.0D+1.0DI ) | pow( 1.0DL+1.0LI, 1.0DL+1.0LI ) | endl;
+
+	//---------------------- Trigonometric ----------------------
+
+	sout | "sin:" | sin( 1.0F ) | sin( 1.0D ) | sin( 1.0L ) | sin( 1.0F+1.0FI ) | sin( 1.0D+1.0DI ) | sin( 1.0DL+1.0LI ) | endl;
+	sout | "cos:" | cos( 1.0F ) | cos( 1.0D ) | cos( 1.0L ) | cos( 1.0F+1.0FI ) | cos( 1.0D+1.0DI ) | cos( 1.0DL+1.0LI ) | endl;
+	sout | "tan:" | tan( 1.0F ) | tan( 1.0D ) | tan( 1.0L ) | tan( 1.0F+1.0FI ) | tan( 1.0D+1.0DI ) | tan( 1.0DL+1.0LI ) | endl;
+	sout | "asin:" | asin( 1.0F ) | asin( 1.0D ) | asin( 1.0L ) | asin( 1.0F+1.0FI ) | asin( 1.0D+1.0DI ) | asin( 1.0DL+1.0LI ) | endl;
+	sout | "acos:" | acos( 1.0F ) | acos( 1.0D ) | acos( 1.0L ) | acos( 1.0F+1.0FI ) | acos( 1.0D+1.0DI ) | acos( 1.0DL+1.0LI ) | endl;
+	sout | "atan:" | atan( 1.0F ) | atan( 1.0D ) | atan( 1.0L ) | atan( 1.0F+1.0FI ) | atan( 1.0D+1.0DI ) | atan( 1.0DL+1.0LI ) | endl;
+	sout | "atan2:" | atan2( 1.0F, 1.0F ) | atan2( 1.0D, 1.0D ) | atan2( 1.0L, 1.0L );
+	sout | "atan:" | atan( 1.0F, 1.0F ) | atan( 1.0D, 1.0D ) | atan( 1.0L, 1.0L );
+
+	//---------------------- Hyperbolic ----------------------
+
+	sout | "sinh:" | sinh( 1.0F ) | sinh( 1.0D ) | sinh( 1.0L ) | sinh( 1.0F+1.0FI ) | sinh( 1.0D+1.0DI ) | sinh( 1.0DL+1.0LI ) | endl;
+	sout | "cosh:" | cosh( 1.0F ) | cosh( 1.0D ) | cosh( 1.0L ) | cosh( 1.0F+1.0FI ) | cosh( 1.0D+1.0DI ) | cosh( 1.0DL+1.0LI ) | endl;
+	sout | "tanh:" | tanh( 1.0F ) | tanh( 1.0D ) | tanh( 1.0L ) | tanh( 1.0F+1.0FI ) | tanh( 1.0D+1.0DI ) | tanh( 1.0DL+1.0LI ) | endl;
+	sout | "acosh:" | acosh( 1.0F ) | acosh( 1.0D ) | acosh( 1.0L ) | acosh( 1.0F+1.0FI ) | acosh( 1.0D+1.0DI ) | acosh( 1.0DL+1.0LI ) | endl;
+	sout | "asinh:" | asinh( 1.0F ) | asinh( 1.0D ) | asinh( 1.0L ) | asinh( 1.0F+1.0FI ) | asinh( 1.0D+1.0DI ) | asinh( 1.0DL+1.0LI ) | endl;
+	sout | "atanh:" | atanh( 1.0F ) | atanh( 1.0D ) | atanh( 1.0L ) | atanh( 1.0F+1.0FI ) | atanh( 1.0D+1.0DI ) | atanh( 1.0DL+1.0LI ) | endl;
+
+	//---------------------- Error / Gamma ----------------------
+
+	sout | "erf:" | erf( 1.0F ) | erf( 1.0D ) | erf( 1.0L ) | endl;
+	sout | "erfc:" | erfc( 1.0F ) | erfc( 1.0D ) | erfc( 1.0L ) | endl;
+	sout | "lgamma:" | lgamma( 4.0F ) | lgamma( 4.0D ) | lgamma( 4.0L ) | endl;
+	int sign;
+	f = lgamma( 4.0F, &sign );
+	sout | "lgamma:" | f | sign;
+	d = lgamma( 4.0D, &sign );
+	sout | d | sign;
+	l = lgamma( 4.0L, &sign );
+	sout | l | sign | endl;
+	sout | "tgamma:" | tgamma( 4.0F ) | tgamma( 4.0D ) | tgamma( 4.0L ) | endl;
+
+	//---------------------- Nearest Integer ----------------------
+
+	sout | "floor:" | floor( 1.2F ) | floor( 1.2D ) | floor( 1.2L ) | endl;
+	sout | "ceil:" | ceil( 1.6F ) | ceil( 1.6D ) | ceil( 1.6L ) | endl;
+	sout | "trunc:" | trunc( 3.5F ) | trunc( 3.5D ) | trunc( 3.5L ) | endl;
+	sout | "rint:" | (float)rint( 1.5F ) | (double)rint( 1.5D ) | (long double)rint( 1.5L ) | endl;
+	sout | "rint:" | (long int)rint( 1.5F ) | (long int)rint( 1.5D ) | (long int)rint( 1.5L ) | endl;
+	sout | "rint:" | (long long int)rint( 1.5F ) | (long long int)rint( 1.5D ) | (long long int)rint( 1.5L ) | endl;
+	sout | "lrint:" | lrint( 1.5F ) | lrint( 1.5D ) | lrint( 1.5L ) | endl;
+	sout | "llrint:" | llrint( 1.5F ) | llrint( 1.5D ) | llrint( 1.5L ) | endl;
+	sout | "nearbyint:" | nearbyint( 3.5F ) | nearbyint( 3.5D ) | nearbyint( 3.5L ) | endl;
+	sout | "round:" | (float)round( 1.5F ) | (double)round( 1.5D ) | (long double)round( 1.5L ) | endl;
+	sout | "round:" | (long int)round( 1.5F ) | (long int)round( 1.5D ) | (long int)round( 1.5L ) | endl;
+	sout | "round:" | (long long int)round( 1.5F ) | (long long int)round( 1.5D ) | (long long int)round( 1.5L ) | endl;
+	sout | "lround:" | lround( 1.5F ) | lround( 1.5D ) | lround( 1.5L ) | endl;
+	sout | "llround:" | llround( 1.5F ) | llround( 1.5D ) | llround( 1.5L ) | endl;
+
+	//---------------------- Manipulation ----------------------
+
+	sout | "copysign:" | copysign( 1.0F, -1.0F ) | copysign( 1.0D, -1.0D ) | copysign( 1.0L, -1.0L ) | endl;
+	int exp;
+	f = frexp( 4.0F, &exp );
+	sout | "frexp:" | f | exp;
+	d = frexp( 4.0D, &exp );
+	sout | d | exp;
+	l = frexp( 4.0L, &exp );
+	sout | l | exp | endl;
+	sout | "ldexp:" | ldexp( 2.0F, 2 ) | ldexp( 2.0D, 2 ) | ldexp( 2.0L, 2 ) | endl;
+	float fi;
+	double di;
+	long double ldi;
+	f = modf( 2.3F, &fi );
+	sout | "modf:" | fi | f;
+	d = modf( 2.3D, &di );
+	sout | di | d;
+	l = modf( 2.3L, &ldi );
+	sout | ldi | l;
+	sout | "nextafter:" | nextafter( 2.0F, 3.0F ) | nextafter( 2.0D, 3.0D ) | nextafter( 2.0L, 3.0L ) | endl;
+	sout | "nexttoward:" | nexttoward( 2.0F, 3.0F ) | nexttoward( 2.0D, 3.0D ) | nexttoward( 2.0L, 3.0L ) | endl;
+
+	sout | "scalbn:" | scalbn( 2.0F, 3 ) | scalbn( 2.0D, 3 ) | scalbn( 2.0L, 3 ) | endl;
+	sout | "scalbln:" | scalbln( 2.0F, 3L ) | scalbln( 2.0D, 3L ) | scalbln( 2.0L, 3L ) | endl;
+} // main
+
+// Local Variables: //
+// tab-width: 4 //
+// compile-command: "cfa math.c" //
+// End: //
Index: src/libcfa/math
===================================================================
--- src/libcfa/math	(revision 845cedcae61ca707834e349e69b07e41e766d7c0)
+++ src/libcfa/math	(revision add371c15f16e8ece211a386eec66e2e03a21afa)
@@ -10,10 +10,144 @@
 // Created On       : Mon Apr 18 23:37:04 2016
 // Last Modified By : Peter A. Buhr
-// Last Modified On : Thu Apr 21 07:56:48 2016
-// Update Count     : 38
+// Last Modified On : Sun Apr 24 12:45:02 2016
+// Update Count     : 59
 // 
 
+extern "C" {
+#include <math.h>										// fpclassify, isfinite, isnormal, isnan, isinf
+} // extern "C"
+
+float fabs( float );
+// extern "C" { double fabs( double ); }
+long double fabs( long double );
+float cabs( float _Complex );
+// extern "C" { double cabs( double _Complex ); }
+long double cabs( long double _Complex );
+
+float ?%?( float, float );
+float fmod( float, float );
+double ?%?( double, double );
+// extern "C" { double fmod( double, double ); }
+long double ?%?( long double, long double );
+long double fmod( long double, long double );
+
+float remainder( float, float );
+// extern "C" { double remainder( double, double ); }
+long double remainder( long double, long double );
+
+[ int, float ] remquo( float, float );
+float remquo( float, float, int * );
+[ int, double ] remquo( double, double );
+// extern "C" { double remquo( double, double, int * ); }
+[ int, long double ] remquo( long double, long double );
+long double remquo( long double, long double, int * );
+
+[ int, float ] div( float, float );						// alternative name for remquo
+float div( float, float, int * );
+[ int, double ] div( double, double );
+// extern "C" { double div( double, double, int * ); }
+[ int, long double ] div( long double, long double );
+long double div( long double, long double, int * );
+
+float fma( float, float, float );
+// extern "C" { double fma( double, double, double ); }
+long double fma( long double, long double, long double );
+
+float fdim( float, float );
+// extern "C" { double fdim( double, double ); }
+long double fdim( long double, long double );
+
+float nan( const char * );
+// extern "C" { double nan( const char * ); }
+long double nan( const char * );
+
+//---------------------- Exponential ----------------------
+
+float exp( float );
+// extern "C" { double exp( double ); }
+long double exp( long double );
+float _Complex exp( float _Complex );
+double _Complex exp( double _Complex );
+long double _Complex exp( long double _Complex );
+
+float exp2( float );
+// extern "C" { double exp2( double ); }
+long double exp2( long double );
+// float _Complex exp2( float _Complex );
+// double _Complex exp2( double _Complex );
+// long double _Complex exp2( long double _Complex );
+
+float expm1( float );
+// extern "C" { double expm1( double ); }
+long double expm1( long double );
+
+float log( float );
+// extern "C" { double log( double ); }
+long double log( long double );
+float _Complex log( float _Complex );
+double _Complex log( double _Complex );
+long double _Complex log( long double _Complex );
+
+float log2( float );
+// extern "C" { double log2( double ); }
+long double log2( long double );
+// float _Complex log2( float _Complex );
+// double _Complex log2( double _Complex );
+// long double _Complex log2( long double _Complex );
+
+float log10( float );
+// extern "C" { double log10( double ); }
+long double log10( long double );
+// float _Complex log10( float _Complex );
+// double _Complex log10( double _Complex );
+// long double _Complex log10( long double _Complex );
+
+float log1p( float );
+// extern "C" { double log1p( double ); }
+long double log1p( long double );
+
+int ilogb( float );
+// extern "C" { int ilogb( double ); }
+int ilogb( long double );
+
+float logb( float );
+// extern "C" { double logb( double ); }
+long double logb( long double );
+
+//---------------------- Power ----------------------
+
+float sqrt( float );
+// extern "C" { double sqrt( double ); }
+long double sqrt( long double );
+float _Complex sqrt( float _Complex );
+double _Complex sqrt( double _Complex );
+long double _Complex sqrt( long double _Complex );
+
+float cbrt( float );
+// extern "C" { double cbrt( double ); }
+long double cbrt( long double );
+
+float hypot( float, float );
+// extern "C" { double hypot( double, double ); }
+long double hypot( long double, long double );
+
+float pow( float, float );
+// extern "C" { double pow( double, double ); }
+long double pow( long double, long double );
+float _Complex pow( float _Complex, float _Complex );
+double _Complex pow( double _Complex, double _Complex );
+long double _Complex pow( long double _Complex, long double _Complex );
+
+//---------------------- Trigonometric ----------------------
+
+float sin( float );
+// extern "C" { double sin( double ); }
+long double sin( long double );
+float _Complex sin( float _Complex );
+double _Complex sin( double _Complex );
+long double _Complex sin( long double _Complex );
+
 float cos( float );
-extern "C" { double cos( double ); }
+// extern "C" { double cos( double ); }
 long double cos( long double );
 float _Complex cos( float _Complex );
@@ -21,6 +155,51 @@
 long double _Complex cos( long double _Complex );
 
+float tan( float );
+// extern "C" { double tan( double ); }
+long double tan( long double );
+float _Complex tan( float _Complex );
+double _Complex tan( double _Complex );
+long double _Complex tan( long double _Complex );
+
+float asin( float );
+// extern "C" { double asin( double ); }
+long double asin( long double );
+float _Complex asin( float _Complex );
+double _Complex asin( double _Complex );
+long double _Complex asin( long double _Complex );
+
+float acos( float );
+// extern "C" { double acos( double ); }
+long double acos( long double );
+float _Complex acos( float _Complex );
+double _Complex acos( double _Complex );
+long double _Complex acos( long double _Complex );
+
+float atan( float );
+// extern "C" { double atan( double ); }
+long double atan( long double );
+float _Complex atan( float _Complex );
+double _Complex atan( double _Complex );
+long double _Complex atan( long double _Complex );
+
+float atan2( float, float );
+double atan2( double, double );
+long double atan2( long double, long double );
+
+float atan( float, float );								// alternative name for atan2
+double atan( double, double );
+long double atan( long double, long double );
+
+//---------------------- Hyperbolic ----------------------
+
+float sinh( float );
+// extern "C" { double sinh( double ); }
+long double sinh( long double );
+float _Complex sinh( float _Complex );
+double _Complex sinh( double _Complex );
+long double _Complex sinh( long double _Complex );
+
 float cosh( float );
-extern "C" { double cosh( double ); }
+// extern "C" { double cosh( double ); }
 long double cosh( long double );
 float _Complex cosh( float _Complex );
@@ -28,13 +207,20 @@
 long double _Complex cosh( long double _Complex );
 
-float acos( float );
-extern "C" { double acos( double ); }
-long double acos( long double );
-float _Complex acos( float _Complex );
-double _Complex acos( double _Complex );
-long double _Complex acos( long double _Complex );
+float tanh( float );
+// extern "C" { double tanh( double ); }
+long double tanh( long double );
+float _Complex tanh( float _Complex );
+double _Complex tanh( double _Complex );
+long double _Complex tanh( long double _Complex );
+
+float asinh( float );
+// extern "C" { double asinh( double ); }
+long double asinh( long double );
+float _Complex asinh( float _Complex );
+double _Complex asinh( double _Complex );
+long double _Complex asinh( long double _Complex );
 
 float acosh( float );
-extern "C" { double acosh( double ); }
+// extern "C" { double acosh( double ); }
 long double acosh( long double );
 float _Complex acosh( float _Complex );
@@ -42,59 +228,6 @@
 long double _Complex acosh( long double _Complex );
 
-float sin( float );
-extern "C" { double sin( double ); }
-long double sin( long double );
-float _Complex sin( float _Complex );
-double _Complex sin( double _Complex );
-long double _Complex sin( long double _Complex );
-
-float sinh( float );
-extern "C" { double sinh( double ); }
-long double sinh( long double );
-float _Complex sinh( float _Complex );
-double _Complex sinh( double _Complex );
-long double _Complex sinh( long double _Complex );
-
-float asin( float );
-extern "C" { double asin( double ); }
-long double asin( long double );
-float _Complex asin( float _Complex );
-double _Complex asin( double _Complex );
-long double _Complex asin( long double _Complex );
-
-float asinh( float );
-extern "C" { double asinh( double ); }
-long double asinh( long double );
-float _Complex asinh( float _Complex );
-double _Complex asinh( double _Complex );
-long double _Complex asinh( long double _Complex );
-
-float tan( float );
-extern "C" { double tan( double ); }
-long double tan( long double );
-float _Complex tan( float _Complex );
-double _Complex tan( double _Complex );
-long double _Complex tan( long double _Complex );
-
-float tanh( float );
-extern "C" { double tanh( double ); }
-long double tanh( long double );
-float _Complex tanh( float _Complex );
-double _Complex tanh( double _Complex );
-long double _Complex tanh( long double _Complex );
-
-float atan( float );
-extern "C" { double atan( double ); }
-long double atan( long double );
-float _Complex atan( float _Complex );
-double _Complex atan( double _Complex );
-long double _Complex atan( long double _Complex );
-
-float _Complex atan( float, float );
-double atan( double, double );
-long double atan( long double, long double );
-
 float atanh( float );
-extern "C" { double atanh( double ); }
+// extern "C" { double atanh( double ); }
 long double atanh( long double );
 float _Complex atanh( float _Complex );
@@ -102,197 +235,117 @@
 long double _Complex atanh( long double _Complex );
 
-float exp( float );
-extern "C" { double exp( double ); }
-long double exp( long double );
-float _Complex exp( float _Complex );
-double _Complex exp( double _Complex );
-long double _Complex exp( long double _Complex );
-
-float log( float );
-extern "C" { double log( double ); }
-long double log( long double );
-float _Complex log( float _Complex );
-double _Complex log( double _Complex );
-long double _Complex log( long double _Complex );
-
-float pow( float, float );
-extern "C" { double pow( double, double ); }
-long double pow( long double, long double );
-float _Complex pow( float _Complex, float _Complex );
-double _Complex pow( double _Complex, double _Complex );
-long double _Complex pow( long double _Complex, long double _Complex );
-
-float sqrt( float );
-extern "C" { double sqrt( double ); }
-long double sqrt( long double );
-float _Complex sqrt( float _Complex );
-double _Complex sqrt( double _Complex );
-long double _Complex sqrt( long double _Complex );
-
-float fabs( float );
-extern "C" { double fabs( double ); }
-long double fabs( long double );
-float cabs( float _Complex );
-extern "C" { double cabs( double _Complex ); }
-long double cabs( long double _Complex );
+//---------------------- Error / Gamma ----------------------
+
+float erf( float );
+// extern "C" { double erf( double ); }
+long double erf( long double );
+// float _Complex erf( float _Complex );
+// double _Complex erf( double _Complex );
+// long double _Complex erf( long double _Complex );
+
+float erfc( float );
+// extern "C" { double erfc( double ); }
+long double erfc( long double );
+// float _Complex erfc( float _Complex );
+// double _Complex erfc( double _Complex );
+// long double _Complex erfc( long double _Complex );
+
+float lgamma( float );
+// extern "C" { double lgamma( double ); }
+long double lgamma( long double );
+float lgamma( float, int * );
+double lgamma( double, int * );
+long double lgamma( long double, int * );
+
+float tgamma( float );
+// extern "C" { double tgamma( double ); }
+long double tgamma( long double );
+
+//---------------------- Nearest Integer ----------------------
 
 float floor( float );
-extern "C" { double floor( double ); }
+// extern "C" { double floor( double ); }
 long double floor( long double );
 
 float ceil( float );
-extern "C" { double ceil( double ); }
+// extern "C" { double ceil( double ); }
 long double ceil( long double );
 
-//---------------------------------------
-
-float cbrt( float );
-extern "C" { double cbrt( double ); }
-long double cbrt( long double );
-
-float copysign( float, float );
-extern "C" { double copysign( double, double ); }
-long double copysign( long double, long double );
-
-float erf( float );
-extern "C" { double erf( double ); }
-long double erf( long double );
-
-float erfc( float );
-extern "C" { double erfc( double ); }
-long double erfc( long double );
-
-float exp2( float );
-extern "C" { double exp2( double ); }
-long double exp2( long double );
-
-float expm1( float );
-extern "C" { double expm1( double ); }
-long double expm1( long double );
-
-float fdim( float, float );
-extern "C" { double fdim( double, double ); }
-long double fdim( long double, long double );
-
-float fma( float, float, float );
-extern "C" { double fma( double, double, double ); }
-long double fma( long double, long double, long double );
-
-float ?%?( float, float );
-float fmod( float, float );
-double ?%?( double, double );
-extern "C" { double fmod( double, double ); }
-long double ?%?( long double, long double );
-long double fmod( long double, long double );
-
-[ int, float ] div( float x );
-float frexp( float, int * );
-[ int, double ] div( double x );
-extern "C" { double frexp( double, int * ); }
-[ int, long double ] div( long double x );
-long double frexp( long double, int * );
-
-float hypot( float, float );
-extern "C" { double hypot( double, double ); }
-long double hypot( long double, long double );
-
-int ilogb( float );
-extern "C" { int ilogb( double ); }
-int ilogb( long double );
-
-float ldexp( float, int );
-extern "C" { double ldexp( double, int ); }
-long double ldexp( long double, int );
-
-float lgamma( float );
-extern "C" { double lgamma( double ); }
-long double lgamma( long double );
-
-float rint( float x );
-long double rint( long double x );
-long int rint( float x );
-long int rint( double x );
-long int rint( long double x );
-long long int rint( float x );
-long long int rint( double x );
-long long int rint( long double x );
+float trunc( float );
+// extern "C" { double trunc( double ); }
+long double trunc( long double );
+
+float rint( float );
+long double rint( long double );
+long int rint( float );
+long int rint( double );
+long int rint( long double );
+long long int rint( float );
+long long int rint( double );
+long long int rint( long double );
 
 long int lrint( float );
-extern "C" { long int lrint( double ); }
+// extern "C" { long int lrint( double ); }
 long int lrint( long double );
 long long int llrint( float );
-extern "C" { long long int llrint( double ); }
+// extern "C" { long long int llrint( double ); }
 long long int llrint( long double );
 
-float round( float x );
-long double round( long double x );
-long int round( float x );
-long int round( double x );
-long int round( long double x );
-long long int round( float x );
-long long int round( double x );
-long long int round( long double x );
+float nearbyint( float );
+// extern "C" { double nearbyint( double ); }
+long double nearbyint( long double );
+
+float round( float );
+long double round( long double );
+long int round( float );
+long int round( double );
+long int round( long double );
+long long int round( float );
+long long int round( double );
+long long int round( long double );
 
 long int lround( float );
-extern "C" { long int lround( double ); }
+// extern "C" { long int lround( double ); }
 long int lround( long double );
 long long int llround( float );
-extern "C" { long long int llround( double ); }
+// extern "C" { long long int llround( double ); }
 long long int llround( long double );
 
-float log10( float );
-extern "C" { double log10( double ); }
-long double log10( long double );
-
-float log1p( float );
-extern "C" { double log1p( double ); }
-long double log1p( long double );
-
-float log2( float );
-extern "C" { double log2( double ); }
-long double log2( long double );
-
-float logb( float );
-extern "C" { double logb( double ); }
-long double logb( long double );
-
-float nearbyint( float );
-extern "C" { double nearbyint( double ); }
-long double nearbyint( long double );
+//---------------------- Manipulation ----------------------
+
+float copysign( float, float );
+// extern "C" { double copysign( double, double ); }
+long double copysign( long double, long double );
+
+float frexp( float, int * );
+// extern "C" { double frexp( double, int * ); }
+long double frexp( long double, int * );
+
+float ldexp( float, int );
+// extern "C" { double ldexp( double, int ); }
+long double ldexp( long double, int );
+
+[ float, float ] modf( float );
+float modf( float, float * );
+[ double, double ] modf( double );
+// extern "C" { double modf( double, double * ); }
+[ long double, long double ] modf( long double );
+long double modf( long double, long double * );
 
 float nextafter( float, float );
-extern "C" { double nextafter( double, double ); }
+// extern "C" { double nextafter( double, double ); }
 long double nextafter( long double, long double );
 
 float nexttoward( float, long double );
-extern "C" { double nexttoward( double, long double ); }
+// extern "C" { double nexttoward( double, long double ); }
 long double nexttoward( long double, long double );
 
-float remainder( float, float );
-extern "C" { double remainder( double, double ); }
-long double remainder( long double, long double );
-
-[ int, float ] remquo( float x, float y );
-float remquo( float, float, int * );
-[ int, double ] remquo( double x, double y );
-extern "C" { double remquo( double, double, int * ); }
-[ int, long double ] remquo( long double x, long double y );
-long double remquo( long double, long double, int * );
-
 float scalbn( float, int );
-extern "C" { double scalbn( double, int ); }
+// extern "C" { double scalbn( double, int ); }
 long double scalbn( long double, int );
 
 float scalbln( float, long int );
-extern "C" { double scalbln( double, long int ); }
+// extern "C" { double scalbln( double, long int ); }
 long double scalbln( long double, long int );
-
-float tgamma( float );
-extern "C" { double tgamma( double ); }
-long double tgamma( long double );
-
-float trunc( float );
-extern "C" { double trunc( double ); }
-long double trunc( long double );
 
 // Local Variables: //
Index: src/libcfa/math.c
===================================================================
--- src/libcfa/math.c	(revision 845cedcae61ca707834e349e69b07e41e766d7c0)
+++ src/libcfa/math.c	(revision add371c15f16e8ece211a386eec66e2e03a21afa)
@@ -1,3 +1,2 @@
-//                               -*- Mode: C -*- 
 // 
 // Cforall Version 1.0.0 Copyright (C) 2016 University of Waterloo
@@ -11,16 +10,121 @@
 // Created On       : Tue Apr 19 22:23:08 2016
 // Last Modified By : Peter A. Buhr
-// Last Modified On : Thu Apr 21 07:56:49 2016
-// Update Count     : 51
+// Last Modified On : Sun Apr 24 08:52:31 2016
+// Update Count     : 75
 // 
 
 #include "math"
 
-//---------------------------------------
-
 extern "C" {
-#include <math.h>
 #include <complex.h>
 } // extern "C"
+
+float fabs( float x ) { return fabsf( x ); }
+long double fabs( long double x ) { return fabsl( x ); }
+float cabs( float _Complex x ) { return cabsf( x ); }
+double cabs( double _Complex x ) { return cabs( x ); }
+long double cabs( long double _Complex x ) { return cabsl( x ); }
+
+float ?%?( float x, float y ) { return fmodf( x, y ); }
+float fmod( float x, float y ) { return fmodf( x, y ); }
+double ?%?( double x, double y ) { return fmod( x, y ); }
+long double ?%?( long double x, long double y ) { return fmodl( x, y ); }
+long double fmod( long double x, long double y ) { return fmodl( x, y ); }
+
+float remainder( float x, float y ) { return remainderf( x, y ); }
+long double remainder( long double x, long double y ) { return remainderl( x, y ); }
+
+// [ int, float ] remquo( float x, float y ) { int quo; x = remquof( x, y, &quo ); return [ quo, x ]; }
+float remquo( float x, float y, int *quo ) { return remquof( x, y, quo ); }
+// [ int, double ] remquo( double x, double y ) { int quo; x = remquo( x, y, &quo ); return [ quo, x ]; }
+// [ int, long double ] remquo( long double x, long double y ) { int quo; x = remquol( x, y, &quo ); return [ quo, x ]; }
+long double remquo( long double x, long double y, int *quo ) { return remquol( x, y, quo ); }
+
+// [ int, float ] div( float x, float y ) { int quo; x = remquof( x, y, &quo ); return [ quo, x ]; }
+float div( float x, float y, int *quo ) { return remquof( x, y, quo ); }
+// [ int, double ] div( double x, double y ) { int quo; x = remquo( x, y, &quo ); return [ quo, x ]; }
+// [ int, long double ] div( long double x, long double y ) { int quo; x = remquol( x, y, &quo ); return [ quo, x ]; }
+long double div( long double x, long double y, int *quo ) { return remquol( x, y, quo ); }
+
+float fma( float x, float y, float z ) { return fmaf( x, y, z ); }
+long double fma( long double x, long double y, long double z ) { return fmal( x, y, z ); }
+
+float fdim( float x, float y ) { return fdimf( x, y ); }
+long double fdim( long double x, long double y ) { return fdiml( x, y ); }
+
+float nan( const char *tag ) { return nanf( tag ); }
+long double nan( const char *tag ) { return nanl( tag ); }
+
+//---------------------- Exponential ----------------------
+
+float exp( float x ) { return expf( x ); }
+long double exp( long double x ) { return expl( x ); }
+float _Complex exp( float _Complex x ) { return cexpf( x ); }
+double _Complex exp( double _Complex x ) { return cexp( x ); }
+long double _Complex exp( long double _Complex x ) { return cexpl( x ); }
+
+float exp2( float x ) { return exp2f( x ); }
+long double exp2( long double x ) { return exp2l( x ); }
+// float _Complex exp2( float _Complex x ) { return cexp2f( x ); }
+// double _Complex exp2( double _Complex x ) { return cexp2( x ); }
+// long double _Complex exp2( long double _Complex x ) { return cexp2l( x ); }
+
+float expm1( float x ) { return expm1f( x ); }
+long double expm1( long double x ) { return expm1l( x ); }
+
+float log( float x ) { return logf( x ); }
+long double log( long double x ) { return logl( x ); }
+float _Complex log( float _Complex x ) { return clogf( x ); }
+double _Complex log( double _Complex x ) { return clog( x ); }
+long double _Complex log( long double _Complex x ) { return clogl( x ); }
+
+float log2( float x ) { return log2f( x ); }
+long double log2( long double x ) { return log2l( x ); }
+// float _Complex log2( float _Complex x ) { return clog2f( x ); }
+// double _Complex log2( double _Complex x ) { return clog2( x ); }
+// long double _Complex log2( long double _Complex x ) { return clog2l( x ); }
+
+float log10( float x ) { return log10f( x ); }
+long double log10( long double x ) { return log10l( x ); }
+// float _Complex log10( float _Complex x ) { return clog10f( x ); }
+// double _Complex log10( double _Complex x ) { return clog10( x ); }
+// long double _Complex log10( long double _Complex x ) { return clog10l( x ); }
+
+float log1p( float x ) { return log1pf( x ); }
+long double log1p( long double x ) { return log1pl( x ); }
+
+int ilogb( float x ) { return ilogbf( x ); }
+int ilogb( long double x ) { return ilogbl( x ); }
+
+float logb( float x ) { return logbf( x ); }
+long double logb( long double x ) { return logbl( x ); }
+
+//---------------------- Power ----------------------
+
+float sqrt( float x ) { return sqrtf( x ); }
+long double sqrt( long double x ) { return sqrtl( x ); }
+float _Complex sqrt( float _Complex x ) { return csqrtf( x ); }
+double _Complex sqrt( double _Complex x ) { return csqrt( x ); }
+long double _Complex sqrt( long double _Complex x ) { return csqrtl( x ); }
+
+float cbrt( float x ) { return cbrtf( x ); }
+long double cbrt( long double x ) { return cbrtl( x ); }
+
+float hypot( float x, float y ) { return hypotf( x, y ); }
+long double hypot( long double x, long double y ) { return hypotl( x, y ); }
+
+float pow( float x, float y ) { return powf( x, y ); }
+long double pow( long double x, long double y ) { return powl( x, y ); }
+float _Complex pow( float _Complex x, float _Complex y ) { return cpowf( x, y ); }
+double _Complex pow( double _Complex x, double _Complex y ) { return cpow( x, y ); }
+long double _Complex pow( long double _Complex x, long double _Complex y ) { return cpowl( x, y ); }
+
+//---------------------- Trigonometric ----------------------
+
+float sin( float x ) { return sinf( x ); }
+long double sin( long double x ) { return sinl( x ); }
+float _Complex sin( float _Complex x ) { return csinf( x ); }
+double _Complex sin( double _Complex x ) { return csin( x ); }
+long double _Complex sin( long double _Complex x ) { return csinl( x ); }
 
 float cos( float x ) { return cosf( x ); }
@@ -30,4 +134,44 @@
 long double _Complex cos( long double _Complex x ) { return ccosl( x ); }
 
+float tan( float x ) { return tanf( x ); }
+long double tan( long double x ) { return tanl( x ); }
+float _Complex tan( float _Complex x ) { return ctanf( x ); }
+double _Complex tan( double _Complex x ) { return ctan( x ); }
+long double _Complex tan( long double _Complex x ) { return ctanl( x ); }
+
+float asin( float x ) { return asinf( x ); }
+long double asin( long double x ) { return asinl( x ); }
+float _Complex asin( float _Complex x ) { return casinf( x ); }
+double _Complex asin( double _Complex x ) { return casin( x ); }
+long double _Complex asin( long double _Complex x ) { return casinl( x ); }
+
+float acos( float x ) { return acosf( x ); }
+long double acos( long double x ) { return acosl( x ); }
+float _Complex acos( float _Complex x ) { return cacosf( x ); }
+double _Complex acos( double _Complex x ) { return cacos( x ); }
+long double _Complex acos( long double _Complex x ) { return cacosl( x ); }
+
+float atan( float x ) { return atanf( x ); }
+long double atan( long double x ) { return atanl( x ); }
+float _Complex atan( float _Complex x ) { return catanf( x ); }
+double _Complex atan( double _Complex x ) { return catan( x ); }
+long double _Complex atan( long double _Complex x ) { return catanl( x ); }
+
+float atan2( float x, float y ) { return atan2f( x, y ); }
+double atan2( double x, double y ) { return atan2( x, y ); }
+long double atan2( long double x, long double y ) { return atan2l( x, y ); }
+
+float atan( float x, float y ) { return atan2f( x, y ); }
+double atan( double x, double y ) { return atan2( x, y ); }
+long double atan( long double x, long double y ) { return atan2l( x, y ); }
+
+//---------------------- Hyperbolic ----------------------
+
+float sinh( float x ) { return sinhf( x ); }
+long double sinh( long double x ) { return sinhl( x ); }
+float _Complex sinh( float _Complex x ) { return csinhf( x ); }
+double _Complex sinh( double _Complex x ) { return csinh( x ); }
+long double _Complex sinh( long double _Complex x ) { return csinhl( x ); }
+
 float cosh( float x ) { return coshf( x ); }
 long double cosh( long double x ) { return coshl( x ); }
@@ -36,9 +180,15 @@
 long double _Complex cosh( long double _Complex x ) { return ccoshl( x ); }
 
-float acos( float x ) { return acosf( x ); }
-long double acos( long double x ) { return acosl( x ); }
-float _Complex acos( float _Complex x ) { return cacosf( x ); }
-double _Complex acos( double _Complex x ) { return cacos( x ); }
-long double _Complex acos( long double _Complex x ) { return cacosl( x ); }
+float tanh( float x ) { return tanhf( x ); }
+long double tanh( long double x ) { return tanhl( x ); }
+float _Complex tanh( float _Complex x ) { return ctanhf( x ); }
+double _Complex tanh( double _Complex x ) { return ctanh( x ); }
+long double _Complex tanh( long double _Complex x ) { return ctanhl( x ); }
+
+float asinh( float x ) { return asinhf( x ); }
+long double asinh( long double x ) { return asinhl( x ); }
+float _Complex asinh( float _Complex x ) { return casinhf( x ); }
+double _Complex asinh( double _Complex x ) { return casinh( x ); }
+long double _Complex asinh( long double _Complex x ) { return casinhl( x ); }
 
 float acosh( float x ) { return acoshf( x ); }
@@ -48,50 +198,4 @@
 long double _Complex acosh( long double _Complex x ) { return cacoshl( x ); }
 
-float sin( float x ) { return sinf( x ); }
-long double sin( long double x ) { return sinl( x ); }
-float _Complex sin( float _Complex x ) { return csinf( x ); }
-double _Complex sin( double _Complex x ) { return csin( x ); }
-long double _Complex sin( long double _Complex x ) { return csinl( x ); }
-
-float sinh( float x ) { return sinhf( x ); }
-long double sinh( long double x ) { return sinhl( x ); }
-float _Complex sinh( float _Complex x ) { return csinhf( x ); }
-double _Complex sinh( double _Complex x ) { return csinh( x ); }
-long double _Complex sinh( long double _Complex x ) { return csinhl( x ); }
-
-float asin( float x ) { return asinf( x ); }
-long double asin( long double x ) { return asinl( x ); }
-float _Complex asin( float _Complex x ) { return casinf( x ); }
-double _Complex asin( double _Complex x ) { return casin( x ); }
-long double _Complex asin( long double _Complex x ) { return casinl( x ); }
-
-float asinh( float x ) { return asinhf( x ); }
-long double asinh( long double x ) { return asinhl( x ); }
-float _Complex asinh( float _Complex x ) { return casinhf( x ); }
-double _Complex asinh( double _Complex x ) { return casinh( x ); }
-long double _Complex asinh( long double _Complex x ) { return casinhl( x ); }
-
-float tan( float x ) { return tanf( x ); }
-long double tan( long double x ) { return tanl( x ); }
-float _Complex tan( float _Complex x ) { return ctanf( x ); }
-double _Complex tan( double _Complex x ) { return ctan( x ); }
-long double _Complex tan( long double _Complex x ) { return ctanl( x ); }
-
-float tanh( float x ) { return tanhf( x ); }
-long double tanh( long double x ) { return tanhl( x ); }
-float _Complex tanh( float _Complex x ) { return ctanhf( x ); }
-double _Complex tanh( double _Complex x ) { return ctanh( x ); }
-long double _Complex tanh( long double _Complex x ) { return ctanhl( x ); }
-
-float atan( float x ) { return atanf( x ); }
-long double atan( long double x ) { return atanl( x ); }
-float _Complex atan( float _Complex x ) { return catanf( x ); }
-double _Complex atan( double _Complex x ) { return catan( x ); }
-long double _Complex atan( long double _Complex x ) { return catanl( x ); }
-
-float atan( float x, float y ) { return atan2f( x, y ); }
-double atan( double x, double y ) { return atan2( x, y ); }
-long double atan( long double x, long double y ) { return atan2l( x, y ); }
-
 float atanh( float x ) { return atanhf( x ); }
 long double atanh( long double x ) { return atanhl( x ); }
@@ -100,32 +204,28 @@
 long double _Complex atanh( long double _Complex x ) { return catanhl( x ); }
 
-float exp( float x ) { return expf( x ); }
-long double exp( long double x ) { return expl( x ); }
-float _Complex exp( float _Complex x ) { return cexpf( x ); }
-double _Complex exp( double _Complex x ) { return cexp( x ); }
-long double _Complex exp( long double _Complex x ) { return cexpl( x ); }
-
-float log( float x ) { return logf( x ); }
-long double log( long double x ) { return logl( x ); }
-float _Complex log( float _Complex x ) { return clogf( x ); }
-double _Complex log( double _Complex x ) { return clog( x ); }
-long double _Complex log( long double _Complex x ) { return clogl( x ); }
-
-float pow( float x, float y ) { return powf( x, y ); }
-long double pow( long double x, long double y ) { return powl( x, y ); }
-float _Complex pow( float _Complex x, float _Complex y ) { return cpowf( x, y ); }
-double _Complex pow( double _Complex x, double _Complex y ) { return cpow( x, y ); }
-long double _Complex pow( long double _Complex x, long double _Complex y ) { return cpowl( x, y ); }
-
-float sqrt( float x ) { return sqrtf( x ); }
-long double sqrt( long double x ) { return sqrtl( x ); }
-float _Complex sqrt( float _Complex x ) { return csqrtf( x ); }
-double _Complex sqrt( double _Complex x ) { return csqrt( x ); }
-long double _Complex sqrt( long double _Complex x ) { return csqrtl( x ); }
-
-float fabs( float x ) { return fabsf( x ); }
-long double fabs( long double x ) { return fabsl( x ); }
-float cabs( float _Complex x ) { return cabsf( x ); }
-long double cabs( long double _Complex x ) { return cabsl( x ); }
+//---------------------- Error / Gamma ----------------------
+
+float erf( float x ) { return erff( x ); }
+long double erf( long double x ) { return erfl( x ); }
+// float _Complex erf( float _Complex x ) { return crflf( x ); }
+// double _Complex erf( double _Complex x ) { return crfl( x ); }
+// long double _Complex erf( long double _Complex x ) { return crfll( x ); }
+
+float erfc( float x ) { return erfcf( x ); }
+long double erfc( long double x ) { return erfcl( x ); }
+// float _Complex erfc( float _Complex x ) { return cerfcf( x ); }
+// double _Complex erfc( double _Complex x ) { return cerfc( x ); }
+// long double _Complex erfc( long double _Complex x ) { return cerfcl( x ); }
+
+float lgamma( float x ) { return lgammaf( x ); }
+long double lgamma( long double x ) { return lgammal( x ); }
+float lgamma( float x, int *sign ) { return lgammaf_r( x, sign ); }
+double lgamma( double x, int *sign ) { return lgamma_r( x, sign ); }
+long double lgamma( long double x, int *sign ) { return lgammal_r( x, sign ); }
+
+float tgamma( float x ) { return tgammaf( x ); }
+long double tgamma( long double x ) { return tgammal( x ); }
+
+//---------------------- Nearest Integer ----------------------
 
 float floor( float x ) { return floorf( x ); }
@@ -135,53 +235,6 @@
 long double ceil( long double x ) { return ceill( x ); }
 
-//---------------------------------------
-
-float cbrt( float x ) { return cbrtf( x ); }
-long double cbrt( long double x ) { return cbrtl( x ); }
-
-float copysign( float x, float y ) { return copysignf( x, y ); }
-long double copysign( long double x, long double y ) { return copysignl( x, y ); }
-
-float erf( float x ) { return erff( x ); }
-long double erf( long double x ) { return erfl( x ); }
-
-float erfc( float x ) { return erfcf( x ); }
-long double erfc( long double x ) { return erfcl( x ); }
-
-float exp2( float x ) { return exp2f( x ); }
-long double exp2( long double x ) { return exp2l( x ); }
-
-float expm1( float x ) { return expm1f( x ); }
-long double expm1( long double x ) { return expm1l( x ); }
-
-float fdim( float x, float y ) { return fdimf( x, y ); }
-long double fdim( long double x, long double y ) { return fdiml( x, y ); }
-
-float fma( float x, float y, float z ) { return fmaf( x, y, z ); }
-long double fma( long double x, long double y, long double z ) { return fmal( x, y, z ); }
-
-float ?%?( float x, float y ) { return fmodf( x, y ); }
-float fmod( float x, float y ) { return fmodf( x, y ); }
-double ?%?( double x, double y ) { return fmod( x, y ); }
-long double ?%?( long double x, long double y ) { return fmodl( x, y ); }
-long double fmod( long double x, long double y ) { return fmodl( x, y ); }
-
-//[ int, float ] div( float x ) { int i; x = frexpf( x, &i ); return [ i, x ]; }
-float frexp( float x, int *ip ) { return frexpf( x, ip ); }
-//[ int, double ] div( double x ) { int i; x = frexp( x, &i ); return [ i, x ]; }
-//[ int, long double ] div( long double x ) { int i; x = frexpl( x, &i ); return [ i, x ]; }
-long double frexp( long double x, int *ip ) { return frexpl( x, ip ); }
-
-float hypot( float x, float y ) { return hypotf( x, y ); }
-long double hypot( long double x, long double y ) { return hypotl( x, y ); }
-
-int ilogb( float x ) { return ilogbf( x ); }
-int ilogb( long double x ) { return ilogbl( x ); }
-
-float ldexp( float x, int exp2 ) { return ldexpf( x, exp2 ); }
-long double ldexp( long double x, int exp2 ) { return ldexpl( x, exp2 ); }
-
-float lgamma( float x ) { return lgammaf( x ); }
-long double lgamma( long double x ) { return lgammal( x ); }
+float trunc( float x ) { return truncf( x ); }
+long double trunc( long double x ) { return truncl( x ); }
 
 float rint( float x ) { return rintf( x ); }
@@ -199,4 +252,7 @@
 long long int llrint( long double x ) { return llrintl( x ); }
 
+float nearbyint( float x ) { return nearbyintf( x ); }
+long double nearbyint( long double x ) { return nearbyintl( x ); }
+
 float round( float x ) { return roundf( x ); }
 long double round( long double x ) { return roundl( x ); }
@@ -213,18 +269,20 @@
 long long int llround( long double x ) { return llroundl( x ); }
 
-float log10( float x ) { return log10f( x ); }
-long double log10( long double x ) { return log10l( x ); }
-
-float log1p( float x ) { return log1pf( x ); }
-long double log1p( long double x ) { return log1pl( x ); }
-
-float log2( float x ) { return log2f( x ); }
-long double log2( long double x ) { return log2l( x ); }
-
-float logb( float x ) { return logbf( x ); }
-long double logb( long double x ) { return logbl( x ); }
-
-float nearbyint( float x ) { return nearbyintf( x ); }
-long double nearbyint( long double x ) { return nearbyintl( x ); }
+//---------------------- Manipulation ----------------------
+
+float copysign( float x, float y ) { return copysignf( x, y ); }
+long double copysign( long double x, long double y ) { return copysignl( x, y ); }
+
+float frexp( float x, int *ip ) { return frexpf( x, ip ); }
+long double frexp( long double x, int *ip ) { return frexpl( x, ip ); }
+
+float ldexp( float x, int exp2 ) { return ldexpf( x, exp2 ); }
+long double ldexp( long double x, int exp2 ) { return ldexpl( x, exp2 ); }
+
+// [ float, float ] modf( float x ) { float i; x = modff( x, &i ); return [ i, x ]; }
+float modf( float x, float *i ) { return modff( x, i ); }
+// [ double, double ] modf( double x ) { double i; x = modf( x, &i ); return [ i, x ]; }
+// [ long double, long double ] modf( long double x ) { long double i; x = modfl( x, &i ); return [ i, x ]; }
+long double modf( long double x, long double *i ) { return modfl( x, i ); }
 
 float nextafter( float x, float y ) { return nextafterf( x, y ); }
@@ -233,13 +291,4 @@
 float nexttoward( float x, long double y ) { return nexttowardf( x, y ); }
 long double nexttoward( long double x, long double y ) { return nexttowardl( x, y ); }
-
-float remainder( float x, float y ) { return remainderf( x, y ); }
-long double remainder( long double x, long double y ) { return remainderl( x, y ); }
-
-//[ int, float ] remquo( float x, float y ) { int quo; x = remquof( x, y, &quo ); return [ quo, x ]; }
-float remquo( float x, float y, int *quo ) { return remquof( x, y, quo ); }
-//[ int, double ] remquo( double x, double y ) { int quo; x = remquo( x, y, &quo ); return [ quo, x ]; }
-//[ int, long double ] remquo( long double x, long double y ) { int quo; x = remquol( x, y, &quo ); return [ quo, x ]; }
-long double remquo( long double x, long double y, int *quo ) { return remquol( x, y, quo ); }
 
 float scalbn( float x, int exp ) { return scalbnf( x, exp ); }
@@ -252,10 +301,4 @@
 long double scalbln( long double x, long int exp ) { return scalblnl( x, exp ); }
 
-float tgamma( float x ) { return tgammaf( x ); }
-long double tgamma( long double x ) { return tgammal( x ); }
-
-float trunc( float x ) { return truncf( x ); }
-long double trunc( long double x ) { return truncl( x ); }
-
 // Local Variables: //
 // mode: c //