MSNumpressTest.cpp File Reference

#include "MSNumpress.hpp"
#include <assert.h>
#include <iostream>
#include <cmath>
#include <cstdlib>
#include <stdio.h>

Include dependency graph for MSNumpressTest.cpp:

Go to the source code of this file.

Functions
void	encodeLinear1 ()
void	encodeLinear ()
void	decodeLinearNice ()
void	decodeLinearNiceLowFP ()
void	decodeLinearWierd ()
void	decodeLinearWierd_llong_overflow ()
void	decodeLinearWierd_int_overflow ()
void	decodeLinearWierd_int_underflow ()
void	decodeLinearCorrupt1 ()
void	decodeLinearCorrupt2 ()
void	optimalLinearFixedPoint ()
void	optimalLinearFixedPointMass ()
void	encodeDecodeLinearStraight ()
void	encodeDecodeSafeStraight ()
void	encodeDecodeSafe ()
void	encodeDecodeLinear ()
void	encodeDecodeLinear5 ()
void	encodeDecodePic ()
void	encodeDecodePic5 ()
void	optimalSlofFixedPoint ()
void	encodeDecodeSlof ()
void	encodeDecodeSlof5 ()
void	testErroneousDecodePic ()
int	main (int argc, const char *argv[])

Variables
double	ENC_TWO_BYTE_FIXED_POINT = 3000.0

Function Documentation

decodeLinearCorrupt1()

void decodeLinearCorrupt1

(

)

Definition at line 279 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeLinear().

Referenced by main().

                            {
    unsigned char encoded[20];
    double decoded[4];
    
    try {
        ms::numpress::MSNumpress::decodeLinear(&encoded[0], 11, &decoded[0]);
        cout << "- fail    decodeLinearCorrupt1: didn't throw exception for corrupt input " << endl << endl;
        assert(0 == 1);
    } catch (const char *err) {
        
    }
    
    try {
        ms::numpress::MSNumpress::decodeLinear(&encoded[0], 14, &decoded[0]);
        cout << "- fail    decodeLinearCorrupt1: didn't throw exception for corrupt input " << endl << endl;
        assert(0 == 1);
    } catch (const char *err) {
        
    }
    
    cout << "+ pass    decodeLinearCorrupt 1 " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

decodeLinearCorrupt2()

void decodeLinearCorrupt2

(

)

Definition at line 303 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeLinear(), ms::numpress::MSNumpress::encodeLinear(), and ms::numpress::MSNumpress::optimalLinearFixedPoint().

Referenced by main().

                            {
    
    double mzs[4];
    
    mzs[0] = 100.0;
    mzs[1] = 200.0;
    mzs[2] = 300.00005;
    mzs[3] = 0.00010;
    
    size_t nMzs = 4;
    unsigned char encoded[28];
    double fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPoint(&mzs[0], nMzs);
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], fixedPoint);
    
    double decoded[4];
    try {
        ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes-1, &decoded[0]);
        cout << "- fail    decodeLinearCorrupt2: didn't throw exception for corrupt input " << endl << endl;
        assert(0 == 1);
    } catch (const char *err) {
        
    }
    
    cout << "+ pass    decodeLinearCorrupt 2 " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

decodeLinearNice()

void decodeLinearNice

(

)

Definition at line 86 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeLinear(), ms::numpress::MSNumpress::encodeLinear(), and ms::numpress::MSNumpress::optimalLinearFixedPoint().

Referenced by main().

                        {
    
    double mzs[4];
    
    mzs[0] = 100.0;
    mzs[1] = 200.0;
    mzs[2] = 300.00005;
    mzs[3] = 400.00010;
    
    size_t nMzs = 4;
    unsigned char encoded[28];
    double fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPoint(&mzs[0], nMzs);
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], fixedPoint);
    
    double decoded[4];
    size_t numDecoded = ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes, &decoded[0]);
    assert(4 == numDecoded);
    assert(abs(100.0 - decoded[0]) < 0.000005);
    assert(abs(200.0 - decoded[1]) < 0.000005);
    assert(abs(300.00005 - decoded[2]) < 0.000005);
    assert(abs(400.00010 - decoded[3]) < 0.000005);
    
    cout << "+ pass    decodeLinearNice " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

decodeLinearNiceLowFP()

void decodeLinearNiceLowFP

(

)

Definition at line 111 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeLinear(), ms::numpress::MSNumpress::encodeLinear(), and ms::numpress::MSNumpress::optimalLinearFixedPointMass().

Referenced by main().

                             {
  
  double mzs[7];
  
  mzs[0] = 100.0;
  mzs[1] = 200.0;
  mzs[2] = 300.00005;
  mzs[3] = 400.00010;
  mzs[4] = 450.00010;
  mzs[5] = 455.00010;
  mzs[6] = 700.00010;
  
  size_t nMzs = 7;
  unsigned char encoded[33]; // max length is 33 bytes

  // check for fixed points
  {
    double fixedPoint;
    fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], nMzs, 0.1);
    assert( abs(5 - fixedPoint) < 0.000005);

    fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], nMzs, 1e-3);
    assert( abs(500 - fixedPoint) < 0.000005);

    fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], nMzs, 1e-5);
    assert( abs(50000 - fixedPoint) < 0.000005);

    fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], nMzs, 1e-7);
    assert( abs(5000000 - fixedPoint) < 0.000005);

    // cannot fulfill accuracy of 1e-8
    fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], nMzs, 1e-8);
    assert( abs(-1 - fixedPoint) < 0.000005);
  }

  {
    double fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], nMzs, 0.001);
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], fixedPoint);

    double decoded[7];
    size_t numDecoded = ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes, &decoded[0]);
    assert(25 == encodedBytes);
    assert(7 == numDecoded);

    assert(abs(100.0 - decoded[0]) < 0.001);
    assert(abs(200.0 - decoded[1]) < 0.001);
    assert(abs(300.00005 - decoded[2]) < 0.001);
    assert(abs(400.00010 - decoded[3]) < 0.001);
  }
  
  double mz_err[5];
  double encodedLength[5];

  // for higher accuracy, we get longer encoded lengths
  mz_err[0] = 0.1;  encodedLength[0] = 22;
  mz_err[1] = 1e-3; encodedLength[1] = 25;
  mz_err[2] = 1e-5; encodedLength[2] = 29;
  mz_err[3] = 1e-6; encodedLength[3] = 30;
  mz_err[4] = 1e-7; encodedLength[4] = 31;

  for (int k = 0; k < 4; k++)
  {
    double fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], nMzs, mz_err[k]);
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], fixedPoint);

    double decoded[7];
    size_t numDecoded = ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes, &decoded[0]);
    assert( encodedLength[k] == encodedBytes);
    assert(7 == numDecoded);
    assert(abs(100.0 - decoded[0]) < mz_err[k]);
    assert(abs(200.0 - decoded[1]) < mz_err[k]);
    assert(abs(300.00005 - decoded[2]) < mz_err[k]);
    assert(abs(400.00010 - decoded[3]) < mz_err[k]);
    assert(abs(450.00010 - decoded[4]) < mz_err[k]);
    assert(abs(455.00010 - decoded[5]) < mz_err[k]);
    assert(abs(700.00010 - decoded[6]) < mz_err[k]);
  }

  cout << "+ pass    decodeLinearNiceFP " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

decodeLinearWierd()

void decodeLinearWierd

(

)

Definition at line 192 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeLinear(), ms::numpress::MSNumpress::encodeLinear(), and ms::numpress::MSNumpress::optimalLinearFixedPoint().

Referenced by main().

                         {
    double mzs[4];
    
    mzs[0] = 100.0;
    mzs[1] = 200.0;
    mzs[2] = 300.00005;
    mzs[3] = 0.00010;
    
    size_t nMzs = 4;
    unsigned char encoded[28];
    double fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPoint(&mzs[0], nMzs);
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], fixedPoint);
    
    double decoded[4];
    size_t numDecoded = ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes, &decoded[0]);
    assert(4 == numDecoded);
    assert(abs(100.0 - decoded[0]) < 0.000005);
    assert(abs(200.0 - decoded[1]) < 0.000005);
    assert(abs(300.00005 - decoded[2]) < 0.000005);
    assert(abs(0.00010 - decoded[3]) < 0.000005);
    
    cout << "+ pass    decodeLinearWierd " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

decodeLinearWierd_int_overflow()

void decodeLinearWierd_int_overflow

(

)

Definition at line 237 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::encodeLinear().

                                      {
    double mzs[4];
    
    mzs[0] = 100.0;
    mzs[1] = 200.0;
    mzs[2] = 30000000.00005;
    mzs[3] = 0.00006;
    
    size_t nMzs = 4;
    unsigned char encoded[28];

    try {
        ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], 1000000000);
        cout << "- fail    test decodeLinearWierd_int_overflow: didn't throw exception for corrupt input " << endl << endl;
        assert(0 == 1);
    } catch (const char *err) {
        assert( std::string(err) == std::string("[MSNumpress::encodeLinear] Cannot encode a number that exceeds the bounds of [-INT_MAX, INT_MAX]."));
    }
    cout << "+ pass    decodeLinearWierd3 " << endl << endl;
}

Here is the call graph for this function:

decodeLinearWierd_int_underflow()

void decodeLinearWierd_int_underflow

(

)

Definition at line 258 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::encodeLinear().

                                       {
    double mzs[4];
    
    mzs[0] = 30000000.00005;
    mzs[1] = 60000000.00005;
    mzs[2] = 30000000.00005;
    mzs[3] = 0.00006;
    
    size_t nMzs = 4;
    unsigned char encoded[28];

    try {
        ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], 1000000000);
        cout << "- fail    test decodeLinearWierd_int_underflow: didn't throw exception for corrupt input " << endl << endl;
        assert(0 == 1);
    } catch (const char *err) {
        assert( std::string(err) == std::string("[MSNumpress::encodeLinear] Cannot encode a number that exceeds the bounds of [-INT_MAX, INT_MAX]."));
    }
    cout << "+ pass    decodeLinearWierd3 " << endl << endl;
}

Here is the call graph for this function:

decodeLinearWierd_llong_overflow()

void decodeLinearWierd_llong_overflow

(

)

Definition at line 216 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::encodeLinear().

                                        {
    double mzs[4];
    
    mzs[0] = 100.0;
    mzs[1] = 200.0;
    mzs[2] = 30000000.00005;
    mzs[3] = 0.000010;
    
    size_t nMzs = 4;
    unsigned char encoded[28];

    try {
        ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], 1000000000000);
        cout << "- fail    test decodeLinearWierd_llong_overflow: didn't throw exception for corrupt input " << endl << endl;
        assert(0 == 1);
    } catch (const char *err) {
        assert( std::string(err) == std::string("[MSNumpress::encodeLinear] Next number overflows LLONG_MAX."));
    }
    cout << "+ pass    decodeLinearWierd_llong_overflow " << endl << endl;
}

Here is the call graph for this function:

encodeDecodeLinear()

void encodeDecodeLinear

(

)

Definition at line 475 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeLinear(), ms::numpress::MSNumpress::encodeLinear(), and ms::numpress::MSNumpress::optimalLinearFixedPoint().

Referenced by main().

                          {
    srand(123459);
    
    size_t n = 1000;
    double mzs[1000];
    mzs[0] = 300 + rand() / double(RAND_MAX);
    for (size_t i=1; i<n; i++) 
        mzs[i] = mzs[i-1] + rand() / double(RAND_MAX);
    
    double fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPoint(&mzs[0], n);
    
    unsigned char encoded[5000];
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], n, &encoded[0], fixedPoint);
    
    double decoded[1000];
    size_t numDecoded = ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes, &decoded[0]);
    
    assert(n == numDecoded);
    
    double m = 0;
    double error;
    double mLim = 0.000005;
    
    for (size_t i=0; i<n; i++) { 
        error = abs(mzs[i] - decoded[i]);
        m = max(m, error);
        if (error >= mLim) {
            cout << "error   " << error << " above limit " << mLim << endl;
            assert(error < mLim);
        }
    }
    cout << "+     size compressed: " << encodedBytes / double(n*8) * 100 << "% " << endl;
    cout << "+           max error: " << m << "  limit: " << mLim << endl;
    cout << "+ pass    encodeDecodeLinear " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeDecodeLinear5()

void encodeDecodeLinear5

(

)

Definition at line 513 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeLinear(), ms::numpress::MSNumpress::encodeLinear(), and ms::numpress::MSNumpress::optimalLinearFixedPoint().

Referenced by main().

                           {
    srand(123662);
    
    size_t n = 1000;
    double mzs[1000];
    mzs[0] = 100 + (rand() % 1000) / 1000.0;
    for (size_t i=1; i<n; i++) 
        mzs[i] = mzs[i-1] + (rand() % 1000) / 1000.0;
        
    size_t encodedBytes;
    size_t numDecoded;
    unsigned char encoded[5000];
    double firstDecoded[1000];
    double decoded[1000];
    
    double fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPoint(&mzs[0], n);
    
    encodedBytes    = ms::numpress::MSNumpress::encodeLinear(&mzs[0], n, &encoded[0], fixedPoint);
    numDecoded      = ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes, &decoded[0]);
    
    for (size_t i=0; i<n; i++) 
        firstDecoded[i] = decoded[i];
    
    for (size_t i=0; i<5; i++) {
        encodedBytes    = ms::numpress::MSNumpress::encodeLinear(&decoded[0], n, &encoded[0], fixedPoint);
        numDecoded      = ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes, &decoded[0]);
    }
    
    double afterFixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPoint(&decoded[0], n);
    
    assert(fixedPoint == afterFixedPoint);
    assert(n == numDecoded);
        
    for (size_t i=0; i<n; i++)
        if (firstDecoded[i] != decoded[i]) {
            cout << endl << firstDecoded[i] << " " << decoded[i] << endl;
            assert(firstDecoded[i] == decoded[i]);
        }
    
    cout << "+ pass    encodeDecodeLinear5 " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeDecodeLinearStraight()

void encodeDecodeLinearStraight

(

)

Definition at line 370 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeLinear(), ms::numpress::MSNumpress::encodeLinear(), and ms::numpress::MSNumpress::optimalLinearFixedPoint().

Referenced by main().

                                  {
    size_t n = 15;
    double mzs[15];
    for (size_t i=0; i<n; i++) 
        mzs[i] = i;
    
    double fixedPoint = ms::numpress::MSNumpress::optimalLinearFixedPoint(&mzs[0], n);
    
    unsigned char encoded[75];
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], n, &encoded[0], fixedPoint);
    
    double decoded[15];
    size_t numDecoded = ms::numpress::MSNumpress::decodeLinear(&encoded[0], encodedBytes, &decoded[0]);
    
    assert(n == numDecoded);
    
    double m = 0;
    double error;
    double mLim = 0.000005;
    
    for (size_t i=0; i<n; i++) { 
        error = abs(mzs[i] - decoded[i]);
        m = max(m, error);
        if (error >= mLim) {
            cout << "error   " << error << " above limit " << mLim << endl;
            assert(error < mLim);
        }
    }
    cout << "+     size compressed: " << encodedBytes / double(n*8) * 100 << "% " << endl;
    cout << "+           max error: " << m << "  limit: " << mLim << endl;
    cout << "+ pass    encodeDecodeLinearStraight " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeDecodePic()

void encodeDecodePic

(

)

Definition at line 557 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodePic(), and ms::numpress::MSNumpress::encodePic().

Referenced by main().

                       {
    srand(123459);
    
    size_t n = 1000;
    double ics[1000];
    for (size_t i=0; i<n; i++) 
        ics[i] = rand() % 1000000;
    
    unsigned char encoded[5000];
    size_t encodedBytes = ms::numpress::MSNumpress::encodePic(&ics[0], n, &encoded[0]);
    
    double decoded[1000];
    size_t numDecoded = ms::numpress::MSNumpress::decodePic(&encoded[0], encodedBytes, &decoded[0]);
    
    assert(n == numDecoded);
    
    for (size_t i=0; i<n; i++) 
        assert(abs(ics[i] - decoded[i]) < 0.5);
    
    cout << "+ pass    encodeDecodePic " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeDecodePic5()

void encodeDecodePic5

(

)

Definition at line 581 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodePic(), and ms::numpress::MSNumpress::encodePic().

Referenced by main().

                        {
    srand(123459);
    
    size_t n = 1000;
    double ics[1000];
    for (size_t i=0; i<n; i++) 
        ics[i] = rand() % 1000000;
    
    ics[1] = 0.0;
    ics[2] = 0.0001;
    ics[3] = 0.0002;
    
    size_t encodedBytes;
    size_t numDecoded;
    unsigned char encoded[5000];
    double firstDecoded[1000];
    double decoded[1000];
    
    encodedBytes    = ms::numpress::MSNumpress::encodePic(&ics[0], n, &encoded[0]);
    numDecoded      = ms::numpress::MSNumpress::decodePic(&encoded[0], encodedBytes, &decoded[0]);
    
    for (size_t i=0; i<n; i++) 
        firstDecoded[i] = decoded[i];
    
    for (size_t i=0; i<5; i++) {
        encodedBytes    = ms::numpress::MSNumpress::encodePic(&decoded[0], n, &encoded[0]);
        numDecoded      = ms::numpress::MSNumpress::decodePic(&encoded[0], encodedBytes, &decoded[0]);
    }
    
    assert(n == numDecoded);
        
    for (size_t i=0; i<n; i++)
        if (firstDecoded[i] != decoded[i]) {
            cout << endl << firstDecoded[i] << " " << decoded[i] << endl;
            assert(firstDecoded[i] == decoded[i]);
        }
    
    cout << "+ pass    encodeDecodePic5 " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeDecodeSafe()

void encodeDecodeSafe

(

)

Definition at line 443 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeSafe(), and ms::numpress::MSNumpress::encodeSafe().

Referenced by main().

                        {
    srand(123459);
    
    double error;
    double eLim = 1.0e-300;
    size_t n = 1000;
    double mzs[1000];
    mzs[0] = 300 + rand() / double(RAND_MAX);
    for (size_t i=1; i<n; i++) 
        mzs[i] = mzs[i-1] + rand() / double(RAND_MAX);
    
    unsigned char encoded[8000];
    size_t encodedBytes = ms::numpress::MSNumpress::encodeSafe(&mzs[0], n, &encoded[0]);
    
    assert(encodedBytes == n*8);
    
    double decoded[1000];
    size_t numDecoded = ms::numpress::MSNumpress::decodeSafe(&encoded[0], encodedBytes, &decoded[0]);
    
    assert(numDecoded == n);
    for (size_t i=0; i<n; i++) { 
        error = abs(mzs[i] - decoded[i]);
        if (error >= eLim) {
            cout << "error   " << error << " is non-zero ( >= " << eLim << " )" << endl;
            assert(error == 0);
        }
    }
    cout << "+ pass    encodeDecodeSafe " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeDecodeSafeStraight()

void encodeDecodeSafeStraight

(

)

Definition at line 405 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeSafe(), and ms::numpress::MSNumpress::encodeSafe().

Referenced by main().

                                {
    double error;
    double eLim = 1.0e-300;
    size_t n = 15;
    double mzs[15];
    for (size_t i=0; i<n; i++) 
        mzs[i] = i+1;
    
    unsigned char encoded[8000];
    size_t encodedBytes = ms::numpress::MSNumpress::encodeSafe(&mzs[0], n, &encoded[0]);
    
    /*
    printf("\n");
    assert(encodedBytes == n*8);
    for (size_t i=0; i<n; i++) { 
        for (size_t j=0; j<8; j++) {
            printf("%x ", encoded[i*8 + j]);
        }
        printf("\n");
    }
    */
    
    double decoded[1000];
    size_t numDecoded = ms::numpress::MSNumpress::decodeSafe(&encoded[0], encodedBytes, &decoded[0]);
    
    assert(numDecoded == n);
    for (size_t i=0; i<n; i++) { 
        error = abs(mzs[i] - decoded[i]);
        if (error >= eLim) {
            cout << "error   " << error << " is non-zero ( >= " << eLim << " )" << endl;
            assert(error == 0);
        }
    }
    cout << "+ pass    encodeDecodeSafeStraight " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeDecodeSlof()

void encodeDecodeSlof

(

)

Definition at line 638 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeSlof(), ms::numpress::MSNumpress::encodeSlof(), and ms::numpress::MSNumpress::optimalSlofFixedPoint().

Referenced by main().

                        {
    srand(123459);
    
    size_t n = 1000;
    double ics[1000];
    for (size_t i=0; i<n; i++) 
        ics[i] = rand() % 1000000;
    
    ics[1] = 0.0;
    ics[2] = 0.0001;
    ics[3] = 0.0002;
    
    double fixedPoint = ms::numpress::MSNumpress::optimalSlofFixedPoint(&ics[0], n);
    
    unsigned char encoded[5000];
    size_t encodedBytes = ms::numpress::MSNumpress::encodeSlof(&ics[0], n, &encoded[0], fixedPoint);
    
    double decoded[1000];
    size_t numDecoded = ms::numpress::MSNumpress::decodeSlof(&encoded[0], encodedBytes, &decoded[0]);
    
    assert(n == numDecoded);
    
    double m = 0;
    double rm = 0;
    double mLim = 0.0005;
    double rmLim = 0.0005;
    double error;
    
    for (size_t i=0; i<n; i++)
        if (ics[i] < 1.0) {
            error = abs(ics[i] - decoded[i]);
            m = max(m, error);
            if (error >= mLim) {
                cout << endl << ics[i] << " " << decoded[i] << endl;
                assert(error < mLim);
            }
        } else {
            error = abs((ics[i] - decoded[i]) / ((ics[i] + decoded[i])/2));
            rm = max(rm, error);
            if (error >= rmLim) {
                cout << endl << ics[i] << " " << decoded[i] << endl;
                assert(error < rmLim);
            }
        }
    cout << "+               max error: " << m << "  limit: " << mLim << endl;
    cout << "+           max rel error: " << rm << "  limit: " << rmLim << endl;
    cout << "+ pass    encodeDecodeSlof " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeDecodeSlof5()

void encodeDecodeSlof5

(

)

Definition at line 689 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodeSlof(), ms::numpress::MSNumpress::encodeSlof(), and ms::numpress::MSNumpress::optimalSlofFixedPoint().

Referenced by main().

                         {
    srand(123459);
    
    size_t n = 1000;
    double ics[1000];
    for (size_t i=0; i<n; i++) 
        ics[i] = rand() % 1000000;
    
    ics[1] = 0.0;
    ics[2] = 0.0001;
    ics[3] = 0.0002;
    
    double fixedPoint = ms::numpress::MSNumpress::optimalSlofFixedPoint(&ics[0], n);
    
    size_t encodedBytes;
    size_t numDecoded;
    unsigned char encoded[5000];
    double firstDecoded[1000];
    double decoded[1000];
    
    encodedBytes    = ms::numpress::MSNumpress::encodeSlof(&ics[0], n, &encoded[0], fixedPoint);
    numDecoded      = ms::numpress::MSNumpress::decodeSlof(&encoded[0], encodedBytes, &decoded[0]);
    
    for (size_t i=0; i<n; i++) 
        firstDecoded[i] = decoded[i];
    
    for (size_t i=0; i<5; i++) {
        encodedBytes    = ms::numpress::MSNumpress::encodeSlof(&decoded[0], n, &encoded[0], fixedPoint);
        numDecoded      = ms::numpress::MSNumpress::decodeSlof(&encoded[0], encodedBytes, &decoded[0]);
    }
    
    double afterFixedPoint = ms::numpress::MSNumpress::optimalSlofFixedPoint(&decoded[0], n);
    
    assert(n == numDecoded);
    assert(fixedPoint == afterFixedPoint);
    
    for (size_t i=0; i<n; i++)
        if (firstDecoded[i] != decoded[i]) {
            cout << endl << firstDecoded[i] << " " << decoded[i] << endl;
            assert(firstDecoded[i] == decoded[i]);
        }
    
    cout << "+ pass    encodeDecodeSlof5 " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeLinear()

void encodeLinear

(

)

Definition at line 63 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::encodeLinear().

Referenced by main().

                    {
    double mzs[4];
    
    mzs[0] = 100.0;
    mzs[1] = 200.0;
    mzs[2] = 300.00005;
    mzs[3] = 400.00010;
    
    size_t nMzs = 4;
    unsigned char encoded[20];
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], 100000.0);
    
    assert(18 == encodedBytes);
    assert(0x80 == encoded[8]);
    assert(0x96 == encoded[9]);
    assert(0x98 == encoded[10]);
    assert(0x00 == encoded[11]);
    assert(0x75 == encoded[16]);
    assert(0x80 == encoded[17]);
    
    cout << "+ pass    encodeLinear " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

encodeLinear1()

void encodeLinear1

(

)

Definition at line 44 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::encodeLinear().

Referenced by main().

                     {

    double mzs[1];
    
    mzs[0] = 100.0;
    
    size_t nMzs = 1;
    unsigned char encoded[12];
    size_t encodedBytes = ms::numpress::MSNumpress::encodeLinear(&mzs[0], nMzs, &encoded[0], 100000.0);
    
    assert(12 == encodedBytes);
    assert(0x80 == encoded[8]);
    assert(0x96 == encoded[9]);
    assert(0x98 == encoded[10]);
    assert(0x00 == encoded[11]);
    
    cout << "+ pass    encodeLinear1 " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

main()

int main	(	int	argc,
		const char *	argv[]
	)

Definition at line 762 of file MSNumpressTest.cpp.

References decodeLinearCorrupt1(), decodeLinearCorrupt2(), decodeLinearNice(), decodeLinearNiceLowFP(), decodeLinearWierd(), encodeDecodeLinear(), encodeDecodeLinear5(), encodeDecodeLinearStraight(), encodeDecodePic(), encodeDecodePic5(), encodeDecodeSafe(), encodeDecodeSafeStraight(), encodeDecodeSlof(), encodeDecodeSlof5(), encodeLinear(), encodeLinear1(), optimalLinearFixedPoint(), optimalLinearFixedPointMass(), optimalSlofFixedPoint(), and testErroneousDecodePic().

                                       {
    optimalLinearFixedPoint();
    optimalLinearFixedPointMass();
    encodeLinear1();
    encodeLinear();
    decodeLinearNice();
    decodeLinearNiceLowFP();
    decodeLinearWierd();
    decodeLinearCorrupt1();
    decodeLinearCorrupt2();
    encodeDecodeLinearStraight();
    encodeDecodeLinear();
    encodeDecodePic();
    encodeDecodeSafeStraight();
    encodeDecodeSafe();
    optimalSlofFixedPoint();
    encodeDecodeSlof();
    encodeDecodeLinear5();
    encodeDecodePic5();
    encodeDecodeSlof5();
    testErroneousDecodePic();
    
    cout << "=== all tests succeeded! ===" << endl;
    return 0;
}

Here is the call graph for this function:

optimalLinearFixedPoint()

void optimalLinearFixedPoint

(

)

Definition at line 331 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::optimalLinearFixedPoint().

Referenced by main().

                               {
    srand(123459);
    
    size_t n = 1000;
    double mzs[1000];
    mzs[0] = 300 + (rand() % 1000) / 1000.0;
    for (size_t i=1; i<n; i++) 
        mzs[i] = mzs[i-1] + (rand() % 1000) / 1000.0;
    
    cout << "+                    max val: " << mzs[n-1] << endl;
    cout << "+          optimal linear fp: " << ms::numpress::MSNumpress::optimalLinearFixedPoint(&mzs[0], n) << endl;
    cout << "+ pass    optimalLinearFixedPoint " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

optimalLinearFixedPointMass()

void optimalLinearFixedPointMass

(

)

Definition at line 346 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::optimalLinearFixedPointMass().

Referenced by main().

                                   {
    srand(123459);
    
    size_t n = 1000;
    double mzs[1000];
    mzs[0] = 300 + (rand() % 1000) / 1000.0;
    for (size_t i=1; i<n; i++) 
        mzs[i] = mzs[i-1] + (rand() % 1000) / 1000.0;
    
    double fixedP; 

    fixedP = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], n, 0.1); 
    assert(abs(5.0 - fixedP) < 0.000005);

    fixedP = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], n, 0.01); 
    assert(abs(50.0 - fixedP) < 0.000005);

    fixedP = ms::numpress::MSNumpress::optimalLinearFixedPointMass(&mzs[0], n, 0.001); 
    assert(abs(500.0 - fixedP) < 0.000005);

    cout << "+          optimal linear fp: " << fixedP << endl;
    cout << "+ pass    optimalLinearFixedPointMass " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

optimalSlofFixedPoint()

void optimalSlofFixedPoint

(

)

Definition at line 623 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::optimalSlofFixedPoint().

Referenced by main().

                             {

    srand(123459);
    
    size_t n = 1000;
    double ics[1000];
    for (size_t i=0; i<n; i++) 
        ics[i] = rand() % 1000000;
        
    cout << "+           optimal slof fp: " << ms::numpress::MSNumpress::optimalSlofFixedPoint(&ics[0], n) << endl;
    cout << "+ pass    optimalSlofFixedPoint " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

testErroneousDecodePic()

void testErroneousDecodePic

(

)

Definition at line 736 of file MSNumpressTest.cpp.

References ms::numpress::MSNumpress::decodePic().

Referenced by main().

                              {
    std::vector<double> result;

    // set data to [  100, 102, 140, 92, 33, 80, 145  ]; // Base64 is "ZGaMXCFQkQ=="
    std::vector<unsigned char> data;
    data.resize(32);
    data[0] = 100;
    data[1] = 102;
    data[2] = 140;
    data[3] = 92;
    data[4] = 33;
    data[5] = 80;
    data[6] = 145;

    try {
        ms::numpress::MSNumpress::decodePic(data, result);
        cout << "- fail    testErroneousDecodePic: didn't throw exception for corrupt input " << endl << endl;
        assert(0 == 1);
    } catch (const char *err) {
        
    }
    
    cout << "+ pass    testErroneousDecodePic " << endl << endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

ENC_TWO_BYTE_FIXED_POINT

double ENC_TWO_BYTE_FIXED_POINT = 3000.0

Definition at line 40 of file MSNumpressTest.cpp.