mailr4812 - /branches/N_state_model/test_suite/unit_tests/_specific_fns/test_n_state_model.py


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

Header


Content

Posted by edward on January 17, 2008 - 11:40:
Author: bugman
Date: Thu Jan 17 11:40:50 2008
New Revision: 4812

URL: http://svn.gna.org/viewcvs/relax?rev=4812&view=rev
Log:
Created a unit test for the N-state model disassemble_param_vector() method.


Modified:
    
branches/N_state_model/test_suite/unit_tests/_specific_fns/test_n_state_model.py

Modified: 
branches/N_state_model/test_suite/unit_tests/_specific_fns/test_n_state_model.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/N_state_model/test_suite/unit_tests/_specific_fns/test_n_state_model.py?rev=4812&r1=4811&r2=4812&view=diff
==============================================================================
--- 
branches/N_state_model/test_suite/unit_tests/_specific_fns/test_n_state_model.py
 (original)
+++ 
branches/N_state_model/test_suite/unit_tests/_specific_fns/test_n_state_model.py
 Thu Jan 17 11:40:50 2008
@@ -61,3 +61,28 @@
         for i in xrange(len(param_vector)):
             self.assertEqual(param_vector[i], vector_true[i])
 
+
+    def test_disassemble_param_vector(self):
+        """Test the operation of the 
specific_fns.n_state_model.disassemble_param_vector() method."""
+
+        # Alias the current data pipe.
+        cdp = relax_data_store[relax_data_store.current_pipe]
+
+        # Set up the initial N, probabilities and Euler angles.
+        cdp.N = 3
+        cdp.probs = [None]*3
+        cdp.alpha = [None]*3
+        cdp.beta = [None]*3
+        cdp.gamma = [None]*3
+
+        # The parameter vector.
+        param_vector = [0.1, 0.3, 0.0, pi/2, 1.0, pi/2, pi, 3*pi/2, pi, 
3*pi/2, 2*pi]
+
+        # Disassemble the parameter vector.
+        self.n_state_model_fns.disassemble_param_vector(param_vector)
+
+        # Check the probabilities.
+        self.assertEqual(cdp.probs, [0.1, 0.3, 0.6])
+        self.assertEqual(cdp.alpha, [0.0, pi/2, pi])
+        self.assertEqual(cdp.beta, [pi/2, pi, 3*pi/2])
+        self.assertEqual(cdp.gamma, [1.0, 3*pi/2, 2*pi])




Related Messages


Powered by MHonArc, Updated Thu Jan 17 12:00:17 2008