Changeset a807206 in sasmodels for sasmodels/models/rectangular_prism.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/rectangular_prism.c
rdeb7ee0 ra807206 1 double form_volume(double a_side, double b2a_ratio, double c2a_ratio);2 double Iq(double q, double sld, double solvent_sld, double a_side,1 double form_volume(double length_a, double b2a_ratio, double c2a_ratio); 2 double Iq(double q, double sld, double solvent_sld, double length_a, 3 3 double b2a_ratio, double c2a_ratio); 4 4 double Iqxy(double qx, double qy, double sld, double solvent_sld, 5 double a_side, double b2a_ratio, double c2a_ratio);5 double length_a, double b2a_ratio, double c2a_ratio); 6 6 7 double form_volume(double a_side, double b2a_ratio, double c2a_ratio)7 double form_volume(double length_a, double b2a_ratio, double c2a_ratio) 8 8 { 9 return a_side * (a_side*b2a_ratio) * (a_side*c2a_ratio);9 return length_a * (length_a*b2a_ratio) * (length_a*c2a_ratio); 10 10 } 11 11 … … 13 13 double sld, 14 14 double solvent_sld, 15 double a_side,15 double length_a, 16 16 double b2a_ratio, 17 17 double c2a_ratio) … … 19 19 double termA, termB, termC; 20 20 21 double b_side = a_side* b2a_ratio;22 double c_side = a_side* c2a_ratio;23 double volume = a_side* b_side * c_side;24 double a_half = 0.5 * a_side;21 double b_side = length_a * b2a_ratio; 22 double c_side = length_a * c2a_ratio; 23 double volume = length_a * b_side * c_side; 24 double a_half = 0.5 * length_a; 25 25 double b_half = 0.5 * b_side; 26 26 double c_half = 0.5 * c_side; … … 92 92 double sld, 93 93 double solvent_sld, 94 double a_side,94 double length_a, 95 95 double b2a_ratio, 96 96 double c2a_ratio) 97 97 { 98 98 double q = sqrt(qx*qx + qy*qy); 99 double intensity = Iq(q, sld, solvent_sld, a_side, b2a_ratio, c2a_ratio);99 double intensity = Iq(q, sld, solvent_sld, length_a, b2a_ratio, c2a_ratio); 100 100 return intensity; 101 101 }
Note: See TracChangeset
for help on using the changeset viewer.