X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?a=blobdiff_plain;ds=sidebyside;f=ginac%2Fsymmetry.h;h=5cd62c713ae295671dd158b115f77ecc7bc9daa7;hb=7d870583a6bf21a2ffb7b6f051b702064623892e;hp=d5a3ca8d21bd77e547693e432ccc4ef81b9e34f3;hpb=86ae5e2dfe39fb17f3d6c52a497a8309c9814111;p=ginac.git diff --git a/ginac/symmetry.h b/ginac/symmetry.h index d5a3ca8d..5cd62c71 100644 --- a/ginac/symmetry.h +++ b/ginac/symmetry.h @@ -3,7 +3,7 @@ * Interface to GiNaC's symmetry definitions. */ /* - * GiNaC Copyright (C) 1999-2003 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2008 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 @@ -17,7 +17,7 @@ * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ #ifndef __GINAC_SYMMETRY_H__ @@ -80,10 +80,13 @@ public: /** Check whether this node actually represents any kind of symmetry. */ bool has_symmetry() const {return type != none || !children.empty(); } + /** Check whether this node involves a cyclic symmetry. */ + bool has_cyclic() const; protected: void do_print(const print_context & c, unsigned level) const; void do_print_tree(const print_tree & c, unsigned level) const; + unsigned calchash() const; // member variables private: @@ -120,13 +123,23 @@ inline symmetry sy_cycl(const symmetry &c1, const symmetry &c2) { return symmetr inline symmetry sy_cycl(const symmetry &c1, const symmetry &c2, const symmetry &c3) { return symmetry(symmetry::cyclic, c1, c2).add(c3); } inline symmetry sy_cycl(const symmetry &c1, const symmetry &c2, const symmetry &c3, const symmetry &c4) { return symmetry(symmetry::cyclic, c1, c2).add(c3).add(c4); } +// These return references to preallocated common symmetries (similar to +// the numeric flyweights). +const symmetry & not_symmetric(); +const symmetry & symmetric2(); +const symmetry & symmetric3(); +const symmetry & symmetric4(); +const symmetry & antisymmetric2(); +const symmetry & antisymmetric3(); +const symmetry & antisymmetric4(); + /** Canonicalize the order of elements of an expression vector, according to * the symmetry properties defined in a symmetry tree. * * @param v Start of expression vector * @param symm Root node of symmetry tree * @return the overall sign introduced by the reordering (+1, -1 or 0) - * or INT_MAX if nothing changed */ + * or numeric_limits::max() if nothing changed */ extern int canonicalize(exvector::iterator v, const symmetry &symm); /** Symmetrize expression over a set of objects (symbols, indices). */ @@ -158,14 +171,6 @@ inline ex symmetrize_cyclic(const ex & e, const exvector & v) return symmetrize(e, v.begin(), v.end()); } -// utility functions - -/** Specialization of is_exactly_a(obj) for symmetry objects. */ -template<> inline bool is_exactly_a(const basic & obj) -{ - return obj.tinfo()==TINFO_symmetry; -} - } // namespace GiNaC #endif // ndef __GINAC_SYMMETRY_H__