MatchedFilterTest.cpp File Reference

#include "MatchedFilter.hpp"
#include "pwiz/utility/misc/unit.hpp"
#include "pwiz/utility/misc/Std.hpp"
#include <complex>
#include <cstring>
#include <typeinfo>

Go to the source code of this file.

Classes
struct	OneMinusAbs

Functions
void	test_SampledData ()
complex< double >	OneMinusAbsComplex (double d)
template<typename Kernel >
void	test_createFilter (const Kernel &f)
template<typename Kernel >
void	test_createFilters (const Kernel &f)
template<typename Kernel >
void	test_compute (const Kernel &f)
template<typename Kernel >
void	test_kernel (const Kernel &kernel)
int	main (int argc, char *argv[])

Variables
ostream *	os_ = 0

---

Function Documentation

void test_SampledData

(

)

Definition at line 38 of file MatchedFilterTest.cpp.

References os_, unit_assert, and x.

Referenced by main().

{
    if (os_) *os_ << "test_SampledData()\n";
    using namespace MatchedFilter;

    SampledData<DxD> sd;
    sd.domain = make_pair(0.0, 2.0);
    sd.samples.push_back(10);
    sd.samples.push_back(11);
    sd.samples.push_back(12);
    sd.samples.push_back(13);
    sd.samples.push_back(14);

    if (os_) *os_ << sd << endl; 

    if (os_) *os_ << "domainWidth: " << sd.domainWidth() << endl;
    unit_assert(sd.domainWidth() == 2.0);

    if (os_) *os_ << "dx: " << sd.dx() << endl;
    unit_assert(sd.dx() == .5); 
    
    for (unsigned int i=0; i<sd.samples.size(); i++)
    {
        if (os_) *os_ << "x(" << i << "): " << sd.x(i) << endl;
        unit_assert(sd.x(i) == i * .5);
    }

    unsigned int count = 0;
    for (double x=-.2; x<2.3; x+=.1, count++)
    {
        //if (os_) *os_ << "sampleIndex(" << x << "): " << sd.sampleIndex(x) << endl;
        unit_assert(sd.sampleIndex(x) == count/5);
        //if (os_) *os_ << "sample(" << x << "): " << sd.sample(x) << endl;
        unit_assert(sd.sample(x) == 10 + count/5);        
    }
}

complex<double> OneMinusAbsComplex

(

double

d

)

Definition at line 86 of file MatchedFilterTest.cpp.

Referenced by main().

{
    return (d>=-1 && d<=1) ? 1 - abs(d) : 0;  
}

template<typename Kernel >

void test_createFilter

(

const Kernel &

f

)

Definition at line 93 of file MatchedFilterTest.cpp.

References pwiz::math::MatchedFilter::details::createFilter(), os_, and unit_assert.

Referenced by test_kernel().

{
    using namespace MatchedFilter;

    if (os_) *os_ << "test_createFilter() " << typeid(f).name() << endl;

    int sampleRadius = 4;
    double dx = .25;
    double shift = 0;

    typedef typename KernelTraits<Kernel>::filter_type filter_type;
    typedef typename KernelTraits<Kernel>::abscissa_type abscissa_type;
    typedef typename KernelTraits<Kernel>::ordinate_type ordinate_type;

    filter_type filter = details::createFilter(f, sampleRadius, dx, shift);

    if (os_)
    {
        copy(filter.begin(), filter.end(), ostream_iterator<ordinate_type>(*os_, " "));
        *os_ << endl;
    }

    unit_assert((int)filter.size() == sampleRadius*2 + 1);
    for (int i=-sampleRadius; i<=sampleRadius; ++i)
        unit_assert(filter[sampleRadius+i] == f(i*dx - shift));
    
    if (os_) *os_ << endl;
}

template<typename Kernel >

void test_createFilters

(

const Kernel &

f

)

Definition at line 124 of file MatchedFilterTest.cpp.

References pwiz::math::MatchedFilter::details::createFilters(), e(), norm(), os_, unit_assert, and unit_assert_equal.

Referenced by test_kernel().

{
    using namespace MatchedFilter;

    if (os_) *os_ << "test_createFilters() " << typeid(f).name() << endl;

    int sampleRadius = 2;
    int subsampleFactor = 4;
    double dx = 1;

    typedef typename KernelTraits<Kernel>::filter_type filter_type;
    typedef typename KernelTraits<Kernel>::ordinate_type ordinate_type;
    vector<filter_type> filters = details::createFilters(f, 
                                                         sampleRadius, 
                                                         subsampleFactor, 
                                                         dx);
    
    // verify filter count
    unit_assert((int)filters.size() == subsampleFactor);

    for (typename vector<filter_type>::const_iterator it=filters.begin(); it!=filters.end(); ++it)
    {
        if (os_)
        {
            copy(it->begin(), it->end(), ostream_iterator<ordinate_type>(*os_, " "));
            *os_ << endl;
        }

        // verify filter size
        unit_assert((int)it->size() == sampleRadius*2 + 1);
    
        // verify filter normalization
        double sum = 0;
        for (typename filter_type::const_iterator jt=it->begin(); jt!=it->end(); ++jt)
            sum += norm(complex<double>(*jt));
        unit_assert_equal(sum, 1, 1e-14);
    }

    if (os_) *os_ << endl;
}

template<typename Kernel >

void test_compute

(

const Kernel &

f

)

Definition at line 167 of file MatchedFilterTest.cpp.

References pwiz::math::MatchedFilter::computeCorrelationData(), e(), os_, and unit_assert.

Referenced by test_kernel().

{
    using namespace MatchedFilter;

    if (os_) *os_ << "test_compute() " << typeid(f).name() << endl;

    typename KernelTraits<Kernel>::sampled_data_type data;
    data.domain = make_pair(0, 10);
    data.samples.resize(11);
    data.samples[5] = 1.;

    if (os_) *os_ << "data: " << data << endl;

    int sampleRadius = 2;
    int sampleFactor = 4;

    typedef typename KernelTraits<Kernel>::correlation_data_type CorrelationData;

    CorrelationData correlationData = 
        computeCorrelationData(data, f, sampleRadius, sampleFactor); 

    if (os_) *os_ << "correlationData: " << correlationData << endl;

    unit_assert(correlationData.samples.size() == 41);
    unit_assert(abs(correlationData.samples[20].dot - 1.) < 1e-12);
}

template<typename Kernel >

void test_kernel

(

const Kernel &

kernel

)

Definition at line 196 of file MatchedFilterTest.cpp.

References os_, test_compute(), test_createFilter(), and test_createFilters().

Referenced by main().

{
    if (os_) *os_ << "***************************************************************\n";
    if (os_) *os_ << "test_kernel() " << typeid(kernel).name() << endl;
    if (os_) *os_ << "***************************************************************\n";
    if (os_) *os_ << endl;

    test_createFilter(kernel);
    test_createFilters(kernel);
    test_compute(kernel);
}

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

Definition at line 209 of file MatchedFilterTest.cpp.

References e(), OneMinusAbsComplex(), os_, TEST_EPILOG, TEST_FAILED, test_kernel(), TEST_PROLOG, and test_SampledData().

{
    TEST_PROLOG(argc, argv)

    try
    {
        if (argc>1 && !strcmp(argv[1],"-v")) os_ = &cout;
        if (os_) *os_ << "MatchedFilterTest\n";

        test_SampledData();
        test_kernel(OneMinusAbs()); 
        test_kernel(&OneMinusAbsComplex); 

    }
    catch (exception& e)
    {
        TEST_FAILED(e.what())
    }
    catch (...)
    {
        TEST_FAILED("Caught unknown exception.")
    }

    TEST_EPILOG
}

---

Variable Documentation

ostream* os_ = 0

Definition at line 35 of file MatchedFilterTest.cpp.