Re: Complex Digamma function in gsl

Re: Complex Digamma function in gsl

[Gerard Jungman]

On Wed, 2006-05-24 at 15:56 +0200, Giorgio.Levy@physics.unige.ch wrote:
> Hi, I am using gsl for simulate the Density of states of superconductors. I
> would like to know if there is a complex version of the digamma function in gsl
> (psi).

Hi. It's best to post directly to gsl-discuss@sources.redhat.com.
Somebody might be able to help you. I am cc'ing the list on this.

Currently there are only real-argument versions of the
various poly-gamma functions, although there is a
complex version of the gamma function.

-- 
Gerard Jungman <jungman@lanl.gov>
Los Alamos National Laboratory

Re: Complex Digamma function in gsl

[Brian Gough]

Gerard Jungman writes:
 > Hi. It's best to post directly to gsl-discuss@sources.redhat.com.
 > Somebody might be able to help you. I am cc'ing the list on this.
 > 
 > Currently there are only real-argument versions of the
 > various poly-gamma functions, although there is a
 > complex version of the gamma function.

Jerry, what's your opinion on the CPSI routine in SLATEC?  It looks
fairly simple.

-- 
Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/

Re: Complex Digamma function in gsl

[Gerard Jungman]

On Sat, 2006-05-27 at 14:42 +0100, Brian Gough wrote:
> 
> Jerry, what's your opinion on the CPSI routine in SLATEC?  It looks
> fairly simple.

I looked at it; it's simple enough, once you get past the
fortran kookiness. One reason I never did a complex version
was that I wanted to also have complex versions for the
general polygamma functions (not just the first), and I
was too tired to do all of that.

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.

-- 
Gerard Jungman <jungman@lanl.gov>
Los Alamos National Laboratory

Re: Complex Digamma function in gsl

[Brian Gough]

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,

Re: Complex Digamma function in gsl

[Gerard Jungman]

Ok, psi(z) is checked into CVS and tests are in place.
Seems to work.

-- 
Gerard Jungman <jungman@lanl.gov>
Los Alamos National Laboratory