Impala Documentation

template<class DstArrayT, class SrcArrayT>

DstArithT Impala::Core::Array::Trait::M9poNw_log< DstArrayT, SrcArrayT >::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
	)			`[inline]`

Actual operation.

Definition at line 88 of file M9poNw.h.

References Impala::Core::Array::Element::E1Log(), Impala::Core::Array::Element::E1Sqrt(), and Impala::Core::Array::Trait::M9poNw_log< DstArrayT, SrcArrayT >::mScale.

00091     {
00092         SrcArithT norm2 = E*E/mScale;
00093         SrcArithT norm3 = E*E*E/mScale;
00094         SrcArithT Nlx = Elx*E - El*Ex;
00095         SrcArithT Nly = Ely*E - El*Ey;
00096         SrcArithT Nlw = Element::E1Sqrt(Nlx*Nlx + Nly*Nly) / norm2;
00097 
00098         SrcArithT A1 = Ellx*E*E - Ell*Ex*E;
00099         SrcArithT A2 = El*El*Ex - Elx*El*E;
00100         SrcArithT A = A1 + 2*A2;
00101         SrcArithT B1 = Elly*E*E - Ell*Ey*E;
00102         SrcArithT B2 = El*El*Ey - Ely*El*E;
00103         SrcArithT B = B1 + 2*B2;
00104 
00105         SrcArithT Nllw = Element::E1Sqrt(A*A + B*B) / norm3;
00106         //return Element::E1Sqrt(Nlw*Nlw + Nllw*Nllw);
00107         return Element::E1Log(Element::E1Sqrt(Nlw*Nlw + Nllw*Nllw) + 1);
00108     }

Here is the call graph for this function: