RPA10Model.py @ 35aface

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload

Last change on this file since 35aface was 35aface, checked in by Jae Cho <jhjcho@…>, 14 years ago
addede new models and attr. non_fittable
Property mode set to `100644`
File size: 8.6 KB

Line
1
2	from sans.models.BaseComponent import BaseComponent
3	from sans.models.RPAModel import RPAModel
4	from copy import deepcopy
5	max_case_n = 10
6	class RPA10Model(BaseComponent):
7	"""
8	This multi-model is based on Parratt formalism and provides the capability
9	of changing the number of layers between 0 and 10.
10	"""
11	def __init__(self, multfactor=1):
12	BaseComponent.__init__(self)
13	"""
14	:param multfactor: number of cases in the model, assumes 0<= case# <=10.
15	"""
16
17	## Setting model name model description
18	self.description=""
19	model = RPAModel()
20	self.model = model
21	self.name = "RPA10Model"
22	self.description=model.description
23	self.case_num = multfactor
24	## Define parameters
25	self.params = {}
26
27	## Parameter details [units, min, max]
28	self.details = {}
29
30	# non-fittable parameters
31	self.non_fittable = model.non_fittable
32
33	# list of function in order of the function number
34	self.fun_list = self._get_func_list()
35	## dispersion
36	self._set_dispersion()
37	## Define parameters
38	self._set_params()
39
40	## Parameter details [units, min, max]
41	self._set_details()
42
43	#list of parameter that can be fitted
44	self._set_fixed_params()
45	self.model.params['lcase_n'] = self.case_num
46
47	## functional multiplicity of the model
48	self.multiplicity = max_case_n
49
50	def _clone(self, obj):
51	"""
52	Internal utility function to copy the internal
53	data members to a fresh copy.
54	"""
55	obj.params = deepcopy(self.params)
56	obj.non_fittable = deepcopy(self.non_fittable)
57	obj.description = deepcopy(self.description)
58	obj.details = deepcopy(self.details)
59	obj.dispersion = deepcopy(self.dispersion)
60	obj.model = self.model.clone()
61
62	return obj
63
64
65	def _set_dispersion(self):
66	"""
67	model dispersions
68	"""
69	##set dispersion from model
70	self.dispersion = {}
71
72
73	def _set_params(self):
74	"""
75	Concatenate the parameters of the model to create
76	this model parameters
77	"""
78	# set the case number first
79	self.model.params['lcase_n'] = self.case_num
80	# rearrange the parameters for the given # of shells
81	for name , value in self.model.params.iteritems():
82	if name == 'lcase_n':
83	continue
84	elif name.lower() == 'scale' or \
85	name.lower() == 'background':
86	self.params[name] = value
87	elif self.case_num <= 1:
88	if name.lower() == 'nc' or \
89	name.lower() == 'phic' or \
90	name.lower() == 'vc' or \
91	name.lower() == 'lc' or \
92	name.lower() == 'bc' or \
93	name.lower() == 'nd' or \
94	name.lower() == 'phid' or \
95	name.lower() == 'vd' or \
96	name.lower() == 'ld' or \
97	name.lower() == 'bd' or \
98	name.lower() == 'kcd':
99	self.params[name] = value
100	elif self.case_num > 1 and self.case_num <= 4:
101	if name.lower() == 'nb' or \
102	name.lower() == 'phib' or \
103	name.lower() == 'vb' or \
104	name.lower() == 'lb' or \
105	name.lower() == 'bb' or \
106	name.lower() == 'nc' or \
107	name.lower() == 'phic' or \
108	name.lower() == 'vc' or \
109	name.lower() == 'lc' or \
110	name.lower() == 'bc' or \
111	name.lower() == 'nd' or \
112	name.lower() == 'phid' or \
113	name.lower() == 'vd' or \
114	name.lower() == 'ld' or \
115	name.lower() == 'bd' or \
116	name.lower() == 'kbc' or \
117	name.lower() == 'kbd' or \
118	name.lower() == 'kcd' :
119	self.params[name] = value
120	else:
121	if name.lower() == 'na' or \
122	name.lower() == 'phia' or \
123	name.lower() == 'va' or \
124	name.lower() == 'la' or \
125	name.lower() == 'ba' or \
126	name.lower() == 'nb' or \
127	name.lower() == 'phib' or \
128	name.lower() == 'vb' or \
129	name.lower() == 'lb' or \
130	name.lower() == 'bb' or \
131	name.lower() == 'nc' or \
132	name.lower() == 'phic' or \
133	name.lower() == 'vc' or \
134	name.lower() == 'lc' or \
135	name.lower() == 'bc' or \
136	name.lower() == 'nd' or \
137	name.lower() == 'phid' or \
138	name.lower() == 'vd' or \
139	name.lower() == 'ld' or \
140	name.lower() == 'bd' or \
141	name.lower() == 'kab' or \
142	name.lower() == 'kac' or \
143	name.lower() == 'kad' or \
144	name.lower() == 'kbc' or \
145	name.lower() == 'kbd' or \
146	name.lower() == 'kcd' :
147	self.params[name] = value
148
149
150
151	def _set_details(self):
152	"""
153	Concatenate details of the original model to create
154	this model details
155	"""
156	for name ,detail in self.model.details.iteritems():
157	if name in self.params.iterkeys():
158	self.details[name]= detail
159
160	def _get_func_list(self):
161	"""
162	Get the list of functions in each cases
163	"""
164	func_list = {}
165	return func_list
166
167	def getProfile(self):
168	"""
169	Get SLD profile
170
171	: return: None, No SLD profile supporting for this model
172	"""
173	return None
174
175	def setParam(self, name, value):
176	"""
177	Set the value of a model parameter
178
179	: param name: name of the parameter
180	: param value: value of the parameter
181	"""
182	# set param to new model
183	self._setParamHelper( name, value)
184
185	self.model.setParam( name, value)
186
187	def _setParamHelper(self, name, value):
188	"""
189	Helper function to setParam
190	"""
191
192	# Look for standard parameter
193	for item in self.params.keys():
194	if item.lower()==name.lower():
195	self.params[item] = value
196	return
197
198	raise ValueError, "Model does not contain parameter %s" % name
199
200
201	def _set_fixed_params(self):
202	"""
203	Fill the self.fixed list with the model fixed list
204	"""
205	pass
206
207	def run(self, x = 0.0):
208	"""
209	Evaluate the model
210
211	:param x: input q, or [q,phi]
212
213	:return: scattering function P(q)
214
215	"""
216
217	return self.model.run(x)
218
219	def runXY(self, x = 0.0):
220	"""
221	Evaluate the model
222
223	: param x: input q-value (float or [float, float] as [qx, qy])
224	: return: scattering function value
225	"""
226
227	return self.model.runXY(x)
228
229	## Now (May27,10) directly uses the model eval function
230	## instead of the for-loop in Base Component.
231	def evalDistribution(self, x = []):
232	"""
233	Evaluate the model in cartesian coordinates
234
235	: param x: input q[], or [qx[], qy[]]
236	: return: scattering function P(q[])
237	"""
238	# set effective radius and scaling factor before run
239	return self.model.evalDistribution(x)
240	def calculate_ER(self):
241	"""
242	"""
243	return self.model.calculate_ER()
244	def set_dispersion(self, parameter, dispersion):
245	"""
246	Set the dispersion object for a model parameter
247
248	: param parameter: name of the parameter [string]
249	:dispersion: dispersion object of type DispersionModel
250	"""
251	pass

Note: See TracBrowser for help on using the repository browser.

SasView

source: sasview/sansmodels/src/sans/models/RPA10Model.py @ 35aface

Download in other formats: