r7516 - /1.3/test_suite/system_tests/scripts/omp_model_free.py -- October 05, 2008 - 16:51

Author: bugman
Date: Sun Oct  5 16:51:58 2008
New Revision: 7516

URL: http://svn.gna.org/viewcvs/relax?rev=7516&view=rev
Log:
Extended the system test of a minimal model-free analysis using the OMP 
relaxation data.

This now includes the basics of the 'final' stage of the 'full_analysis.py' 
script.


Modified:
    1.3/test_suite/system_tests/scripts/omp_model_free.py

Modified: 1.3/test_suite/system_tests/scripts/omp_model_free.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/test_suite/system_tests/scripts/omp_model_free.py?rev=7516&r1=7515&r2=7516&view=diff
==============================================================================
--- 1.3/test_suite/system_tests/scripts/omp_model_free.py (original)
+++ 1.3/test_suite/system_tests/scripts/omp_model_free.py Sun Oct  5 16:51:58 
2008
@@ -25,6 +25,9 @@
 
 # The optimisation technique.
 MIN_ALGOR = 'simplex'
+
+# Number of Monte Carlo simulations.
+MC_NUM = 3
 
 
 class Main:
@@ -87,6 +90,30 @@
         self.model_selection(pipe='aic')
 
 
+        # Final stage.
+        ##############
+
+        # Unfix all parameters (to switch to the global models).
+        fix('all', fixed=False)
+
+        # Model selection between MI to MV.
+        self.model_selection(pipe='final')
+
+        # Fix the diffusion tensor.
+        fix('diff')
+
+        # Monte Carlo simulations.
+        monte_carlo.setup(number=MC_NUM)
+        monte_carlo.create_data()
+        monte_carlo.initial_values()
+        minimise(MIN_ALGOR)
+        eliminate()
+        monte_carlo.error_analysis()
+
+        # Write the final results.
+        results.write(file='devnull', force=True)
+
+
     def model_selection(self, pipe=None):
         """Model selection function."""