X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?a=blobdiff_plain;f=ginac%2Finifcns.h;h=11850c34896a563d935b52807e1a4a372ab4fa07;hb=69f078195a11c1a9610a6577ce2b9a2ce81e8912;hp=314228f8033972bbeee3944f6fbc84064129361a;hpb=c5ca06e3a25226684028dec4bd8cba0833998be6;p=ginac.git diff --git a/ginac/inifcns.h b/ginac/inifcns.h index 314228f8..11850c34 100644 --- a/ginac/inifcns.h +++ b/ginac/inifcns.h @@ -3,7 +3,7 @@ * Interface to GiNaC's initially known functions. */ /* - * GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2003 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 @@ -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,39 +132,7 @@ 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); - -/** 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()); -} - -/** Symmetrize expression by cyclic permuation over a set of objects - * (symbols, indices). */ -ex symmetrize_cyclic(const ex & e, exvector::const_iterator first, exvector::const_iterator last); - -/** Symmetrize expression by cyclic permutation over a set of objects - * (symbols, indices). */ -inline ex symmetrize_cyclic(const ex & e, const exvector & v) -{ - return symmetrize(e, v.begin(), v.end()); -} +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)