|| || || || || || || ||

---

RecGauss.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_RecGauss_h
#define Impala_Core_Array_RecGauss_h

#include "Core/Array/MakeGaussIIR1d.h"
#include "Core/Array/RecConvKernelSeparated.h"
#include "Core/Array/ConvKernel1d.h"
#include "Core/Array/ConvKernelSeparated.h"
#include "Core/Array/MakeFromData.h"

namespace Impala
{
namespace Core
{
namespace Array
{


template<class DstArrayT, class SrcArrayT>
inline void
RecGauss(DstArrayT*& dst, SrcArrayT* src, double sx, double sy, int dx, int dy,
         int recurOrder)
{
    bool hasDer = ((dx != 0) || (dy != 0));
    DstArrayT* tmp = 0;

    Array2dScalarReal64* gx = MakeGaussIIR1d(sx, dx, recurOrder);
    Array2dScalarReal64* gy = MakeGaussIIR1d(sy, dy, recurOrder);
    if (hasDer)
        RecConvKernelSeparated(tmp, src, gx, gy);
    else
        RecConvKernelSeparated(dst, src, gx, gy);
    delete gx;
    delete gy;

    if (!hasDer)
        return;

    double d1[3] = {-0.5, 0, 0.5};
    double d2[3] = {1, -2, 1};
    Array2dScalarReal64* kerx = 0;
    Array2dScalarReal64* kery = 0;

    if (dx == 1)
        kerx = MakeFromData<Array2dScalarReal64>(d1, 3, 1, 0, 0, true);
//      kerx = MakeFromData<Array2dScalarReal64>(d1, 3, 1);
    if (dx == 2)
        kerx = MakeFromData<Array2dScalarReal64>(d2, 3, 1, 0, 0, true);
//      kerx = MakeFromData<Array2dScalarReal64>(d2, 3, 1);
    if (dy == 1)
        kery = MakeFromData<Array2dScalarReal64>(d1, 3, 1, 0, 0, true);
//      kery = MakeFromData<Array2dScalarReal64>(d1, 3, 1);
    if (dy == 2)
        kery = MakeFromData<Array2dScalarReal64>(d2, 3, 1, 0, 0, true);
//      kery = MakeFromData<Array2dScalarReal64>(d2, 3, 1);

    if ((kerx != 0) && (kery != 0))
    {
        std::cout << "doing it" << std::endl;
        ConvKernelSeparated(dst, tmp, kerx, kery);
    }
    else
    {
        if (kerx != 0)
            ConvKernel1d(dst, tmp, 1, kerx);
        else
            ConvKernel1d(dst, tmp, 2, kery);
    }
    delete kerx;
    delete kery;
    delete tmp;
}

} // namespace Array
} // namespace Core
} // namespace Impala

#endif

---