Changeset ff7119b in sasmodels for sasmodels/models/cylinder.c
- Timestamp:
- Aug 26, 2014 10:27:06 PM (10 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:
- 5d4777d
- Parents:
- a7684e5
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/cylinder.c
ra7684e5 rff7119b 14 14 real length) 15 15 { 16 const real h = REAL(0.5)*length;16 const real halflength = REAL(0.5)*length; 17 17 real summ = REAL(0.0); 18 // real lower=0, upper=M_PI_2; 18 19 for (int i=0; i<76 ;i++) { 19 //const real zi = ( Gauss76Z[i]*(uplim-lolim) + uplim + lolim )/2.0; 20 // translate a point in [-1,1] to a point in [lower,upper] 21 //const real zi = ( Gauss76Z[i]*(upper-lower) + upper + lower )/2.0; 20 22 const real zi = REAL(0.5)*(Gauss76Z[i]*M_PI_2 + M_PI_2); 21 summ += Gauss76Wt[i] * CylKernel(q, radius, h , zi);23 summ += Gauss76Wt[i] * CylKernel(q, radius, halflength, zi); 22 24 } 23 //const real form = (uplim-lolim)/2.0*summ; 25 // translate dx in [-1,1] to dx in [lower,upper] 26 //const real form = (upper-lower)/2.0*summ; 24 27 const real form = summ * M_PI_4; 25 28
Note: See TracChangeset
for help on using the changeset viewer.