source: sasview/park_integration/test/testfitting.py @ be72312c

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 be72312c was 4dd63eb, checked in by Gervaise Alina <gervyh@…>, 16 years ago

change on set_model

  • Property mode set to 100644
File size: 3.4 KB
Line 
1"""
2    Unit tests for fitting module
3"""
4import unittest
5from sans.guitools.plottables import Theory1D
6from sans.guitools.plottables import Data1D
7
8import math
9class testFitModule(unittest.TestCase):
10    """ test fitting """
11    def testLoader(self):
12        """
13            test module Load
14        """
15        from sans.fit.Loader import Load
16        load= Load()
17       
18        load.set_filename("testdata_line.txt")
19        self.assertEqual(load.get_filename(),"testdata_line.txt")
20        load.set_values()
21        x=[]
22        y=[]
23        dx=[]
24        dy=[]
25       
26        x,y,dx,dy = load.get_values()
27       
28        # test that values have been loaded
29        self.assertNotEqual(x, None)
30        self.assertNotEqual(y, [])
31        self.assertNotEqual(dy, None)
32        self.assertEqual(len(x),len(y))
33        self.assertEqual(len(dy),len(y))
34       
35        # test data the two plottables contained values loaded
36        data1 = Theory1D(x=[], y=[], dy=None)
37        data2 = Data1D(x=[], y=[],dx=None, dy=None)
38        data1.name = "data1"
39        data2.name = "data2"
40        load.load_data(data1)
41        load.load_data(data2)
42       
43       
44        for i in range(len(x)):
45            self.assertEqual(data2.x[i],x[i])
46            self.assertEqual(data1.y[i],y[i])
47            self.assertEqual(data2.y[i],y[i])
48            self.assertEqual(data1.dx[i],dx[i])
49            self.assertEqual(data2.dy[i],dy[i])
50            self.assertEqual(data1.x[i],data2.x[i])
51            self.assertEqual(data2.y[i],data2.y[i])
52           
53           
54    def testfit_1Data_1Model(self):
55        """ test fitting for one data and one model park vs scipy"""
56        #load data
57        from sans.fit.Loader import Load
58        load= Load()
59        load.set_filename("testdata_line.txt")
60        load.set_values()
61        data1 = Data1D(x=[], y=[],dx=None, dy=None)
62        load.load_data(data1)
63       
64        #Importing the Fit module
65        from sans.fit.Fitting import Fit
66        fitter= Fit('scipy')
67       
68        # Receives the type of model for the fitting
69        from sans.guitools.LineModel import LineModel
70        model1  = LineModel()
71        model2  = LineModel()
72       
73        #Do the fit SCIPY
74        fitter.set_data(data1,1)
75        fitter.set_model(model1,"M1",1,{'A':2,'B':1})
76       
77        chisqr1, out1, cov1=fitter.fit()
78        """ testing SCIPy results"""
79        self.assert_(math.fabs(out1[1]-2.5)/math.sqrt(cov1[1][1]) < 2)
80        self.assert_(math.fabs(out1[0]-4.0)/math.sqrt(cov1[0][0]) < 2)
81        self.assert_(chisqr1/len(data1.x) < 2)
82        # PARK
83        fitter= Fit('park')
84       
85        #Do the fit
86        fitter.set_data(data1,1)
87        fitter.set_model(model2,"M1",1,{'A':2,'B':1})
88       
89        chisqr2, out2, cov2=fitter.fit(None,None)
90       
91        self.assert_(math.fabs(out2[1]-2.5)/math.sqrt(cov2[1][1]) < 2)
92        self.assert_(math.fabs(out2[0]-4.0)/math.sqrt(cov2[0][0]) < 2)
93        self.assert_(chisqr2/len(data1.x) < 2)
94        print "scipy",chisqr1, out1, cov1
95        print "park",chisqr2, out2, cov2
96        self.assertAlmostEquals(out1[1], out2[1],0)
97        self.assertAlmostEquals(out1[0], out2[0],0)
98        self.assertAlmostEquals(cov1[0][0], cov2[0][0],1)
99        self.assertAlmostEquals(cov1[1][1], cov2[1][1],1)
100        self.assertAlmostEquals(chisqr1, chisqr2)
101       
102     
Note: See TracBrowser for help on using the repository browser.