r27221 - /trunk/test_suite/system_tests/relax_disp.py -- January 19, 2015 - 13:48

Author: tlinnet
Date: Mon Jan 19 13:48:41 2015
New Revision: 27221

URL: http://svn.gna.org/viewcvs/relax?rev=27221&view=rev
Log:
Change to systemtest Relax_disp.x_test_task_7882_kex_conf()

This is just a temporary systemtest, to check for local minima.

This is method in regression book of Graphpad: 
http://www.graphpad.com/faq/file/Prism4RegressionBook.pdf
Page: 109-111.

Task #7882 (https://gna.org/task/?7882): Implement Monte-Carlo simulation, 
where errors are generated with width of standard deviation or residuals.

Modified:
    trunk/test_suite/system_tests/relax_disp.py

Modified: trunk/test_suite/system_tests/relax_disp.py
URL: 
http://svn.gna.org/viewcvs/relax/trunk/test_suite/system_tests/relax_disp.py?rev=27221&r1=27220&r2=27221&view=diff
==============================================================================
--- trunk/test_suite/system_tests/relax_disp.py (original)
+++ trunk/test_suite/system_tests/relax_disp.py Mon Jan 19 13:48:41 2015
@@ -8621,7 +8621,7 @@
 
         # Now generate a list of kex values to try.
         kex_cur = cur_spin.kex
-        kex_list = linspace(kex_cur - 1500, kex_cur + 3000, 1000)
+        kex_list = linspace(kex_cur - 1500, kex_cur + 3000, 200)
 
         chi2_list = []
 
@@ -8653,11 +8653,11 @@
         print(kex_sel)
         print(chi2_sel)
 
-        if False:
+        if True:
             import matplotlib.pyplot as plt
 
-            plt.plot(kex_sel, chi2_sel)
-            plt.plot(SSbest_kex, SSbest_fit)
+            plt.plot(kex_sel, chi2_sel, "bo")
+            plt.plot(SSbest_kex, SSbest_fit, "g*")
             plt.show()