mailr14914 - /1.3/test_suite/system_tests/structure.py


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

Header


Content

Posted by edward on October 25, 2011 - 14:58:
Author: bugman
Date: Tue Oct 25 14:58:38 2011
New Revision: 14914

URL: http://svn.gna.org/viewcvs/relax?rev=14914&view=rev
Log:
Expanded the Structure.test_displacement() system test to test the XML save 
files.


Modified:
    1.3/test_suite/system_tests/structure.py

Modified: 1.3/test_suite/system_tests/structure.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/test_suite/system_tests/structure.py?rev=14914&r1=14913&r2=14914&view=diff
==============================================================================
--- 1.3/test_suite/system_tests/structure.py (original)
+++ 1.3/test_suite/system_tests/structure.py Tue Oct 25 14:58:38 2011
@@ -23,6 +23,7 @@
 # Python module imports.
 from numpy import float64, zeros
 from os import sep
+from tempfile import mktemp
 
 # relax module imports.
 from base_classes import SystemTestCase
@@ -60,18 +61,35 @@
         path = status.install_path + 
sep+'test_suite'+sep+'shared_data'+sep+'structures'
 
         # Load the file as two models.
-        self.interpreter.structure.read_pdb('Ap4Aase_res1-12.pdb', dir=path, 
set_model_num=0)
         self.interpreter.structure.read_pdb('Ap4Aase_res1-12.pdb', dir=path, 
set_model_num=1)
+        self.interpreter.structure.read_pdb('Ap4Aase_res1-12.pdb', dir=path, 
set_model_num=2)
 
         # A rotation.
         R = zeros((3, 3), float64)
         euler_to_R_zyz(1.3, 0.4, 4.5, R)
 
         # Rotate the second model.
-        self.interpreter.structure.rotate(R, model=1)
+        self.interpreter.structure.rotate(R, model=2)
 
         # Calculate the displacement.
         self.interpreter.structure.displacement()
+
+        # Shift a third structure back using the calculated displacement.
+        self.interpreter.structure.read_pdb('Ap4Aase_res1-12.pdb', dir=path, 
set_model_num=3)
+        self.interpreter.structure.rotate(R, model=3)
+
+        # Save the results.
+        self.tmpfile = mktemp()
+        self.interpreter.state.save(self.tmpfile, dir=None, force=True)
+
+        # Reset relax.
+        self.interpreter.reset()
+
+        # Load the results.
+        self.interpreter.state.load(self.tmpfile)
+
+        # Test the loading.
+        self.assert_(hasattr(cdp.structure, 'displacements'))
 
 
     def test_load_internal_results(self):




Related Messages


Powered by MHonArc, Updated Tue Oct 25 16:40:02 2011