From d42631ab885d479c03a62a7a797701a3f40801a9 Mon Sep 17 00:00:00 2001
From: Erich Eckner <git@eckner.net>
Date: Fri, 8 Apr 2016 15:53:47 +0200
Subject: perl-pdl neu

---
 perl-pdl/gsl_sf_ellint.pd.diff | 38 ++++++++++++++++++++++++++++++++++++++
 1 file changed, 38 insertions(+)
 create mode 100644 perl-pdl/gsl_sf_ellint.pd.diff

(limited to 'perl-pdl/gsl_sf_ellint.pd.diff')

diff --git a/perl-pdl/gsl_sf_ellint.pd.diff b/perl-pdl/gsl_sf_ellint.pd.diff
new file mode 100644
index 00000000..9203d9df
--- /dev/null
+++ b/perl-pdl/gsl_sf_ellint.pd.diff
@@ -0,0 +1,38 @@
+--- gsl_sf_ellint.pd_orig	2015-11-15 10:46:20 +1100
++++ gsl_sf_ellint.pd	2015-11-15 10:54:47 +1100
+@@ -82,6 +82,10 @@
+        Doc =>'Legendre form of incomplete elliptic integrals P(phi,k,n) = Integral[(1 + n Sin[t]^2)^(-1)/Sqrt[1 - k^2 Sin[t]^2], {t, 0, phi}]'
+       );
+ 
++my $v = `gsl-config --version`;
++
++if($v < 2.0) {
++
+ pp_def('gsl_sf_ellint_D',
+        GenericTypes => [D],
+        Pars=>'double phi(); double k(); double n();
+@@ -95,6 +99,24 @@
+        Doc =>'Legendre form of incomplete elliptic integrals D(phi,k,n)'
+       );
+ 
++}
++else {
++
++pp_def('gsl_sf_ellint_D',
++       GenericTypes => [D],
++       Pars=>'double phi(); double k();
++              double [o]y(); double [o]e()',
++       Code =>'
++gsl_sf_result r;
++GSLERR(gsl_sf_ellint_D_e,($phi(),$k(),GSL_PREC_DOUBLE,&r))
++$y() = r.val;
++$e() = r.err; 
++',
++       Doc =>'Legendre form of incomplete elliptic integrals D(phi,k)'
++      );
++
++}
++
+ pp_def('gsl_sf_ellint_RC',
+        GenericTypes => [D],
+        Pars=>'double x(); double yy(); double [o]y(); double [o]e()',
-- 
cgit v1.2.3-54-g00ecf