00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 #ifndef ORBIT_SET_H_
00034 #define ORBIT_SET_H_
00035
00036 #include <permlib/transversal/orbit.h>
00037
00038 namespace permlib {
00039
00041 template<class PERM,class DOMAIN>
00042 class OrbitSet : public Orbit<PERM,DOMAIN> {
00043 public:
00044 virtual bool contains(const DOMAIN& val) const;
00045
00047 bool empty() const { return m_orbitSet.empty(); }
00048
00050
00056 template<class Action>
00057 void orbit(const DOMAIN& beta, const PERMlist &generators, Action a);
00058
00060 ulong size() const {return m_orbitSet.size(); }
00061
00062 virtual const DOMAIN& element() const;
00063 protected:
00065 std::set<DOMAIN> m_orbitSet;
00066
00067 virtual bool foundOrbitElement(const DOMAIN& alpha, const DOMAIN& alpha_p, const PERMptr& p);
00068 };
00069
00070 template <class PERM,class DOMAIN>
00071 inline bool OrbitSet<PERM,DOMAIN>::foundOrbitElement(const DOMAIN& alpha, const DOMAIN& alpha_p, const PERMptr& p) {
00072 if (m_orbitSet.insert(alpha_p).second) {
00073 DEBUG(std::cout << " o " << alpha_p << " @ " << (*p) << std::endl;)
00074 return true;
00075 }
00076 return false;
00077 }
00078
00079 template <class PERM,class DOMAIN>
00080 inline bool OrbitSet<PERM,DOMAIN>::contains(const DOMAIN& val) const {
00081 return m_orbitSet.find(val) != m_orbitSet.end();
00082 }
00083
00084 template <class PERM,class DOMAIN>
00085 template<class Action>
00086 inline void OrbitSet<PERM,DOMAIN>::orbit(const DOMAIN& beta, const PERMlist &generators, Action a) {
00087 std::list<DOMAIN> orbitList;
00088 Orbit<PERM,DOMAIN>::orbit(beta, generators, a, orbitList);
00089 }
00090
00091 template <class PERM,class DOMAIN>
00092 inline const DOMAIN& OrbitSet<PERM,DOMAIN>::element() const {
00093 return *(m_orbitSet.begin());
00094 }
00095
00096 }
00097
00098 #endif // -- ORBIT_SET_H_