star-sf

Set of surface and volume scattering functions
git clone git://git.meso-star.fr/star-sf.git
Log | Files | Refs | README | LICENSE

commit ed136ddade56f9404af56bc12a60f3e0dd94af6b
parent f0830178e6791f9e21d5c4833ecd410839a1d668
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Mon,  8 Feb 2021 10:18:11 +0100

Remove dead code from the RDGFA phase function

Diffstat:
Msrc/ssf_phase_rdgfa.c | 29-----------------------------
1 file changed, 0 insertions(+), 29 deletions(-)

diff --git a/src/ssf_phase_rdgfa.c b/src/ssf_phase_rdgfa.c @@ -101,35 +101,6 @@ compute_interval_angular_range } } -#if 0 -/* Not normalized */ -static INLINE double -eval__(struct rdgfa* rdgfa, const double theta) -{ - /* Input arguments */ - const double lambda = rdgfa->wavelength; - const double Df = rdgfa->fractal_dimension; - const double Rg = rdgfa->gyration_radius; - const double Rg2 = Rg*Rg; - - /* TODO precompute constants */ - const double k = 2*PI / lambda; /* [nm^-1] */ - const double k2 = k*k; - - const double g = pow(1 + 4*k2*Rg2/(3*Df), -Df*0.5); - - const double q = 4*PI / lambda * sin(theta*0.5); - const double qRg = q * Rg; - const double qRg2 = qRg*qRg; - const double f = qRg2 < 1.5*Df - ? exp(-1.0/3.0 * qRg2) - : pow(3.0*Df/(2.0*EXP1) * 1.0/qRg2, Df*0.5); - - ASSERT(d3_is_normalized(wo) && d3_is_normalized(wi)); - return f; -} -#endif - static INLINE double eval_f(struct rdgfa* rdgfa, const double theta) {