Changeset 34f3ad0 in sasview for pr_inversion/src/sans/pr/distance_explorer.py
- Timestamp:
- Apr 27, 2012 9:07:00 AM (12 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:
- 8a621ac
- Parents:
- 1f9f3c8a
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
pr_inversion/src/sans/pr/distance_explorer.py
rd9621b17 r34f3ad0 16 16 over that range. 17 17 """ 18 import sys 18 19 19 import numpy20 import math21 import sys22 20 23 21 class Results: … … 40 38 self.bck = [] 41 39 self.d_max = [] 42 ## List of errors found during the last exploration 40 ## List of errors found during the last exploration 43 41 self.errors = [] 42 44 43 45 44 class DistExplorer(object): … … 55 54 56 55 """ 57 self.pr_state 56 self.pr_state = pr_state 58 57 self._default_min = 0.8 * self.pr_state.d_max 59 58 self._default_max = 1.2 * self.pr_state.d_max 60 61 59 62 60 def __call__(self, dmin=None, dmax=None, npts=10): … … 80 78 81 79 # Loop over d_max values 82 for i in range(npts): 83 d = dmin + i * (dmax - dmin) /(npts-1.0)80 for i in range(npts): 81 d = dmin + i * (dmax - dmin) / (npts - 1.0) 84 82 self.pr_state.d_max = d 85 83 try: 86 out, cov = self.pr_state.invert(self.pr_state.nfunc) 84 out, cov = self.pr_state.invert(self.pr_state.nfunc) 87 85 88 86 # Store results … … 100 98 results.pos.append(pos) 101 99 results.pos_err.append(pos_err) 102 results.osc.append(osc) 100 results.osc.append(osc) 103 101 except: 104 102 # This inversion failed, skip this D_max value … … 108 106 109 107 return results 110
Note: See TracChangeset
for help on using the changeset viewer.