]> www.ginac.de Git - ginac.git/blobdiff - ginac/inifcns_gamma.cpp
- most functions were not static, but they should have been
[ginac.git] / ginac / inifcns_gamma.cpp
index 34a390ed23e23e1b412609399f926328596570cf..fd3c32efef5359475862b5f413c6b543ac97d572 100644 (file)
@@ -30,6 +30,8 @@
 #include "power.h"
 #include "symbol.h"
 
+namespace GiNaC {
+
 //////////
 // gamma function
 //////////
@@ -39,7 +41,7 @@
  *  evaluation some day...
  *
  *  @exception fail_numeric("complex_infinity") or something similar... */
-ex gamma_eval(ex const & x)
+static ex gamma_eval(ex const & x)
 {
     if ( x.info(info_flags::numeric) ) {
 
@@ -74,7 +76,7 @@ ex gamma_eval(ex const & x)
     return gamma(x).hold();
 }    
     
-ex gamma_evalf(ex const & x)
+static ex gamma_evalf(ex const & x)
 {
     BEGIN_TYPECHECK
         TYPECHECK(x,numeric)
@@ -83,16 +85,16 @@ ex gamma_evalf(ex const & x)
     return gamma(ex_to_numeric(x));
 }
 
-ex gamma_diff(ex const & x, unsigned diff_param)
+static ex gamma_diff(ex const & x, unsigned diff_param)
 {
     ASSERT(diff_param==0);
 
-    return power(x, -1);       //!!
+    return power(x, -1);       // FIXME
 }
 
-ex gamma_series(ex const & x, symbol const & s, ex const & point, int order)
+static ex gamma_series(ex const & x, symbol const & s, ex const & point, int order)
 {
-       //!! Only handle one special case for now...
+       // FIXME: Only handle one special case for now...
        if (x.is_equal(s) && point.is_zero()) {
                ex e = 1 / s - EulerGamma + s * (power(Pi, 2) / 12 + power(EulerGamma, 2) / 2) + Order(power(s, 2));
                return e.series(s, point, order);
@@ -101,3 +103,5 @@ ex gamma_series(ex const & x, symbol const & s, ex const & point, int order)
 }
 
 REGISTER_FUNCTION(gamma, gamma_eval, gamma_evalf, gamma_diff, gamma_series);
+
+} // namespace GiNaC