compiled newlib, need to hook up the system calls

rocklibc/src/math/isnan.c

@@ -0,0 +1,123 @@
+/* isnan( double )
+
+   This file is part of the Public Domain C Library (PDCLib).
+   Permission is granted to use, modify, and / or redistribute at will.
+*/
+
+#include <math.h>
+
+#ifndef REGTEST
+
+#include "pdclib/_PDCLIB_print.h"
+
+#include <string.h>
+
+int _PDCLIB_isnand( double x )
+{
+    _PDCLIB_bigint_digit_t data[ sizeof( double ) / sizeof( _PDCLIB_bigint_digit_t ) ];
+    size_t size;
+    int exp;
+
+    memcpy( data, &x, sizeof( double ) );
+    exp = _PDCLIB_DBL_EXP( data );
+    size = _PDCLIB_DBL_SIZE( data );
+
+    if ( exp == ( ( _PDCLIB_DBL_MAX_EXP - 1 ) + _PDCLIB_DBL_MAX_EXP ) )
+    {
+        _PDCLIB_bigint_digit_t digit = 0;
+        size_t i;
+
+        for ( i = 0; i < size; ++i )
+        {
+            digit |= data[i];
+        }
+
+        return digit != 0;
+    }
+
+    return 0;
+}
+
+int _PDCLIB_isnanf( float x )
+{
+    _PDCLIB_bigint_digit_t data[ sizeof( float ) / sizeof( _PDCLIB_bigint_digit_t ) ];
+    size_t size;
+    int exp;
+
+    memcpy( data, &x, sizeof( float ) );
+    exp = _PDCLIB_FLT_EXP( data );
+    size = _PDCLIB_FLT_SIZE( data );
+
+    if ( exp == ( ( _PDCLIB_FLT_MAX_EXP - 1 ) + _PDCLIB_FLT_MAX_EXP ) )
+    {
+        _PDCLIB_bigint_digit_t digit = 0;
+        size_t i;
+
+        for ( i = 0; i < size; ++i )
+        {
+            digit |= data[i];
+        }
+
+        return digit != 0;
+    }
+
+    return 0;
+}
+
+int _PDCLIB_isnanl( long double x )
+{
+    _PDCLIB_bigint_digit_t data[ sizeof( long double ) / sizeof( _PDCLIB_bigint_digit_t ) ];
+    size_t size;
+    int exp;
+
+    memcpy( data, &x, sizeof( long double ) );
+    exp = _PDCLIB_LDBL_EXP( data );
+    size = _PDCLIB_LDBL_SIZE( data );
+
+    if ( exp == ( ( _PDCLIB_LDBL_MAX_EXP - 1 ) + _PDCLIB_LDBL_MAX_EXP ) )
+    {
+        _PDCLIB_bigint_digit_t digit = 0;
+        size_t i;
+
+        for ( i = 0; i < size; ++i )
+        {
+            digit |= data[i];
+        }
+
+        return digit != 0;
+    }
+
+    return 0;
+}
+
+#endif
+
+#ifdef TEST
+
+#include "_PDCLIB_test.h"
+
+int main( void )
+{
+    float f;
+    double d;
+    long double ld;
+
+    f = 1e100;
+    TESTCASE( isnan( f ) == 0 );
+    f = 0.0 / 0.0;
+    TESTCASE( isnan( f ) == 1 );
+
+    d = 1e500;
+    TESTCASE( isnan( d ) == 0 );
+    d = 0.0 / 0.0;
+    TESTCASE( isnan( d ) == 1 );
+
+    ld = 1e500;
+    TESTCASE( isnan( ld ) == 0 );
+    ld = 0.0 / 0.0;
+    TESTCASE( isnan( ld ) == 1 );
+
+    return TEST_RESULTS;
+}
+
+#endif