Changeset 64f0c5d in sasview for sansmodels/src/sans/models


Ignore:
Timestamp:
Jul 18, 2011 7:40:10 PM (13 years ago)
Author:
Jae Cho <jhjcho@…>
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:
458c930
Parents:
d9341f2
Message:

updated utests according to the correction and changes in models

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sansmodels/src/sans/models/test/utest_other_models.py

    r7ef319e r64f0c5d  
    2222            - hardSphereStructure 
    2323             
     24    @ Note: We don't use matrix for 2D anymore so testEval2D can be ignored (JC)      
    2425             
    2526""" 
     
    4849        self.x_array = self.comp.evalDistribution(self.x) 
    4950        self.y_array = self.comp.evalDistribution(self.y) 
    50         qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    51         qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    52         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     51         
     52        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
     53        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
     54        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    5355         
    5456    def test1D(self): 
     
    6769    def testEval_2D(self): 
    6870        """ Test 2D model for a coreshell  with evalDistribution""" 
    69         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    70         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    71         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    72         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    73          
    74     def testCriticalPoint(self): 
    75         """ Test coreshell at the critical point""" 
    76         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     71        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     72        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     73 
     74    # No more singular point     
     75    #def testCriticalPoint(self): 
     76    #    """ Test coreshell at the critical point""" 
     77    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    7778   
    7879class TestMultiShellModel(unittest.TestCase): 
     
    9899        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    99100        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    100         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     101        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    101102         
    102103    def test1D(self): 
     
    115116    def testEval_2D(self): 
    116117        """ Test 2D model for a MultiShell  with evalDistribution""" 
    117         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    118         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    119         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    120         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    121      
    122     def testCriticalPoint(self): 
    123         """ Test multishell at the critical point""" 
    124         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     118        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     119        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     120 
     121    # No more singular point 
     122    #def testCriticalPoint(self): 
     123    #    """ Test multishell at the critical point""" 
     124    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    125125     
    126126class TestVesicleModel(unittest.TestCase): 
     
    142142        self.x_array = self.comp.evalDistribution(self.x) 
    143143        self.y_array = self.comp.evalDistribution(self.y) 
    144         qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    145         qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    146         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     144        #qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
     145        #qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
     146        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    147147         
    148148    def test1D(self): 
     
    161161    def testEval_2D(self): 
    162162        """ Test 2D model for a Vesicle with evalDistribution""" 
    163         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    164         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    165         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    166         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     163        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     164        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     165 
    167166     
    168167    def testCriticalPoint(self): 
     
    191190        self.x_array = self.comp.evalDistribution(self.x) 
    192191        self.y_array = self.comp.evalDistribution(self.y) 
    193         qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    194         qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    195         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     192        #qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
     193        #qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
     194        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    196195         
    197196    def test1D(self): 
     
    210209    def testEval_2D(self): 
    211210        """ Test 2D model for a BinaryHS with evalDistribution""" 
    212         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    213         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    214         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    215         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    216          
    217     def testCriticalPoint(self): 
    218         """ Test BinaryHS at the critical point""" 
    219         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     211        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     212        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     213 
     214    # No more singular point     
     215    #def testCriticalPoint(self): 
     216    #    """ Test BinaryHS at the critical point""" 
     217    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    220218        
    221219      
     
    243241        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    244242        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    245         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     243        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    246244         
    247245        self.q = 0.001 
     
    266264    def testEval_2D(self): 
    267265        """ Test 2D model for a CoreShellCylinder with evalDistribution""" 
    268         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    269         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    270         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    271         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    272   
    273     def testCriticalPoint(self): 
    274         """ Test CoreShellCylinder at the critical point""" 
    275         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     266        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     267        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     268 
     269    # No more singular point 
     270    #def testCriticalPoint(self): 
     271    #    """ Test CoreShellCylinder at the critical point""" 
     272    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    276273  
    277274  
     
    287284        self.comp.setParam("radius",30) 
    288285        self.comp.setParam("length", 400) 
    289         self.comp.setParam("contrast",5.3e-6 ) 
     286        self.comp.setParam("sldCyl",6.3e-6 ) 
     287        self.comp.setParam("sldSolv",1e-6 ) 
    290288        self.comp.setParam("Background", 0.01) 
    291289         
     
    296294        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    297295        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    298         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     296        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    299297         
    300298        self.q = 0.001 
     
    319317    def testEval_2D(self): 
    320318        """ Test 2D model for a HollowCylinder with evalDistribution""" 
    321         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    322         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    323         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    324         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    325          
    326     def testCriticalPoint(self): 
    327         """ Test HollowCylinder at the critical point""" 
    328         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     319        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     320        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     321 
     322    # No more singular point     
     323    #def testCriticalPoint(self): 
     324    #    """ Test HollowCylinder at the critical point""" 
     325    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    329326         
    330327         
     
    337334        #Give initial value to model 
    338335        self.comp.setParam("scale", 1.0) 
    339         self.comp.setParam("contrast",5.3e-6 ) 
     336        self.comp.setParam("sldSolv",6.3e-6 ) 
     337        self.comp.setParam("sldCyl",1e-6 ) 
    340338        self.comp.setParam("kuhn_length",100) 
    341339        self.comp.setParam("length", 1000) 
     
    350348        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    351349        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    352         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     350        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    353351         
    354352        self.q = 0.001 
     
    372370    def testEval_2D(self): 
    373371        """ Test 2D model for a FlexibleCylinder Model with evalDistribution""" 
    374         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    375         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    376         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    377         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    378     
    379     def testCriticalPoint(self): 
    380         """ Test FlexibleCylinder at the critical point""" 
    381         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     372        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     373        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     374 
     375    # No more singular point 
     376    #def testCriticalPoint(self): 
     377    #    """ Test FlexibleCylinder at the critical point""" 
     378    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    382379   
    383380               
     
    409406        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    410407        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    411         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     408        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    412409         
    413410        self.q = 0.001 
     
    432429    def testEval_2D(self): 
    433430        """ Test 2D model for a StackedDisks Model with evalDistribution""" 
    434         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    435         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    436         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    437         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    438      
    439     def testCriticalPoint(self): 
    440         """ Test StackedDisks at the critical point""" 
    441         self.assert_(numpy.isfinite(self.comp.run(0.0)))  
     431        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     432        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     433 
     434    # No more singular point 
     435    #def testCriticalPoint(self): 
     436    #    """ Test StackedDisks at the critical point""" 
     437    #    self.assert_(numpy.isfinite(self.comp.run(0.0)))  
    442438    
    443439   
     
    453449        self.comp.setParam("short_b", 75) 
    454450        self.comp.setParam("long_c",400 ) 
    455         self.comp.setParam("contrast", 5.3e-006 ) 
     451        self.comp.setParam("sldPipe", 6.3e-006 ) 
     452        self.comp.setParam("sldSolv", 1e-006 ) 
    456453        self.comp.setParam("scale",1.0 ) 
    457454         
     
    463460        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    464461        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    465         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     462        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    466463         
    467464        self.q = 0.001 
     
    487484    def testEval_2D(self): 
    488485        """ Test 2D model for a Parallelepiped Model with evalDistribution""" 
    489         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    490         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    491         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    492         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     486        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     487        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     488 
    493489    
    494490    def testCriticalPoint(self): 
     
    506502        self.comp.setParam("r_ratio",1.5)  
    507503        self.comp.setParam("length",400.0) 
    508         self.comp.setParam("contrast",5.3e-006) 
     504        self.comp.setParam("sldCyl",4e-006) 
     505        self.comp.setParam("sldSolv",1e-006) 
    509506        self.comp.setParam("background",0.0) 
    510507        self.comp.setParam("cyl_theta",0.0) 
     
    519516        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    520517        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    521         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     518        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    522519         
    523520        self.q = 0.001 
     
    528525    def test1D(self): 
    529526        """ Test 1D model for a EllipticalCylinder Model""" 
    530         self.assertAlmostEqual(self.comp.run(0.001),2108.32, 1) 
     527        self.assertAlmostEqual(self.comp.run(0.001), 675.504402, 4) 
    531528         
    532529    def test1D_2(self): 
     
    542539    def testEval_2D(self): 
    543540        """ Test 2D model for a EllipticalCylinder with evalDistribution""" 
    544         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    545         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    546         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    547         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     541        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     542        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     543 
    548544     
    549545    def testCriticalPoint(self): 
     
    561557        self.comp.setParam("radius_a",20.0) 
    562558        self.comp.setParam("radius_b",400.0) 
    563         self.comp.setParam("contrast",3e-006) 
     559        self.comp.setParam("sldEll",4e-006) 
     560        self.comp.setParam("sldSolv",1e-006) 
    564561        self.comp.setParam("background",0.0) 
    565562        self.comp.setParam("axis_theta",1.57) 
     
    573570        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    574571        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    575         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     572        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    576573         
    577574        self.q = 0.001 
     
    595592    def testEval_2D(self): 
    596593        """ Test 2D model for a Ellipsoid Model with evalDistribution""" 
    597         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    598         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    599         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    600         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     594        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     595        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     596 
    601597     
    602598    def testCriticalPoint(self): 
     
    616612        self.comp.setParam("equat_shell",250.0) 
    617613        self.comp.setParam("polar_shell", 30.0) 
    618         self.comp.setParam("contrast",1e-006) 
     614        self.comp.setParam("sld_shell",1e-006) 
     615        self.comp.setParam("sld_core",2e-006) 
    619616        self.comp.setParam("sld_solvent",6.3e-006) 
    620617        self.comp.setParam("background",0.001) 
     
    629626        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    630627        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    631         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     628        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    632629         
    633630        self.q = 0.001 
     
    638635    def test1D(self): 
    639636        """ Test 1D model for a CoreShellEllipsoid Model""" 
    640         self.assertAlmostEqual(self.comp.run(1.0), 0.002578943, 4) 
     637        self.assertAlmostEqual(self.comp.run(1.0), 0.001894, 4) 
    641638         
    642639    def test1D_2(self): 
     
    652649    def testEval_2D(self): 
    653650        """ Test 2D model for a CoreShellEllipsoid with evalDistribution""" 
    654         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    655         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    656         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    657         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     651        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     652        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     653 
    658654     
    659655    def testCriticalPoint(self): 
     
    671667        self.comp.setParam("semi_axisB", 100.0) 
    672668        self.comp.setParam("semi_axisC",400.0 ) 
    673         self.comp.setParam("contrast",5.3e-6) 
     669        self.comp.setParam("sldSolv",6.3e-6) 
     670        self.comp.setParam("sldEll",1e-6) 
    674671        self.comp.setParam("background",0.0) 
    675672        self.comp.setParam("axis_theta", 1.0) 
     
    684681        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    685682        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    686         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     683        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    687684         
    688685        self.q = 0.001 
     
    708705    def testEval_2D(self): 
    709706        """ Test 2D model for a TriaxialEllipsoid Model with evalDistribution""" 
    710         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    711         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    712         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    713         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     707        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     708        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     709 
    714710     
    715711    def testCriticalPoint(self): 
     
    736732        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    737733        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    738         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     734        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    739735         
    740736    def test1D(self): 
    741737        """ Test 1D model for a Lamellar Model""" 
    742         self.assertAlmostEquals(self.comp.run(0.001), 902120) 
     738        self.assertAlmostEquals(self.comp.run(0.001), 882289.54309, 3) 
    743739         
    744740    def test1D_2(self): 
    745741        """ Test 2D model for a Lamellar Model""" 
    746         self.assertAlmostEquals(self.comp.run([0.001, 1.3]),902120)      
     742        self.assertAlmostEquals(self.comp.run([0.001, 1.3]),882289.54309, 3)      
    747743     
    748744    def testEval_1D(self): 
     
    753749    def testEval_2D(self): 
    754750        """ Test 2D model for a Lamellar Model with evalDistribution""" 
    755         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    756         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    757         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    758         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8)         
    759     
    760     def testCriticalPoint(self): 
    761         """ Test Lamellar at the critical point""" 
    762         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     751        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     752        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     753        
     754    # No more singular point 
     755    #def testCriticalPoint(self): 
     756    #    """ Test Lamellar at the critical point""" 
     757    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    763758      
    764759class TestLamellarFFHGModel(unittest.TestCase): 
     
    783778        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    784779        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    785         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     780        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    786781         
    787782         
    788783    def test1D(self): 
    789784        """ Test 1D model for a LamellarFFHG Model""" 
    790         self.assertAlmostEquals(self.comp.run(0.001),653144) 
     785        self.assertAlmostEquals(self.comp.run(0.001),653143.9209, 3) 
    791786         
    792787    def test1D_2(self): 
    793788        """ Test 2D model for a LamellarFFHG Model""" 
    794         self.assertAlmostEquals(self.comp.run([0.001, 1.3]),653144) 
     789        self.assertAlmostEquals(self.comp.run([0.001, 1.3]),653143.9209, 3) 
    795790     
    796791    def testEval_1D(self): 
     
    801796    def testEval_2D(self): 
    802797        """ Test 2D model for a LamellarFFHG Model with evalDistribution""" 
    803         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    804         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    805         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    806         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8)         
    807      
    808     def testCriticalPoint(self): 
    809         """ Test LamellarFFHG at the critical point""" 
    810         self.assert_(numpy.isfinite(self.comp.run(0.0)))               
     798        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     799        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     800       
     801    # No more singular point 
     802    #def testCriticalPoint(self): 
     803    #    """ Test LamellarFFHG at the critical point""" 
     804    #    self.assert_(numpy.isfinite(self.comp.run(0.0)))               
    811805     
    812806class TestLamellarPSModel(unittest.TestCase): 
     
    819813        self.comp.setParam("spacing",400.0) 
    820814        self.comp.setParam("delta",30.0) 
    821         self.comp.setParam("contrast",5.3e-006) 
     815        self.comp.setParam("sld_bi",6.3e-006) 
     816        self.comp.setParam("sld_sol",1e-006) 
    822817        self.comp.setParam("n_plates",20.0)  
    823818        self.comp.setParam("caille", 0.1) 
     
    831826        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    832827        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    833         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     828        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    834829         
    835830         
    836831    def test1D(self): 
    837832        """ Test 1D model for a LamellarPS Model""" 
    838         self.assertAlmostEquals(self.comp.run(0.001), 28526.8, 1) 
     833        self.assertAlmostEquals(self.comp.run(0.001), 27899.30836, 1) 
    839834         
    840835    def test1D_2(self): 
    841836        """ Test 2D model for a LamellarPS Model""" 
    842         self.assertAlmostEquals(self.comp.run([0.001, 1.3]),28526.8, 1)  
     837        self.assertAlmostEquals(self.comp.run([0.001, 1.3]),27899.30836, 1)  
    843838     
    844839    def testEval_1D(self): 
     
    849844    def testEval_2D(self): 
    850845        """ Test 2D model for a LamellarPS Model with evalDistribution""" 
    851         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    852         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    853         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    854         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8)         
    855              
    856     def testCriticalPoint(self): 
    857         """ Test LamellarPS at the critical point""" 
    858         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     846        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     847        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     848       
     849    # No more singular point         
     850    #def testCriticalPoint(self): 
     851    #    """ Test LamellarPS at the critical point""" 
     852    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    859853      
    860854class TestLamellarPSHGModel(unittest.TestCase): 
     
    882876        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    883877        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    884         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     878        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    885879         
    886880    def test1D(self): 
    887881        """ Test 1D model for a LamellarPSHG Model""" 
    888         self.assertAlmostEquals(self.comp.run(0.001),6.83139e6, 3) 
     882        self.assertAlmostEquals(self.comp.run(0.001),6831387.29466, 3) 
    889883         
    890884    def test1D_2(self): 
    891885        """ Test 2D model for a LamellarPSHG Model""" 
    892         self.assertAlmostEquals(self.comp.run([0.001, 1.3]),6.83139e6,3) 
     886        self.assertAlmostEquals(self.comp.run([0.001, 1.3]),6831387.29466,3) 
    893887         
    894888    def testEval_1D(self): 
     
    899893    def testEval_2D(self): 
    900894        """ Test 2D model for a LamellarPSHG Model with evalDistribution""" 
    901         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    902         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    903         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    904         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8)         
    905      
    906     def testCriticalPoint(self): 
    907         """ Test LamellarPSHG at the critical point""" 
    908         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     895        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     896        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     897        
     898    # No more singular point 
     899    #def testCriticalPoint(self): 
     900    #    """ Test LamellarPSHG at the critical point""" 
     901    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    909902      
    910903class TestSquareWellStructure(unittest.TestCase): 
     
    926919        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    927920        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    928         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     921        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    929922         
    930923         
     
    944937    def testEval_2D(self): 
    945938        """ Test 2D model for a SquareWellStructure with evalDistribution""" 
    946         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    947         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    948         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    949         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    950      
    951     def testCriticalPoint(self): 
    952         """ Test SquareWellStructure at the critical point""" 
    953         self.assert_(numpy.isfinite(self.comp.run(0.0)))                      
     939        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     940        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     941 
     942    # No more singular point 
     943    #def testCriticalPoint(self): 
     944    #    """ Test SquareWellStructure at the critical point""" 
     945    #    self.assert_(numpy.isfinite(self.comp.run(0.0)))                      
    954946      
    955947class TestHardsphereStructure(unittest.TestCase): 
     
    968960        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    969961        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    970         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     962        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    971963         
    972964         
     
    986978    def testEval_2D(self): 
    987979        """ Test 2D model for a HardsphereStructure with evalDistribution""" 
    988         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    989         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    990         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    991         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    992      
    993     def testCriticalPoint(self): 
    994         """ Test HardsphereStructure at the critical point""" 
    995         self.assert_(numpy.isfinite(self.comp.run(0.0)))                   
     980        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     981        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     982 
     983    # No more singular point 
     984    #def testCriticalPoint(self): 
     985    #    """ Test HardsphereStructure at the critical point""" 
     986    #    self.assert_(numpy.isfinite(self.comp.run(0.0)))                   
    996987     
    997988class TestStickyHSStructure(unittest.TestCase): 
     
    10131004        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    10141005        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1015         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1006        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    10161007         
    10171008    def test1D(self): 
     
    10301021    def testEval_2D(self): 
    10311022        """ Test 2D model for a StickyHSStructure with evalDistribution""" 
    1032         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1033         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1034         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1035         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    1036                 
    1037     def testCriticalPoint(self): 
    1038         """ Test StickyHSStructure at the critical point""" 
    1039         self.assert_(numpy.isfinite(self.comp.run(0.0)))  
     1023        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1024        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1025 
     1026    # No more singular point            
     1027    #def testCriticalPoint(self): 
     1028    #    """ Test StickyHSStructure at the critical point""" 
     1029    #    self.assert_(numpy.isfinite(self.comp.run(0.0)))  
    10401030         
    10411031class TestHayterMSAStructure(unittest.TestCase): 
     
    10591049        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    10601050        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1061         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1051        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    10621052 
    10631053    def test1D(self): 
     
    10761066    def testEval_2D(self): 
    10771067        """ Test 2D model for a HayterMSAStructure with evalDistribution""" 
    1078         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1079         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1080         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1081         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1068        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1069        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1070 
    10821071     
    10831072    def testCriticalPoint(self): 
     
    11091098        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    11101099        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1111         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
    1112          
     1100        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
     1101 
    11131102    def test1D(self): 
    11141103        """ Test 1D model for a BEPolyelectrolyte""" 
     
    11261115    def testEval_2D(self): 
    11271116        """ Test 2D model for a BEPolyelectrolyte with evalDistribution""" 
    1128         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1129         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1130         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1131         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1117        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1118        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1119 
    11321120    
    11331121    def testCriticalPoint(self): 
     
    11531141        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    11541142        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1155         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1143        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    11561144         
    11571145    def test1D(self): 
    11581146        """ Test 1D model for a DABModel""" 
    1159         self.assertAlmostEqual(self.comp.run(0.001),10.9681, 3) 
     1147        self.assertAlmostEqual(self.comp.run(0.001),637957.9047, 3) 
    11601148         
    11611149    def test1D_2(self): 
    11621150        """ Test 2D model for a DABModel""" 
    1163         self.assertAlmostEqual(self.comp.run([0.001, 1.3]),10.9681, 3) 
     1151        self.assertAlmostEqual(self.comp.run([0.001, 1.3]),637957.90473, 3) 
    11641152          
    11651153    def testEval_1D(self): 
     
    11701158    def testEval_2D(self): 
    11711159        """ Test 2D model for a DABModel with evalDistribution""" 
    1172         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1173         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1174         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1175         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1160        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1161        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1162 
    11761163     
    11771164    def testCriticalPoint(self): 
     
    11961183        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    11971184        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1198         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1185        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    11991186         
    12001187    def test1D(self): 
     
    12131200    def testEval_2D(self): 
    12141201        """ Test 2D model for a GuinierModel with evalDistribution""" 
    1215         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1216         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1217         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1218         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1202        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1203        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1204 
    12191205     
    12201206    def testCriticalPoint(self): 
     
    12401226        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    12411227        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1242         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1228        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    12431229         
    12441230    def test1D(self): 
     
    12571243    def testEval_2D(self): 
    12581244        """ Test 2D model for a DebyeModel with evalDistribution""" 
    1259         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1260         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1261         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1262         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1245        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1246        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1247 
    12631248     
    12641249    def testCriticalPoint(self): 
     
    12831268        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    12841269        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1285         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1270        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    12861271         
    12871272    def test1D(self): 
     
    13001285    def testEval_2D(self): 
    13011286        """ Test 2D model for a PorodModel with evalDistribution""" 
    1302         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1303         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1304         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1305         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1287        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1288        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1289 
    13061290     
    13071291    def testCreaticalPoint(self): 
     
    13281312        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    13291313        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1330         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1314        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    13311315         
    13321316    def test1D(self): 
     
    13451329    def testEval_2D(self): 
    13461330        """ Test 2D model for a PeakGauss Model with evalDistribution""" 
    1347         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1348         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1349         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1350         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1331        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1332        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
    13511333     
    13521334    def testCriticalPoint(self): 
     
    13721354        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    13731355        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1374         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1356        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    13751357         
    13761358    def test1D(self): 
     
    13891371    def testEval_2D(self): 
    13901372        """ Test 2D model for a PeakLorentz Model with evalDistribution""" 
    1391         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1392         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1393         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1394         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1373        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1374        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1375 
    13951376         
    13961377    def testCriticalPoint(self): 
     
    13981379        self.comp.setParam('B', 0.0) 
    13991380        self.assertRaises(ZeroDivisionError, self.comp.run, 10) 
    1400         self.assert_(numpy.isfinite(self.comp.run(0.0)))   
     1381        #self.assert_(numpy.isfinite(self.comp.run(0.0)))   
    14011382         
    14021383        
     
    14101391        self.comp.setParam('radius', 5.0) 
    14111392        self.comp.setParam('fractal_dim', 2.0) 
    1412         self.comp.setParam('corr_length',100.0) 
    1413         self.comp.setParam('block_sld', 2.0e-6) 
    1414         self.comp.setParam('solvent_sld', 6.35e-6) 
     1393        self.comp.setParam('cor_length',100.0) 
     1394        self.comp.setParam('sldBlock', 2.0e-6) 
     1395        self.comp.setParam('sldSolv', 6.35e-6) 
    14151396        self.comp.setParam('background',0.0) 
    14161397         
     
    14221403        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    14231404        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1424         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1405        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    14251406         
    14261407    def test1D(self): 
     
    14391420    def testEval_2D(self): 
    14401421        """ Test 2D model for a Fractal Abs Model with evalDistribution""" 
    1441         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1442         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1443         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1444         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    1445      
    1446     def testCriticalPoint(self): 
    1447         """ Test Fractal Abs at the critical point""" 
    1448         self.assert_(numpy.isfinite(self.comp.run(0.0)))     
     1422        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1423        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1424 
     1425    # No more singular point 
     1426    #def testCriticalPoint(self): 
     1427    #    """ Test Fractal Abs at the critical point""" 
     1428    #    self.assert_(numpy.isfinite(self.comp.run(0.0)))     
    14491429    
    14501430class TestFractalModel(unittest.TestCase): 
     
    14571437        self.comp.setParam('radius', 5.0) 
    14581438        self.comp.setParam('fractal_dim', 2.0) 
    1459         self.comp.setParam('corr_length',100.0) 
    1460         self.comp.setParam('block_sld', 2.0e-6) 
    1461         self.comp.setParam('solvent_sld', 6.35e-6) 
     1439        self.comp.setParam('cor_length',100.0) 
     1440        self.comp.setParam('sldBlock', 2.0e-6) 
     1441        self.comp.setParam('sldSolv', 6.35e-6) 
    14621442        self.comp.setParam('background',0.0) 
    14631443         
     
    14691449        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    14701450        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1471         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1451        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    14721452         
    14731453    def test1D(self): 
     
    14861466    def testEval_2D(self): 
    14871467        """ Test 2D model for a Fractal Model with evalDistribution""" 
    1488         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1489         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1490         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1491         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
    1492          
    1493     def testCriticalPoint(self): 
    1494         """ Test Fractal at the critical point""" 
    1495         self.assert_(numpy.isfinite(self.comp.run(0.0))) 
     1468        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1469        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1470 
     1471    # No more singular point     
     1472    #def testCriticalPoint(self): 
     1473    #    """ Test Fractal at the critical point""" 
     1474    #    self.assert_(numpy.isfinite(self.comp.run(0.0))) 
    14961475     
    14971476class TestLorentzModel(unittest.TestCase): 
     
    15121491        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    15131492        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1514         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1493        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    15151494         
    15161495    def test1D(self): 
     
    15291508    def testEval_2D(self): 
    15301509        """ Test 2D model for a Lorentz Model with evalDistribution""" 
    1531         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1532         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1533         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1534         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1510        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1511        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1512 
    15351513           
    15361514    def testCriticalPoint(self): 
     
    15551533        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    15561534        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1557         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1535        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    15581536         
    15591537    def test1D(self): 
     
    15721550    def testEval_2D(self): 
    15731551        """ Test 2D model for a PowerLawAbs Model with evalDistribution""" 
    1574         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1575         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1576         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1577         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1552        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1553        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1554 
    15781555     
    15791556    def testCriticalPoint(self): 
     
    15981575        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    15991576        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1600         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1577        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    16011578         
    16021579    def test1D(self): 
     
    16151592    def testEval_2D(self): 
    16161593        """ Test 2D model for a PowerLawModel with evalDistribution""" 
    1617         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1618         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1619         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1620         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1594        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1595        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1596 
    16211597     
    16221598    def testCriticalPoint(self): 
     
    16421618        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    16431619        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1644         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1620        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    16451621         
    16461622    def test1D(self): 
     
    16591635    def testEval_2D(self): 
    16601636        """ Test 2D model for a TeubnerStrey  with evalDistribution""" 
    1661         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1662         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1663         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1664         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1637        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0],8) 
     1638        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1639 
    16651640     
    16661641    def testCriticalPoint(self): 
     
    16841659        qx_prime = numpy.reshape(self.x, [1,len(self.x)]) 
    16851660        qy_prime = numpy.reshape(self.y, [len(self.y),1]) 
    1686         self.xy_matrix = self.comp.evalDistribution([qx_prime, qy_prime]) 
     1661        self.xy_matrix = self.comp.evalDistribution([self.x, self.y]) 
    16871662         
    16881663    def test1D(self): 
     
    16971672    def testEval_2D(self): 
    16981673        """ Test 2D model for a Line  with evalDistribution""" 
    1699         self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0][0],8) 
    1700         self.assertAlmostEquals(self.comp.runXY([0.4,1.57]),self.xy_matrix[1][0], 8) 
    1701         self.assertAlmostEquals(self.comp.runXY([1.3, 0.5]),self.xy_matrix[0][1],8) 
    1702         self.assertAlmostEquals(self.comp.runXY([1.3, 1.57]),self.xy_matrix[1][1],8) 
     1674        self.assertAlmostEquals(self.comp.runXY([0.4, 0.5]),self.xy_matrix[0], 8) 
     1675        self.assertAlmostEquals(self.comp.runXY([1.3,1.57]),self.xy_matrix[1], 8) 
     1676 
    17031677     
    17041678    def testCriticalPoint(self): 
Note: See TracChangeset for help on using the changeset viewer.