X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?a=blobdiff_plain;f=ginac%2Finifcns.h;h=3d170ff8df5435d332971e9731e3795c3cbaefc4;hb=f1637dceefe0fbf79308958f8bc28dedd38aab85;hp=cff2bb2484309200a3ff777665e1c5b6dbb90979;hpb=8d08b902f1bb734fe91b5aa29c66c59788464be7;p=ginac.git diff --git a/ginac/inifcns.h b/ginac/inifcns.h index cff2bb24..3d170ff8 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 @@ -26,12 +26,16 @@ #include "function.h" #include "ex.h" -#ifndef NO_NAMESPACE_GINAC namespace GiNaC { -#endif // ndef NO_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) @@ -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(const ex & p1) { - return function(function_index_zeta1, p1); + return function(function_index_zeta1, 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); + 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) +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(const ex & p1) { - return function(function_index_psi1, p1); + return function(function_index_psi1, 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); + return function(function_index_psi2, p1, p2); } - + /** Factorial function. */ DECLARE_FUNCTION_1P(factorial) @@ -128,15 +133,30 @@ DECLARE_FUNCTION_2P(Derivative) ex lsolve(const ex &eqns, const ex &symbols); -ex ncpower(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) { - return is_ex_the_function(e, Order); + return is_ex_the_function(e, Order); } -#ifndef NO_NAMESPACE_GINAC } // namespace GiNaC -#endif // ndef NO_NAMESPACE_GINAC #endif // ndef __GINAC_INIFCNS_H__