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/GraphvizFormatting.h>
#include <gtsam/inference/FactorGraph.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 std::function<bool(const Factor* /*factor*/, double /*whitenedError*/, size_t /*index*/)>&
        printCondition = [](const Factor *,double, size_t) {return true;}) const;

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


    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;

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

    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;

#ifdef 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