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

---

M9N4poFSph.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Trait_M9N4poFSph_h
#define Impala_Core_Array_Trait_M9N4poFSph_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 M9N4poFSph
{
public:
    typedef Pattern::TagTransInVar TransVarianceCategory;

    typedef Pattern::TagCallValue CallCategory;

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

    M9N4poFSph()
    {
    }

    void
    DoIt(DstStorT* FThetax, DstStorT* FThetay, DstStorT* FPhix, DstStorT* FPhiy,
         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 I1 = Element::E1Sqrt(R*R + G*G + B*B);
        SrcArithT I2 = Element::E1Sqrt(R*R + G*G);
        SrcArithT I1I2 = I1*I2;

        *FThetax = (I2 > 0) ? (R*Gx - G*Rx) / I2 : I2;
        *FThetay = (I2 > 0) ? (R*Gy - G*Ry) / I2 : I2;
        *FPhix = (I1I2 > 0) ? (G*(B*Gx - G*Bx) + R*(B*Rx - R*Bx)) / (I1I2) : 0;
        *FPhiy = (I1I2 > 0) ? (G*(B*Gy - G*By) + R*(B*Ry - R*By)) / (I1I2) : 0;
    }
};

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

#endif

---