X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fsymmetry.h;h=9232180fbeb2837d4f593897fdf748a7ef9d81bb;hp=d5a3ca8d21bd77e547693e432ccc4ef81b9e34f3;hb=7769ce49235ed6510785baa0a29801e3a59b563e;hpb=86ae5e2dfe39fb17f3d6c52a497a8309c9814111 diff --git a/ginac/symmetry.h b/ginac/symmetry.h index d5a3ca8d..9232180f 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-2010 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,18 +17,18 @@ * * 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__ -#define __GINAC_SYMMETRY_H__ - -#include +#ifndef GINAC_SYMMETRY_H +#define GINAC_SYMMETRY_H #include "ex.h" +#include "archive.h" -namespace GiNaC { +#include +namespace GiNaC { class sy_is_less; class sy_swap; @@ -80,10 +80,17 @@ 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; + /** Save (a.k.a. serialize) object into archive. */ + void archive(archive_node& n) const; + /** Read (a.k.a. deserialize) object from archive. */ + void read_archive(const archive_node& n, lst& syms); 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: @@ -96,6 +103,7 @@ private: /** Vector of child nodes. */ exvector children; }; +GINAC_DECLARE_UNARCHIVER(symmetry); // global functions @@ -120,13 +128,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 +176,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__ +#endif // ndef GINAC_SYMMETRY_H