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

---

PatM6PixOp.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Pattern_PatM6PixOp_h
#define Impala_Core_Array_Pattern_PatM6PixOp_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
PatM6PixOp(DstArrayT*& dst,
           SrcArrayT* s1, SrcArrayT* s2, SrcArrayT* s3,
           SrcArrayT* s4, SrcArrayT* s5, SrcArrayT* s6, MpoT& mpo)
{
    if (dst == 0)
        dst = ArrayClone<DstArrayT>(s1);

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

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

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

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

#endif

---