On 16 October 2015 at 17:42, Sze Chan <samuelsw.c...@mail.utoronto.ca> wrote: > Hello Edward, > > The problem to load into @N* appears when I use my full data, but using a > truncated version does not cause a problem. > > I attached, the truncated data, full data, and the PDB I've been using (its a > truncated version of a PDB from the RCSB) in the bug report.
Cheers! I've create a system test and now fixed the problem in relax. See https://gna.org/bugs/?23933#comment3 . The problem was introduced in Nov. 2014. It might take me a while to release a new version of minfx, bmrblib, and relax to have all these issues resolved. Anyway, you can ignore the comparison to 'None' warnings, as these are harmless for now. And they are fixed in minfx, bmrblib, and relax and so will be part of relax 4.0.1. With the fixes, the relax log will look like: """ relax> pipe.create(pipe_name='mf', pipe_type='mf', bundle=None) relax> structure.read_pdb(file='LARA_N_term_no_helixFH_reg.pdb', dir='/data/relax/relax-trunk/test_suite/shared_data/model_free/bug_23933_relax_data_read_ids', read_mol=None, set_mol_name=None, read_model=None, set_model_num=None, alt_loc=None, verbosity=1, merge=False) Internal relax PDB parser. Opening the file '/data/relax/relax-trunk/test_suite/shared_data/model_free/bug_23933_relax_data_read_ids/LARA_N_term_no_helixFH_reg.pdb' for reading. Adding molecule 'LARA_N_term_no_helixFH_reg_mol24' (from the original molecule number 24). relax> structure.load_spins(spin_id='@N', from_mols=None, mol_name_target=None, ave_pos=True) Adding the following spins to the relax data store. # mol_name res_num res_name spin_num spin_name LARA_N_term_no_helixFH_reg_mol24 329 GLN 1 N LARA_N_term_no_helixFH_reg_mol24 330 GLN 10 N LARA_N_term_no_helixFH_reg_mol24 331 SER 19 N relax> structure.load_spins(spin_id='@H', from_mols=None, mol_name_target=None, ave_pos=True) Adding the following spins to the relax data store. # mol_name res_num res_name spin_num spin_name LARA_N_term_no_helixFH_reg_mol24 330 GLN 0 H LARA_N_term_no_helixFH_reg_mol24 331 SER 0 H relax> structure.load_spins(spin_id='@NE1', from_mols=None, mol_name_target=None, ave_pos=True) Adding the following spins to the relax data store. RelaxWarning: No spins matching the '@NE1' ID string could be found. relax> structure.load_spins(spin_id='@HE1', from_mols=None, mol_name_target=None, ave_pos=True) Adding the following spins to the relax data store. RelaxWarning: No spins matching the '@HE1' ID string could be found. relax> relax_data.read(ri_id='R1_600', ri_type='R1', frq=600402816.0, file='r1_600.txt', dir='/data/relax/relax-trunk/test_suite/shared_data/model_free/bug_23933_relax_data_read_ids', spin_id_col=None, mol_name_col=None, res_num_col=1, res_name_col=None, spin_num_col=None, spin_name_col=None, data_col=2, error_col=3, sep=None, spin_id=None) Opening the file '/data/relax/relax-trunk/test_suite/shared_data/model_free/bug_23933_relax_data_read_ids/r1_600.txt' for reading. Traceback (most recent call last): File "/data/relax/relax-trunk/test_suite/system_tests/model_free.py", line 499, in test_bug_23933_relax_data_read_ids self.interpreter.relax_data.read(ri_id='R1_600', ri_type='R1', frq=600402816.0, file='r1_600.txt', dir=path, res_num_col=1, data_col=2, error_col=3) File "/data/relax/relax-trunk/prompt/uf_objects.py", line 225, in __call__ self._backend(*new_args, **uf_kargs) File "/data/relax/relax-trunk/pipe_control/relax_data.py", line 925, in read pack_data(ri_id, ri_type, frq, values, errors, mol_names=mol_names, res_nums=res_nums, res_names=res_names, spin_nums=spin_nums, spin_names=spin_names, spin_id=spin_id) File "/data/relax/relax-trunk/pipe_control/relax_data.py", line 773, in pack_data raise RelaxMultiSpinIDError(spin_ids[i], new_ids) RelaxMultiSpinIDError: RelaxError: The spin ID '#LARA_N_term_no_helixFH_reg_mol24:331' corresponds to multiple spins, including '#LARA_N_term_no_helixFH_reg_mol24:331@N' and '#LARA_N_term_no_helixFH_reg_mol24:331@H'. """ The bug that was hidden from you in the GUI will now stop relax and pop up a RelaxError dialog. Therefore you don't need the new relax version or use a checked out copy of the source code repository (unless you'd like the new code to remove all of the numpy FutureWarning messages). You just need to set the spin_id argument to '@N'. This allows relax to understand that you have nitrogen and not proton relaxation data - relax will accept both and it is currently ambiguous if it is 'N' or 'H' data. Anyway, thank you for the bug report and attached files. This has allowed me to fix both the numpy FutureWarning and missing 'ids' variable problems. If you see any other strange warnings or errors, it would be appreciated if you could report them so that they can be resolved for future relax versions. Cheers, Edward _______________________________________________ relax (http://www.nmr-relax.com) This is the relax-users mailing list relax-users@gna.org To unsubscribe from this list, get a password reminder, or change your subscription options, visit the list information page at https://mail.gna.org/listinfo/relax-users