X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Finifcns_gamma.cpp;h=f1a2455a10bef1d7c979105336bdecab79eec7ce;hp=c9adc19967ad37a4da71f91b1929103a0cd58586;hb=619d77d2676f7f1a562fb9fefc0ba6754fe2d750;hpb=d67dadd063fbae8e9a64560d2ea97c7af0248203

diff --git a/ginac/inifcns_gamma.cpp b/ginac/inifcns_gamma.cpp
index c9adc199..f1a2455a 100644
--- a/ginac/inifcns_gamma.cpp
+++ b/ginac/inifcns_gamma.cpp
@@ -4,7 +4,7 @@
  *  some related stuff. */
 
 /*
- *  GiNaC Copyright (C) 1999-2005 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2007 Johannes Gutenberg University Mainz, Germany
  *
  *  This program is free software; you can redistribute it and/or modify
  *  it under the terms of the GNU General Public License as published by
@@ -54,8 +54,7 @@ static ex lgamma_evalf(const ex & x)
 
 
 /** Evaluation of lgamma(x), the natural logarithm of the Gamma function.
- *  Knows about integer arguments and that's it.  Somebody ought to provide
- *  some good numerical evaluation some day...
+ *  Handles integer arguments as a special case.
  *
  *  @exception GiNaC::pole_error("lgamma_eval(): logarithmic pole",0) */
 static ex lgamma_eval(const ex & x)
@@ -69,7 +68,8 @@ static ex lgamma_eval(const ex & x)
 			else
 				throw (pole_error("lgamma_eval(): logarithmic pole",0));
 		}
-		//  lgamma_evalf should be called here once it becomes available
+		if (!ex_to<numeric>(x).is_rational())
+			return lgamma(ex_to<numeric>(x));
 	}
 	
 	return lgamma(x).hold();
@@ -165,7 +165,8 @@ static ex tgamma_eval(const ex & x)
 				return (pow(*_num_2_p, n).div(doublefactorial(n.mul(*_num2_p).sub(*_num1_p))))*sqrt(Pi);
 			}
 		}
-		//  tgamma_evalf should be called here once it becomes available
+		if (!ex_to<numeric>(x).is_rational())
+			return tgamma(ex_to<numeric>(x));
 	}
 	
 	return tgamma(x).hold();
@@ -262,7 +263,8 @@ static ex beta_eval(const ex & x, const ex & y)
 		    (nx+ny).is_integer() &&
 		   !(nx+ny).is_positive())
 			 return _ex0;
-		// beta_evalf should be called here once it becomes available
+		if (!ex_to<numeric>(x).is_rational() || !ex_to<numeric>(x).is_rational())
+			return evalf(beta(x, y).hold());
 	}
 	
 	return beta(x,y).hold();