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

---

M9poANG2w.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Trait_M9poANG2w_h
#define Impala_Core_Array_Trait_M9poANG2w_h

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

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


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

    typedef Pattern::TagCallValue CallCategory;

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

    M9poANG2w()
    {
    }

    DstArithT
    DoIt(const SrcArithT& R, const SrcArithT& G, const SrcArithT& B,
         const SrcArithT& Rx, const SrcArithT& Gx, const SrcArithT& Bx,
         const SrcArithT& Ry, const SrcArithT& Gy, const SrcArithT& By)
    {
        SrcArithT I = (R*R+G*G+B*B);
        SrcArithT fw2 = Rx*Rx+Ry*Ry+Gx*Gx+Gy*Gy+Bx*Bx+By*By;
        SrcArithT A1 = R*Rx+G*Gx+B*Bx;
        SrcArithT A2 = R*Ry+G*Gy+B*By;
        SrcArithT f_r = (A1*A1+A2*A2)/I;
        SrcArithT f_ang = Element::E1Sqrt(fw2-f_r);
        SrcArithT f_ang2 = f_ang/Element::E1Sqrt(I);
        return f_ang2;
    }

};

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

#endif

---