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

---

BpoTestEqual.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Trait_BpoTestEqual_h
#define Impala_Core_Array_Trait_BpoTestEqual_h

#include "Core/Array/Pattern/Categories.h"
#include "Core/Array/Element/E1Cast.h"
#include "Core/Array/Element/E1Abs.h"
#include "Core/Array/Element/E1Sum.h"

namespace Impala
{
namespace Core
{
namespace Array
{
namespace Trait
{

using namespace Impala::Core::Array::Element;


template<class DstArrayT, class Src1ArrayT, class Src2ArrayT>
class BpoTestEqual
{
public:
    typedef Pattern::TagTransInVar TransVarianceCategory;

    typedef Pattern::TagCallValue CallCategory;

    typedef typename DstArrayT::ArithType DstArithT;
    typedef typename Src1ArrayT::ArithType Src1ArithT;
    typedef typename Src2ArrayT::ArithType Src2ArithT;

    BpoTestEqual(Real64 epsilon)
    { 
        mEpsilon = epsilon;
        mNr = 0;
        mNrDiff = 0;
        mSumSrc1 = 0;
        mSumSrc2 = 0;
        mSumDiff = 0;
    }

    DstArithT
    DoIt(const Src1ArithT& x, const Src2ArithT& y)
    {
        mSumSrc1 += E1Sum(E1Abs(x));
        mSumSrc2 += E1Sum(E1Abs(y));
        mNr++;
        Real64 d = E1Sum(E1Abs(x - y));
        if (d > mEpsilon) {
            mNrDiff++;
            mSumDiff += d;
            return E1Cast(d, DstArithT());
        }
        return E1Cast(0, DstArithT());
    }

    //private:
    Real64 mEpsilon;
    Real64 mSumSrc1;
    Real64 mSumSrc2;
    Real64 mSumDiff;
    Int32 mNr;
    Int32 mNrDiff;
};

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

#endif

---