Changeset f95556f in sasmodels
- Timestamp:
- Aug 1, 2016 8:54:35 PM (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:
- d119f34
- Parents:
- 33875e3
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/hollow_cylinder.c
r58210db rf95556f 6 6 double solvent_sld, double theta, double phi); 7 7 8 #define INVALID(v) (v.core_radius >= v.radius || v.radius >= v.length)8 #define INVALID(v) (v.core_radius >= v.radius) 9 9 10 10 // From Igor library … … 29 29 double q, double core_radius, double radius, double length, double dum) 30 30 { 31 //Note: lim_{r -> r_c} psi = J0(core_radius*qs)32 31 const double qs = q*sqrt(1.0-dum*dum); 33 32 const double lam1 = sas_J1c(radius*qs); 34 33 const double lam2 = sas_J1c(core_radius*qs); 35 34 const double gamma_sq = square(core_radius/radius); 35 //Note: lim_{r -> r_c} psi = J0(core_radius*qs) 36 36 const double psi = (lam1 - gamma_sq*lam2)/(1.0 - gamma_sq); //SRK 10/19/00 37 37 const double t2 = sinc(q*length*dum/2.0); 38 39 38 return square(psi*t2); 40 39 }
Note: See TracChangeset
for help on using the changeset viewer.