mailr19988 - /branches/relax_disp/specific_analyses/relax_disp/__init__.py


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

Header


Content

Posted by edward on June 09, 2013 - 13:38:
Author: bugman
Date: Sun Jun  9 13:38:17 2013
New Revision: 19988

URL: http://svn.gna.org/viewcvs/relax?rev=19988&view=rev
Log:
Added a printout at the end of the optimisation of the final dispersion 
parameter values.


Modified:
    branches/relax_disp/specific_analyses/relax_disp/__init__.py

Modified: branches/relax_disp/specific_analyses/relax_disp/__init__.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/specific_analyses/relax_disp/__init__.py?rev=19988&r1=19987&r2=19988&view=diff
==============================================================================
--- branches/relax_disp/specific_analyses/relax_disp/__init__.py (original)
+++ branches/relax_disp/specific_analyses/relax_disp/__init__.py Sun Jun  9 
13:38:17 2013
@@ -57,8 +57,8 @@
 from pipe_control.result_files import add_result_file
 from specific_analyses.api_base import API_base
 from specific_analyses.api_common import API_common
-from specific_analyses.relax_disp.disp_data import average_intensity, 
find_intensity_keys, loop_cluster, loop_frq, loop_frq_point, 
loop_frq_point_key, loop_frq_point_time, loop_point, loop_time, relax_time, 
return_cpmg_frqs, return_index_from_disp_point, return_index_from_frq, 
return_key_from_disp_point_index, return_param_key_from_data, 
return_r2eff_arrays, return_spin_lock_nu1, spin_ids_to_containers
-from specific_analyses.relax_disp.parameters import assemble_param_vector, 
assemble_scaling_matrix, disassemble_param_vector, linear_constraints, 
param_index_to_param_info, param_num
+from specific_analyses.relax_disp.disp_data import average_intensity, 
find_intensity_keys, loop_cluster, loop_frq, loop_frq_point, 
loop_frq_point_key, loop_frq_point_time, loop_point, loop_time, relax_time, 
return_cpmg_frqs, return_index_from_disp_point, return_index_from_frq, 
return_key_from_disp_point_index, return_param_key_from_data, 
return_r2eff_arrays, return_spin_lock_nu1, return_value_from_frq_index, 
spin_ids_to_containers
+from specific_analyses.relax_disp.parameters import assemble_param_vector, 
assemble_scaling_matrix, disassemble_param_vector, linear_constraints, 
loop_parameters, param_index_to_param_info, param_num
 from specific_analyses.relax_disp.variables import CPMG_EXP, FIXED_TIME_EXP, 
MODEL_LIST_FULL, MODEL_LM63, MODEL_CR72, MODEL_M61, MODEL_M61B, MODEL_NOREX, 
MODEL_R2EFF, R1RHO_EXP, VAR_TIME_EXP
 from target_functions.relax_disp import Dispersion
 from user_functions.data import Uf_tables; uf_tables = Uf_tables()
@@ -1432,6 +1432,21 @@
                         # Store the back-calculated data.
                         spin.r2eff_bc[key] = model.back_calc[spin_index, 
frq_index, disp_pt_index]
 
+            # Optimisation printout.
+            if verbosity:
+                print("\nOptimised parameter values:")
+                for param_name, param_index, spin_index, frq_index in 
loop_parameters(spins=spins):
+                    # The parameter with additional details.
+                    param_text = param_name
+                    if param_name == 'r2':
+                        frq = return_value_from_frq_index(frq_index)
+                        if frq:
+                            param_text += " (%.3f MHz)" % (frq / 1e6) 
+                    param_text += ":"
+
+                    # The printout.
+                    print("%-20s %20.15g" % (param_text, 
param_vector[param_index]))
+
 
     def model_desc(self, model_info):
         """Return a description of the model.




Related Messages


Powered by MHonArc, Updated Sun Jun 09 14:00:02 2013