someff Namespace Reference

Functions
complex_t	ffSphere (cvector_t q, double R)

Detailed Description

Some form factor functions.

Function Documentation

ffSphere()

complex_t someff::ffSphere	(	cvector_t	q,
		double	R
	)

Returns the form factor of a sphere of radius R.

Used by the hard sphere and by several soft sphere classes.

Definition at line 21 of file SomeFormFactors.cpp.

{
    complex_t q1 =
        sqrt(q.x() * q.x() + q.y() * q.y() + q.z() * q.z()); // NO sesquilinear dot product!
    complex_t qR = q1 * R;
 
    if (std::abs(qR) < 1e-4) { // relative error is O(qR^4) with small prefactor
        // expand sin(qR)-qR*cos(qR) up to qR^5
        return 4 * M_PI / 3 * pow(R, 3) * (1. - 0.1 * pow(qR, 2));
    }
    return 4 * M_PI * pow(q1, -3) * (sin(qR) - qR * cos(qR));
}

References M_PI, BasicVector3D< T >::x(), BasicVector3D< T >::y(), and BasicVector3D< T >::z().

Referenced by FormFactorFullSphere::evaluate_for_q(), FormFactorSphereGaussianRadius::evaluate_for_q(), and FormFactorSphereLogNormalRadius::evaluate_for_q().

Here is the call graph for this function: