mailRe: [bug #9238] Error in setting the heteronuclear frequency in 'maths_fns/jw_mapping.py'


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by Edward d'Auvergne on May 30, 2007 - 17:05:
That one definitely looks like a fatal bug!  I have no idea how I let
that one get through.  Thanks for catching that one - this makes quite
a difference for reduced spectral density mapping.  The bug is a
little more problematic though.  The multiplication of 2pi is required
to convert to the correct units:

frq = frq * 2 * pi

And the ratio of gyromagnetic ratios is required as well:

frqX = frq * self.data.gx / self.data.gh

So as you can see, the RSDM results from relax will be rubbish!  I'll
make your changes to the code, fix the test-suite, and will hopefully
soon be able to release a new relax 1.2 version.

Thanks,

Edward



On 5/30/07, anonymous <NO-REPLY.INVALID-ADDRESS@xxxxxxx> wrote:

URL:
  <http://gna.org/bugs/?9238>

                 Summary: Error in setting the heteronuclear frequency in
'maths_fns/jw_mapping.py'
                 Project: relax
            Submitted by: None
            Submitted on: Wednesday 05/30/2007 at 02:10 CEST
                Category: relax's source code
                Severity: 4 - Important
                Priority: 5 - Normal
                  Status: None
                 Privacy: Public
             Assigned to: None
         Originator Name: Sébastien Morin
        Originator Email: sebastien.morin.1@xxxxxxxxx
             Open/Closed: Open
         Discussion Lock: Any
                 Release: Repository: 1.2 line
        Operating System: GNU/Linux

    _______________________________________________________

Details:

In the file 'maths_fns/jw_mapping.py' of the 1.2 repository (revision 3296),
there is an error in setting the heteronuclear frequency.

Line 52 is : frqX = 2.0 * pi * frq

It is erroneous since 'frqX' is not about transforming 'frq' to angular
frequency, but calculating the heteronuclear frequency corresponding to the
entered proton frequency...

Thus, it should be something like :

frqX = frq * self.data.gx / self.data.gh

or, more explicitely,

frqX = frq * 0.101584088

Thanks


Sébastien  :)




    _______________________________________________________

Reply to this item at:

  <http://gna.org/bugs/?9238>

_______________________________________________
  Message sent via/by Gna!
  http://gna.org/


_______________________________________________
relax (http://nmr-relax.com)

This is the relax-devel mailing list
relax-devel@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-devel




Related Messages


Powered by MHonArc, Updated Wed May 30 17:20:19 2007