static double form_volume(double length_a, double b2a_ratio, double c2a_ratio) { double length_b = length_a * b2a_ratio; double length_c = length_a * c2a_ratio; double vol_shell = 2.0 * (length_a*length_b + length_a*length_c + length_b*length_c); return vol_shell; } static void Fq(double q, double *F1, double *F2, double sld, double solvent_sld, double length_a, double b2a_ratio, double c2a_ratio) { const double length_b = length_a * b2a_ratio; const double length_c = length_a * c2a_ratio; const double a_half = 0.5 * length_a; const double b_half = 0.5 * length_b; const double c_half = 0.5 * length_c; //Integration limits to use in Gaussian quadrature const double v1a = 0.0; const double v1b = M_PI_2; //theta integration limits const double v2a = 0.0; const double v2b = M_PI_2; //phi integration limits double outer_sum_F1 = 0.0; double outer_sum_F2 = 0.0; for(int i=0; i