iterative.cpp

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
 
 * -------------------------------------------------------------------------- */
 
#include <gtsam/linear/iterative-inl.h>
#include <gtsam/base/Vector.h>
#include <gtsam/base/Matrix.h>
#include <gtsam/linear/GaussianFactorGraph.h>
#include <gtsam/linear/IterativeSolver.h>
 
#include <iostream>
 
using namespace std;
 
namespace gtsam {
 
  /* ************************************************************************* */
  void System::print(const string& s) const {
    cout << s << endl;
    gtsam::print(A_, "A");
    gtsam::print(b_, "b");
  }
 
  /* ************************************************************************* */
 
  Vector steepestDescent(const System& Ab, const Vector& x,
      const ConjugateGradientParameters & parameters) {
    return conjugateGradients<System, Vector, Vector>(Ab, x, parameters, true);
  }
 
  Vector conjugateGradientDescent(const System& Ab, const Vector& x,
      const ConjugateGradientParameters & parameters) {
    return conjugateGradients<System, Vector, Vector>(Ab, x, parameters);
  }
 
  /* ************************************************************************* */
  Vector steepestDescent(const Matrix& A, const Vector& b, const Vector& x,
      const ConjugateGradientParameters & parameters) {
    System Ab(A, b);
    return conjugateGradients<System, Vector, Vector>(Ab, x, parameters, true);
  }
 
  Vector conjugateGradientDescent(const Matrix& A, const Vector& b,
      const Vector& x, const ConjugateGradientParameters & parameters) {
    System Ab(A, b);
    return conjugateGradients<System, Vector, Vector>(Ab, x, parameters);
  }
 
  /* ************************************************************************* */
  VectorValues steepestDescent(const GaussianFactorGraph& fg,
      const VectorValues& x, const ConjugateGradientParameters & parameters) {
    return conjugateGradients<GaussianFactorGraph, VectorValues, Errors>(
        fg, x, parameters, true);
  }
 
  VectorValues conjugateGradientDescent(const GaussianFactorGraph& fg,
      const VectorValues& x, const ConjugateGradientParameters & parameters) {
    return conjugateGradients<GaussianFactorGraph, VectorValues, Errors>(
        fg, x, parameters);
  }
 
/* ************************************************************************* */
 
} // namespace gtsam