mailr24318 - in /branches/frame_order_cleanup/test_suite/system_tests: ./ scripts/frame_order/cam/


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by edward on June 25, 2014 - 16:53:
Author: bugman
Date: Wed Jun 25 16:53:16 2014
New Revision: 24318

URL: http://svn.gna.org/viewcvs/relax?rev=24318&view=rev
Log:
Created 3 frame order system tests for the new large angle double rotor CaM 
synthetic data.

These are the Frame_order.test_cam_double_rotor_large_angle,
Frame_order.test_cam_double_rotor_large_angle_rdc, and
Frame_order.test_cam_double_rotor_large_angle_pcs system tests.


Added:
    
branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor_large_angle.py
      - copied, changed from r24317, 
branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor.py
Modified:
    branches/frame_order_cleanup/test_suite/system_tests/frame_order.py

Modified: branches/frame_order_cleanup/test_suite/system_tests/frame_order.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_cleanup/test_suite/system_tests/frame_order.py?rev=24318&r1=24317&r2=24318&view=diff
==============================================================================
--- branches/frame_order_cleanup/test_suite/system_tests/frame_order.py 
(original)
+++ branches/frame_order_cleanup/test_suite/system_tests/frame_order.py Wed 
Jun 25 16:53:16 2014
@@ -267,6 +267,33 @@
         self.check_chi2(0.080669539157379247)
 
 
+    def test_cam_double_rotor_large_angle(self):
+        """Test the double rotor frame order model of CaM."""
+
+        # The flags, execute the script, and then check the chi2 value.
+        self.flags()
+        
self.interpreter.run(script_file=self.cam_path+'double_rotor_large_angle.py')
+        self.check_chi2(0.081769384900455433)
+
+
+    def test_cam_double_rotor_large_angle_pcs(self):
+        """Test the double rotor frame order model of CaM (with only PCS 
data)."""
+
+        # The flags, execute the script, and then check the chi2 value.
+        self.flags(rdc=False)
+        
self.interpreter.run(script_file=self.cam_path+'double_rotor_large_angle.py')
+        self.check_chi2(0.00016331682903726827)
+
+
+    def test_cam_double_rotor_large_angle_rdc(self):
+        """Test the double rotor frame order model of CaM (with only RDC 
data)."""
+
+        # The flags, execute the script, and then check the chi2 value.
+        self.flags(pcs=False)
+        
self.interpreter.run(script_file=self.cam_path+'double_rotor_large_angle.py')
+        self.check_chi2(0.080669539157379247)
+
+
     def test_cam_free_rotor(self):
         """Test the free rotor frame order model of CaM."""
 

Copied: 
branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor_large_angle.py
 (from r24317, 
branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor.py)
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor_large_angle.py?p2=branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor_large_angle.py&p1=branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor.py&r1=24317&r2=24318&rev=24318&view=diff
==============================================================================
--- 
branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor.py
        (original)
+++ 
branches/frame_order_cleanup/test_suite/system_tests/scripts/frame_order/cam/double_rotor_large_angle.py
    Wed Jun 25 16:53:16 2014
@@ -39,10 +39,10 @@
     # The model parameters.
     AXIS_THETA = 1.494291741547518
     AXIS_PHI = 2.525044022476957
-    CONE_SIGMA_MAX = 10.5 / 2.0 / 360.0 * 2.0 * pi
+    CONE_SIGMA_MAX = 55.0 / 2.0 / 360.0 * 2.0 * pi
     AXIS_THETA2 = 2.30381499622381
     AXIS_PHI2 = -2.249696457768556
-    CONE_SIGMA_MAX2 = 11.5 / 2.0 / 360.0 * 2.0 * pi
+    CONE_SIGMA_MAX2 = 85.0 / 2.0 / 360.0 * 2.0 * pi
 
     # The pivot points.
     PIVOT = array([41.739, 6.03, -0.764], float32)




Related Messages


Powered by MHonArc, Updated Wed Jun 25 17:20:03 2014