X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Finifcns_zeta.cpp;h=287059835180a06dea6a2763b5d77ae11074cda3;hp=0f069089b2785f6890da63113066a11ce53504f6;hb=7cf5e1f6346223de683c9e7986d098fd042d691c;hpb=4e3a4ac2bcb0837611ea31bc8fc05d84a20c33ac diff --git a/ginac/inifcns_zeta.cpp b/ginac/inifcns_zeta.cpp index 0f069089..28705983 100644 --- a/ginac/inifcns_zeta.cpp +++ b/ginac/inifcns_zeta.cpp @@ -38,35 +38,40 @@ namespace GiNaC { static ex zeta1_evalf(const ex & x) { - BEGIN_TYPECHECK - TYPECHECK(x,numeric) - END_TYPECHECK(zeta(x)) - - return zeta(ex_to_numeric(x)); + if (is_exactly_a(x)) { + try { + return zeta(ex_to(x)); + } catch (const dunno &e) { } + } + + return zeta(x).hold(); } static ex zeta1_eval(const ex & x) { if (x.info(info_flags::numeric)) { - numeric y = ex_to_numeric(x); + numeric y = ex_to(x); // trap integer arguments: if (y.is_integer()) { if (y.is_zero()) - return -_ex1_2(); - if (x.is_equal(_ex1())) + return _ex_1_2; + if (x.is_equal(_ex1)) throw(std::domain_error("zeta(1): infinity")); if (x.info(info_flags::posint)) { if (x.info(info_flags::odd)) return zeta(x).hold(); else - return abs(bernoulli(y))*pow(Pi,x)*pow(_num2(),y-_num1())/factorial(y); + return abs(bernoulli(y))*pow(Pi,x)*pow(_num2,y-_num1)/factorial(y); } else { if (x.info(info_flags::odd)) - return -bernoulli(_num1()-y)/(_num1()-y); + return -bernoulli(_num1-y)/(_num1-y); else - return _num0(); + return _num0; } } + // zeta(float) + if (x.info(info_flags::numeric) && !x.info(info_flags::crational)) + return zeta1_evalf(x); } return zeta(x).hold(); } @@ -75,7 +80,7 @@ static ex zeta1_deriv(const ex & x, unsigned deriv_param) { GINAC_ASSERT(deriv_param==0); - return zeta(_ex1(), x); + return zeta(_ex1, x); } const unsigned function_index_zeta1 = @@ -83,6 +88,7 @@ const unsigned function_index_zeta1 = eval_func(zeta1_eval). evalf_func(zeta1_evalf). derivative_func(zeta1_deriv). + latex_name("\\zeta"). overloaded(2)); ////////// @@ -116,6 +122,7 @@ const unsigned function_index_zeta2 = function::register_new(function_options("zeta"). eval_func(zeta2_eval). derivative_func(zeta2_deriv). + latex_name("\\zeta"). overloaded(2)); } // namespace GiNaC