Domain.h

Go to the documentation of this file.

/*
 * Domain.h
 * @brief Domain restriction constraint
 * @date Feb 13, 2012
 * @author Frank Dellaert
 */
 
#pragma once
 
#include <gtsam/discrete/DiscreteKey.h>
#include <gtsam_unstable/discrete/Constraint.h>
#include <optional>
 
namespace gtsam {
 
class GTSAM_UNSTABLE_EXPORT Domain : public Constraint {
  size_t cardinality_;       
  std::set<size_t> values_;  
 
 public:
  typedef std::shared_ptr<Domain> shared_ptr;
 
  // Constructor on Discrete Key initializes an "all-allowed" domain
  Domain(const DiscreteKey& dkey)
      : Constraint(dkey.first), cardinality_(dkey.second) {
    for (size_t v = 0; v < cardinality_; v++) values_.insert(v);
  }
 
  // Constructor on Discrete Key with single allowed value
  // Consider SingleValue constraint
  Domain(const DiscreteKey& dkey, size_t v)
      : Constraint(dkey.first), cardinality_(dkey.second) {
    values_.insert(v);
  }
 
  Key key() const { return keys_[0]; }
 
  // The associated discrete key
  DiscreteKey discreteKey() const { return DiscreteKey(key(), cardinality_); }
 
  void insert(size_t value) { values_.insert(value); }
 
  void erase(size_t value) { values_.erase(value); }
 
  size_t nrValues() const { return values_.size(); }
 
  bool isSingleton() const { return nrValues() == 1; }
 
  size_t firstValue() const { return *values_.begin(); }
 
  // print
  void print(const std::string& s = "", const KeyFormatter& formatter =
                                            DefaultKeyFormatter) const override;
 
  bool equals(const DiscreteFactor& other, double tol) const override {
    if (!dynamic_cast<const Domain*>(&other))
      return false;
    else {
      const Domain& f(static_cast<const Domain&>(other));
      return (cardinality_ == f.cardinality_) && (values_ == f.values_);
    }
  }
 
  AlgebraicDecisionTree<Key> errorTree() const override {
    throw std::runtime_error("Domain::error not implemented");
  }
 
  // Return concise string representation, mostly to debug arc consistency.
  // Converts from base 0 to base1.
  std::string base1Str() const;
 
  // Check whether domain cotains a specific value.
  bool contains(size_t value) const { return values_.count(value) > 0; }
 
  double operator()(const DiscreteValues& values) const override;
 
  DecisionTreeFactor toDecisionTreeFactor() const override;
 
  DecisionTreeFactor operator*(const DecisionTreeFactor& f) const override;
 
  /*
   * Ensure Arc-consistency by checking every possible value of domain j.
   * @param j domain to be checked
   * @param (in/out) domains all domains, but only domains->at(j) will be
   * checked.
   * @return true if domains->at(j) was changed, false otherwise.
   */
  bool ensureArcConsistency(Key j, Domains* domains) const override;
 
  std::optional<Domain> checkAllDiff(const KeyVector keys,
                                       const Domains& domains) const;
 
  Constraint::shared_ptr partiallyApply(const DiscreteValues& values) const override;
 
  Constraint::shared_ptr partiallyApply(const Domains& domains) const override;
};
 
}  // namespace gtsam