X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?a=blobdiff_plain;ds=sidebyside;f=ginac%2Finifcns.h;h=e0ddc8bbe6ef0d885cdf614a49f757b3a1771ab9;hb=24f75d316ba00dcdf32d1623ade26b0c15a01dff;hp=6a7efbe610a24930690b1462eb42f02427585f9a;hpb=0117bd6ef4af029934703940d59e1c70866937b0;p=ginac.git diff --git a/ginac/inifcns.h b/ginac/inifcns.h index 6a7efbe6..e0ddc8bb 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,38 +132,9 @@ 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()); -} - -/** Symmetrize expression over a list of objects (symbols, indices). */ -ex symmetrize(const ex & e, const lst & l); - -/** 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()); -} - -/** Antisymmetrize expression over a list of objects (symbols, indices). */ -ex antisymmetrize(const ex & e, const lst & l); +ex lsolve(const ex &eqns, const ex &symbols, unsigned options = determinant_algo::automatic); +/** 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);