X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?a=blobdiff_plain;ds=sidebyside;f=ginac%2Finifcns.h;h=0defba9fbca56f41354f7c0abec922d4f175dea7;hb=9546e7613b679b28659b0774d2398ad92db26b61;hp=0ccde8de774a9d2900581d905908fc17b82325cc;hpb=8d16780dbd9f4da9397e638aca213745589818c0;p=ginac.git diff --git a/ginac/inifcns.h b/ginac/inifcns.h index 0ccde8de..0defba9f 100644 --- a/ginac/inifcns.h +++ b/ginac/inifcns.h @@ -3,7 +3,7 @@ * Interface to GiNaC's initially known functions. */ /* - * GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2001 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 @@ -23,11 +23,20 @@ #ifndef __GINAC_INIFCNS_H__ #define __GINAC_INIFCNS_H__ -#include -#include +#include "function.h" +#include "ex.h" namespace GiNaC { +/** Absolute value. */ +DECLARE_FUNCTION_1P(abs) + +/** Complex sign. */ +DECLARE_FUNCTION_1P(csgn) + +/** Eta function: log(a*b) == log(a) + log(b) + eta(a, b). */ +DECLARE_FUNCTION_2P(eta) + /** Sine. */ DECLARE_FUNCTION_1P(sin) @@ -79,12 +88,41 @@ DECLARE_FUNCTION_1P(Li2) /** Trilogarithm. */ DECLARE_FUNCTION_1P(Li3) -/** Gamma function. */ -DECLARE_FUNCTION_1P(gamma) +// overloading at work: we cannot use the macros here +/** Riemann's Zeta-function. */ +extern const unsigned function_index_zeta1; +template +inline function zeta(const T1 & p1) { + return function(function_index_zeta1, ex(p1)); +} +/** Derivatives of Riemann's Zeta-function. */ +extern const unsigned function_index_zeta2; +template +inline function zeta(const T1 & p1, const T2 & p2) { + return function(function_index_zeta2, ex(p1), ex(p2)); +} + +/** Gamma-function. */ +DECLARE_FUNCTION_1P(lgamma) +DECLARE_FUNCTION_1P(tgamma) -/** Psi function (aka polygamma-function) */ -DECLARE_FUNCTION_2P(psi) - +/** Beta-function. */ +DECLARE_FUNCTION_2P(beta) + +// overloading at work: we cannot use the macros here +/** Psi-function (aka digamma-function). */ +extern const unsigned function_index_psi1; +template +inline function psi(const T1 & p1) { + return function(function_index_psi1, ex(p1)); +} +/** Derivatives of Psi-function (aka polygamma-functions). */ +extern const unsigned function_index_psi2; +template +inline function psi(const T1 & p1, const T2 & p2) { + return function(function_index_psi2, ex(p1), ex(p2)); +} + /** Factorial function. */ DECLARE_FUNCTION_1P(factorial) @@ -94,11 +132,10 @@ DECLARE_FUNCTION_2P(binomial) /** Order term function (for truncated power series). */ DECLARE_FUNCTION_1P(Order) -ex lsolve(ex const &eqns, ex const &symbols); - -ex ncpower(ex const &basis, unsigned exponent); +ex lsolve(const ex &eqns, const ex &symbols); -inline bool is_order_function(ex const & e) +/** Check whether a function is the Order (O(n)) function. */ +inline bool is_order_function(const ex & e) { return is_ex_the_function(e, Order); }