Impala Documentation

template<class DstArrayT, class SrcArrayT>

DstArithT Impala::Core::Array::Trait::M9poNw< 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 41 of file M9poNw.h.

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

00044     {
00045         SrcArithT norm2 = E*E/mScale;
00046         SrcArithT norm3 = E*E*E/mScale;
00047         SrcArithT Nlx = Elx*E - El*Ex;
00048         SrcArithT Nly = Ely*E - El*Ey;
00049         SrcArithT Nlw = Element::E1Sqrt(Nlx*Nlx + Nly*Nly) / norm2;
00050 
00051         SrcArithT A1 = Ellx*E*E - Ell*Ex*E;
00052         SrcArithT A2 = El*El*Ex - Elx*El*E;
00053         SrcArithT A = A1 + 2*A2;
00054         SrcArithT B1 = Elly*E*E - Ell*Ey*E;
00055         SrcArithT B2 = El*El*Ey - Ely*El*E;
00056         SrcArithT B = B1 + 2*B2;
00057 
00058         SrcArithT Nllw = Element::E1Sqrt(A*A + B*B) / norm3;
00059         return Element::E1Sqrt(Nlw*Nlw + Nllw*Nllw);
00060         //return Element::E1Log(Element::E1Sqrt(Nlw*Nlw + Nllw*Nllw) + 1);
00061     }

Here is the call graph for this function: