Changeset f980d4a in sasview


Ignore:
Timestamp:
Nov 13, 2007 9:39:05 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:
883a2ef
Parents:
5a2070e
Message:

Updated testing of 2D sim (made sure orientation was right)

Location:
sansrealspace/src/realspace/test
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • sansrealspace/src/realspace/test/cylinder_d=0.5_Iq2D.txt

    ra2c1196 rf980d4a  
    1 PARS: {'scale': 1.0, 'length': 100, 'radius': 20, 'background': 0.0, 'cyl_theta': 1.0, 'contrast': 1.0, 'cyl_phi': 1.0} 
     1PARS: {'scale': 1.0, 'length': 100, 'radius': 20, 'background': 0.0, 'cyl_theta': 1.5707963267948966, 'contrast': 1.0, 'cyl_phi': 1.5707963267948966} 
    22<q>  <ana>  <sim>  <err> 
    3 0.001  1.25586e+013  1.25584e+013  5.26233e+008 
    4 0.00929167  1.1909e+013  1.19068e+013  1.04923e+010 
    5 0.0175833  1.0341e+013  1.03408e+013  3.11388e+010 
    6 0.025875  8.16205e+012  8.15049e+012  5.18067e+010 
    7 0.0341667  5.77306e+012  5.75279e+012  6.23249e+010 
    8 0.0424583  3.56425e+012  3.55305e+012  5.89665e+010 
    9 0.05075  1.82415e+012  1.84081e+012  4.53723e+010 
    10 0.0590417  6.86403e+011  6.74985e+011  2.6955e+010 
    11 0.0673333  1.26379e+011  1.248e+011  1.07889e+010 
    12 0.075625  1.69416e+009  2.89547e+009  1.50778e+009 
    13 0.0839167  1.17674e+011  1.07376e+011  8.66466e+009 
    14 0.0922083  2.93262e+011  2.90816e+011  1.4123e+010 
    15 0.1005  4.06932e+011  3.98585e+011  1.69585e+010 
    16 0.108792  4.1248e+011  4.30879e+011  1.82588e+010 
    17 0.117083  3.26735e+011  3.32331e+011  1.66224e+010 
    18 0.125375  2.00372e+011  1.9155e+011  1.2709e+010 
    19 0.133667  8.63351e+010  9.06573e+010  8.67162e+009 
    20 0.141958  1.76833e+010  1.4993e+010  3.5005e+009 
    21 0.15025  2.29778e+008  2.08421e+008  4.0582e+008 
    22 0.158542  1.83678e+010  1.58227e+010  3.54108e+009 
    23 0.166833  4.78555e+010  4.61547e+010  6.02634e+009 
    24 0.175125  6.83196e+010  6.41105e+010  7.15624e+009 
    25 0.183417  7.03999e+010  7.16712e+010  7.5763e+009 
    26 0.191708  5.61199e+010  5.95249e+010  6.90467e+009 
    27 0.2  3.43293e+010  4.68345e+010  6.07661e+009 
     30.001  1.25586e+013  1.25584e+013  5.25262e+008 
     40.00929167  1.1909e+013  1.19058e+013  1.05064e+010 
     50.0175833  1.0341e+013  1.03345e+013  3.12179e+010 
     60.025875  8.16205e+012  8.14926e+012  5.17962e+010 
     70.0341667  5.77306e+012  5.81045e+012  6.21145e+010 
     80.0424583  3.56425e+012  3.55223e+012  5.90114e+010 
     90.05075  1.82415e+012  1.80424e+012  4.48859e+010 
     100.0590417  6.86403e+011  6.93632e+011  2.71769e+010 
     110.0673333  1.26379e+011  1.2136e+011  1.06077e+010 
     120.075625  1.69416e+009  1.50569e+009  1.08773e+009 
     130.0839167  1.17674e+011  1.24732e+011  9.32509e+009 
     140.0922083  2.93262e+011  2.7841e+011  1.38067e+010 
     150.1005  4.06932e+011  4.09882e+011  1.72092e+010 
     160.108792  4.1248e+011  4.37694e+011  1.84369e+010 
     170.117083  3.26735e+011  3.07555e+011  1.59575e+010 
     180.125375  2.00372e+011  1.98719e+011  1.30112e+010 
     190.133667  8.63351e+010  8.31996e+010  8.33716e+009 
     200.141958  1.76833e+010  1.4757e+010  3.47028e+009 
     210.15025  2.29778e+008  5.88028e+008  6.83052e+008 
     220.158542  1.83678e+010  1.6151e+010  3.56511e+009 
     230.166833  4.78555e+010  5.04122e+010  6.33775e+009 
     240.175125  6.83196e+010  7.36705e+010  7.69267e+009 
     250.183417  7.03999e+010  7.70008e+010  7.86228e+009 
     260.191708  5.61199e+010  5.68128e+010  6.74312e+009 
     270.2  3.43293e+010  3.12432e+010  4.97778e+009 
  • sansrealspace/src/realspace/test/early_test.py

    ra2c1196 rf980d4a  
    209209     
    210210    # Along Y 
    211     #ana.setParam('cyl_theta', 1.57) 
    212     #ana.setParam('cyl_phi', 1.57) 
     211    ana.setParam('cyl_theta', 1.57) 
     212    ana.setParam('cyl_phi', 1.57) 
    213213     
    214214    # Along Z 
    215     ana.setParam('cyl_theta', 0) 
    216     ana.setParam('cyl_phi', 0) 
     215    #ana.setParam('cyl_theta', 0) 
     216    #ana.setParam('cyl_phi', 0) 
    217217     
    218218    model = VolumeCanvas.VolumeCanvas()     
     
    226226     
    227227    # Along Y 
    228     #model.setParam('%s.orientation' % handle, [0,0,0]) 
     228    model.setParam('%s.orientation' % handle, [0,0,0]) 
    229229     
    230230    # Along Z 
    231     model.setParam('%s.orientation' % handle, [1.57,0,0]) 
    232      
    233      
    234     ana = ana.runXY([0.1, 0.01]) 
    235     sim = model.getIq2D(0.1, 0.01) 
     231    #model.setParam('%s.orientation' % handle, [1.57,0,0]) 
     232     
     233     
    236234    print model.npts 
    237     print ana, sim, sim/ana, ana/sim 
     235    for i in range(40): 
     236        qmax = 0.5 
     237        anaX = ana.runXY([qmax*i/40.0, 0.0]) 
     238        simX = model.getIq2D(qmax*i/40.0, 0.0) 
     239         
     240        anaY = ana.runXY([0, qmax*i/40.0]) 
     241        simY = model.getIq2D(0, qmax*i/40.0) 
     242        print anaX, simX, simX/anaX, '|', anaY, simY, simY/anaY 
    238243     
    239244def test_7(): 
     
    284289        print ana_val, sim_val, sim_val/ana_val, err, (sim_val-ana_val)/err 
    285290     
    286  
    287291     
    288292     
    289293if __name__ == "__main__": 
    290     test_7() 
     294    test_6() 
  • sansrealspace/src/realspace/test/sim_validation.py

    ra2c1196 rf980d4a  
    109109        self.ana.setParam('radius', radius) 
    110110        self.ana.setParam('length', length) 
     111        self.ana.setParam('cyl_theta', math.pi/2.0) 
     112        self.ana.setParam('cyl_phi', math.pi/2.0) 
    111113        self.create() 
    112114        
     
    141143        self.ana.setParam('radius_b', radius_b) 
    142144        #self.ana.setParam('radius', radius_a) 
     145         
     146        # Default orientation is there=1.57, phi=0 
     147        # Radius_a is along the x direction 
    143148 
    144149        self.create() 
     
    311316     
    312317    # 2D: Density=5, 71.2 secs for 50 points 
    313     vali = CoreShellValidator(radius = 15, thickness=5, density = 5.0) 
     318    #vali = CoreShellValidator(radius = 15, thickness=5, density = 5.0) 
    314319    #validate_model(vali, q_min=0.001, q_max=1, n_q=50) 
    315     validate_model_2D(vali, q_min=0.001, q_max=1, phi=1.0, n_q=50) 
     320    #validate_model_2D(vali, q_min=0.001, q_max=1, phi=1.0, n_q=50) 
    316321 
    317322    # 2D: Density=2, 11.1 secs for 25 points 
     
    325330    #vali = CylinderValidator(radius = 20, length=100, density = 0.1) 
    326331    #validate_model(vali, q_min=0.001, q_max=0.5, n_q=25) 
    327     #vali = CylinderValidator(radius = 20, length=100, density = 0.5) 
    328     #validate_model_2D(vali, q_min=0.001, q_max=0.2, phi=1.0, n_q=25) 
     332    vali = CylinderValidator(radius = 20, length=100, density = 0.5) 
     333    validate_model_2D(vali, q_min=0.001, q_max=0.2, phi=1.0, n_q=25) 
    329334     
    330335    # 2D: Density=0.5, 2.26 secs for 25 points 
  • sansrealspace/src/realspace/test/utest_oriented.py

    ra2c1196 rf980d4a  
    187187        self.ana.setParam('radius_a', radius_a) 
    188188        self.ana.setParam('radius_b', radius_b) 
    189  
     189        # Default orientation is there=1.57, phi=0 
     190        # Radius_a is along the x direction 
    190191        
    191192        canvas = VolumeCanvas.VolumeCanvas() 
Note: See TracChangeset for help on using the changeset viewer.