r8993 - /branches/relax_disp/maths_fns/dispersion.c -- March 24, 2009 - 20:45

Author: semor
Date: Tue Mar 24 20:45:38 2009
New Revision: 8993

URL: http://svn.gna.org/viewcvs/relax?rev=8993&view=rev
Log:
Implemented the equation for back-calculation of R2eff for fast-exchange CPMG 
relaxation dispersion.


Modified:
    branches/relax_disp/maths_fns/dispersion.c

Modified: branches/relax_disp/maths_fns/dispersion.c
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/maths_fns/dispersion.c?rev=8993&r1=8992&r2=8993&view=diff
==============================================================================
--- branches/relax_disp/maths_fns/dispersion.c (original)
+++ branches/relax_disp/maths_fns/dispersion.c Tue Mar 24 20:45:38 2009
@@ -36,20 +36,14 @@
        */
 
     /* Declarations */
-    double Rx, I0;
+    double R2, Rex, kex;
     int i;
 
 
     /* Loop over the time points */
-    /* for (i = 0; i < num_times; i++) { */
     for (i = 0; i < num_times; i++) {
-        /* Zero Rx value */
-        if (params[0] == 0.0)
-            back_calc[i] = 0.0;
-
         /* Back calculate */
-        else
-            back_calc[i] = params[1] * exp(-cpmg_frqs[i] * params[0]);
+        back_calc[i] = params[0] + params[1] * (1 - 2 * tanh(params[2] / (2 
* 4 * cpmg_frqs[i])) * ((4 * cpmg_frqs[i] ) / params[2]));
 
     }
 }