Hi,
we were thinking about next necessary changes in the CST branch.
Now
we need to break through the problem, which implies the fact we split
the relaxation rate calculation into contributions of individual
interactions. Each of them has its own data class to store its
parameters (so far called data[i][j], where the [i] was a spin index
and [j] was the interaction index).
It seems to us, that it the best way to deal with it is to edit
functions: func_mf, func_local_tm , func_diff, func_all and their
equivalents for a first and second derivatives (defined in mf.py file).
Within
these functions the calculations of direction cosines, diffusion tensor
weight calculations, components of the spectral densities and so on are
performed. All these must be calculated for each interaction
separately, because each interaction has its own data storage (which
replaced previously used one data class container for the whole IS spin
system).
Instead of:
# Direction cosine calculations.
if self.diff_data.calc_di:
self.diff_data.calc_di(data, self.diff_data)
# Diffusion tensor weight calculations.
self.diff_data.calc_ci(data, self.diff_data)
# Diffusion tensor correlation times.
self.diff_data.calc_ti(data, self.diff_data)
# Calculate the components of the spectral densities.
if data.calc_jw_comps:
data.calc_jw_comps(data, params)
# Calculate the R1, R2, and sigma_noe values.
data.ri_prime = data.create_ri_prime(data)
we would suggest to introduce a loop over interacations
for j in xrange(self.num_interactions[i]):
# Direction cosine calculations.
if self.diff_data.calc_di:
self.diff_data.calc_di(data[j], self.diff_data)
# Diffusion tensor weight calculations.
self.diff_data.calc_ci(data[j], self.diff_data)
# Diffusion tensor correlation times.
self.diff_data.calc_ti(data[j], self.diff_data)
# Calculate the components of the spectral densities.
if data.calc_jw_comps:
data.calc_jw_comps(data[j], params)
# Calculate the R1, R2, and sigma_noe components.
data.ri_prime = data.create_ri_prime(data[j])
It
must be accompanied in the next step by a change the ri_prime function
so that it just calculate only a product of specific interaction
constant and corresponding linear combination of spectral densities.
While the final sumation over all interactions should be done in a
separate step.
Moreover it will be also necessary to distinguish within the
function setup_equation the type of equation used for contribution of
individual interactions according to their type.
Best
Pavel
On 10 September 2010 15:43, Edward d'Auvergne
<edward@xxxxxxxxxxxxx> wrote:
Hi Pavel,
I missed it in the patches, but there were tab characters '\t' causing
problems. These are now fixed. relax requires that a tab is replaced
by 4 spaces. I have also added you to the copyright notices
(http://svn.gna.org/viewcvs/relax?view=rev&revision=11543) of the
files you have modified.
Regards,
Edward
On 10 September 2010 15:36, Edward d'Auvergne <
edward@xxxxxxxxxxxxx> wrote:
> Hi,
>
> I've carefully checked the patches and committed them with the
> messages you wrote. Sorry again for the delays. It should be faster
> now that I am no longer in the tropical wilderness of far north
> Australia.
>
> Regards,
>
> Edward
>
>
> On 6 September 2010 13:23, Edward d'Auvergne <
edward@xxxxxxxxxxxxx> wrote:
>> Hi Pavel,
>>
>> Sorry for the delayed response. I was at the ICMRBS conference in
>> Australia and then travelled through the tropical north end of
>> Australia afterwards. I came back yesterday out of the remote
>> wilderness and can soon start looking at this patches.
>>
>> Regards,
>>
>> Edward
>>
>>
>> On 31 August 2010 18:00, Pavel Kaderavek <
pavel.kaderavek@xxxxxxxxx> wrote:
>>> Hi,
>>>
>>> some time ago, we submitted two patches regarding CST branch. We are not
>>> sure if we should wait for some additional comment from your side, or we can
>>> continue with introducing further changes in the code.
>>> Next step would be a splitting of the relaxation equation so that
>>> contribution to the relaxation due to the individual types of interaction
>>> (dipole-dipole, CSA) can be calculated separately.
>>>
>>>
>>> Regars,
>>>
>>> Pavel, Petr
>>>
>>
>
Re: CST branch