source: branches/OptimisedGrowerTesting/src/ATrous/HaarFilter.cc @ 1441

#include <duchamp/ATrous/HaarFilter.hh>
#include <duchamp/ATrous/Filter.hh>

namespace duchamp {

    HaarFilter::HaarFilter():
    Filter()
    {
        const unsigned int filtersize=3;
        double filter[filtersize] = {0., 1./2., 1./2.};
        this->filter1D = std::vector<double>(filter,filter+filtersize);
        this->name = "B3 spline function";
        
        const unsigned int size1D=6;
        this->maxNumScales[0] = size1D;
        double sigmaFactors1D[size1D+1] = 
            {1.00000000000,7.07167810e-1,5.00000000e-1,3.53553391e-1,
             2.50000000e-1,1.76776695e-1,1.25000000e-1};

        const unsigned int size2D=6;
        this->maxNumScales[1] = size2D;
        double sigmaFactors2D[size2D+1] = 
            {1.00000000000,4.33012702e-1,2.16506351e-1,1.08253175e-1,
             5.41265877e-2,2.70632939e-2,1.35316469e-2};

        const unsigned int size3D=8;
        this->maxNumScales[2] = size3D;
        double sigmaFactors3D[size3D+1] = 
            {1.00000000000,9.35414347e-1,3.30718914e-1,1.16926793e-1,
             4.13398642e-2,1.46158492e-2,5.16748303e-3,1.82698115e-3,
             6.45935379e-44};

        this->sigmaFactors[0] = new std::vector<double>(sigmaFactors1D,sigmaFactors1D+this->maxNumScales[0]);
        this->sigmaFactors[1] = new std::vector<double>(sigmaFactors2D,sigmaFactors2D+this->maxNumScales[1]);
        this->sigmaFactors[2] = new std::vector<double>(sigmaFactors3D,sigmaFactors3D+this->maxNumScales[2]);

    }

    HaarFilter::HaarFilter(const HaarFilter& other):
    Filter(other)
    {
    }

    HaarFilter& HaarFilter::operator= (const HaarFilter& other)
    {
        if(this == &other) return *this;
        ((Filter &) *this) = other;
        return *this;
    }

}