Changeset 2e2452b in sasview for sansmodels/src/sans/models/c_smearer/smearer.cpp
- Timestamp:
- Dec 9, 2011 5:15:40 PM (13 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:
- 2371363
- Parents:
- 0e8b249
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sansmodels/src/sans/models/c_smearer/smearer.cpp
r09e89b7 r2e2452b 169 169 // Find Weights 170 170 // Find q where the resolution smearing calculation of I(q) occurs 171 double q, q_min, q_max, q_0 ;171 double q, q_min, q_max, q_0=0.0; 172 172 get_bin_range(i, &q, &q_min, &q_max); 173 173 // Block q becomes <=0 … … 175 175 continue; 176 176 } 177 bool last_qpoint = true;178 177 // Find q[0] value to normalize the weight later, 179 178 // otherwise, we will have a precision problem. … … 182 181 } 183 182 // Loop over all qj-values 184 bool first_w = true;185 183 for(int j=0; j<nbins; j++) { 186 184 double q_j, q_high, q_low; … … 263 261 // Later, it will be normalized again. 264 262 265 double q_shift_min = q - width;266 267 double u = (q_j * q_j - (q_shift) * (q_shift));268 263 // The fabs below are not necessary but in case: the weight should never be imaginary. 269 264 // At the edge of each sub_width. weight += u(at q_high bin) - u(0), where u(0) = 0, … … 292 287 293 288 // Loop over all q-values 294 double step = (qmax-qmin)/((double)nbins-1.0);295 289 double q, q_min, q_max; 296 290 double q_j, q_jmax, q_jmin;
Note: See TracChangeset
for help on using the changeset viewer.