Hi again, I just ran a test with the minfx library instead of the "minimise" directory, even though the test-suite fails for the two last model-free tests. As expected, the script failed when it reached the minimise() function... =========================== relax> minimise(*args=('newton',), run='final', func_tol=1e-25, max_iterations=10000000, constraints=1, scaling=1, print_flag=1) Simulation 1 Traceback (most recent call last): File "/home/semor/packages/relax-1.2/relax", line 454, in ? Relax() File "/home/semor/packages/relax-1.2/relax", line 166, in __init__ self.interpreter.run() File "/home/semor/packages/relax-1.2/prompt/interpreter.py", line 216, in run run_script(intro=self.relax.intro_string, local=self.local, script_file=self.relax.script_file, quit=1) File "/home/semor/packages/relax-1.2/prompt/interpreter.py", line 391, in run_script console.interact(intro, local, script_file, quit) File "/home/semor/packages/relax-1.2/prompt/interpreter.py", line 343, in interact_script execfile(script_file, local) File "./full_analysis_final.py", line 587, in ? Main(self.relax) File "./full_analysis_final.py", line 313, in __init__ minimise(MIN_ALGOR, run='final') File "/home/semor/packages/relax-1.2/prompt/minimisation.py", line 365, in minimise self.relax.generic.minimise.minimise(run=run, min_algor=min_algor, min_options=min_options, func_tol=func_tol, grad_tol=grad_tol, max_iterations=max_iterations, constraints=constraints, scaling=scaling, print_flag=print_flag) File "/home/semor/packages/relax-1.2/generic_fns/minimise.py", line 135, in minimise minimise(run=run, min_algor=min_algor, min_options=min_options, func_tol=func_tol, grad_tol=grad_tol, max_iterations=max_iterations, constraints=constraints, scaling=scaling, print_flag=print_flag-1, sim_index=i) File "/home/semor/packages/relax-1.2/specific_fns/model_free.py", line 2372, in minimise results = generic_minimise(func=self.mf.func, dfunc=self.mf.dfunc, d2func=self.mf.d2func, args=(), x0=self.param_vector, min_algor=min_algor, min_options=min_options, func_tol=func_tol, grad_tol=grad_tol, maxiter=max_iterations, A=A, b=b, full_output=1, print_flag=print_flag) File "/home/semor/packages/relax-1.3/minfx/generic.py", line 392, in generic_minimise results = method_of_multipliers(func=func, dfunc=dfunc, d2func=d2func, args=args, x0=x0, min_options=min_options, A=A, b=b, l=l, u=u, c=c, dc=dc, d2c=d2c, func_tol=func_tol, grad_tol=grad_tol, maxiter=maxiter, full_output=full_output, print_flag=print_flag) File "/home/semor/packages/relax-1.3/minfx/method_of_multipliers.py", line 143, in method_of_multipliers results = min.minimise() File "/home/semor/packages/relax-1.3/minfx/method_of_multipliers.py", line 397, in minimise results = self.generic_minimise(func=self.func_LA, dfunc=self.func_dLA, d2func=self.func_d2LA, args=self.args, x0=self.xk, min_algor=self.min_algor, min_options=self.min_options, func_tol=None, grad_tol=self.tk, maxiter=maxiter, full_output=1, print_flag=sub_print_flag, print_prefix="\t") File "/home/semor/packages/relax-1.3/minfx/generic.py", line 323, in generic_minimise results = newton(func=func, dfunc=dfunc, d2func=d2func, args=args, x0=x0, min_options=min_options, func_tol=func_tol, grad_tol=grad_tol, maxiter=maxiter, full_output=full_output, print_flag=print_flag, print_prefix=print_prefix) File "/home/semor/packages/relax-1.3/minfx/newton.py", line 44, in newton results = min.minimise() File "/home/semor/packages/relax-1.3/minfx/base_classes.py", line 233, in minimise self.new_param_func() File "/home/semor/packages/relax-1.3/minfx/newton.py", line 163, in new_param_func self.pk = self.get_pk() File "/home/semor/packages/relax-1.3/minfx/base_classes.py", line 642, in gmw return gmw(self.dfk, self.d2fk, self.I, self.n, self.mach_acc, self.print_prefix, self.print_flag, return_matrix) File "/home/semor/packages/relax-1.3/minfx/hessian_mods/gmw81.py", line 150, in gmw y = solve(L, dfk) File "/usr/lib/python2.4/site-packages/numpy/linalg/linalg.py", line 179, in solve t, result_t = _commonType(a, b) File "/usr/lib/python2.4/site-packages/numpy/linalg/linalg.py", line 71, in _commonType if issubclass(a.dtype.type, inexact): AttributeError: dtype =========================== I'm not familiar with the minfx package, but this error as probably nothing to do with the constraints problem... More to come... Stay tuned... Séb Sébastien Morin wrote:
|