Changeset 62827da in sasview for sansmodels

Timestamp:

Dec 18, 2012 9:12:03 AM (12 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:

5b07138

Parents:

6550b64

Message:

updating tests

Location:

sansmodels/test

Files:

• utest_dispersity.py (modified) (3 diffs)
• utest_models.py (modified) (3 diffs)
• utest_models_array.py (modified) (1 diff)
• utest_other_dispersity.py (modified) (1 diff)
• validate_2D_3axes_shape_model.py (modified) (1 diff)
• validate_2D_elliptical_cyl.py (modified) (2 diffs)
• validate_2D_model.py (modified) (2 diffs)

sansmodels/test/utest_dispersity.py

@@ -29 +29 @@
         self.model.setParam('background', 0.0)
         self.model.setParam('cyl_theta', 0.0)
-        self.model.setParam('cyl_phi', 0.0)
+        self.model.setParam('cyl_phi', 90.0)
         
     def test_simple(self):
@@ -134 +134 @@
         self.model.setParam('background', 0.0)
         self.model.setParam('axis_theta', 0.0)
-        self.model.setParam('axis_phi', 0.0)
+        self.model.setParam('axis_phi', 90.0)
         
     def test_simple(self):
@@ -282 +282 @@
         self.model.setParam('sldSolv', 1.e-6)
         self.model.setParam('background', 0.0)
-        self.model.setParam('axis_theta', 89.95445)
+        self.model.setParam('axis_theta', 0.0)
         self.model.setParam('axis_phi', 0.0)
         

sansmodels/test/utest_models.py

@@ -60 +60 @@
     def test2D(self):
         """ Test 2D model of a cylinder """ 
-        self.comp.setParam('cyl_theta', 1*180/math.pi)
-        self.comp.setParam('cyl_phi', 1*180/math.pi)
+        self.comp.setParam('cyl_theta', 10.0)
+        self.comp.setParam('cyl_phi', 10.0)
         self.assertAlmostEqual(self.comp.run([0.2, 2.5]), 
-                               0.038176446608393366, 4) 
+                               0.038176446608393366, 2) 
  
 class TestGaussian(unittest.TestCase):
@@ -96 +96 @@
         self.assertAlmostEqual(self.gauss.run([math.sqrt(8.0), math.pi/4.0]), 2.0*2.0*value*value, 5)
         
-class TestStarPolymer(unittest.TestCase):
-    """Unit tests for Star Polymer"""
-    def setUp(self):
-        from sans.models.StarPolymer import StarPolymer
-        self.star = StarPolymer()
-
-    def test1D(self):
-        self.star.setParam('scale', 3.0)
-        self.star.setParam('arms', 5.0)
-        self.star.setParam('R2', 100.0)
-        self.star.setParam('background', 1.0)
-        self.assertAlmostEqual(self.star.run(0.01), 3.9900214150405686)
-
-    def test2D(self):
-        x = 0.009539392014169456
-        y = 0.003
-        # together these have modulus 1
-        self.star.setParam('scale', 3.0)
-        self.star.setParam('arms', 5.0)
-        self.star.setParam('R2', 100.0)
-        self.star.setParam('background', 1.0)
-        self.assertAlmostEqual(self.star.runXY([x, y]), 3.9900214150405686)
-
-    def test2Dphi(self):
-        self.star.setParam('scale', 3.0)
-        self.star.setParam('arms', 5.0)
-        self.star.setParam('R2', 100.0)
-        self.star.setParam('background', 1.0)
-        
-        self.assertAlmostEqual(self.star.run([0.01, math.pi/4.0]), 3.9900214150405686)
-
-
-
 class TestLorentzian(unittest.TestCase):
     """Unit tests for Lorentzian function"""
@@ -134 +101 @@
     def setUp(self):
         from sans.models.Lorentzian import Lorentzian
-        self.lor = Lorentzian()
+        self.lor= Lorentzian()
         
     def test1D(self):

sansmodels/test/utest_models_array.py

@@ -69 +69 @@
     def test1D_2(self):
         """ Test 2D model of a cylinder """ 
-        self.comp.setParam('cyl_theta', 1.0 * 180/3.1415)
-        self.comp.setParam('cyl_phi', 1.0 * 180/3.1415)
+        self.comp.setParam('cyl_theta', 10.0)
+        self.comp.setParam('cyl_phi', 10.0)
         self.assertAlmostEqual(self.comp.run([0.2, 2.5]), 
-                               0.038176446608393366, 3)
+                               0.038176446608393366, 2)
         
     def test1D_3(self):

sansmodels/test/utest_other_dispersity.py

@@ -21 +21 @@
         self.model.setParam('background', 0.0)
         self.model.setParam('cyl_theta', 0.0)
-        self.model.setParam('cyl_phi', 0.0)
+        self.model.setParam('cyl_phi', 90.0)
         
     def test_simple(self):

sansmodels/test/validate_2D_3axes_shape_model.py

@@ -103 +103 @@
                             
                         # sin() is due to having not uniform bin number density wrt the q plane.
-                        sum += model.run([q, 0])*math.sin(theta*math.pi/180.0)
-                        weight += math.sin(theta*math.pi/180.0) 
+                        sum += model.run([q, 0])*math.fabs(math.cos(theta*math.pi/180.0))
+                        weight += math.fabs(math.cos(theta*math.pi/180.0) )
 
             value = sum/weight #*math.pi/2.0

sansmodels/test/validate_2D_elliptical_cyl.py

@@ -87 +87 @@
                             pass
                     else:
-                        sum += math.sin(theta)*model.run([q, 0])
+                        sum += math.fabs(math.cos(theta))*model.run([q, 0])
          
         value = sum/npts/npts/npts
@@ -205 +205 @@
                         print "ERROR", q, theta, 180 * 2.0 / npts * j
                 else:
-                    sum += math.sin(theta)*model.run([q, 0])
+                    sum += math.fabs(math.cos(theta))*model.run([q, 0])
         
         value = sum/npts/npts*math.pi/2.0

sansmodels/test/validate_2D_model.py

@@ -88 +88 @@
                     if str(model.run([q, 0])).count("INF")>0:                        
                         print "ERROR", q, theta, 180.0 * 2.0 / npts * j
-                    sum += math.sin(theta*math.pi/180.0)*model.run([q, 0])
+                    sum += math.fabs(math.cos(theta*math.pi/180.0))*model.run([q, 0])
                     #sum += model.run([q, 0])
             
@@ -112 +112 @@
         model = model_class()
         model.setParam('scale', 1.0)
-        model.setParam('contrast', 1.0)
+        #model.setParam('contrast', 1.0)
         
         theta_label = 'cyl_theta'