PolyBoRi
RelatedTermsBase.h
Go to the documentation of this file.
1 // -*- c++ -*-
2 //*****************************************************************************
14 //*****************************************************************************
15 
16 #ifndef polybori_groebner_RelatedTermsBase_h_
17 #define polybori_groebner_RelatedTermsBase_h_
18 
19 #include "minimal_elements.h"
21 #include "BoundedDivisorsOf.h"
22 
23 // include basic definitions
24 #include "groebner_defs.h"
25 
27 
33 
34 public:
37  m_lead(ring), m_ignorable(ring), m_ext_prod_terms(ring),
38  m_related(ring), m_divisors_of(ring, ring) { }
39 
41  RelatedTermsBase(const Monomial& lead, const MonomialSet& terms,
42  const MonomialSet& ignorable, const MonomialSet& significant):
43  m_lead(lead), m_ignorable(ignorable), m_ext_prod_terms(lead.ring()),
44  m_related(lead.ring()), m_divisors_of(lead, significant) {
45 
46  init(significant, unrelated_terms(significant));
47  }
48 
50  const MonomialSet& related_terms() const { return m_related; }
51 
53  const BoundedDivisorsOf& divisors_of() const { return m_divisors_of; }
54 
56  MonomialSet factors(const MonomialSet& terms) const {
57  return minimal_elements_divided(critical_terms_base(terms), m_lead,
58  m_ext_prod_terms);
59  }
60 
63  return mod_mon_set(m_related.intersect(terms), m_ignorable);
64  }
65 
67  MonomialSet unrelated_terms(const MonomialSet& terms) const {
68  return std::accumulate(m_lead.begin(), m_lead.end(), terms,
69  member_function_operator(&MonomialSet::subset0) );
70  }
71 
72 private:
74  void init(const MonomialSet& significant, const MonomialSet& unrelated) {
75  PBORI_ASSERT(!m_ignorable.owns(m_lead) && !significant.ownsOne());
76  m_related = significant.diff(unrelated);
77  m_ext_prod_terms = m_ignorable.existAbstract(m_lead).diff(unrelated);
78  }
79 
80  Monomial m_lead;
81  MonomialSet m_ignorable;
82  MonomialSet m_ext_prod_terms; // Contains terms of @c m_ignorable related to
83  // @c m_lead with its variables substituted by 1
84  MonomialSet m_related;
85  BoundedDivisorsOf m_divisors_of;
86 };
87 
89 
90 #endif /* polybori_groebner_RelatedTermsBase_h_ */