NonlinearFactorGraph.h

Go to the documentation of this file.

/* ----------------------------------------------------------------------------
 
 * GTSAM Copyright 2010, Georgia Tech Research Corporation,
 * Atlanta, Georgia 30332-0415
 * All Rights Reserved
 * Authors: Frank Dellaert, et al. (see THANKS for the full author list)
 
 * See LICENSE for the license information
 
 * -------------------------------------------------------------------------- */

 
// \callgraph
 
#pragma once
 
#include <gtsam/geometry/Point2.h>
#include <gtsam/nonlinear/NonlinearFactor.h>
#include <gtsam/nonlinear/NoiseModelFactorN.h>
#include <gtsam/nonlinear/GraphvizFormatting.h>
#include <gtsam/inference/FactorGraph.h>
#include <gtsam/inference/FactorErrorPredicate.h>
#include <gtsam/nonlinear/PriorFactor.h>
 
#include <memory>
#include <functional>
 
namespace gtsam {
 
  // Forward declarations
  class Values;
  class Ordering;
  class GaussianFactorGraph;
  class SymbolicFactorGraph;
  template<typename T>
  class Expression;
  template<typename T>
  class ExpressionFactor;

  class GTSAM_EXPORT NonlinearFactorGraph: public FactorGraph<NonlinearFactor> {
 
  public:
 
    typedef FactorGraph<NonlinearFactor> Base;
    typedef NonlinearFactorGraph This;
    typedef std::shared_ptr<This> shared_ptr;


    NonlinearFactorGraph() {}

    template<typename ITERATOR>
    NonlinearFactorGraph(ITERATOR firstFactor, ITERATOR lastFactor) : Base(firstFactor, lastFactor) {}

    template<class CONTAINER>
    explicit NonlinearFactorGraph(const CONTAINER& factors) : Base(factors) {}

    template<class DERIVEDFACTOR>
    NonlinearFactorGraph(const FactorGraph<DERIVEDFACTOR>& graph) : Base(graph) {}


    void print(
        const std::string& str = "NonlinearFactorGraph: ",
        const KeyFormatter& keyFormatter = DefaultKeyFormatter) const override;

    void printErrors(const Values& values, const std::string& str = "NonlinearFactorGraph: ",
      const KeyFormatter& keyFormatter = DefaultKeyFormatter,
      const FactorErrorPredicate&
        printCondition = FactorErrorPredicate{
            [](const Factor*, double, size_t) { return true; }}) const;

    bool equals(const NonlinearFactorGraph& other, double tol = 1e-9) const;


    virtual double error(const Values& values) const;

    double probPrime(const Values& values) const;

    std::shared_ptr<SymbolicFactorGraph> symbolic() const;

    Ordering orderingCOLAMD() const;

    Ordering orderingCOLAMDConstrained(const FastMap<Key, int>& constraints) const;

    virtual std::shared_ptr<GaussianFactorGraph> linearize(const Values& linearizationPoint) const;

    virtual std::shared_ptr<const NonlinearFactorGraph> cloneShared() const {
      return std::make_shared<NonlinearFactorGraph>(*this);
    }

    typedef std::function<void(const std::shared_ptr<HessianFactor>& hessianFactor)> Dampen;

    std::shared_ptr<HessianFactor> linearizeToHessianFactor(
        const Values& values, const Dampen& dampen = nullptr) const;

    std::shared_ptr<HessianFactor> linearizeToHessianFactor(
        const Values& values, const Ordering& ordering, const Dampen& dampen = nullptr) const;

    Values updateCholesky(const Values& values,
                          const Dampen& dampen = nullptr) const;

    Values updateCholesky(const Values& values, const Ordering& ordering,
                          const Dampen& dampen = nullptr) const;

    NonlinearFactorGraph clone() const;

    NonlinearFactorGraph rekey(const std::map<Key,Key>& rekey_mapping) const;

    template<typename T>
    void addExpressionFactor(const SharedNoiseModel& R, const T& z,
                             const Expression<T>& h) {
      this->emplace_shared<ExpressionFactor<T>>(R, z, h);
    }

    template<typename T>
    void addPrior(Key key, const T& prior,
                  const SharedNoiseModel& model = nullptr) {
      emplace_shared<PriorFactor<T>>(key, prior, model);
    }

    template<typename T>
    void addPrior(Key key, const T& prior, const Matrix& covariance) {
      emplace_shared<PriorFactor<T>>(key, prior, covariance);
    }

 
    using FactorGraph::dot;
    using FactorGraph::saveGraph;

    void dot(std::ostream& os, const Values& values,
             const KeyFormatter& keyFormatter = DefaultKeyFormatter,
             const GraphvizFormatting& writer = GraphvizFormatting()) const;

    std::string dot(
        const Values& values,
        const KeyFormatter& keyFormatter = DefaultKeyFormatter,
        const GraphvizFormatting& writer = GraphvizFormatting()) const;

    void saveGraph(
        const std::string& filename, const Values& values,
        const KeyFormatter& keyFormatter = DefaultKeyFormatter,
        const GraphvizFormatting& writer = GraphvizFormatting()) const;
 
   private:

    std::shared_ptr<HessianFactor> linearizeToHessianFactor(
        const Values& values, const Scatter& scatter, const Dampen& dampen = nullptr) const;
 
#if GTSAM_ENABLE_BOOST_SERIALIZATION
    friend class boost::serialization::access;
    template<class ARCHIVE>
    void serialize(ARCHIVE & ar, const unsigned int /*version*/) {
      ar & boost::serialization::make_nvp("NonlinearFactorGraph",
                boost::serialization::base_object<Base>(*this));
    }
#endif
  };

template<>
struct traits<NonlinearFactorGraph> : public Testable<NonlinearFactorGraph> {
};
 
} //\ namespace gtsam