Changeset a807206 in sasmodels for sasmodels/models/parallelepiped.c
- Timestamp:
- Oct 1, 2016 12:42:06 AM (8 years ago)
- Branches:
- master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- caddb14, 5031ca3
- Parents:
- 2222134
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/parallelepiped.c
rdeb7ee0 ra807206 1 double form_volume(double a_side, double b_side, double c_side);2 double Iq(double q, double sld, double solvent_sld, double a_side, double b_side, double c_side);1 double form_volume(double length_a, double length_b, double length_c); 2 double Iq(double q, double sld, double solvent_sld, double length_a, double length_b, double length_c); 3 3 double Iqxy(double qx, double qy, double sld, double solvent_sld, 4 double a_side, double b_side, double c_side, double theta, double phi, double psi);4 double length_a, double length_b, double length_c, double theta, double phi, double psi); 5 5 6 6 // From Igor library … … 32 32 33 33 34 double form_volume(double a_side, double b_side, double c_side)34 double form_volume(double length_a, double length_b, double length_c) 35 35 { 36 return a_side * b_side * c_side;36 return length_a * length_b * length_c; 37 37 } 38 38 … … 41 41 double sld, 42 42 double solvent_sld, 43 double a_side,44 double b_side,45 double c_side)43 double length_a, 44 double length_b, 45 double length_c) 46 46 { 47 47 double tmp1, tmp2; 48 48 49 double mu = q * b_side;49 double mu = q * length_b; 50 50 51 51 // Scale sides by B 52 double a_scaled = a_side / b_side;53 double c_scaled = c_side / b_side;52 double a_scaled = length_a / length_b; 53 double c_scaled = length_c / length_b; 54 54 55 55 //Order of integration … … 102 102 } 103 103 104 const double vd = (sld-solvent_sld) * form_volume( a_side, b_side, c_side);104 const double vd = (sld-solvent_sld) * form_volume(length_a, length_b, length_c); 105 105 106 106 // convert from [1e-12 A-1] to [cm-1] and 0.5 factor for outer integral … … 113 113 double sld, 114 114 double solvent_sld, 115 double a_side,116 double b_side,117 double c_side,115 double length_a, 116 double length_b, 117 double length_c, 118 118 double theta, 119 119 double phi, … … 161 161 162 162 // Call the IGOR library function to get the kernel 163 double form = _pkernel( q* a_side, q*b_side, q*c_side, cos_val_a, cos_val_b, cos_val_c);163 double form = _pkernel( q*length_a, q*length_b, q*length_c, cos_val_a, cos_val_b, cos_val_c); 164 164 165 165 // Multiply by contrast^2 166 const double vd = (sld - solvent_sld) * form_volume( a_side, b_side, c_side);166 const double vd = (sld - solvent_sld) * form_volume(length_a, length_b, length_c); 167 167 return 1.0e-4 * vd * vd * form; 168 168 }
Note: See TracChangeset
for help on using the changeset viewer.