Changeset 0d6e865 in sasmodels for sasmodels/models/capped_cylinder.c
- Timestamp:
- Oct 11, 2016 1:42:00 PM (8 years ago)
- Branches:
- master, core_shell_microgels, costrafo411, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- 1fdb555
- Parents:
- 19e7ca7
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/capped_cylinder.c
r2222134 r0d6e865 116 116 // Compute angle alpha between q and the cylinder axis 117 117 double sn, cn; 118 SINCOS( theta*M_PI_180, sn, cn);118 SINCOS(phi*M_PI_180, sn, cn); 119 119 const double q = sqrt(qx*qx+qy*qy); 120 const double cos_val = cn*cos(phi*M_PI_180)*(qx/q) + sn*(qy/q);120 const double cos_val = (q==0. ? 1.0 : (cn*qx + sn*qy)*sin(theta*M_PI_180)/q); 121 121 const double alpha = acos(cos_val); // rod angle relative to q 122 122
Note: See TracChangeset
for help on using the changeset viewer.