Complex Digamma function in gsl

Brian Gough bjg@network-theory.co.uk
Thu Jun 1 00:07:00 GMT 2006


Gerard Jungman writes:
 > But since there is some demand, I have now implemented
 > psi(z) for complex z. It is checked into CVS, but it
 > is not tested yet. The problem is that I cannot get
 > test_sf.c to build; it fails to find some symbols
 > from the library for complex math, declared in 
 > gsl/complex/gsl_complex_math.h.
 > 
 > Brian, can you look at this and figure out what
 > is going on? If I can get the tests to build, I
 > can wrap this up.

Here are the changes to get it to link (basically just to add
../complex/libgslcomplex.la to the list of libraries)

-- 
Brian Gough

Index: Makefile.am
===================================================================
RCS file: /cvs/gsl/gsl/specfunc/Makefile.am,v
retrieving revision 1.99
diff -u -r1.99 Makefile.am
--- Makefile.am	18 Apr 2006 17:59:46 -0000	1.99
+++ Makefile.am	31 May 2006 11:41:48 -0000
@@ -12,7 +12,7 @@
 
 check_PROGRAMS = test
 
-test_LDADD = libgslspecfunc.la ../eigen/libgsleigen.la ../linalg/libgsllinalg.la  ../sort/libgslsort.la ../matrix/libgslmatrix.la ../vector/libgslvector.la ../blas/libgslblas.la ../cblas/libgslcblas.la ../block/libgslblock.la ../ieee-utils/libgslieeeutils.la  ../err/libgslerr.la ../test/libgsltest.la ../sys/libgslsys.la ../utils/libutils.la
+test_LDADD = libgslspecfunc.la ../eigen/libgsleigen.la ../linalg/libgsllinalg.la  ../sort/libgslsort.la ../matrix/libgslmatrix.la ../vector/libgslvector.la ../blas/libgslblas.la ../cblas/libgslcblas.la ../block/libgslblock.la ../complex/libgslcomplex.la ../ieee-utils/libgslieeeutils.la  ../err/libgslerr.la ../test/libgsltest.la ../sys/libgslsys.la ../utils/libutils.la
 
 test_SOURCES = test_sf.c test_sf.h test_airy.c test_bessel.c test_coulomb.c test_dilog.c test_gamma.c test_hyperg.c test_legendre.c test_mathieu.c
   
Index: gsl_sf_psi.h
===================================================================
RCS file: /cvs/gsl/gsl/specfunc/gsl_sf_psi.h,v
retrieving revision 1.21
diff -u -r1.21 gsl_sf_psi.h
--- gsl_sf_psi.h	30 May 2006 23:11:01 -0000	1.21
+++ gsl_sf_psi.h	31 May 2006 11:41:48 -0000
@@ -73,7 +73,7 @@
  *
  * exceptions: GSL_EDOM
  */
-int gsl_sf_psi_complex_e(
+int gsl_sf_complex_psi_e(
   const double x,
   const double y,
   gsl_sf_result * result_re,
Index: psi.c
===================================================================
RCS file: /cvs/gsl/gsl/specfunc/psi.c,v
retrieving revision 1.41
diff -u -r1.41 psi.c
--- psi.c	30 May 2006 23:11:01 -0000	1.41
+++ psi.c	31 May 2006 11:41:48 -0000
@@ -784,7 +784,7 @@
 
 
 int
-gsl_sf_psi_complex_e(
+gsl_sf_complex_psi_e(
   double x,
   double y,
   gsl_sf_result * result_re,



More information about the Gsl-discuss mailing list