Correct dependancies

1 files changed, 21 insertions, 2 deletions

diff --git a/src/elementaryFunctions/log1p/dlog1ps.c b/src/elementaryFunctions/log1p/dlog1ps.c
index 741b1844..8275ea57 100644
--- a/src/elementaryFunctions/log1p/dlog1ps.c
+++ b/src/elementaryFunctions/log1p/dlog1ps.c
@@ -10,9 +10,28 @@
  *
  */
 
-#include <math.h>
 #include "log1p.h"
+#include "log.h"
 
 double		dlog1ps(double in) {
-  return (log1p(in));
+  static double A = -1.0/3.0;
+  static double B = 0.5;
+
+  if(in < -1)
+    {/* got NaN */
+      return (in - in) / (in - in); /* NaN */
+    }
+  else if(A <= in && in <= B)
+    {/* use the function log((1+g)/(1-g)) with g = x/(x + 2) */
+      double Temp = 0;
+      Temp = in / ( in + 2);
+      /*
+return lnp1m1(Temp);
+      */
+      return 0;
+    }
+  else
+    {/* use the standard formula */
+      return dlogs(in + 1);
+    }
 }