mailr26866 - /trunk/test_suite/system_tests/scripts/relax_fit_exp_3param_inv_neg.py


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

Header


Content

Posted by edward on November 30, 2014 - 00:51:
Author: bugman
Date: Sun Nov 30 00:51:14 2014
New Revision: 26866

URL: http://svn.gna.org/viewcvs/relax?rev=26866&view=rev
Log:
Updated Séb's relax_fit_exp_3param_inv_neg.py system test script to work with 
the current relax design.


Modified:
    trunk/test_suite/system_tests/scripts/relax_fit_exp_3param_inv_neg.py

Modified: 
trunk/test_suite/system_tests/scripts/relax_fit_exp_3param_inv_neg.py
URL: 
http://svn.gna.org/viewcvs/relax/trunk/test_suite/system_tests/scripts/relax_fit_exp_3param_inv_neg.py?rev=26866&r1=26865&r2=26866&view=diff
==============================================================================
--- trunk/test_suite/system_tests/scripts/relax_fit_exp_3param_inv_neg.py     
  (original)
+++ trunk/test_suite/system_tests/scripts/relax_fit_exp_3param_inv_neg.py     
  Sun Nov 30 00:51:14 2014
@@ -5,7 +5,7 @@
 import sys
 
 # relax module imports.
-from data import Relax_data_store; ds = Relax_data_store()
+from data_store import Relax_data_store; ds = Relax_data_store()
 from status import Status; status = Status()
 
 
@@ -18,7 +18,7 @@
 pipe.create('rx', 'relax_fit')
 
 # The path to the data files.
-data_path = status.install_path + 
sep+'test_suite'+sep+'shared_data'+sep+'curve_fitting'
+data_path = status.install_path + 
sep+'test_suite'+sep+'shared_data'+sep+'curve_fitting'+sep+'inversion_recovery'
 
 # Load the sequence.
 sequence.read('Ap4Aase.seq', dir=status.install_path + 
sep+'test_suite'+sep+'shared_data', res_num_col=1, res_name_col=2)
@@ -75,22 +75,22 @@
             spectrum.delete(names[i])
 
 # Deselect unresolved spins.
-deselect.read(file='unresolved', dir=data_path, res_num_col=1)
+deselect.read(file='unresolved', dir=data_path+sep+'..', res_num_col=1)
 
 # Set the relaxation curve type.
-relax_fit.select_model('exp_2param_neg')
+relax_fit.select_model('inv')
 
 # Grid search.
-grid_search(inc=11)
+minimise.grid_search(inc=11)
 
 # Minimise.
-minimise('simplex', constraints=False)
+minimise.execute('newton', constraints=False)
 
 # Monte Carlo simulations.
 monte_carlo.setup(number=3)
 monte_carlo.create_data()
 monte_carlo.initial_values()
-minimise('simplex', constraints=False)
+minimise.execute('newton', constraints=False)
 monte_carlo.error_analysis()
 
 # Save the relaxation rates.
@@ -103,8 +103,8 @@
 grace.write(y_data_type='chi2', file='devnull', force=True)    # Minimised 
chi-squared value.
 grace.write(y_data_type='i0', file='devnull', force=True)    # Initial peak 
intensity.
 grace.write(y_data_type='rx', file='devnull', force=True)    # Relaxation 
rate.
-grace.write(x_data_type='relax_times', y_data_type='intensities', 
file='devnull', force=True)    # Average peak intensities.
-grace.write(x_data_type='relax_times', y_data_type='intensities', norm=True, 
file='devnull', force=True)    # Average peak intensities (normalised).
+grace.write(x_data_type='relax_times', y_data_type='peak_intensity', 
file='devnull', force=True)    # Average peak intensities.
+grace.write(x_data_type='relax_times', y_data_type='peak_intensity', 
norm=True, file='devnull', force=True)    # Average peak intensities 
(normalised).
 
 # Save the program state.
 state.save('devnull', force=True)




Related Messages


Powered by MHonArc, Updated Sun Nov 30 01:00:02 2014