Changeset 2cb89e7 in sasview for sansmodels/src/sans/models/c_extensions/parallelepiped.c
- Timestamp:
- Aug 16, 2009 10:30:27 AM (15 years ago)
- Branches:
- master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload
- Children:
- 7d11b81
- Parents:
- b0d0723
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sansmodels/src/sans/models/c_extensions/parallelepiped.c
rd5bd424 r2cb89e7 22 22 // Fill paramater array 23 23 dp[0] = pars->scale; 24 dp[1] = pars->short_ edgeA;25 dp[2] = pars->long er_edgeB;26 dp[3] = pars->long uest_edgeC;24 dp[1] = pars->short_a; 25 dp[2] = pars->long_b; 26 dp[3] = pars->longer_c; 27 27 dp[4] = pars->contrast; 28 28 dp[5] = pars->background; … … 38 38 39 39 //handle arg=0 separately, as sin(t)/t -> 1 as t->0 40 argA = a*ala ;41 argB = b*alb ;42 argC = c*alc ;40 argA = a*ala/2; 41 argB = b*alb/2; 42 argC = c*alc/2; 43 43 if(argA==0.0) { 44 44 tmp1 = 1.0; … … 105 105 double pi = 4.0*atan(1.0); 106 106 107 edgeA = pars->short_ edgeA;108 edgeB = pars->long er_edgeB;109 edgeC = pars->long uest_edgeC;107 edgeA = pars->short_a; 108 edgeB = pars->long_b; 109 edgeC = pars->longer_c; 110 110 111 111 … … 152 152 //normalize by cylinder volume 153 153 //NOTE that for this (Fournet) definition of the integral, one must MULTIPLY by Vparallel 154 vol = 8*edgeA* edgeB * edgeC;154 vol = edgeA* edgeB * edgeC; 155 155 answer *= vol; 156 156
Note: See TracChangeset
for help on using the changeset viewer.