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

---

M9poCw.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Trait_M9poCw_h
#define Impala_Core_Array_Trait_M9poCw_h

#include "Core/Array/Pattern/Categories.h"
#include "Core/Array/Element/E1Sqrt.h"
#include "Core/Array/Element/E1Log.h"

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


template<class DstArrayT, class SrcArrayT>
class M9poCw
{
public:
    typedef Pattern::TagTransInVar TransVarianceCategory;

    typedef Pattern::TagCallValue CallCategory;

    typedef typename DstArrayT::ArithType DstArithT;
    typedef typename SrcArrayT::ArithType SrcArithT;

    M9poCw(SrcArithT scale)
    {
        mScale = scale;
    }

    DstArithT
    DoIt(const SrcArithT& E, const SrcArithT& El, const SrcArithT& Ell,
         const SrcArithT& Ex, const SrcArithT& Elx, const SrcArithT& Ellx,
         const SrcArithT& Ey, const SrcArithT& Ely, const SrcArithT& Elly)
    {
        SrcArithT Clx_a = Elx*E - El*Ex;
        SrcArithT Cly_a = Ely*E - El*Ey;
        SrcArithT Clw_a = Clx_a*Clx_a + Cly_a*Cly_a;
        SrcArithT Cllx_a = Ellx*E - Ell*Ex;
        SrcArithT Clly_a = Elly*E - Ell*Ey;
        SrcArithT Cllw_a = Cllx_a*Cllx_a + Clly_a*Clly_a;
        SrcArithT Cw_a = Element::E1Sqrt(Clw_a + Cllw_a);
        SrcArithT norm = E*E/mScale;
        return Cw_a/norm;
        //return Element::E1Log(Cw_a/norm + 1);
        //return Element::E1Log(100000000*Cw_a+1) - Element::E1Log(100000000*norm+1); 
}

private:
    SrcArithT mScale;
};

template<class DstArrayT, class SrcArrayT>
class M9poCw_log
{
public:
    typedef Pattern::TagTransInVar TransVarianceCategory;

    typedef Pattern::TagCallValue CallCategory;

    typedef typename DstArrayT::ArithType DstArithT;
    typedef typename SrcArrayT::ArithType SrcArithT;

    M9poCw_log(SrcArithT scale)
    {
        mScale = scale;
    }

    DstArithT
    DoIt(const SrcArithT& E, const SrcArithT& El, const SrcArithT& Ell,
         const SrcArithT& Ex, const SrcArithT& Elx, const SrcArithT& Ellx,
         const SrcArithT& Ey, const SrcArithT& Ely, const SrcArithT& Elly)
    {
        SrcArithT Clx_a = Elx*E - El*Ex;
        SrcArithT Cly_a = Ely*E - El*Ey;
        SrcArithT Clw_a = Clx_a*Clx_a + Cly_a*Cly_a;
        SrcArithT Cllx_a = Ellx*E - Ell*Ex;
        SrcArithT Clly_a = Elly*E - Ell*Ey;
        SrcArithT Cllw_a = Cllx_a*Cllx_a + Clly_a*Clly_a;
        SrcArithT Cw_a = Element::E1Sqrt(Clw_a + Cllw_a);
        SrcArithT norm = E*E/mScale;
        //return Cw_a/norm;
        return Element::E1Log(Cw_a/norm + 1);
        //return Element::E1Log(100000000*Cw_a+1) - Element::E1Log(100000000*norm+1); 
}

private:
    SrcArithT mScale;
};

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

#endif

---