Changeset 197260f in sasview

Timestamp:

Mar 24, 2016 3:37:55 AM (9 years ago)

Author:

wojciech

Branches:

master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload

Children:

273a6d7

Parents:

75898d4

Message:

Libfunc cleaned up

File:

• src/sas/sascalc/calculator/c_extensions/libfunc.c (modified) (6 diffs)

src/sas/sascalc/calculator/c_extensions/libfunc.c

@@ -236 +236 @@
 
 /**
-   Implements eq 6.2.5 (small gamma) of Numerical Recipes in C, essentially
-   the incomplete gamma function multiplied by the gamma function.
-   Required for implementation of fast error function (erf)
-**/
-
-/*
+Wojtek's comment Mar 22 2016: The remaing code can mostly likely be deleated
+Keeping it in order to check if it is not breaking anything
+**/
+
+
 #define ITMAX 100
 #define EPS 3.0e-7
@@ -250 +249 @@
   float sum,del,ap;
 
-  *gln = gamln(a);
+  *gln = lgamma(a);
   if(x <= 0.0) {
     if (x < 0.0) printf("Error: x less than 0 in routine gser");
@@ -275 +274 @@
 
 }
-*/
 
 /**
@@ -281 +279 @@
    representation
 **/
-/*
 void gcf(float *gammcf, float a, float x, float *gln) {
   int i;
   float an,b,c,d,del,h;
 
-  *gln = gamln(a);
+  *gln = lgamma(a);
   b = x+1.0-a;
   c = 1.0/FPMIN;
@@ -307 +304 @@
   return;
 }
-*/
 
 /**
    Represents incomplete error function, P(a,x)
 **/
-/*
 float gammp(float a, float x) {
   float gamser,gammcf,gln;
@@ -324 +319 @@
   }
 }
-*/
 /**
     Implementation of the error function, erf(x)
 **/
-/*
 float erff(float x) {
   return x < 0.0 ? -gammp(0.5,x*x) : gammp(0.5,x*x);
 }
-*/