Hi Edward. ./relax -s Relax_disp.test_bug_21665_cpmg_two_fields_two_delaytimes_fail_relax_disp is now passing the test until now. But I am unsure how to handle: ./relax -s Relax_disp.test_bug_21665_cpmg_two_fields_two_delaytimes_fail_calc File "/sbinlab2/tlinnet/software/NMR-relax/relax_trunk/specific_analyses/relax_disp/checks.py", line 112, in check_exp_type_fixed_time raise RelaxError("The experiment '%s' is not of the fixed relaxation time period data type." % exp_type) RelaxError: RelaxError: The experiment 'SQ CPMG' is not of the fixed relaxation time period data type. ------ The code is: def check_exp_type_fixed_time(): """Check that only fixed time experiment types have been set up. @raises RelaxError: If exponential curves are present. """ # Get the times. times = get_times() # Loop over all experiment types. for exp_type in times: if len(times[exp_type]) > 1: raise RelaxError("The experiment '%s' is not of the fixed relaxation time period data type." % exp_type) And return {'SQ CPMG': [0.040000000000000001, 0.059999999999999998]} Should I modify the code to loop over all spectra IDs, and count the number of times insted for each spectrum? ---------- Forwarded message ---------- From: <tlinnet@xxxxxxxxxxxxx> Date: 2014-02-24 18:28 GMT+01:00 Subject: r22294 - in /trunk/specific_analyses/relax_disp: api.py disp_data.py To: relax-commits@xxxxxxx Author: tlinnet Date: Mon Feb 24 18:28:11 2014 New Revision: 22294 URL: http://svn.gna.org/viewcvs/relax?rev=22294&view=rev Log: Replacing looping over time points from cdp.relax_time_list to loop_time(frq=frq). Regarding bug #21665, (https://gna.org/bugs/?21665) - Running a CPMG analysis with two fields at two delay times. loop_time has been modified to accept spectrometer frequency as input to restrict the looping. Modified: trunk/specific_analyses/relax_disp/api.py trunk/specific_analyses/relax_disp/disp_data.py Modified: trunk/specific_analyses/relax_disp/api.py URL: http://svn.gna.org/viewcvs/relax/trunk/specific_analyses/relax_disp/api.py?rev=22294&r1=22293&r2=22294&view=diff ============================================================================== --- trunk/specific_analyses/relax_disp/api.py (original) +++ trunk/specific_analyses/relax_disp/api.py Mon Feb 24 18:28:11 2014 @@ -149,7 +149,7 @@ values = [] errors = [] times = [] - for time in cdp.relax_time_list: + for time in loop_time(frq=frq): # The data. values.append(average_intensity(spin=spin, exp_type=exp_type, frq=frq, offset=offset, point=point, time=time)) errors.append(average_intensity(spin=spin, exp_type=exp_type, frq=frq, offset=offset, point=point, time=time, error=True)) @@ -378,7 +378,7 @@ values = [] errors = [] times = [] - for time in cdp.relax_time_list: + for time in loop_time(frq=frq): values.append(average_intensity(spin=spin, exp_type=exp_type, frq=frq, offset=offset, point=point, time=time, sim_index=sim_index)) errors.append(average_intensity(spin=spin, exp_type=exp_type, frq=frq, offset=offset, point=point, time=time, error=True)) times.append(time) @@ -1373,7 +1373,7 @@ # Generate the data structure to return. errors = [] - for time in cdp.relax_time_list: + for time in loop_time(frq=frq): errors.append(average_intensity(spin=spin, exp_type=exp_type, frq=frq, offset=offset, point=point, time=time, error=True)) # All other models. @@ -1597,7 +1597,7 @@ # Loop over each time point. ti = 0 - for time in loop_time(): + for time in loop_time(frq=frq): # Get the intensity keys. int_keys = find_intensity_keys(exp_type=exp_type, frq=frq, point=point, time=time) Modified: trunk/specific_analyses/relax_disp/disp_data.py URL: http://svn.gna.org/viewcvs/relax/trunk/specific_analyses/relax_disp/disp_data.py?rev=22294&r1=22293&r2=22294&view=diff ============================================================================== --- trunk/specific_analyses/relax_disp/disp_data.py (original) +++ trunk/specific_analyses/relax_disp/disp_data.py Mon Feb 24 18:28:11 2014 @@ -1930,7 +1930,7 @@ set_labels.append("Spin %s" % id) # Loop over the relaxation time periods. - for time in cdp.relax_time_list: + for time in loop_time(frq=frq): # The key. keys = find_intensity_keys(exp_type=exp_type, frq=frq, point=point, time=time) _______________________________________________ relax (http://www.nmr-relax.com) This is the relax-commits mailing list relax-commits@xxxxxxx 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-commits