Changeset 786 for trunk/source/conv_eden_ioniz.cpp

Timestamp:

01/16/07 10:31:41 (2 years ago)

Author:

peter

Message:

source/service.cpp

• Bug-fix. Fix problem in e2() with negative return values for very large tau.
  

source/thirdparty_bevington.cpp
source/conv_eden_ioniz.cpp
source/thirdparty.cpp
source/thirdparty.h

• Replace Bevington linfit routine with open source version.
  
• Bug-fix. Replace call to y0() -> bessel_y0(), y1() -> bessel_y1() in bessel_yn().
  

Files:

• trunk/source/conv_eden_ioniz.cpp (modified) (1 diff)

trunk/source/conv_eden_ioniz.cpp

@@ -263 +263 @@
                                                 if( nEval>3 && lgHistUpdate )
                                                 {
-                                                        double y_intercept , stdx , stdy ,reg_coef,
-                                                                slope_save;
+                                                        double y_intercept, stdx, stdy, slope_save;
                                                         slope_save = slope_ls;
                                                         /* enough data to do least squares */
-                                                        if( linfit(eden_assumed_history, 
-                                                                eden_error_history, 
-                                                                stdarray, 
-                                                                nEval, 
-                                                                /* this says same weight for all */
-                                                                0, 
-                                                                &y_intercept, 
-                                                                &stdy, 
-                                                                &slope_ls, 
-                                                                &stdx, 
-                                                                &reg_coef) )
+                                                        if( linfit(nEval,
+                                                                   eden_assumed_history, 
+                                                                   eden_error_history, 
+                                                                   y_intercept, 
+                                                                   stdy, 
+                                                                   slope_ls, 
+                                                                   stdx) )
                                                         {
                                                                 int i;