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

---

PatGenConv2dSep.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Pattern_PatGenConv2dSep_h
#define Impala_Core_Array_Pattern_PatGenConv2dSep_h

#include "Core/Array/CheckBorderSize.h"
#include "Core/Array/Pattern/PatSetBorder.h"
#include "Core/Array/Pattern/FuncGenConv2dSep.h"
#include "Util/TimePlot.h"

#ifdef PX_HORUS_USED
#include "Core/Array/Pattern/PxArrayFunc.h"
#include "Core/Array/Pattern/PxStateTrans.h"
#endif

namespace Impala
{
namespace Core
{
namespace Array
{
namespace Pattern
{


template<class DstArrayT, class SrcArrayT, class KerArrayT,
         class PixOpT, class RedOpT>
inline void
PatGenConv2dSep(DstArrayT*& dst, SrcArrayT* src, KerArrayT* ker1,
                KerArrayT* ker2, PixOpT& pixOp, RedOpT& redOp,
                int vType, Util::TimePlot* plotter)
{
    if (vType == 6) // the default has not been implemented yet
        vType = 0;

    int numX = ArrayCW(ker1) / 2;
    int numY = ArrayCW(ker2) / 2;
    SrcArrayT* checkedSrc = CheckBorderSize(src, numX, numY);

    if ((dst == 0) || (dst == checkedSrc))
        dst = ArrayClone<DstArrayT>(checkedSrc);

#ifdef PX_HORUS_USED
    if (!PxRunParallel()) {                         // run sequential
#endif
        if (plotter)
            plotter->SplitTime();
        if ((vType == 0) || (vType == 1) || (vType == 2))
        {
            PatSetBorder(checkedSrc, numX, 0);
            if (plotter)
                plotter->SplitTime();
            KerArrayT* scratch = ArrayClone<KerArrayT>(checkedSrc);
            if (plotter)
                plotter->SplitTime();
            FuncGenConv2dSepDispatch_H(scratch, checkedSrc,
                                         ker1, pixOp, redOp, vType);
            if (plotter)
                plotter->SplitTime();
            PatSetBorder(scratch, 0, numY);
            if (plotter)
                plotter->SplitTime();
            FuncGenConv2dSepDispatch_V(dst, scratch,
                                         ker2, pixOp, redOp, vType);
            delete scratch;
        }
        else
        {
            PatSetBorder(checkedSrc, numX, numY);
            if (plotter)
                plotter->SplitTime();
            if (plotter)
                plotter->SplitTime();
            if (plotter)
                plotter->SplitTime();
            if (plotter)
                plotter->SplitTime();
            //FuncGenConv2dSepDispatch(dst, checkedSrc, ker1,
            //                           ker2, pixOp, redOp, vType);
        }

#ifdef PX_HORUS_USED
    } else {                                        // run parallel
        PxArrayPreStateTrans(checkedSrc, PAR_PART, STRONG);
        PxArrayPreStateTrans(dst, PAR_PART, WEAK);
        if ((vType == 0) || (vType == 1) || (vType == 2)) {
            PatSetBorder(PxArrayPD(checkedSrc), numX, 0);
            KerArrayT* scratch =
                        ArrayClone<KerArrayT>(PxArrayPD(checkedSrc));
            FuncGenConv2dSepDispatch_H(scratch, PxArrayPD(checkedSrc),
                                         ker1, pixOp, redOp, vType);
            PatSetBorder(scratch, 0, numY);
            FuncGenConv2dSepDispatch_V(PxArrayPD(dst), scratch,
                                         ker2, pixOp, redOp, vType);
            delete scratch;
        } else {
            PatSetBorder(PxArrayPD(checkedSrc), numX, numY);
            //FuncGenConv2dSepDispatch(PxArrayPD(dst),
            //                           PxArrayPD(checkedSrc), ker1,
            //                           ker2, pixOp, redOp, vType);
        }
        if (PxArrayPD(checkedSrc) != PxArrayPD(src))
            PxArrayDeletePD(checkedSrc);
        PxArrayPostStateTrans(dst);

        if (!PxRunLazyParallel()) {
            PxArrayForceNonDistributed(src);
            PxArrayForceNonDistributed(dst);
        }
    }
#endif

    if (checkedSrc != src)
        delete checkedSrc;
}

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

#endif

---