testBuiltins.cpp File Reference

#include "header.h"
#include "DiagonalMsg.h"
#include "OneToAllMsg.h"
#include "../scheduling/Clock.h"
#include "Arith.h"
#include "TableBase.h"
#include "Table.h"
#include <queue>
#include "../shell/Shell.h"

Include dependency graph for testBuiltins.cpp:

Go to the source code of this file.

Functions
void	testArith ()
void	testBuiltins ()
void	testBuiltinsProcess ()
void	testFibonacci ()
void	testGetMsg ()
void	testMpiBuiltins ()
void	testMpiFibonacci ()
void	testNSDF ()
void	testStats ()
void	testTable ()
void	testUtilsForCompareXplot ()
void	testUtilsForLoadXplot ()

Function Documentation

void testArith

(

)

Definition at line 23 of file testBuiltins.cpp.

References Arith::arg1(), Arith::arg2(), Element::data(), Id::destroy(), Arith::getOutput(), Arith::initCinfo(), Id::nextId(), and Arith::process().

Referenced by testBuiltins().

{
    Id a1id = Id::nextId();
    unsigned int size = 10;
    Element* a1 = new GlobalDataElement(
                    a1id, Arith::initCinfo(), "a1", size );

    Eref a1_0( a1, 0 );
    Eref a1_1( a1, 1 );

    Arith* data1_0 = reinterpret_cast< Arith* >( a1->data( 0 ) );

    data1_0->arg1( 1 );
    data1_0->arg2( 0 );

    ProcInfo p;
    data1_0->process( a1_0, &p );

    assert( data1_0->getOutput() == 1 );

    data1_0->arg1( 1 );
    data1_0->arg2( 2 );

    data1_0->process( a1_0, &p );

    assert( data1_0->getOutput() == 3 );

    a1id.destroy();

    cout << "." << flush;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void testBuiltins

(

)

Definition at line 424 of file testBuiltins.cpp.

References testArith(), testNSDF(), and testTable().

Referenced by nonMpiTests().

{
    testArith();
    testTable();
        testNSDF();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void testBuiltinsProcess

(

)

Definition at line 431 of file testBuiltins.cpp.

References testGetMsg(), and testStats().

Referenced by processTests().

{
//  testFibonacci(); Nov 2013: Waiting till we have the MsgObjects fixed.
    testGetMsg();
    testStats();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void testFibonacci

(

)

This test uses the Diagonal Msg and summing in the Arith element to generate a Fibonacci series.

Definition at line 59 of file testBuiltins.cpp.

References Finfo::addMsg(), Arith::arg1(), Arith::arg2(), Eref::data(), Element::data(), Id::destroy(), Shell::doSetClock(), Shell::doStart(), Eref::element(), Id::eref(), Cinfo::findFinfo(), Arith::getOutput(), Arith::initCinfo(), Clock::initCinfo(), Msg::mid(), Id::nextId(), Shell::numNodes(), and DiagonalMsg::setStride().

{
    if ( Shell::numNodes() > 1 )
        return;
    unsigned int numFib = 20;
    Id a1id = Id::nextId();
    Element* a1 = new LocalDataElement(
                    a1id, Arith::initCinfo(), "a1", numFib );

    Arith* data = reinterpret_cast< Arith* >( a1->data( 0 ) );
    if ( data ) {
        data->arg1( 0 );
        data->arg2( 1 );
    }

    const Finfo* outFinfo = Arith::initCinfo()->findFinfo( "output" );
    const Finfo* arg1Finfo = Arith::initCinfo()->findFinfo( "arg1" );
    const Finfo* arg2Finfo = Arith::initCinfo()->findFinfo( "arg2" );
    const Finfo* procFinfo = Arith::initCinfo()->findFinfo( "process" );
    DiagonalMsg* dm1 = new DiagonalMsg( a1, a1, 0 );
    bool ret = outFinfo->addMsg( arg1Finfo, dm1->mid(), a1 );
    assert( ret );
    dm1->setStride( 1 );

    DiagonalMsg* dm2 = new DiagonalMsg( a1, a1, 0 );
    ret = outFinfo->addMsg( arg2Finfo, dm2->mid(), a1 );
    assert( ret );
    dm1->setStride( 2 );

    /*
    bool ret = DiagonalMsg::add( a1, "output", a1, "arg1", 1 );
    assert( ret );
    ret = DiagonalMsg::add( a1, "output", a1, "arg2", 2 );
    assert( ret );
    */

    Shell* shell = reinterpret_cast< Shell* >( Id().eref().data() );
    shell->doSetClock( 0, 1.0 );
    Eref clocker = Id( 1 ).eref();

    const Finfo* proc0Finfo = Clock::initCinfo()->findFinfo( "process0" );
    OneToAllMsg* otam = new OneToAllMsg( clocker, a1, 0 );
    ret = proc0Finfo->addMsg( procFinfo, otam->mid(), clocker.element() );

    assert( ret );

    shell->doStart( numFib );
    unsigned int f1 = 1;
    unsigned int f2 = 0;
    for ( unsigned int i = 0; i < numFib; ++i ) {
        // if ( a1->dataHandler()->isDataHere( i ) ) {
            Arith* data = reinterpret_cast< Arith* >( a1->data( i ) );
            // cout << Shell::myNode() << ": i = " << i << ", " << data->getOutput() << ", " << f1 << endl;
            assert( data->getOutput() == f1 );
        // }
        unsigned int temp = f1;
        f1 = temp + f2;
        f2 = temp;
    }

    a1id.destroy();
    cout << "." << flush;
}

Here is the call graph for this function:

void testGetMsg

(

)

Tests capacity to send a request for a field value to an object

Definition at line 272 of file testBuiltins.cpp.

References Eref::data(), Shell::doAddMsg(), Shell::doCopy(), Shell::doCreate(), Shell::doDelete(), Shell::doReinit(), Shell::doSetClock(), Shell::doStart(), doubleEq(), Shell::doUseClock(), ObjId::eref(), Id::eref(), Field< A >::get(), LookupField< L, A >::get(), numEntries, Eref::objId(), and SetGet1< A >::set().

Referenced by testBuiltinsProcess().

{
    Shell* shell = reinterpret_cast< Shell* >( Id().eref().data() );
    ObjId tabid = shell->doCreate( "Table", ObjId(), "tab", 1 );
    assert( tabid != ObjId() );
    ObjId arithid = shell->doCreate( "Arith", ObjId(), "arith", 1 );
    assert( arithid != ObjId() );
    // Table* t = reinterpret_cast< Table* >( tabid.eref().data() );
    ObjId ret = shell->doAddMsg( "Single",
        tabid.eref().objId(), "requestOut",
        arithid.eref().objId(), "getOutputValue" );
    assert( ret != ObjId() );
    ret = shell->doAddMsg( "Single", arithid.eref().objId(), "output",
        arithid.eref().objId(), "arg1" );
    assert( ret != ObjId() );
    shell->doSetClock( 0, 1 );
    shell->doSetClock( 1, 1 );
    shell->doUseClock( "/arith", "process", 0 );
    shell->doUseClock( "/tab", "process", 1 );
    unsigned int numEntries = Field< unsigned int >::get( tabid, "size" );
    assert( numEntries == 0 );
    Id clockId( 1 );
    // clockId.element()->digestMessages();
    // tabid.element()->digestMessages();
    // arithid.element()->digestMessages();
    shell->doReinit();
    numEntries = Field< unsigned int >::get( tabid, "size" );
    assert( numEntries == 1 ); // One for reinit call.
    SetGet1< double >::set( arithid, "arg1", 0.0 );
    SetGet1< double >::set( arithid, "arg2", 2.0 );
    shell->doStart( 100 );

    numEntries = Field< unsigned int >::get( tabid, "size" );
    assert( numEntries == 101 ); // One for reinit call, 100 for process.
    vector< double > temp = Field< vector< double > >::get( tabid, "vector" );

    for ( unsigned int i = 0; i < 100; ++i ) {
        double ret = LookupField< unsigned int, double >::get( tabid, "y", i );
        assert( doubleEq( ret, 2 * i ) );
        assert( doubleEq( temp[i], 2 * i ) );
    }

    // Here we check using the 'get' message with multiple targets
    Id arith2 = shell->doCopy( arithid, ObjId(), "arith2", 1, false, false);
    shell->doUseClock( "/arith2", "process", 0 );
    ret = shell->doAddMsg( "Single",
        tabid.eref().objId(), "requestOut",
        arith2.eref().objId(), "getOutputValue" );
    shell->doReinit();
    SetGet1< double >::set( arithid, "arg1", 0.0 );
    SetGet1< double >::set( arithid, "arg2", 2.0 );
    SetGet1< double >::set( arith2, "arg1", 10.0 );
    SetGet1< double >::set( arith2, "arg2", 12.0 );
    shell->doStart( 100 );

    numEntries = Field< unsigned int >::get( tabid, "size" );
    // One for reinit call, 100 for process, and there are two targets.
    assert( numEntries == 202 );
    temp = Field< vector< double > >::get( tabid, "vector" );

    for ( unsigned int i = 1; i < 100; ++i ) {
        assert( doubleEq( temp[2 * i], 2 * i ) );
        assert( doubleEq( temp[2 * i + 1], 10 + 12 * i ) );
    }


    // Perhaps I should do another test without reinit.
    /*
    SetGet2< string, string >::set(
        tabid.eref(), "xplot", "testfile", "testplot" );
    tabentry.destroy();
        */
    shell->doDelete( arithid );
    shell->doDelete( arith2 );
    shell->doDelete( tabid );
    cout << "." << flush;

}

Here is the call graph for this function:

Here is the caller graph for this function:

void testMpiBuiltins

(

)

Definition at line 438 of file testBuiltins.cpp.

Referenced by mpiTests().

{
    // testMpiFibonacci();
    cout << "." << flush;
}

Here is the caller graph for this function:

void testMpiFibonacci

(

)

This test uses the Diagonal Msg and summing in the Arith element to generate a Fibonacci series.

Definition at line 127 of file testBuiltins.cpp.

References Eref::data(), Id::destroy(), Shell::doAddMsg(), Shell::doCreate(), Shell::doSetClock(), Shell::doStart(), doubleEq(), Shell::doUseClock(), Id::eref(), Field< A >::getVec(), SetGet1< A >::set(), and Field< A >::set().

{
    unsigned int numFib = 20;

    // Id a1id = Id::nextId();
    Shell* shell = reinterpret_cast< Shell* >( Id().eref().data() );

    Id a1id = shell->doCreate( "Arith", ObjId(), "a1", numFib );
    SetGet1< double >::set( a1id, "arg1", 0 );
    SetGet1< double >::set( a1id, "arg2", 1 );

    /*
    Arith* data = reinterpret_cast< Arith* >( a1->dataHandler()->data( 0 ) );

    if ( data ) {
        data->arg1( 0 );
        data->arg2( 1 );
    }
    */

    ObjId mid1 = shell->doAddMsg( "Diagonal",
        ObjId( a1id, 0 ), "output", ObjId( a1id, 0 ), "arg1" );
    bool ret = Field< int >::set( mid1, "stride", 1 );
    assert( ret );

    ObjId mid2 = shell->doAddMsg( "Diagonal",
        ObjId( a1id, 0 ), "output", ObjId( a1id, 0 ), "arg2" );
    ret = Field< int >::set( mid2, "stride", 2 );
    assert( ret );

    /*
    bool ret = DiagonalMsg::add( a1, "output", a1, "arg1", 1 );
    assert( ret );
    ret = DiagonalMsg::add( a1, "output", a1, "arg2", 2 );
    assert( ret );
    */

    shell->doSetClock( 0, 1.0 );
    shell->doUseClock( "/a1", "process", 0 );

    shell->doStart( numFib );

    vector< double > retVec;
    Field< double >::getVec( a1id, "outputValue", retVec );
    assert( retVec.size() == numFib );

    unsigned int f1 = 1;
    unsigned int f2 = 0;
    for ( unsigned int i = 0; i < numFib; ++i ) {
        /*
        if ( a1->dataHandler()->isDataHere( i ) ) {
            Arith* data = reinterpret_cast< Arith* >( a1->dataHandler()->data( i ) );
            // cout << Shell::myNode() << ": i = " << i << ", " << data->getOutput() << ", " << f1 << endl;
            assert( data->getOutput() == f1 );
        }
        */
        assert( doubleEq( retVec[i], f1 ) );
        unsigned int temp = f1;
        f1 = temp + f2;
        f2 = temp;
    }

    a1id.destroy();
    cout << "." << flush;
}

Here is the call graph for this function:

void testNSDF

(

)

Definition at line 97 of file testNSDF.cpp.

References testCreateStringDataset().

Referenced by testBuiltins().

{
    testCreateStringDataset();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void testStats

(

)

Definition at line 353 of file testBuiltins.cpp.

References Eref::data(), Shell::doAddMsg(), Shell::doCreate(), Shell::doDelete(), Shell::doReinit(), Shell::doSetClock(), Shell::doStart(), doubleApprox(), Shell::doUseClock(), Id::eref(), Field< A >::get(), PI, and Field< A >::set().

Referenced by testBuiltinsProcess().

{
    Shell* shell = reinterpret_cast< Shell* >( Id().eref().data() );
    unsigned int size = 1000;
    vector< double > sinewave( size, 0 );
    for ( unsigned int i = 0; i < size; ++i ) {
        sinewave[i] = sin( 2 * i * PI / static_cast< double >( size ) );
    }
    ObjId tabid = shell->doCreate( "StimulusTable", ObjId(), "tab", 1 );
    assert( tabid != ObjId() );
    Field< vector< double > >::set( tabid, "vector", sinewave );
    Field< double >::set( tabid, "startTime", 0 );
    Field< double >::set( tabid, "stopTime", size );
    Field< double >::set( tabid, "loopTime", size );
    Field< double >::set( tabid, "stepSize", 0 );
    Field< double >::set( tabid, "stepPosition", 0 );
    Field< bool >::set( tabid, "doLoop", true );

    ObjId stat1 = shell->doCreate( "Stats", ObjId(), "stat1", 1 );
    Field< unsigned int >::set( stat1, "windowLength", size );

    ObjId mid = shell->doAddMsg( "Single", tabid, "output",
                    stat1, "input" );

    shell->doUseClock( "/tab", "process", 0 );
    shell->doUseClock( "/stat1", "process", 1 );
    shell->doSetClock( 0, 1 );
    shell->doSetClock( 1, 1 );

    shell->doReinit();
    shell->doStart( size/2 );
    double mean = Field< double >::get( stat1, "mean" );
    double sdev = Field< double >::get( stat1, "sdev" );
    double sum = Field< double >::get( stat1, "sum" );
    unsigned int num = Field< unsigned int >::get( stat1, "num" );
    double wmean = Field< double >::get( stat1, "wmean" );
    double wsdev = Field< double >::get( stat1, "wsdev" );
    double wsum = Field< double >::get( stat1, "wsum" );
    unsigned int wnum = Field< unsigned int >::get( stat1, "wnum" );
    assert( doubleApprox( mean, 2.0 / PI ) );
    assert( doubleApprox( sdev, 0.3077627 ) );
    assert( doubleApprox( sum, size / PI ) );
    assert( doubleApprox( num, size/2 ) );
    assert( doubleApprox( wmean, 2.0 / PI ) );
    // assert( doubleApprox( wsdev, 1.0/sqrt( 2.0 ) ) );
    assert( doubleApprox( wsum, size / PI ) );
    assert( wnum == size/2 );

    shell->doStart( size );
    mean = Field< double >::get( stat1, "mean" );
    sdev = Field< double >::get( stat1, "sdev" );
    sum = Field< double >::get( stat1, "sum" );
    num = Field< unsigned int >::get( stat1, "num" );
    wmean = Field< double >::get( stat1, "wmean" );
    wsdev = Field< double >::get( stat1, "wsdev" );
    wsum = Field< double >::get( stat1, "wsum" );
    wnum = Field< unsigned int >::get( stat1, "wnum" );
    assert( doubleApprox( mean, 2.0 / (PI * 3.0) ) );
    assert( doubleApprox( sdev, 0.6745136 ) );
    assert( doubleApprox( sum, size / PI ) );
    assert( doubleApprox( num, 3 * size/2 ) );
    assert( doubleApprox( wmean, 0 ) );
    assert( doubleApprox( wsdev, 1.0/sqrt( 2.0 ) ) );
    assert( doubleApprox( wsum, 0 ) );
    assert( wnum == size );

    shell->doDelete( tabid );
    shell->doDelete( stat1 );
    cout << "." << flush;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void testTable

(

)

Definition at line 245 of file testBuiltins.cpp.

References Neutral::children(), Eref::data(), Shell::doCreate(), Shell::doDelete(), doubleEq(), ObjId::eref(), Id::eref(), LookupField< L, A >::get(), Table::input(), testUtilsForCompareXplot(), and testUtilsForLoadXplot().

Referenced by testBuiltins().

{
    testUtilsForLoadXplot();
    testUtilsForCompareXplot();
    Shell* shell = reinterpret_cast< Shell* >( Id().eref().data() );
    vector< Id > ret;
    Neutral::children( Id().eref(), ret );
    ObjId tabid = shell->doCreate( "Table", ObjId(), "tab", 1 );
    assert( tabid != ObjId() );
    Table* t = reinterpret_cast< Table* >( tabid.eref().data() );
    for ( unsigned int i = 0; i < 100; ++i ) {
        t->input( sqrt((double) i ) );
    }
    vector< double > values = Field< vector< double > >::get( tabid, "vector");
    assert( values.size() == 100 );
    for ( unsigned int i = 0; i < 100; ++i ) {
        double ret = LookupField< unsigned int, double >::get( tabid, "y", i );
        assert( doubleEq( values[i] , sqrt((double) i ) ) );
        assert( doubleEq( ret , sqrt((double) i ) ) );
    }
    shell->doDelete( tabid );
    cout << "." << flush;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void testUtilsForCompareXplot

(

)

Definition at line 213 of file testBuiltins.cpp.

References doubleEq(), getRMS(), getRMSDiff(), and getRMSRatio().

Referenced by testTable().

{
    double getRMSDiff( const vector< double >& v1, const vector< double >& v2 );
    double getRMS( const vector< double >& v );

    double getRMSRatio( const vector< double >& v1, const vector< double >& v2 );

    vector< double > v1;
    vector< double > v2;
    v1.push_back( 0 );
    v1.push_back( 1 );
    v1.push_back( 2 );

    v2.push_back( 0 );
    v2.push_back( 1 );
    v2.push_back( 2 );

    double r1 = sqrt( 5.0 / 3.0 );
    double r2 = sqrt( 1.0 / 3.0 );

    assert( doubleEq( getRMS( v1 ), r1 ) );
    assert( doubleEq( getRMS( v2 ), r1 ) );
    assert( doubleEq( getRMSDiff( v1, v2 ), 0 ) );
    assert( doubleEq( getRMSRatio( v1, v2 ), 0 ) );

    v2[2] = 3;
    assert( doubleEq( getRMS( v2 ), sqrt( 10.0/3.0 ) ) );
    assert( doubleEq( getRMSDiff( v1, v2 ), r2 ) );
    assert( doubleEq( getRMSRatio( v1, v2 ), r2 / ( sqrt( 10.0/3.0 ) + r1 ) ) );
    cout << "." << flush;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void testUtilsForLoadXplot

(

)

Definition at line 193 of file testBuiltins.cpp.

References doubleEq(), getYcolumn(), and isNamedPlot().

Referenced by testTable().

{
    bool isNamedPlot( const string& line, const string& plotname );
    double getYcolumn( const string& line );

    assert( isNamedPlot( "/plotname foo", "foo" ) );
    assert( !isNamedPlot( "/plotname foo", "bar" ) );
    assert( !isNamedPlot( "/newplot", "bar" ) );
    assert( !isNamedPlot( "", "bar" ) );
    assert( !isNamedPlot( "1234.56", "bar" ) );

    assert( doubleEq( getYcolumn( "123.456" ), 123.456 ) );
    assert( doubleEq( getYcolumn( "987  123.456" ), 123.456 ) );
    assert( doubleEq( getYcolumn( "987 23.456" ), 23.456 ) );
    assert( doubleEq( getYcolumn( "987   3.456" ), 3.456 ) );
    assert( doubleEq( getYcolumn( "987   0.456" ), 0.456 ) );
    assert( doubleEq( getYcolumn( "987.6     0.456  1111.1" ), 987.6 ) );
    cout << "." << flush;
}

Here is the call graph for this function:

Here is the caller graph for this function: