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

---

PatMPixOp.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Pattern_PatMPixOp_h
#define Impala_Core_Array_Pattern_PatMPixOp_h

#include "Core/Array/Array2dTem.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
PatMPixOp(DstArrayT*& dst,
          const std::vector<SrcArrayT*>& srcList,
          MpoT& mpo)
{
    int length = srcList.size();
    if (length < 1)
        return;
    if (dst == 0)
        dst = ArrayClone<DstArrayT>(srcList[0]);

#ifdef PX_HORUS_USED
    if (!PxRunParallel()) {                         // run sequential
#endif
        FuncMpo(dst, srcList, mpo);

#ifdef PX_HORUS_USED
    } else {                                        // run parallel
        int i;
        for(i=0 ; i<length ; ++i)
            PxArrayPreStateTrans(srcList[i],  PAR_PART, STRONG);

        std::vector<SrcArrayT*> pdList;
        pdList.resize(length);
        for(i=0 ; i<length ; ++i)
            pdList[i] = PxArrayPD(srcList[i]);

        FuncMpo(PxArrayPD(dst), pdList, mpo);
        PxArrayPostStateTrans(dst);

        if (!PxRunLazyParallel()) {
            for(i=0 ; i<length ; ++i)
                PxArrayForceNonDistributed(srcList[i]);
            PxArrayForceNonDistributed(dst);
        }
    }
#endif
}

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

#endif

---