1 Patch taken from the Sage packaging system.
3 --- lcalc-1.23/include/Lgamma.h 2012-08-08 23:21:55.000000000 +0200
4 +++ lcalc-1.23/include/Lgamma.h 2014-05-18 21:15:27.786889718 +0200
6 //n=0 should just give log_GAMMA(z)... thus making log_GAMMA
7 //code obsolete. But leave log_GAMMA intact anyways.
9 -precise(ttype) log_GAMMA (ttype z,int n=0)
10 +precise(ttype) log_GAMMA (ttype z,int n)
13 precise(ttype) log_G,r,r2,y;
15 //value exp_w which holds exp(-w)
16 //computes G(z,w), so there's an extra w^(-z) factor.
17 template <class ttype>
18 -Complex inc_GAMMA (ttype z,ttype w, const char *method="temme", ttype exp_w = 0, bool recycle=false)
19 +Complex inc_GAMMA (ttype z,ttype w, const char *method, ttype exp_w, bool recycle)
26 template <class ttype>
27 -ttype cfrac_GAMMA (ttype z,ttype w, ttype exp_w=0, bool recycle=false) //computes G(z,w) via continued fraction
28 +ttype cfrac_GAMMA (ttype z,ttype w, ttype exp_w, bool recycle) //computes G(z,w) via continued fraction
35 template <class ttype>
36 -ttype asympt_GAMMA (ttype z,ttype w, ttype exp_w = 0, bool recycle=false) //computes G(z,w) via asymptotic series
37 +ttype asympt_GAMMA (ttype z,ttype w, ttype exp_w, bool recycle) //computes G(z,w) via asymptotic series
40 if(my_verbose>3) cout << "called asympt_GAMMA("<<z<<","<<w<<")"<< endl;
44 template <class ttype>
45 -ttype comp_inc_GAMMA (ttype z,ttype w,ttype exp_w = 0, bool recycle=false) //computes g(z,w)
46 +ttype comp_inc_GAMMA (ttype z,ttype w,ttype exp_w, bool recycle) //computes g(z,w)
53 template <class ttype>
54 -Complex gamma_sum(Complex s, int what_type, ttype *coeff, int N, Double g, Complex l, Double Q, Long Period, Complex delta=1, const char *method="temme")
55 +Complex gamma_sum(Complex s, int what_type, ttype *coeff, int N, Double g, Complex l, Double Q, Long Period, Complex delta, const char *method)