Changeset 4dd63eb in sasview for park_integration/test
- Timestamp:
- Jul 10, 2008 4:05:09 PM (17 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:
- 49b3ddd
- Parents:
- 83ca047
- Location:
- park_integration/test
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
park_integration/test/constrainttestpark.py
r792db7d5 r4dd63eb 27 27 #Importing the Fit module 28 28 from sans.fit.Fitting import Fit 29 fitter= Fit( )29 fitter= Fit('park') 30 30 # Receives the type of model for the fitting 31 31 from sans.guitools.LineModel import LineModel 32 32 model1 = LineModel() 33 33 model2 = LineModel() 34 #set engine for scipy 35 fitter.fit_engine('park') 36 engine = fitter.returnEngine() 34 37 35 #Do the fit 38 engine.set_param( model1,"M1", {'A':2,'B':3})39 engine.set_model(model1,1)40 engine.set_data(data1,1)41 #print "m1 param b", M1.B42 engine.set_param( model2,"M2", {'A':'M1.A','B':'M1.B'})43 engine.set_model(model2,2)44 36 45 engine.set_data(data2,2) 37 fitter.set_model(model1,"M1",1, {'A':1,'B':2}) 38 fitter.set_data(data1,1) 39 40 fitter.set_model(model2,"M2",2, {'A':'M1.A','B':'M1.B'}) 41 fitter.set_data(data2,2) 46 42 47 43 48 chisqr2, out2, cov2= engine.fit({'A':2,'B':1},None,None)44 chisqr2, out2, cov2= fitter.fit() 49 45 print "chisqr2",chisqr2 50 46 print "out2", out2 -
park_integration/test/testfitting.py
r792db7d5 r4dd63eb 5 5 from sans.guitools.plottables import Theory1D 6 6 from sans.guitools.plottables import Data1D 7 from sans.fit.ScipyFitting import Parameter 7 8 8 import math 9 9 class testFitModule(unittest.TestCase): … … 64 64 #Importing the Fit module 65 65 from sans.fit.Fitting import Fit 66 fitter= Fit() 67 fitter.fit_engine('scipy') 68 engine = fitter.returnEngine() 66 fitter= Fit('scipy') 67 69 68 # Receives the type of model for the fitting 70 69 from sans.guitools.LineModel import LineModel … … 73 72 74 73 #Do the fit SCIPY 75 engine.set_data(data1,1) 76 #engine.set_param( model1,"M1", {'A':2,'B':4}) 77 engine.set_model(model1,1) 74 fitter.set_data(data1,1) 75 fitter.set_model(model1,"M1",1,{'A':2,'B':1}) 78 76 79 chisqr1, out1, cov1= engine.fit({'A':2,'B':1},None,None)77 chisqr1, out1, cov1=fitter.fit() 80 78 """ testing SCIPy results""" 81 79 self.assert_(math.fabs(out1[1]-2.5)/math.sqrt(cov1[1][1]) < 2) … … 83 81 self.assert_(chisqr1/len(data1.x) < 2) 84 82 # PARK 85 fitter= Fit() 86 fitter.fit_engine('park') 87 engine = fitter.returnEngine() 83 fitter= Fit('park') 88 84 89 85 #Do the fit 90 engine.set_data(data1,1) 91 engine.set_param( model2,"M1", {'A':2,'B':1}) 92 engine.set_model(model2,1) 86 fitter.set_data(data1,1) 87 fitter.set_model(model2,"M1",1,{'A':2,'B':1}) 93 88 94 chisqr2, out2, cov2= engine.fit({'A':2,'B':1},None,None)89 chisqr2, out2, cov2=fitter.fit(None,None) 95 90 96 91 self.assert_(math.fabs(out2[1]-2.5)/math.sqrt(cov2[1][1]) < 2) -
park_integration/test/testpark.py
r792db7d5 r4dd63eb 33 33 #Importing the Fit module 34 34 from sans.fit.Fitting import Fit 35 fitter= Fit( )35 fitter= Fit('park') 36 36 # Receives the type of model for the fitting 37 37 from sans.guitools.LineModel import LineModel 38 38 model1 = LineModel() 39 39 model2 = LineModel() 40 #set engine for scipy 41 fitter.fit_engine('park') 42 engine = fitter.returnEngine() 40 43 41 #Do the fit 44 engine.set_param( model1,"M1", {'A':2.5,'B':4}) 45 engine.set_model(model1,1) 46 engine.set_data(data1,1) 42 fitter.set_model(model1,"M1",1, {'A':2.5,'B':4}) 43 fitter.set_data(data1,1) 47 44 48 engine.set_param( model2,"M2", {'A':2,'B':4}) 49 engine.set_model(model2,2) 50 engine.set_data(data2,2) 45 fitter.set_model(model2,"M2",2, {'A':2,'B':3}) 46 fitter.set_data(data2,2) 51 47 52 chisqr1, out1, cov1= engine.fit({'A':2,'B':1},None,None)48 chisqr1, out1, cov1= fitter.fit() 53 49 54 50 self.assert_(math.fabs(out1[1]-2.5)/math.sqrt(cov1[1][1]) < 2) … … 63 59 64 60 65 engine.set_data(data3,1)66 chisqr2, out2, cov2= engine.fit({'A':2,'B':1},None,None)61 fitter.set_data(data3,1) 62 chisqr2, out2, cov2= fitter.fit(None,None) 67 63 self.assert_(math.fabs(out2[1]-2.5)/math.sqrt(cov2[1][1]) < 2) 68 64 print math.fabs(out2[0]-4.0)/math.sqrt(cov2[0][0]) … … 75 71 #self.assert_(chisqr2/len(data2.x) < 2) 76 72 73 fitter.remove_Fit_Problem(2) 77 74 78 79 engine.remove_Fit_Problem(2) 80 chisqr3, out3, cov3= engine.fit({'A':2,'B':1},None,None) 75 chisqr3, out3, cov3= fitter.fit() 81 76 #print "park",chisqr3, out3, cov3 82 77 self.assert_(math.fabs(out1[1]-2.5)/math.sqrt(cov1[1][1]) < 2)
Note: See TracChangeset
for help on using the changeset viewer.