SavitzkyGolaySmootherTest.cpp File Reference

#include "SavitzkyGolaySmoother.hpp"
#include "pwiz/utility/misc/unit.hpp"
#include "pwiz/utility/misc/Std.hpp"

Go to the source code of this file.

Classes
struct	TestData

Functions
vector< double >	parseDoubleArray (const string &doubleArray)
void	test ()
int	main (int argc, char *argv[])

Variables
ostream *	os_ = 0
const TestData	testData []
const size_t	testDataSize = sizeof(testData) / sizeof(TestData)

---

Function Documentation

vector<double> parseDoubleArray

(

const string &

doubleArray

)

Definition at line 162 of file SavitzkyGolaySmootherTest.cpp.

References max().

{
    vector<double> doubleVector;
    vector<string> tokens;
    bal::split(tokens, doubleArray, bal::is_space(), bal::token_compress_on);
    if (!tokens.empty() && !tokens[0].empty())
        for (size_t i=0; i < tokens.size(); ++i)
            doubleVector.push_back(max(0.0, lexical_cast<double>(tokens[i])));
    return doubleVector;
}

void test

(

)

Definition at line 174 of file SavitzkyGolaySmootherTest.cpp.

References e(), parseDoubleArray(), TestData::polynomialOrder, testData, testDataSize, unit_assert, unit_assert_equal, TestData::windowSize, TestData::xRaw, TestData::xSmoothed, TestData::yRaw, and TestData::ySmoothed.

{
    // test that invalid value exceptions are thrown

    // order too low
    /*unit_assert_throws_what(SavitzkyGolaySmoother(1, 15), runtime_error, \
        "[SavitzkyGolaySmoother::ctor()] Invalid value for polynomial order; valid range is [2, 20]");

    // order too high
    unit_assert_throws_what(SavitzkyGolaySmoother(21, 15), runtime_error, \
        "[SavitzkyGolaySmoother::ctor()] Invalid value for polynomial order; valid range is [2, 20]");

    // window size too small
    unit_assert_throws_what(SavitzkyGolaySmoother(2, 3), runtime_error, \
        "[SavitzkyGolaySmoother::ctor()] Invalid value for window size; value must be odd and in range [5, infinity)");

    // window size isn't odd
    unit_assert_throws_what(SavitzkyGolaySmoother(2, 6), runtime_error, \
        "[SavitzkyGolaySmoother::ctor()] Invalid value for window size; value must be odd and in range [5, infinity)");
*/
    SavitzkyGolaySmoother(2, 100001); // window size is valid up to numeric limits

    /*if (os_)
    {
        *os_ << "Unsmoothed data (" << testY.size() << "):\t";
        copy(testY.begin(), testY.end(), ostream_iterator<double>(*os_, "\t"));
        *os_ << endl;
    }*/

    /*if (os_)
    {
        *os_ << "Smoothed data (" << smoothData.size() << "):\t";
        copy(smoothData.begin(), smoothData.end(), ostream_iterator<double>(*os_, "\t"));
        *os_ << endl;
    }*/

    for (size_t i=0; i < testDataSize; ++i)
    {
        const TestData& data = testData[i];
        
        vector<double> xRaw = parseDoubleArray(data.xRaw);
        vector<double> yRaw = parseDoubleArray(data.yRaw);
        vector<double> targetSmoothedX = parseDoubleArray(data.xSmoothed);
        vector<double> targetSmoothedY = parseDoubleArray(data.ySmoothed);

        // sanity checks
        unit_assert(xRaw.size() == yRaw.size());
        unit_assert(targetSmoothedX.size() == targetSmoothedY.size());

        SavitzkyGolaySmoother smoother(data.polynomialOrder, data.windowSize);
        vector<double> smoothedX, smoothedY;
        smoother.smooth(xRaw, yRaw, smoothedX, smoothedY);

        unit_assert(smoothedX.size() == targetSmoothedX.size());
        unit_assert(smoothedY.size() == targetSmoothedY.size());

        for (size_t j=0; j < smoothedY.size(); ++j)
            unit_assert_equal(smoothedY[j], targetSmoothedY[j], 1e-5);
    }
}

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

Definition at line 236 of file SavitzkyGolaySmootherTest.cpp.

References e(), os_, test(), TEST_EPILOG, TEST_FAILED, and TEST_PROLOG.

{
    TEST_PROLOG(argc, argv)

    try
    {
        if (argc>1 && !strcmp(argv[1],"-v")) os_ = &cout;
        test();
    }
    catch (exception& e)
    {
        TEST_FAILED(e.what())
    }
    catch (...)
    {
        TEST_FAILED("Caught unknown exception.")
    }

    TEST_EPILOG
}

---

Variable Documentation

ostream* os_ = 0

Definition at line 32 of file SavitzkyGolaySmootherTest.cpp.

const TestData testData[]

Definition at line 44 of file SavitzkyGolaySmootherTest.cpp.

const size_t testDataSize = sizeof(testData) / sizeof(TestData)

Definition at line 159 of file SavitzkyGolaySmootherTest.cpp.