URL:
<http://gna.org/bugs/?22033>
Summary: Not possible to minimise with "Levenburg-Marquardt".
Project: relax
Submitted by: tlinnet
Submitted on: Sun 11 May 2014 06:37:03 PM UTC
Category: relax's source code
Specific analysis category: Relaxation curve-fitting
Priority: 3 - Low
Severity: 3 - Normal
Status: None
Assigned to: None
Originator Name:
Originator Email:
Open/Closed: Open
Release: Repository: trunk
Discussion Lock: Any
Operating System: All systems
_______________________________________________________
Details:
I tried to minimise with LM:
minimise(min_algor='LM', func_tol=ds.set_func_tol, max_iter=ds.set_max_iter,
constraints=True, scaling=True, verbosity=ds.verbosity)
relax> minimise(min_algor='LM', line_search=None, hessian_mod=None,
hessian_type=None, func_tol=1e-25, grad_tol=None, max_iter=10000000,
constraints=True, scaling=True, verbosity=1)
Fitting to the spin block [':1@N']
----------------------------------
Logarithmic barrier function
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
k: 0 xk: [ 2.504123959, 2.402856606, 0.5,
0.975, 0.20008] fk: 1.76614984356
Entering sub-algorithm.
Traceback (most recent call last):
File "/sbinlab2/software/NMR-relax/relax_trunk/multi/processor.py", line
494, in run
self.callback.init_master(self)
File "/sbinlab2/software/NMR-relax/relax_trunk/multi/__init__.py", line 318,
in default_init_master
self.master.run()
File "/sbinlab2/software/NMR-relax/relax_trunk/relax.py", line 194, in run
self.interpreter.run(self.script_file)
File "/sbinlab2/software/NMR-relax/relax_trunk/prompt/interpreter.py", line
275, in run
return run_script(intro=self.__intro_string, local=locals(),
script_file=script_file, show_script=self.__show_script,
raise_relax_error=self.__raise_relax_error)
File "/sbinlab2/software/NMR-relax/relax_trunk/prompt/interpreter.py", line
569, in run_script
return console.interact(intro, local, script_file,
show_script=show_script, raise_relax_error=raise_relax_error)
File "/sbinlab2/software/NMR-relax/relax_trunk/prompt/interpreter.py", line
468, in interact_script
exec_script(script_file, local)
File "/sbinlab2/software/NMR-relax/relax_trunk/prompt/interpreter.py", line
347, in exec_script
runpy.run_module(module, globals)
File "/usr/lib64/python2.6/runpy.py", line 140, in run_module
fname, loader, pkg_name)
File "/usr/lib64/python2.6/runpy.py", line 34, in _run_code
exec code in run_globals
File
"/sbinlab2/tlinnet/phd/relax/chi2_problems/CR72/10_dw_error_Levenburg_Marquardt/CR72_10_dw_error_10_dw_error_Levenburg_Marquardt.py",
line 476, in <module>
minimise(min_algor='LM', func_tol=ds.set_func_tol,
max_iter=ds.set_max_iter, constraints=True, scaling=True,
verbosity=ds.verbosity)
File "/sbinlab2/software/NMR-relax/relax_trunk/prompt/uf_objects.py", line
223, in __call__
self._backend(*new_args, **uf_kargs)
File "/sbinlab2/software/NMR-relax/relax_trunk/pipe_control/minimise.py",
line 238, in minimise
processor.run_queue()
File "/sbinlab2/software/NMR-relax/relax_trunk/multi/uni_processor.py", line
184, in run_queue
command.run(self, completed)
File
"/sbinlab2/software/NMR-relax/relax_trunk/specific_analyses/relax_disp/optimisation.py",
line 755, in run
results = generic_minimise(func=model.func, args=(), x0=self.param_vector,
min_algor=self.min_algor, min_options=self.min_options,
func_tol=self.func_tol, grad_tol=self.grad_tol, maxiter=self.max_iterations,
A=self.A, b=self.b, full_output=True, print_flag=self.verbosity)
File "/usr/lib/python2.6/site-packages/minfx/generic.py", line 399, in
generic_minimise
results = log_barrier_function(func=func, dfunc=dfunc, d2func=d2func,
args=args, x0=x0, min_options=min_options, A=A, b=b, func_tol=func_tol,
grad_tol=grad_tol, maxiter=maxiter, full_output=full_output,
print_flag=print_flag)
File "/usr/lib/python2.6/site-packages/minfx/log_barrier_function.py", line
96, in log_barrier_function
results = min.minimise()
File "/usr/lib/python2.6/site-packages/minfx/log_barrier_function.py", line
264, in minimise
results = self.generic_minimise(func=self.func_log, dfunc=self.func_dlog,
d2func=self.func_d2log, args=self.args, x0=self.xk, min_algor=self.min_algor,
min_options=self.min_options, func_tol=self.func_tol, grad_tol=self.grad_tol,
maxiter=maxiter, full_output=1, print_flag=self.print_flag,
print_prefix="\t")
File "/usr/lib/python2.6/site-packages/minfx/generic.py", line 387, in
generic_minimise
results = levenberg_marquardt(chi2_func=func, dchi2_func=dfunc,
dfunc=min_options[0], errors=min_options[1], args=args, x0=x0,
func_tol=func_tol, grad_tol=grad_tol, maxiter=maxiter,
full_output=full_output, print_flag=print_flag, print_prefix=print_prefix)
IndexError: tuple index out of range
_______________________________________________________
Reply to this item at:
<http://gna.org/bugs/?22033>
_______________________________________________
Message sent via/by Gna!
http://gna.org/
[bug #22033] Not possible to minimise with "Levenburg-Marquardt".