#include "number.h"
#include "packedintarray.h"
#include <math.h>
#include <sys/time.h>

#include <bu/sio.h>
using namespace Bu;

void packedtest1()
{
    println("-==- Packed Int Test -==-");

    PackedIntArray a(4);
    a.append( 3 );
    a.append( 9 );
    a.append( 5 );

    println("%1").arg( a.toString() );
    println("%1").arg( a.toBitString() );
    println("%1").arg( PackedIntArray(4, 10).toString() );

    PackedIntArray b(5);
    for( int j = 0; j < 16; j++ )
    {
        b.append( 21 );
        if( b[j] != 21 )
        {
            println("Error at position %1").arg( j );
            println("Raw: %1 (%2)").arg( b.toBitString() ).arg( b.toString() );
        }
    }
}

void numbertest1()
{
    println("-==- Number test -==-");

    Number a("1000902491523000321");
    Number b("3004392012498000700");

    println("%1 + %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a + b).toString() );
    
    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );

    a = "-872";
    b = "123";

    println("%1 + %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a + b).toString() );

    a = "728";
    b = "-51";

    println("%1 + %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a + b).toString() );

    a = "44";
    b = "-55";

    println("%1 + %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a + b).toString() );
    
    a = "44";
    b = "-66";

    println("%1 + %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a + b).toString() );
    
    a = "44";
    b = "-66";

    println("%1 - %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a - b).toString() );
    
    a = "44";
    b = "66";

    println("%1 - %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a - b).toString() );

    a = "7814";
    b = "24";

    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );
    
    a = "12345";
    b = "678";

    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );
 
    a = "3592846";
    b = "944634757";

    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );
    
    a = "3592846";
    b = "";

    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );
    
    a = "123";
    b = "-45";

    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );
    
    a = "-123";
    b = "45";

    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );
    
    a = "-123";
    b = "-45";

    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );
    
    a = "123";
    b = "45";

    println("%1 / %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a / b).toString() );
    
    a = "12345";
    b = "45";

    println("%1 / %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a / b).toString() );

    a = "3007103450821050020096034077958224700";
    b = "898239467";

    println("%1 / %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a / b).toString() );
    println("%1 %% %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a % b).toString() );
    
    a = "983429807324875233421784598754987439873472349875329853298732";
    b = "18446744073709551615";

    println("%1 / %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a / b).toString() );
    println("%1 %% %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a % b).toString() );
}

#define compcheck( anum, op, bnum ) \
    a = #anum; b = #bnum; \
    println("%4: %1 " #op " %2 = %3").arg( a.toString() ).arg( b.toString() ). \
        arg( a op b ).arg( ((a op b) == (anum op bnum)) ? "pass" : "fail" )
    
void numbertestcomp()
{
    Number a, b;

    println("-==- Greater Than -==-");
    compcheck( 5, >, 10 );
    compcheck( 10, >, 5 );
    compcheck( 5, >, 5 );
    compcheck( 7, >, 5 );
    compcheck( 5, >, 7 );
    compcheck( 123, >, 122 );
    compcheck( 123, >, 123 );
    compcheck( 123, >, 124 );
    compcheck( -123, >, 122 );
    compcheck( -123, >, -122 );
    compcheck( -122, >, -123 );
    compcheck( 123, >, -122 );
   
    println("-==- Less Than -==-");
    compcheck( 5, <, 10 );
    compcheck( 10, <, 5 );
    compcheck( 5, <, 5 );
    compcheck( 7, <, 5 );
    compcheck( 5, <, 7 );
    compcheck( 123, <, 122 );
    compcheck( 123, <, 123 );
    compcheck( 123, <, 124 );
    compcheck( -123, <, 122 );
    compcheck( -123, <, -122 );
    compcheck( -122, <, -123 );
    compcheck( 123, <, -122 );
    
    println("-==- Greater Than or Equal To -==-");
    compcheck( 5, >=, 10 );
    compcheck( 10, >=, 5 );
    compcheck( 5, >=, 5 );
    compcheck( 7, >=, 5 );
    compcheck( 5, >=, 7 );
    compcheck( 123, >=, 122 );
    compcheck( 123, >=, 123 );
    compcheck( 123, >=, 124 );
    compcheck( -123, >=, 122 );
    compcheck( -123, >=, -122 );
    compcheck( -122, >=, -123 );
    compcheck( 123, >=, -122 );
   
    println("-==- Less Than or Equal To -==-");
    compcheck( 5, <=, 10 );
    compcheck( 10, <=, 5 );
    compcheck( 5, <=, 5 );
    compcheck( 7, <=, 5 );
    compcheck( 5, <=, 7 );
    compcheck( 123, <=, 122 );
    compcheck( 123, <=, 123 );
    compcheck( 123, <=, 124 );
    compcheck( -123, <=, 122 );
    compcheck( -123, <=, -122 );
    compcheck( -122, <=, -123 );
    compcheck( 123, <=, -122 );
}

int getHob( int x )
{
    for( int j = 31; j >= 0; j-- )
    {
        if( x&(1<<j) )
        {
            return j+1;
        }
    }
    return 0;
}

int getOrd( int x )
{
    return (int)(log(x)*0xb.8aa3b295c17fp-3+1.0);
}

double getTime()
{
    timeval t;
    gettimeofday( &t, NULL );
    return (double)t.tv_sec + (double)(t.tv_usec*.000001);
}

void ordertest()
{
    for( int j = 2; j < 64; j++ )
    {
        if( getOrd( j ) != getHob(j-1) )
            println("Failure on radix %1").arg( j );
    }

    double dStart, dEnd;

    dStart = getTime();
    for( int k = 0; k < 100000; k++ )
    {
        for( int j = 2; j < 32; j++ )
        {
            getOrd( j );
        }
    }
    dEnd = getTime();

    println("getOrd: %1 in %2 (%3 cps)").arg( 30*100000 ).arg( dEnd-dStart ).
        arg( (30*100000)/(dEnd-dStart) );
    
    dStart = getTime();
    for( int k = 0; k < 100000; k++ )
    {
        for( int j = 2; j < 32; j++ )
        {
            getHob( j-1 );
        }
    }
    dEnd = getTime();

    println("getHob: %1 in %2 (%3 cps)").arg( 30*100000 ).arg( dEnd-dStart ).
        arg( (30*100000)/(dEnd-dStart) );
}

void radixtest()
{
    Number a( 16 ), b( 16 );

    a = "f8a72bce3";
    b = "9ea8cb3";
    println("%1 + %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a + b).toString() );
    println("%1 - %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a - b).toString() );
    println("%1 / %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a / b).toString() );
    println("%1 * %2 = %3").
        arg( a.toString() ).
        arg( b.toString() ).
        arg( (a * b).toString() );
}

int main( int , char *[] )
{
    println("CliC");

//    packedtest1();
//    numbertest1();
//    numbertestcomp();
    radixtest();

//    ordertest();

    return 0;
}