00001 #ifndef Impala_Core_Array_Trait_M9poJOSAw_h 00002 #define Impala_Core_Array_Trait_M9poJOSAw_h 00003 00004 #include "Core/Array/Pattern/Categories.h" 00005 #include "Core/Array/Element/E1Sqrt.h" 00006 00007 namespace Impala 00008 { 00009 namespace Core 00010 { 00011 namespace Array 00012 { 00013 namespace Trait 00014 { 00015 00016 00019 template<class DstArrayT, class SrcArrayT> 00020 class M9poJOSAw 00021 { 00022 public: 00024 typedef Pattern::TagTransInVar TransVarianceCategory; 00025 00027 typedef Pattern::TagCallValue CallCategory; 00028 00029 typedef typename DstArrayT::ArithType DstArithT; 00030 typedef typename SrcArrayT::ArithType SrcArithT; 00031 00033 M9poJOSAw() 00034 { 00035 } 00036 00038 DstArithT 00039 DoIt(const SrcArithT& R, const SrcArithT& G, const SrcArithT& B, 00040 const SrcArithT& Rx, const SrcArithT& Gx, const SrcArithT& Bx, 00041 const SrcArithT& Ry, const SrcArithT& Gy, const SrcArithT& By) 00042 { 00043 SrcArithT rrx = (G*Rx-Gx*R)/(G*G); 00044 SrcArithT rry = (G*Ry-Gy*R)/(G*G); 00045 SrcArithT ggx = (B*Gx-Bx*G)/(B*B); 00046 SrcArithT ggy = (B*Gy-By*G)/(B*B); 00047 SrcArithT bbx = (R*Bx-Rx*B)/(R*R); 00048 SrcArithT bby = (R*By-Ry*B)/(R*R); 00049 00050 return Element::E1Sqrt(rrx*rrx+rry*rry+ggx*ggx+ggy*ggy+bbx*bbx+bby*bby); 00051 } 00052 00053 }; 00054 00055 } // namespace Trait 00056 } // namespace Array 00057 } // namespace Core 00058 } // namespace Impala 00059 00060 #endif