r4896 - /branches/N_state_model/test_suite/system_tests/n_state_model.py -- January 28, 2008 - 17:11

Author: bugman
Date: Mon Jan 28 17:11:14 2008
New Revision: 4896

URL: http://svn.gna.org/viewcvs/relax?rev=4896&view=rev
Log:
Updated the 5-state xz plane system test to check the optimised parameter 
values.


Modified:
    branches/N_state_model/test_suite/system_tests/n_state_model.py

Modified: branches/N_state_model/test_suite/system_tests/n_state_model.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/N_state_model/test_suite/system_tests/n_state_model.py?rev=4896&r1=4895&r2=4896&view=diff
==============================================================================
--- branches/N_state_model/test_suite/system_tests/n_state_model.py (original)
+++ branches/N_state_model/test_suite/system_tests/n_state_model.py Mon Jan 
28 17:11:14 2008
@@ -21,7 +21,7 @@
 
###############################################################################
 
 # Python module imports.
-from math import sqrt
+from math import pi, sqrt
 from unittest import TestCase
 
 # relax module imports.
@@ -51,3 +51,34 @@
 
         # Execute the script.
         
self.relax.interpreter.run(script_file='test_suite/system_tests/scripts/5_state_xz.py')
+
+        # Alias the current data pipe.
+        cdp = relax_data_store[relax_data_store.current_pipe]
+
+        # Test the optimised probabilities.
+        self.assertAlmostEqual(cdp.probs[0], 0.2)
+        self.assertAlmostEqual(cdp.probs[1], 0.2)
+        self.assertAlmostEqual(cdp.probs[2], 0.2)
+        self.assertAlmostEqual(cdp.probs[3], 0.2)
+        self.assertAlmostEqual(cdp.probs[4], 0.2)
+
+        # Test the optimised alpha Euler angles.
+        self.assertAlmostEqual(cdp.alpha[0], 0.0)
+        self.assertAlmostEqual(cdp.alpha[1], 0.0)
+        self.assertAlmostEqual(cdp.alpha[2], 0.0)
+        self.assertAlmostEqual(cdp.alpha[3], 0.0)
+        self.assertAlmostEqual(cdp.alpha[4], 0.0)
+
+        # Test the optimised beta Euler angles.
+        self.assertAlmostEqual(cdp.beta[0], pi/4)
+        self.assertAlmostEqual(cdp.beta[1], pi/8)
+        self.assertAlmostEqual(cdp.beta[2], 0.0)
+        self.assertAlmostEqual(cdp.beta[3], -pi/8)
+        self.assertAlmostEqual(cdp.beta[4], -pi/4)
+
+        # Test the optimised gamma Euler angles.
+        self.assertAlmostEqual(cdp.gamma[0], 0.0)
+        self.assertAlmostEqual(cdp.gamma[1], 0.0)
+        self.assertAlmostEqual(cdp.gamma[2], 0.0)
+        self.assertAlmostEqual(cdp.gamma[3], 0.0)
+        self.assertAlmostEqual(cdp.gamma[4], 0.0)