r13064 - /branches/inversion-recovery/maths_fns/inversion-recovery.c -- June 16, 2011 - 12:07

Author: semor
Date: Thu Jun 16 12:07:47 2011
New Revision: 13064

URL: http://svn.gna.org/viewcvs/relax?rev=13064&view=rev
Log:

Modified the inversion-recovery c code to include the inversion-recovery 
equation.


Modified:
    branches/inversion-recovery/maths_fns/inversion-recovery.c

Modified: branches/inversion-recovery/maths_fns/inversion-recovery.c
URL: 
http://svn.gna.org/viewcvs/relax/branches/inversion-recovery/maths_fns/inversion-recovery.c?rev=13064&r1=13063&r2=13064&view=diff
==============================================================================
--- branches/inversion-recovery/maths_fns/inversion-recovery.c (original)
+++ branches/inversion-recovery/maths_fns/inversion-recovery.c Thu Jun 16 
12:07:47 2011
@@ -1,5 +1,6 @@
 /*
  * Copyright (C) 2006 Edward d'Auvergne
+ * Copyright (C) 2011 Sebastien Morin
  *
  * This file is part of the program relax.
  *
@@ -27,12 +28,12 @@
 
 
 
-void exponential(double *params, double *relax_times, double *back_calc, int 
num_times) {
-       /* Function to back calculate the peak intensities.
-       */
+void inversion-recovery(double *params, double *relax_times, double 
*back_calc, int num_times) {
+        /* Function to back calculate the peak intensities.
+        */
 
     /* Declarations */
-    double Rx, I0;
+    double Rx, I0, Iinf;
     int i;
 
 
@@ -45,7 +46,7 @@
 
         /* Back calculate */
         else
-            back_calc[i] = params[1] * exp(-relax_times[i] * params[0]);
+            back_calc[i] = params[2] - (params[2] - params[1]) * 
exp(-relax_times[i] * params[0]);
 
     }
 }