mailr2882 - /1.3/sample_scripts/full_analysis.py


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

Header


Content

Posted by edward . dauvergne on November 23, 2006 - 09:08:
Author: bugman
Date: Thu Nov 23 09:07:41 2006
New Revision: 2882

URL: http://svn.gna.org/viewcvs/relax?rev=2882&view=rev
Log:
Improvements to the printed out convergence test messages.


Modified:
    1.3/sample_scripts/full_analysis.py

Modified: 1.3/sample_scripts/full_analysis.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/sample_scripts/full_analysis.py?rev=2882&r1=2881&r2=2882&view=diff
==============================================================================
--- 1.3/sample_scripts/full_analysis.py (original)
+++ 1.3/sample_scripts/full_analysis.py Thu Nov 23 09:07:41 2006
@@ -322,20 +322,20 @@
         # Chi-squared test.
         ###################
 
-        print "# Chi-squared test.\n"
-        print "chi2 (k-1): " + `self.relax.data.chi2['previous']`
-        print "chi2 (k):   " + `self.relax.data.chi2[run]`
+        print "Chi-squared test:"
+        print "    chi2 (k-1): " + `self.relax.data.chi2['previous']`
+        print "    chi2 (k):   " + `self.relax.data.chi2[run]`
         if self.relax.data.chi2['previous'] == self.relax.data.chi2[run]:
-            print "The chi-squared value has converged.\n"
-        else:
-            print "The chi-squared value has not converged.\n"
+            print "    The chi-squared value has converged.\n"
+        else:
+            print "    The chi-squared value has not converged.\n"
             chi2_converged = 0
 
 
         # Identical model-free model test.
         ##################################
 
-        print "# Identical model-free models test."
+        print "Identical model-free models test:"
 
         # Create a string representation of the model-free models of the 
previous run.
         prev_models = ''
@@ -349,16 +349,16 @@
 
         # The test.
         if prev_models == curr_models:
-            print "The model-free models have converged.\n"
-        else:
-            print "The model-free models have not converged.\n"
+            print "    The model-free models have converged.\n"
+        else:
+            print "    The model-free models have not converged.\n"
             models_converged = 0
 
 
         # Identical parameter value test.
         #################################
 
-        print "# Identical parameter test."
+        print "Identical parameter test:"
 
         # Only run the tests if the model-free models have converged.
         if models_converged:
@@ -378,10 +378,10 @@
 
                 # Test if not identical.
                 if prev_val != curr_val:
-                    print "Parameter:   " + param
-                    print "Value (k-1): " + `prev_val`
-                    print "Value (k):   " + `curr_val`
-                    print "The diffusion parameters have not converged.\n"
+                    print "    Parameter:   " + param
+                    print "    Value (k-1): " + `prev_val`
+                    print "    Value (k):   " + `curr_val`
+                    print "    The diffusion parameters have not 
converged.\n"
                     params_converged = 0
 
             # Skip the rest if the diffusion tensor parameters have not 
converged.
@@ -400,24 +400,28 @@
 
                         # Test if not identical.
                         if prev_val != curr_val:
-                            print "Spin system: " + 
`self.relax.data.res[run][i].num` + ' ' + self.relax.data.res[run][i].name
-                            print "Parameter:   " + 
self.relax.data.res[run][i].params[j]
-                            print "Value (k-1): " + `prev_val`
-                            print "Value (k):   " + `curr_val`
-                            print "The model-free parameters have not 
converged.\n"
+                            print "    Spin system: " + 
`self.relax.data.res[run][i].num` + ' ' + self.relax.data.res[run][i].name
+                            print "    Parameter:   " + 
self.relax.data.res[run][i].params[j]
+                            print "    Value (k-1): " + `prev_val`
+                            print "    Value (k):   " + `curr_val`
+                            print "    The model-free parameters have not 
converged.\n"
                             params_converged = 0
                             break
 
         # The model-free models haven't converged hence the parameter values 
haven't converged.
         else:
-            print "The model-free models haven't converged hence the 
parameters haven't converged."
+            print "    The model-free models haven't converged hence the 
parameters haven't converged.\n"
             params_converged = 0
+
+        # Print out.
+        if params_converged:
+            print "    The diffusion tensor and model-free parameters have 
converged.\n"
 
 
         # Final print out.
         ##################
 
-        print "\n# Convergence:"
+        print "\nConvergence:"
         if chi2_converged and models_converged and params_converged:
             print "    [ Yes ]"
         else:




Related Messages


Powered by MHonArc, Updated Thu Nov 23 09:20:06 2006