Changeset 02cb31a in sasview for sansrealspace/src/realspace/test/early_test.py

Timestamp:

Sep 27, 2007 11:38:08 AM (17 years ago)

Author:

Mathieu Doucet <doucetm@…>

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:

1b0707e1

Parents:

2bb0b26

Message:

Added some 2D simulation tests.

File:

• sansrealspace/src/realspace/test/early_test.py (modified) (1 diff)

sansrealspace/src/realspace/test/early_test.py

@@ -170 +170 @@
         print "q=%g  sim=%g  ana=%g   ratio=%g" % (q, sim_1, ana_1, sim_1/ana_1)
 
+def test_5():
+    from sans.models.SphereModel import SphereModel
+    model = VolumeCanvas.VolumeCanvas()
+    
+    handle = model.add('sphere')
+    
+    radius = 10
+    density = .1
+    
+    ana = SphereModel()
+    ana.setParam('scale', 1.0)
+    ana.setParam('contrast', 1.0)
+    ana.setParam('background', 0.0)
+    ana.setParam('radius', radius)
+    
+    model.setParam('lores_density', density)
+    model.setParam('%s.radius' % handle, radius)
+    model.setParam('scale' , 1.0)
+    model.setParam('%s.contrast' % handle, 1.0)
+    model.setParam('background' , 0.0)
+    
+    ana = ana.runXY([0.1, 0.1])
+    sim = model.getIq2D(0.1, 0.1)
+    print ana, sim, sim/ana, ana/sim
+
+def test_6():
+    from sans.models.CylinderModel import CylinderModel
+    radius = 5
+    length = 40
+    density = 20
+    
+    ana = CylinderModel()
+    ana.setParam('scale', 1.0)
+    ana.setParam('contrast', 1.0)
+    ana.setParam('background', 0.0)
+    ana.setParam('radius', radius)
+    ana.setParam('length', length)
+    
+    # Along Y
+    #ana.setParam('cyl_theta', 1.57)
+    #ana.setParam('cyl_phi', 1.57)
+    
+    # Along Z
+    ana.setParam('cyl_theta', 0)
+    ana.setParam('cyl_phi', 0)
+    
+    model = VolumeCanvas.VolumeCanvas()    
+    handle = model.add('cylinder')
+    model.setParam('lores_density', density)
+    model.setParam('%s.radius' % handle, radius)
+    model.setParam('%s.length' % handle, length)
+    model.setParam('scale' , 1.0)
+    model.setParam('%s.contrast' % handle, 1.0)
+    model.setParam('background' , 0.0)
+    
+    # Along Y
+    #model.setParam('%s.orientation' % handle, [0,0,0])
+    
+    # Along Z
+    model.setParam('%s.orientation' % handle, [1.57,0,0])
+    
+    
+    ana = ana.runXY([0.1, 0.01])
+    sim = math.pi/2*model.getIq2D(0.1, 0.01)
+    print model.npts
+    print ana, sim, sim/ana, ana/sim
+    
+
     
 if __name__ == "__main__":
-    test_1()
+    test_6()