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

---

PatM3PixOp.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Pattern_PatM3PixOp_h
#define Impala_Core_Array_Pattern_PatM3PixOp_h

#include "Core/Array/Pattern/FuncMpo.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 MpoT>
inline void
PatM3PixOp(DstArrayT*& dst,
           SrcArrayT* s1, SrcArrayT* s2, SrcArrayT* s3, MpoT& mpo)
{
    if (dst == 0)
        dst = ArrayClone<DstArrayT>(s1);

#ifdef PX_HORUS_USED
    if (!PxRunParallel()) {                         // run sequential
#endif
        FuncM3po(dst, s1, s2, s3, mpo);

#ifdef PX_HORUS_USED
    } else {                                        // run parallel
        PxArrayPreStateTrans(s1,  PAR_PART, STRONG);
        PxArrayPreStateTrans(s2,  PAR_PART, STRONG);
        PxArrayPreStateTrans(s3,  PAR_PART, STRONG);
        PxArrayPreStateTrans(dst, PAR_PART, WEAK);
        FuncM3po(PxArrayPD(dst),
                 PxArrayPD(s1), PxArrayPD(s2), PxArrayPD(s3), mpo);
        PxArrayPostStateTrans(dst);

        if (!PxRunLazyParallel()) {
            PxArrayForceNonDistributed(s1);
            PxArrayForceNonDistributed(s2);
            PxArrayForceNonDistributed(s3);
            PxArrayForceNonDistributed(dst);
        }
    }
#endif
}

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

#endif

---