Changes

Jump to navigation Jump to search

Tutorial for model free SBiNLab

3,302 bytes added, 11:07, 14 October 2017
Number of Monte-Carlo simulations. The protocol will find optimum parameter values in this protocol, but error
estimation will not be very reliable. Standard is 500.
 
We use [http://www.dayid.org/comp/tm.html tmux] to make a terminal-session, we can get back to,
if our own terminal connection get closed.
 
* start a new session: '''tmux'''
* re-attach a detached session: '''tmux attach'''
Run with
<source lang="bash">
# Make terminal-session
tmux
 
relax 04_run_default_with_tolerance_lim.py -t 04_run_default_with_tolerance_lim.log
</source>
<source lang="bash">
less +F 04_run_default_with_tolerance_lim.log
</source>
 
* To scroll up and down, use keyboard: '''Ctrl+c'''
* To return to follow mode, use keyboard: '''Shift+f'''
* To exit, use keyboard: '''Ctrl+c''' and then: '''q'''
 
=== Try normal run with MC 20 ===
We have 4 CPU on our lab computers.<br>
Let us assign 1 to a run normal settings, and only MC=20.
 
Copy '''04_run_default_with_tolerance_lim.py''' to '''05_run_def_MC20.py''' and modify last lines:
 
<source lang="python">
# The number of Monte Carlo simulations to be used for error analysis at the end of the analysis.
#MC_NUM = 500
MC_NUM = 20
 
# The diffusion model. Standard is 'Fully automated', which means: DIFF_MODEL=['local_tm', 'sphere', 'prolate', 'oblate', 'ellipsoid', 'final']
# 'local_tm', 'sphere', ''prolate', 'oblate', 'ellipsoid', or 'final'
#DIFF_MODEL = 'local_tm'
DIFF_MODEL = ['local_tm', 'sphere', 'prolate', 'oblate', 'ellipsoid', 'final']
 
# The maximum number of iterations for the global iteration. Set to None, then the algorithm iterates until convergence.
MAX_ITER = None
 
# Automatic looping over all rounds until convergence (must be a boolean value of True or False). Standard is: conv_loop=True : and cannot be changed in the GUI.
CONV_LOOP = True
 
# Change some minimise opt params.
# This goes into: minimise.execute(self.min_algor, func_tol=self.opt_func_tol, max_iter=self.opt_max_iterations)
#####################
#dAuvergne_protocol.opt_func_tol = 1e-5 # Standard: opt_func_tol = 1e-25
#dAuvergne_protocol.opt_max_iterations = 1000 # Standard: opt_max_iterations = int(1e7)
#dAuvergne_protocol.opt_func_tol = 1e-10 # Standard: opt_func_tol = 1e-25
#dAuvergne_protocol.opt_max_iterations = int(1e5) # Standard: opt_max_iterations = int(1e7)
 
#####################################
 
# The results dir.
var = 'result_05'
results_dir = os.getcwd() + os.sep + var
 
# Save the state before running. Open and check in GUI!
state.save(state=var+'_ini.bz2', dir=results_dir, force=True)
 
# To check in GUI
# relax -g
# File -> Open relax state
# In folder "result_03" open "result_03_ini.bz2"
# View -> Data pipe editor
# Right click on pipe, and select "Associate with a new auto-analysis"
 
dAuvergne_protocol(pipe_name=name, pipe_bundle=bundle_name, results_dir=results_dir, diff_model=DIFF_MODEL, mf_models=MF_MODELS, local_tm_models=LOCAL_TM_MODELS, grid_inc=GRID_INC, min_algor=MIN_ALGOR, mc_sim_num=MC_NUM, max_iter=MAX_ITER, conv_loop=CONV_LOOP)
</source>
 
* MC_NUM = 20
Number of Monte-Carlo simulations. The protocol will find optimum parameter values in this protocol, but error
estimation will not be very reliable. Standard is 500.
 
We use [http://www.dayid.org/comp/tm.html tmux] to make a terminal-session, we can get back to,
if our own terminal connection get closed.
 
* start a new session: '''tmux'''
* re-attach a detached session: '''tmux attach'''
 
Run with
<source lang="bash">
# Make terminal-session
tmux
 
relax 05_run_def_MC20.py -t 05_run_def_MC20.log
</source>
 
You can then in another terminal follow the logfile by
<source lang="bash">
less +F 05_run_def_MC20.log
</source>
Trusted, Bureaucrats
1,382

edits

Navigation menu