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 692d726 was
af4af5a,
checked in by Gervaise Alina <gervyh@…>, 16 years ago
|
make modification on testing
|
-
Property mode set to
100644
|
File size:
1.7 KB
|
Line | |
---|
1 | """ |
---|
2 | Unit tests for fitting module |
---|
3 | """ |
---|
4 | import unittest |
---|
5 | from sans.guitools.plottables import Theory1D |
---|
6 | from sans.guitools.plottables import Data1D |
---|
7 | from sans.fit.AbstractFitEngine import Data, Model,FitData1D |
---|
8 | import math |
---|
9 | from sans.fit.Fitting import Fit |
---|
10 | from DataLoader.loader import Loader |
---|
11 | |
---|
12 | class testFitModule(unittest.TestCase): |
---|
13 | """ test fitting """ |
---|
14 | def test_park(self): |
---|
15 | """ Simple cylinder model fit (park) """ |
---|
16 | |
---|
17 | out=Loader().load("cyl_400_20.txt") |
---|
18 | data1 = Data1D(x=out.x, y=out.y, dx=out.dx, dy=out.y) |
---|
19 | |
---|
20 | fitter = Fit('park') |
---|
21 | # Receives the type of model for the fitting |
---|
22 | from sans.models.CylinderModel import CylinderModel |
---|
23 | model1 = CylinderModel() |
---|
24 | model1.setParam('contrast', 1) |
---|
25 | data1.smearer = None |
---|
26 | data = FitData1D(data1) |
---|
27 | model = Model(model1) |
---|
28 | |
---|
29 | pars1 =['length','radius','scale'] |
---|
30 | fitter.set_data(data,1) |
---|
31 | model.set( scale=1e-10 ) |
---|
32 | fitter.set_model(model,1,pars1) |
---|
33 | fitter.select_problem_for_fit(Uid=1,value=1) |
---|
34 | result1 = fitter.fit() |
---|
35 | |
---|
36 | self.assert_(result1) |
---|
37 | self.assertTrue(len(result1.pvec)>0 or len(result1.pvec)==0 ) |
---|
38 | self.assertTrue(len(result1.stderr)> 0 or len(result1.stderr)==0) |
---|
39 | |
---|
40 | #print result1.pvec[0]-400.0, result1.pvec[0] |
---|
41 | #print math.fabs(result1.pvec[0]-400.0)/3.0 |
---|
42 | self.assertTrue( math.fabs(result1.pvec[0]-400.0)/3.0 < result1.stderr[0] ) |
---|
43 | self.assertTrue( math.fabs(result1.pvec[1]-20.0)/3.0 < result1.stderr[1] ) |
---|
44 | self.assertTrue( math.fabs(result1.pvec[2]-9.0e-12)/3.0 < result1.stderr[2] ) |
---|
45 | self.assertTrue( result1.fitness < 1.0 ) |
---|
46 | |
---|
Note: See
TracBrowser
for help on using the repository browser.