PolyhedralEdge Class Reference

Collaboration diagram for PolyhedralEdge:

Public Member Functions
	PolyhedralEdge (const kvector_t _Vlow, const kvector_t _Vhig)
kvector_t	E () const
kvector_t	R () const
complex_t	qE (cvector_t q) const
complex_t	qR (cvector_t q) const
complex_t	contrib (int m, cvector_t qpa, complex_t qrperp) const

Private Attributes
kvector_t	m_E
kvector_t	m_R

Detailed Description

One edge of a polygon, for form factor computation.

Definition at line 24 of file PolyhedralComponents.h.

Constructor & Destructor Documentation

PolyhedralEdge()

PolyhedralEdge::PolyhedralEdge	(	const kvector_t	_Vlow,
		const kvector_t	_Vhig
	)

Definition at line 31 of file PolyhedralComponents.cpp.

    : m_E((_Vhig - _Vlow) / 2), m_R((_Vhig + _Vlow) / 2)
{
    if (m_E.mag2() == 0)
        throw std::invalid_argument("At least one edge has zero length");
};

References m_E, and BasicVector3D< T >::mag2().

Here is the call graph for this function:

Member Function Documentation

E()

kvector_t PolyhedralEdge::E ( ) const

inline

Definition at line 29 of file PolyhedralComponents.h.

{ return m_E; }

References m_E.

Referenced by PolyhedralFace::edge_sum_ff(), and PolyhedralFace::ff_n_core().

R()

kvector_t PolyhedralEdge::R ( ) const

inline

Definition at line 30 of file PolyhedralComponents.h.

{ return m_R; }

References m_R.

qE()

complex_t PolyhedralEdge::qE ( cvector_t  q ) const

inline

Definition at line 31 of file PolyhedralComponents.h.

{ return m_E.dot(q); }

References BasicVector3D< T >::dot(), and m_E.

Referenced by contrib(), and PolyhedralFace::edge_sum_ff().

Here is the call graph for this function:

qR()

complex_t PolyhedralEdge::qR ( cvector_t  q ) const

inline

Definition at line 32 of file PolyhedralComponents.h.

{ return m_R.dot(q); }

References BasicVector3D< T >::dot(), and m_R.

Referenced by PolyhedralFace::edge_sum_ff().

Here is the call graph for this function:

contrib()

complex_t PolyhedralEdge::contrib	(	int	m,
		cvector_t	qpa,
		complex_t	qrperp
	)		const

Returns sum_l=0^M/2 u^2l v^(M-2l) / (2l+1)!(M-2l)! - vperp^M/M!

Definition at line 40 of file PolyhedralComponents.cpp.

{
    complex_t u = qE(qpa);
    complex_t v2 = m_R.dot(qpa);
    complex_t v1 = qrperp;
    complex_t v = v2 + v1;
#ifdef POLYHEDRAL_DIAGNOSTIC
    if (diagnosis.debmsg >= 5)
        std::cout << std::scientific << std::showpos << std::setprecision(16) << "contrib: u=" << u
                  << " v1=" << v1 << " v2=" << v2 << "\n";
#endif
    if (v == 0.) { // only 2l=M contributes
        if (M & 1) // M is odd
            return 0.;
        else
            return ReciprocalFactorialArray[M] * (pow(u, M) / (M + 1.) - pow(v1, M));
    }
    complex_t ret = 0;
    // the l=0 term, minus (qperp.R)^M, which cancels under the sum over E*contrib()
    if (v1 == 0.) {
        ret = ReciprocalFactorialArray[M] * pow(v2, M);
    } else if (v2 == 0.) {
        ; // leave ret=0
    } else {
        // binomial expansion
        for (int mm = 1; mm <= M; ++mm) {
            complex_t term = ReciprocalFactorialArray[mm] * ReciprocalFactorialArray[M - mm]
                             * pow(v2, mm) * pow(v1, M - mm);
            ret += term;
#ifdef POLYHEDRAL_DIAGNOSTIC
            if (diagnosis.debmsg >= 6)
                std::cout << "contrib mm=" << mm << " t=" << term << " s=" << ret << "\n";
#endif
        }
    }
    if (u == 0.)
        return ret;
    for (int l = 1; l <= M / 2; ++l) {
        complex_t term = ReciprocalFactorialArray[M - 2 * l] * ReciprocalFactorialArray[2 * l + 1]
                         * pow(u, 2 * l) * pow(v, M - 2 * l);
        ret += term;
#ifdef POLYHEDRAL_DIAGNOSTIC
        if (diagnosis.debmsg >= 6)
            std::cout << "contrib l=" << l << " t=" << term << " s=" << ret << "\n";
#endif
    }
    return ret;
}

References BasicVector3D< T >::dot(), m_R, qE(), anonymous_namespace{PolyhedralComponents.cpp}::ReciprocalFactorialArray, and StringUtils::scientific().

Referenced by PolyhedralFace::ff_n_core().

Here is the call graph for this function:

Member Data Documentation

m_E

kvector_t PolyhedralEdge::m_E

private

vector pointing from mid of edge to upper vertex

Definition at line 37 of file PolyhedralComponents.h.

Referenced by E(), PolyhedralEdge(), and qE().

m_R

kvector_t PolyhedralEdge::m_R

private

position vector of edge midpoint

Definition at line 38 of file PolyhedralComponents.h.

Referenced by contrib(), qR(), and R().

---

The documentation for this class was generated from the following files:

• /home/www/ba/Sample/HardParticle/PolyhedralComponents.h
• /home/www/ba/Sample/HardParticle/PolyhedralComponents.cpp