X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Finifcns.h;h=0defba9fbca56f41354f7c0abec922d4f175dea7;hp=96c256077f92f0964c4f7fa8360f407e61e8405d;hb=c366645d84ba4cfdd5b98785a9c5e66e0aea3e62;hpb=9961dfdd7a383f09b1040e195a8817a85d945e58 diff --git a/ginac/inifcns.h b/ginac/inifcns.h index 96c25607..0defba9f 100644 --- a/ginac/inifcns.h +++ b/ginac/inifcns.h @@ -91,13 +91,15 @@ DECLARE_FUNCTION_1P(Li3) // overloading at work: we cannot use the macros here /** Riemann's Zeta-function. */ extern const unsigned function_index_zeta1; -inline function zeta(const ex & p1) { - return function(function_index_zeta1, p1); +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; -inline function zeta(const ex & p1, const ex & p2) { - return function(function_index_zeta2, p1, p2); +template +inline function zeta(const T1 & p1, const T2 & p2) { + return function(function_index_zeta2, ex(p1), ex(p2)); } /** Gamma-function. */ @@ -110,13 +112,15 @@ DECLARE_FUNCTION_2P(beta) // overloading at work: we cannot use the macros here /** Psi-function (aka digamma-function). */ extern const unsigned function_index_psi1; -inline function psi(const ex & p1) { - return function(function_index_psi1, p1); +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; -inline function psi(const ex & p1, const ex & p2) { - return function(function_index_psi2, p1, p2); +template +inline function psi(const T1 & p1, const T2 & p2) { + return function(function_index_psi2, ex(p1), ex(p2)); } /** Factorial function. */ @@ -128,32 +132,8 @@ DECLARE_FUNCTION_2P(binomial) /** Order term function (for truncated power series). */ DECLARE_FUNCTION_1P(Order) -/** Inert partial differentiation operator. */ -DECLARE_FUNCTION_2P(Derivative) - ex lsolve(const ex &eqns, const ex &symbols); -/** Power of non-commutative basis. */ -ex ncpow(const ex & basis, unsigned exponent); - -/** Symmetrize expression over a set of objects (symbols, indices). */ -ex symmetrize(const ex & e, exvector::const_iterator first, exvector::const_iterator last); - -/** Symmetrize expression over a set of objects (symbols, indices). */ -inline ex symmetrize(const ex & e, const exvector & v) -{ - return symmetrize(e, v.begin(), v.end()); -} - -/** Antisymmetrize expression over a set of objects (symbols, indices). */ -ex antisymmetrize(const ex & e, exvector::const_iterator first, exvector::const_iterator last); - -/** Antisymmetrize expression over a set of objects (symbols, indices). */ -inline ex antisymmetrize(const ex & e, const exvector & v) -{ - return antisymmetrize(e, v.begin(), v.end()); -} - /** Check whether a function is the Order (O(n)) function. */ inline bool is_order_function(const ex & e) {