X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Finifcns.h;h=6a7efbe610a24930690b1462eb42f02427585f9a;hp=644492d852dcb42f9f2c697a86409f7b8e8a505f;hb=0117bd6ef4af029934703940d59e1c70866937b0;hpb=15d4b353c85f8815a95d97ab977c9ca48e155574 diff --git a/ginac/inifcns.h b/ginac/inifcns.h index 644492d8..6a7efbe6 100644 --- a/ginac/inifcns.h +++ b/ginac/inifcns.h @@ -3,7 +3,7 @@ * Interface to GiNaC's initially known functions. */ /* - * GiNaC Copyright (C) 1999-2000 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,15 +23,19 @@ #ifndef __GINAC_INIFCNS_H__ #define __GINAC_INIFCNS_H__ -#include -#include +#include "function.h" +#include "ex.h" -#ifndef NO_GINAC_NAMESPACE namespace GiNaC { -#endif // ndef NO_GINAC_NAMESPACE /** 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) @@ -84,36 +88,37 @@ DECLARE_FUNCTION_1P(Li2) /** Trilogarithm. */ DECLARE_FUNCTION_1P(Li3) -// overloading at work: we cannot use the macros +// overloading at work: we cannot use the macros here /** Riemann's Zeta-function. */ extern const unsigned function_index_zeta1; -inline function zeta(ex const & p1) { - return function(function_index_zeta1, p1); +inline function zeta(const ex & p1) { + return function(function_index_zeta1, p1); } /** Derivatives of Riemann's Zeta-function. */ extern const unsigned function_index_zeta2; -inline function zeta(ex const & p1, ex const & p2) { - return function(function_index_zeta2, p1, p2); +inline function zeta(const ex & p1, const ex & p2) { + return function(function_index_zeta2, p1, p2); } /** Gamma-function. */ -DECLARE_FUNCTION_1P(gamma) +DECLARE_FUNCTION_1P(lgamma) +DECLARE_FUNCTION_1P(tgamma) /** Beta-function. */ DECLARE_FUNCTION_2P(beta) -// overloading at work: we cannot use the macros +// overloading at work: we cannot use the macros here /** Psi-function (aka digamma-function). */ extern const unsigned function_index_psi1; -inline function psi(ex const & p1) { - return function(function_index_psi1, p1); +inline function psi(const ex & p1) { + return function(function_index_psi1, p1); } /** Derivatives of Psi-function (aka polygamma-functions). */ extern const unsigned function_index_psi2; -inline function psi(ex const & p1, ex const & p2) { - return function(function_index_psi2, p1, p2); +inline function psi(const ex & p1, const ex & p2) { + return function(function_index_psi2, p1, p2); } - + /** Factorial function. */ DECLARE_FUNCTION_1P(factorial) @@ -123,17 +128,43 @@ DECLARE_FUNCTION_2P(binomial) /** Order term function (for truncated power series). */ DECLARE_FUNCTION_1P(Order) -ex lsolve(ex const &eqns, ex const &symbols); +/** 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()); +} -ex ncpower(ex const &basis, unsigned exponent); +/** Antisymmetrize expression over a list of objects (symbols, indices). */ +ex antisymmetrize(const ex & e, const lst & l); -inline bool is_order_function(ex const & e) +inline bool is_order_function(const ex & e) { - return is_ex_the_function(e, Order); + return is_ex_the_function(e, Order); } -#ifndef NO_GINAC_NAMESPACE } // namespace GiNaC -#endif // ndef NO_GINAC_NAMESPACE #endif // ndef __GINAC_INIFCNS_H__