* Re: Complex Digamma function in gsl
[not found] <1148479001.4474661950e5c@webmail.physics.unige.ch>
@ 2006-05-26 8:29 ` Gerard Jungman
2006-05-30 10:55 ` Brian Gough
0 siblings, 1 reply; 5+ messages in thread
From: Gerard Jungman @ 2006-05-26 8:29 UTC (permalink / raw)
To: Giorgio.Levy; +Cc: gsl-discuss
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
^ permalink raw reply [flat|nested] 5+ messages in thread
* Re: Complex Digamma function in gsl
2006-05-26 8:29 ` Complex Digamma function in gsl Gerard Jungman
@ 2006-05-30 10:55 ` Brian Gough
2006-05-31 11:44 ` Gerard Jungman
0 siblings, 1 reply; 5+ messages in thread
From: Brian Gough @ 2006-05-30 10:55 UTC (permalink / raw)
To: Gerard Jungman; +Cc: Giorgio.Levy, gsl-discuss
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/
^ permalink raw reply [flat|nested] 5+ messages in thread
* Re: Complex Digamma function in gsl
2006-05-30 10:55 ` Brian Gough
@ 2006-05-31 11:44 ` Gerard Jungman
2006-06-01 0:07 ` Brian Gough
0 siblings, 1 reply; 5+ messages in thread
From: Gerard Jungman @ 2006-05-31 11:44 UTC (permalink / raw)
To: Brian Gough; +Cc: Giorgio.Levy, gsl-discuss
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
^ permalink raw reply [flat|nested] 5+ messages in thread
* Re: Complex Digamma function in gsl
2006-05-31 11:44 ` Gerard Jungman
@ 2006-06-01 0:07 ` Brian Gough
2006-06-01 12:25 ` Gerard Jungman
0 siblings, 1 reply; 5+ messages in thread
From: Brian Gough @ 2006-06-01 0:07 UTC (permalink / raw)
To: Gerard Jungman; +Cc: gsl-discuss
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,
^ permalink raw reply [flat|nested] 5+ messages in thread
* Re: Complex Digamma function in gsl
2006-06-01 0:07 ` Brian Gough
@ 2006-06-01 12:25 ` Gerard Jungman
0 siblings, 0 replies; 5+ messages in thread
From: Gerard Jungman @ 2006-06-01 12:25 UTC (permalink / raw)
To: Brian Gough; +Cc: gsl-discuss
Ok, psi(z) is checked into CVS and tests are in place.
Seems to work.
--
Gerard Jungman <jungman@lanl.gov>
Los Alamos National Laboratory
^ permalink raw reply [flat|nested] 5+ messages in thread
end of thread, other threads:[~2006-06-01 0:07 UTC | newest]
Thread overview: 5+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
[not found] <1148479001.4474661950e5c@webmail.physics.unige.ch>
2006-05-26 8:29 ` Complex Digamma function in gsl Gerard Jungman
2006-05-30 10:55 ` Brian Gough
2006-05-31 11:44 ` Gerard Jungman
2006-06-01 0:07 ` Brian Gough
2006-06-01 12:25 ` Gerard Jungman
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).