Hi Edward,
I am running "dauvergne_protocol" script for my protein. The
"local_tm" and "sphere" models calculations proceeded fine, but with
prolate/oblate/ellipsoid diffusion tensors an exception is thrown (the
output is below). I have run the same script successfully before. The
difference is that I have excluded several mobile residues and using
"ave=False" in "structure.vectors" (the pdb file is an NMR ensemble,
20 models). If necessary, I can send the input files separately. Relax
revision is 14809.
Thank you,
Vitaly
Traceback (most recent call last):
File "/home/vitalyv/relax/multi/uni_processor.py", line 136, in run
self.callback.init_master(self)
File "/home/vitalyv/relax/multi/processor.py", line 263, in
default_init_master
self.master.run()
File "/home/vitalyv/relax/relax.py", line 174, in run
self.interpreter.run(self.script_file)
File "/home/vitalyv/relax/prompt/interpreter.py", line 300, in run
return run_script(intro=self.__intro_string, local=locals(),
script_file=script_file, quit=self.__quit_flag,
show_script=self.__show_script,
raise_relax_error=self.__raise_relax_error)
File "/home/vitalyv/relax/prompt/interpreter.py", line 610, in run_script
return console.interact(intro, local, script_file, quit,
show_script=show_script, raise_relax_error=raise_relax_error)
File "/home/vitalyv/relax/prompt/interpreter.py", line 495, in
interact_script
exec_script(script_file, local)
File "/home/vitalyv/relax/prompt/interpreter.py", line 383, in exec_script
runpy.run_module(module, globals)
File "/usr/lib/python2.7/runpy.py", line 180, in run_module
fname, loader, pkg_name)
File "/usr/lib/python2.7/runpy.py", line 72, in _run_code
exec code in run_globals
File "/home/vitalyv/relax/wtpln/dauvergne_protocol.py", line 218, in
<module>
dAuvergne_protocol(pipe_name=name, 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,
conv_loop=CONV_LOOP)
File "/home/vitalyv/relax/auto_analyses/dauvergne_protocol.py", line
223, in __init__
self.execute()
File "/home/vitalyv/relax/auto_analyses/dauvergne_protocol.py", line
596, in execute
self.interpreter.grid_search(inc=inc)
File "/home/vitalyv/relax/prompt/minimisation.py", line 106, in
grid_search
minimise.grid_search(lower=lower, upper=upper, inc=inc,
constraints=constraints, verbosity=verbosity)
File "/home/vitalyv/relax/generic_fns/minimise.py", line 146, in
grid_search
processor_box.processor.run_queue()
File "/home/vitalyv/relax/multi/uni_processor.py", line 149, in run_queue
command.run(self, completed)
File
"/home/vitalyv/relax/specific_fns/model_free/multi_processor_commands.py",
line 144, in run
raise Capturing_exception(rank=processor.rank(),
name=processor.get_name())
Capturing_exception:
------------------------------------------------------------------------------------------------------------------------
File
"/home/vitalyv/relax/specific_fns/model_free/multi_processor_commands.py",
line 132, in run
results = self.optimise()
File
"/home/vitalyv/relax/specific_fns/model_free/multi_processor_commands.py",
line 185, in optimise
results = grid_point_array(func=self.mf.func, args=(),
points=self.opt_params.subdivision,
verbosity=self.opt_params.verbosity)
File "/usr/lib/python2.7/site-packages/minfx/grid.py", line 282, in
grid_point_array
f = func(*(points[k],)+args)
File "/home/vitalyv/relax/maths_fns/mf.py", line 517, in func_diff
self.diff_data.calc_ci(data, self.diff_data)
File "/home/vitalyv/relax/maths_fns/weights.py", line 80, in
calc_spheroid_ci
data.ci[0] = 0.25 * data.three_dz2_one**2
Nested Exception from sub processor
Exception type: ValueError
Message: setting an array element with a sequence.