Subversion Repositories openrisc_2011-10-31

/*

===============================================================================

This C source file is part of TestFloat, Release 2a, a package of programs

for testing the correctness of floating-point arithmetic complying to the

IEC/IEEE Standard for Floating-Point.

Written by John R. Hauser.  More information is available through the Web

page `http://HTTP.CS.Berkeley.EDU/~jhauser/arithmetic/TestFloat.html'.

THIS SOFTWARE IS DISTRIBUTED AS IS, FOR FREE.  Although reasonable effort

has been made to avoid it, THIS SOFTWARE MAY CONTAIN FAULTS THAT WILL AT

TIMES RESULT IN INCORRECT BEHAVIOR.  USE OF THIS SOFTWARE IS RESTRICTED TO

PERSONS AND ORGANIZATIONS WHO CAN AND WILL TAKE FULL RESPONSIBILITY FOR ANY

AND ALL LOSSES, COSTS, OR OTHER PROBLEMS ARISING FROM ITS USE.

Derivative works are acceptable, even for commercial purposes, so long as

(1) they include prominent notice that the work is derivative, and (2) they

include prominent notice akin to these four paragraphs for those parts of

this code that are retained.

Modified for use with or1ksim's testsuite.

Contributor Julius Baxter <julius.baxter@orsoc.se>

===============================================================================

*/

                                                                 /*

#include <stdlib.h>

#include <stdio.h>

                                                                 */

#include "support.h" // OR1k support C library

#include "milieu.h"

#include "softfloat.h"

#include "testCases.h"

#include "writeHex.h"

#include "testLoops.h"

volatile flag stop = FALSE;

char *trueName, *testName;

flag forever, errorStop;

uint32 maxErrorCount = 0;

flag checkNaNs = FALSE;

int8 *trueFlagsPtr;

int8 ( *testFlagsFunctionPtr )( void );

char *functionName;

char *roundingPrecisionName, *roundingModeName, *tininessModeName;

flag anyErrors = FALSE;

void

writeFunctionName( /* FILE *stream */ )

{

    printf( functionName/*, stream */ );

    if ( roundingModeName ) {

        if ( roundingPrecisionName ) {

            printf( ", precision "/*, stream */ );

            printf( roundingPrecisionName/*, stream */ );

        }

        printf( ", rounding "/*, stream */ );

        printf( roundingModeName/*, stream */ );

        if ( tininessModeName ) {

            printf( ", tininess "/*, stream */ );

            printf( tininessModeName/*, stream */ );

            printf( " rounding"/*, stream */ );

        }

    }

}

void

exitWithStatus( void )

{

  exit( anyErrors ? /*EXIT_FAILURE */ 0x1 : /* EXIT_SUCCESS */ 0 );

}

static uint32 tenthousandsCount, errorCount = 0;

static void

writeTestsTotal( void )

{

    if ( forever ) {

        printf( "Unbounded tests.\n"/*, stderr */ );

    }

    else {

        printf( "\r%d tests total.\n", testCases_total );

    }

}

static void

writeTestsPerformed( int16 count )

{

    if ( tenthousandsCount ) {

        printf(

             "\r%d%04d tests performed", tenthousandsCount, count );

    }

    else {

        printf(  "\r%d tests performed", count );

    }

    if ( errorCount ) {

        printf(

            "; %d error%s found.\n",

            errorCount,

            ( errorCount == 1 ) ? "" : "s"

        );

    }

    else {

        printf( ".\n"/*, stderr */ );

        printf( "No errors found in "/*, stdout*/ );

        writeFunctionName( /*stdout */ );

        printf( ".\n"/*, stdout*/ );

        //fflush( stdout );

    }

}

static void

checkEarlyExit( void )

{

    ++tenthousandsCount;

    if ( stop ) {

        writeTestsPerformed( 0 );

        exitWithStatus();

    }

    printf(  "\r%3d0000", tenthousandsCount );

}

static void

writeErrorFound( int16 count )

{

    putchar( '\r'/*, stderr */ );

    if ( errorCount == 1 ) {

        printf( "Errors found in "/*, stdout*/ );

        writeFunctionName( /*stdout*/ );

        printf( ":\n"/*, stdout*/ );

    }

    if ( stop ) {

        writeTestsPerformed( count );

        exitWithStatus();

    }

    anyErrors = TRUE;

}

INLINE void

writeInput_a_int32( void )

{

    writeHex_bits32( testCases_a_int32/*, stdout*/ );

}

#ifdef BITS64

INLINE void

writeInput_a_int64( void )

{

    writeHex_bits64( testCases_a_int64/*, stdout*/ );

}

#endif

INLINE void

writeInput_a_float32( void )

{

    writeHex_float32( testCases_a_float32/*, stdout*/ );

}

static void

writeInputs_ab_float32( void )

{

    writeHex_float32( testCases_a_float32/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    writeHex_float32( testCases_b_float32/*, stdout*/ );

}

INLINE void

writeInput_a_float64( void )

{

    writeHex_float64( testCases_a_float64/*, stdout*/ );

}

static void

writeInputs_ab_float64( void )

{

    writeHex_float64( testCases_a_float64/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    writeHex_float64( testCases_b_float64/*, stdout*/ );

}

#ifdef FLOATX80

INLINE void

writeInput_a_floatx80( void )

{

    writeHex_floatx80( testCases_a_floatx80/*, stdout*/ );

}

static void

writeInputs_ab_floatx80( void )

{

    writeHex_floatx80( testCases_a_floatx80/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    writeHex_floatx80( testCases_b_floatx80/*, stdout*/ );

}

#endif

#ifdef FLOAT128

INLINE void

writeInput_a_float128( void )

{

    writeHex_float128( testCases_a_float128/*, stdout*/ );

}

static void

writeInputs_ab_float128( void )

{

    writeHex_float128( testCases_a_float128/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    writeHex_float128( testCases_b_float128/*, stdout*/ );

}

#endif

static void

 writeOutputs_z_flag(

     flag trueZ, uint8 trueFlags, flag testZ, uint8 testFlags )

{

    printf( trueName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_flag( trueZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( trueFlags/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    printf( testName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_flag( testZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( testFlags/*, stdout*/ );

    putchar( '\n'/*, stdout*/ );

}

static void

 writeOutputs_z_int32(

     int32 trueZ, uint8 trueFlags, int32 testZ, uint8 testFlags )

{

    printf( trueName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_bits32( trueZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( trueFlags/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    printf( testName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_bits32( testZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( testFlags/*, stdout*/ );

    putchar( '\n'/*, stdout*/ );

}

#ifdef BITS64

static void

 writeOutputs_z_int64(

     int64 trueZ, uint8 trueFlags, int64 testZ, uint8 testFlags )

{

    printf( trueName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_bits64( trueZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( trueFlags/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    printf( testName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_bits64( testZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( testFlags/*, stdout*/ );

    putchar( '\n'/*, stdout*/ );

}

#endif

static void

 writeOutputs_z_float32(

     float32 trueZ, uint8 trueFlags, float32 testZ, uint8 testFlags )

{

    printf( trueName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_float32( trueZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( trueFlags/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    printf( testName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_float32( testZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( testFlags/*, stdout*/ );

    putchar( '\n'/*, stdout*/ );

}

static void

 writeOutputs_z_float64(

     float64 trueZ, uint8 trueFlags, float64 testZ, uint8 testFlags )

{

    printf( trueName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_float64( trueZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( trueFlags/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    printf( testName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_float64( testZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( testFlags/*, stdout*/ );

    putchar( '\n'/*, stdout*/ );

}

#ifdef FLOATX80

static void

 writeOutputs_z_floatx80(

     floatx80 trueZ, uint8 trueFlags, floatx80 testZ, uint8 testFlags )

{

    printf( trueName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_floatx80( trueZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( trueFlags/*, stdout*/ );

    printf( "  "/*, stdout*/ );

    printf( testName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_floatx80( testZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( testFlags/*, stdout*/ );

    putchar( '\n'/*, stdout*/ );

}

#endif

#ifdef FLOAT128

static void

 writeOutputs_z_float128(

     float128 trueZ, uint8 trueFlags, float128 testZ, uint8 testFlags )

{

    printf( trueName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_float128( trueZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( trueFlags/*, stdout*/ );

    printf( "\n\t"/*, stdout*/ );

    printf( testName/*, stdout*/ );

    printf( ": "/*, stdout*/ );

    writeHex_float128( testZ/*, stdout*/ );

    putchar( ' '/*, stdout*/ );

    writeHex_float_flags( testFlags/*, stdout*/ );

    putchar( '\n'/*, stdout*/ );

}

#endif

INLINE flag float32_isNaN( float32 a )

{

    return 0x7F800000 < ( a & 0x7FFFFFFF );

}

#ifdef BITS64

INLINE flag float64_same( float64 a, float64 b )

{

    return a == b;

}

INLINE flag float64_isNaN( float64 a )

{

    return LIT64( 0x7FF0000000000000 ) < ( a & LIT64( 0x7FFFFFFFFFFFFFFF ) );

}

#else

INLINE flag float64_same( float64 a, float64 b )

{

    return ( a.high == b.high ) && ( a.low == b.low );

}

INLINE flag float64_isNaN( float64 a )

{

    bits32 absAHigh;

    absAHigh = a.high & 0x7FFFFFFF;

    return

        ( 0x7FF00000 < absAHigh ) || ( ( absAHigh == 0x7FF00000 ) && a.low );

}

#endif

#ifdef FLOATX80

INLINE flag floatx80_same( floatx80 a, floatx80 b )

{

    return ( a.high == b.high ) && ( a.low == b.low );

}

INLINE flag floatx80_isNaN( floatx80 a )

{

    return ( ( a.high & 0x7FFF ) == 0x7FFF ) && a.low;

}

#endif

#ifdef FLOAT128

INLINE flag float128_same( float128 a, float128 b )

{

    return ( a.high == b.high ) && ( a.low == b.low );

}

INLINE flag float128_isNaN( float128 a )

{

    bits64 absAHigh;

    absAHigh = a.high & LIT64( 0x7FFFFFFFFFFFFFFF );

    return

           ( LIT64( 0x7FFF000000000000 ) < absAHigh )

        || ( ( absAHigh == LIT64( 0x7FFF000000000000 ) ) && a.low );

}

#endif

void

 test_a_int32_z_float32(

     float32 trueFunction( int32 ), float32 testFunction( int32 ) )

{

    int16 count;

    float32 trueZ, testZ;

    uint8 trueFlags, testFlags;

    errorCount = 0;

    tenthousandsCount = 0;

    count = 10000;

    testCases_initSequence( testCases_sequence_a_int32 );

    writeTestsTotal();

    while ( ! testCases_done || forever ) {

        testCases_next();

        *trueFlagsPtr = 0;

        trueZ = trueFunction( testCases_a_int32 );

        trueFlags = *trueFlagsPtr;

        (void) testFlagsFunctionPtr();

        testZ = testFunction( testCases_a_int32 );

        testFlags = testFlagsFunctionPtr();

        --count;

        if ( count == 0 ) {

            checkEarlyExit();

            count = 10000;

        }

        if ( ( trueZ != testZ ) || ( trueFlags != testFlags ) ) {

            if (    ! checkNaNs

                 && float32_isNaN( trueZ )

                 && float32_isNaN( testZ )

                 && ! float32_is_signaling_nan( testZ )

                 && ( trueFlags == testFlags )

               ) {

                /* no problem */

            }

            else {

                ++errorCount;

                writeErrorFound( 10000 - count );

                writeInput_a_int32();

                printf( "  "/*, stdout*/ );

                writeOutputs_z_float32( trueZ, trueFlags, testZ, testFlags );

                //fflush( stdout );

                if ( errorCount == maxErrorCount ) goto exit;

            }

        }

    }

 exit:

    writeTestsPerformed( 10000 - count );

}

void

 test_a_int32_z_float64(

     float64 trueFunction( int32 ), float64 testFunction( int32 ) )

{

    int16 count;

    float64 trueZ, testZ;

    uint8 trueFlags, testFlags;

    errorCount = 0;

    tenthousandsCount = 0;

    count = 10000;

    testCases_initSequence( testCases_sequence_a_int32 );

    writeTestsTotal();

    while ( ! testCases_done || forever ) {

        testCases_next();

        *trueFlagsPtr = 0;

        trueZ = trueFunction( testCases_a_int32 );

        trueFlags = *trueFlagsPtr;

        (void) testFlagsFunctionPtr();

        testZ = testFunction( testCases_a_int32 );

        testFlags = testFlagsFunctionPtr();

        --count;

        if ( count == 0 ) {

            checkEarlyExit();

            count = 10000;

        }

        if ( ! float64_same( trueZ, testZ ) || ( trueFlags != testFlags ) ) {

            if (    ! checkNaNs

                 && float64_isNaN( trueZ )

                 && float64_isNaN( testZ )

                 && ! float64_is_signaling_nan( testZ )

                 && ( trueFlags == testFlags )

               ) {

                /* no problem */

            }

            else {

                ++errorCount;

                writeErrorFound( 10000 - count );

                writeInput_a_int32();

                printf( "  "/*, stdout*/ );

                writeOutputs_z_float64( trueZ, trueFlags, testZ, testFlags );

                //fflush( stdout );

                if ( errorCount == maxErrorCount ) goto exit;

            }

        }

    }

 exit:

    writeTestsPerformed( 10000 - count );

}

#ifdef FLOATX80

void

 test_a_int32_z_floatx80(

     floatx80 trueFunction( int32 ), floatx80 testFunction( int32 ) )

{

    int16 count;

    floatx80 trueZ, testZ;

    uint8 trueFlags, testFlags;

    errorCount = 0;

    tenthousandsCount = 0;

    count = 10000;

    testCases_initSequence( testCases_sequence_a_int32 );

    writeTestsTotal();

    while ( ! testCases_done || forever ) {

        testCases_next();

        *trueFlagsPtr = 0;

        trueZ = trueFunction( testCases_a_int32 );

        trueFlags = *trueFlagsPtr;

        (void) testFlagsFunctionPtr();

        testZ = testFunction( testCases_a_int32 );

        testFlags = testFlagsFunctionPtr();

        --count;

        if ( count == 0 ) {

            checkEarlyExit();

            count = 10000;

        }

        if ( ! floatx80_same( trueZ, testZ ) || ( trueFlags != testFlags ) ) {

            if (    ! checkNaNs

                 && floatx80_isNaN( trueZ )

                 && floatx80_isNaN( testZ )

                 && ! floatx80_is_signaling_nan( testZ )

                 && ( trueFlags == testFlags )

               ) {

                /* no problem */

            }

            else {

                ++errorCount;

                writeErrorFound( 10000 - count );

                writeInput_a_int32();

                printf( "  "/*, stdout*/ );

                writeOutputs_z_floatx80( trueZ, trueFlags, testZ, testFlags );

                //fflush( stdout );

                if ( errorCount == maxErrorCount ) goto exit;

            }

        }

    }

 exit:

    writeTestsPerformed( 10000 - count );

}

#endif

#ifdef FLOAT128

void

 test_a_int32_z_float128(

     float128 trueFunction( int32 ), float128 testFunction( int32 ) )

{

    int16 count;

    float128 trueZ, testZ;

    uint8 trueFlags, testFlags;

    errorCount = 0;

    tenthousandsCount = 0;

    count = 10000;

    testCases_initSequence( testCases_sequence_a_int32 );

    writeTestsTotal();

    while ( ! testCases_done || forever ) {

        testCases_next();

        *trueFlagsPtr = 0;

        trueZ = trueFunction( testCases_a_int32 );

        trueFlags = *trueFlagsPtr;

        (void) testFlagsFunctionPtr();

        testZ = testFunction( testCases_a_int32 );

        testFlags = testFlagsFunctionPtr();

        --count;

        if ( count == 0 ) {

            checkEarlyExit();

            count = 10000;

        }

        if ( ! float128_same( trueZ, testZ ) || ( trueFlags != testFlags ) ) {

            if (    ! checkNaNs

                 && float128_isNaN( trueZ )

                 && float128_isNaN( testZ )

                 && ! float128_is_signaling_nan( testZ )

                 && ( trueFlags == testFlags )

               ) {

                /* no problem */

            }