Changes in sasmodels/models/parallelepiped.c [d605080:9b7b23f] in sasmodels
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/parallelepiped.c
rd605080 r9b7b23f 1 double form_volume(double length_a, double length_b, double length_c); 2 double Iq(double q, double sld, double solvent_sld, 3 double length_a, double length_b, double length_c); 4 double Iqxy(double qx, double qy, double sld, double solvent_sld, 5 double length_a, double length_b, double length_c, 6 double theta, double phi, double psi); 7 8 double form_volume(double length_a, double length_b, double length_c) 1 static double 2 form_volume(double length_a, double length_b, double length_c) 9 3 { 10 4 return length_a * length_b * length_c; … … 12 6 13 7 14 double Iq(double q, 8 static double 9 Iq(double q, 15 10 double sld, 16 11 double solvent_sld, … … 20 15 { 21 16 const double mu = 0.5 * q * length_b; 22 17 23 18 // Scale sides by B 24 19 const double a_scaled = length_a / length_b; 25 20 const double c_scaled = length_c / length_b; 26 21 27 22 // outer integral (with gauss points), integration limits = 0, 1 28 23 double outer_total = 0; //initialize integral … … 57 52 58 53 59 double Iqxy(double qx, double qy, 54 static double 55 Iqxy(double qa, double qb, double qc, 60 56 double sld, 61 57 double solvent_sld, 62 58 double length_a, 63 59 double length_b, 64 double length_c, 65 double theta, 66 double phi, 67 double psi) 60 double length_c) 68 61 { 69 double q, xhat, yhat, zhat; 70 ORIENT_ASYMMETRIC(qx, qy, theta, phi, psi, q, xhat, yhat, zhat); 71 72 const double siA = sas_sinx_x(0.5*length_a*q*xhat); 73 const double siB = sas_sinx_x(0.5*length_b*q*yhat); 74 const double siC = sas_sinx_x(0.5*length_c*q*zhat); 62 const double siA = sas_sinx_x(0.5*length_a*qa); 63 const double siB = sas_sinx_x(0.5*length_b*qb); 64 const double siC = sas_sinx_x(0.5*length_c*qc); 75 65 const double V = form_volume(length_a, length_b, length_c); 76 66 const double drho = (sld - solvent_sld);
Note: See TracChangeset
for help on using the changeset viewer.