FixedLagSmoother.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_unstable/dllexport.h>
#include <gtsam/inference/Key.h>
#include <gtsam/nonlinear/NonlinearFactorGraph.h>
#include <gtsam/nonlinear/Values.h>
 
#include <map>
#include <vector>
 
namespace gtsam {
 
class GTSAM_UNSTABLE_EXPORT FixedLagSmoother {
 
public:

  typedef boost::shared_ptr<FixedLagSmoother> shared_ptr;

  typedef std::map<Key, double> KeyTimestampMap;
  typedef std::multimap<double, Key> TimestampKeyMap;

  // TODO: Think of some more things to put here
  struct Result {
    size_t iterations; 
    size_t intermediateSteps; 
    size_t nonlinearVariables; 
    size_t linearVariables; 
    double error; 
    Result() : iterations(0), intermediateSteps(0), nonlinearVariables(0), linearVariables(0), error(0) {};

    size_t getIterations() const { return iterations; }
    size_t getIntermediateSteps() const { return intermediateSteps; }
    size_t getNonlinearVariables() const { return nonlinearVariables; }
    size_t getLinearVariables() const { return linearVariables; }
    double getError() const { return error; }
    void print() const;
  };

  FixedLagSmoother(double smootherLag = 0.0) : smootherLag_(smootherLag) { }

  virtual ~FixedLagSmoother() { }

  virtual void print(
      const std::string& s = "FixedLagSmoother:\n",
      const KeyFormatter& keyFormatter = DefaultKeyFormatter) const;

  virtual bool equals(const FixedLagSmoother& rhs, double tol = 1e-9) const;

  double smootherLag() const {
    return smootherLag_;
  }

  double& smootherLag() {
    return smootherLag_;
  }

  const KeyTimestampMap& timestamps() const {
    return keyTimestampMap_;
  }

  virtual Result update(const NonlinearFactorGraph& newFactors = NonlinearFactorGraph(),
                        const Values& newTheta = Values(),
                        const KeyTimestampMap& timestamps = KeyTimestampMap(),
                        const FactorIndices& factorsToRemove = FactorIndices()) = 0;

  virtual Values calculateEstimate() const  = 0;
 
 
protected:

  double smootherLag_;

  TimestampKeyMap timestampKeyMap_;
  KeyTimestampMap keyTimestampMap_;

  void updateKeyTimestampMap(const KeyTimestampMap& newTimestamps);

  void eraseKeyTimestampMap(const KeyVector& keys);

  double getCurrentTimestamp() const;

  KeyVector findKeysBefore(double timestamp) const;

  KeyVector findKeysAfter(double timestamp) const;
 
}; // FixedLagSmoother

typedef FixedLagSmoother::KeyTimestampMap FixedLagSmootherKeyTimestampMap;
typedef FixedLagSmootherKeyTimestampMap::value_type FixedLagSmootherKeyTimestampMapValue;
typedef FixedLagSmoother::Result FixedLagSmootherResult;
 
}