No commit message

etheory · Jul 17, 2011 · eb4ea4f · eb4ea4f
1 parent 962ffde
commit eb4ea4f
Show file tree

Hide file tree

Showing 3 changed files with 21 additions and 3 deletions.
diff --git a/fastapprox/configure.ac.local b/fastapprox/configure.ac.local
@@ -1,4 +1,6 @@
 dnl -- include additional autoconf commands here
 dnl -- do not include AC_OUTPUT, this is called for you
 
-AC_CHECK_HEADERS([emmintrin.h boost/math/special_functions/digamma.hpp])
+AC_CHECK_HEADERS([emmintrin.h boost/math/special_functions/digamma.hpp gsl/gsl_sf_lambert.h])
+
+PKG_CHECK_MODULES([GSL], [gsl])
diff --git a/fastapprox/tests/Makefile.am.local b/fastapprox/tests/Makefile.am.local
@@ -77,7 +77,7 @@ testfastlambertwcpp_SOURCES =		\
   testfastlambertwcpp.cc
 
 testfastlambertw_LDADD =		\
-  -lm
+  -lm ${GSL_LIBS}
 
 testfastlambertwcpp_LDADD =		\
   -lm

diff --git a/fastapprox/tests/testfastlambertw.c b/fastapprox/tests/testfastlambertw.c
@@ -9,6 +9,11 @@
 
 #include "testmacros.h"
 
+#include "../src/config.h"
+#ifdef HAVE_GSL_GSL_SF_LAMBERT_H 
+#include <gsl/gsl_sf_lambert.h>
+#endif
+
 static inline float
 lambertwrange (void)
 {
@@ -28,7 +33,7 @@ lambertwf (float x)
   float w = (x < 5) ? 0 : log (x) - log (log (x)) + log (log (x)) / log (x);
   unsigned int n;
 
-  for (n = 0; n < 10; ++n)
+  for (n = 0; n < 20; ++n)
     {
       w = (w * w + exp (-w) * x) / (1.0 + w);
     }
@@ -39,6 +44,10 @@ lambertwf (float x)
 test_scalar (fastlambertw, lambertwf, lambertwrange (), 1e-4f, 100000000)
 test_scalar (fasterlambertw, lambertwf, lambertwrange (), 1e-2f, 100000000)
 
+#ifdef HAVE_GSL_GSL_SF_LAMBERT_H 
+test_scalar (gsl_sf_lambert_W0, lambertwf, lambertwrange (), 1e-2f, 1000000)
+#endif
+
 test_vector (vfastlambertw, lambertwf, lambertwrange (), 1e-4f, 100000000)
 test_vector (vfasterlambertw, lambertwf, lambertwrange (), 1e-2f, 100000000)
 
@@ -69,11 +78,18 @@ main (int   argc,
   test_fasterlambertw ();
   test_vfastlambertw ();
   test_vfasterlambertw ();
+#ifdef HAVE_GSL_GSL_SF_LAMBERT_H 
+  test_gsl_sf_lambert_W0 ();
+#endif
 
   time_fastlambertw ();
   time_fasterlambertw ();
   time_vfastlambertw ();
   time_vfasterlambertw ();
 
+#ifdef HAVE_GSL_GSL_SF_LAMBERT_H 
+  time_gsl_sf_lambert_W0 ();
+#endif
+
   return 0;
 }