2010-10-14 12:54:38 +08:00
|
|
|
/* ----------------------------------------------------------------------------
|
|
|
|
|
|
|
|
* 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
|
|
|
|
|
|
|
|
* -------------------------------------------------------------------------- */
|
|
|
|
|
2011-10-14 11:23:14 +08:00
|
|
|
/**
|
|
|
|
* @file testGraph.cpp
|
|
|
|
* @date Jan 12, 2010
|
|
|
|
* @author nikai
|
|
|
|
* @brief unit test for graph-inl.h
|
2010-01-13 00:12:25 +08:00
|
|
|
*/
|
|
|
|
|
|
|
|
#include <iostream>
|
|
|
|
#include <boost/shared_ptr.hpp>
|
2010-01-19 03:11:22 +08:00
|
|
|
#include <boost/assign/std/list.hpp> // for operator +=
|
|
|
|
using namespace boost::assign;
|
2010-01-13 00:12:25 +08:00
|
|
|
|
2010-10-26 04:10:33 +08:00
|
|
|
#include <CppUnitLite/TestHarness.h>
|
2010-01-13 00:12:25 +08:00
|
|
|
|
2010-01-18 13:38:53 +08:00
|
|
|
// TODO: DANGEROUS, create shared pointers
|
|
|
|
#define GTSAM_MAGIC_GAUSSIAN 3
|
|
|
|
|
2010-08-20 01:23:19 +08:00
|
|
|
#include <gtsam/slam/pose2SLAM.h>
|
|
|
|
#include <gtsam/inference/graph-inl.h>
|
2010-01-13 00:12:25 +08:00
|
|
|
|
|
|
|
using namespace std;
|
|
|
|
using namespace gtsam;
|
|
|
|
|
2012-02-21 08:53:35 +08:00
|
|
|
Key kx(size_t i) { return Symbol('x',i); }
|
|
|
|
|
2010-01-19 03:11:22 +08:00
|
|
|
/* ************************************************************************* */
|
|
|
|
// x1 -> x2
|
|
|
|
// -> x3 -> x4
|
|
|
|
// -> x5
|
|
|
|
TEST ( Ordering, predecessorMap2Keys ) {
|
2012-02-21 08:53:35 +08:00
|
|
|
PredecessorMap<Key> p_map;
|
|
|
|
p_map.insert(kx(1),kx(1));
|
|
|
|
p_map.insert(kx(2),kx(1));
|
|
|
|
p_map.insert(kx(3),kx(1));
|
|
|
|
p_map.insert(kx(4),kx(3));
|
|
|
|
p_map.insert(kx(5),kx(1));
|
2012-02-07 12:58:11 +08:00
|
|
|
|
2012-02-21 08:53:35 +08:00
|
|
|
list<Key> expected;
|
|
|
|
expected += kx(4),kx(5),kx(3),kx(2),kx(1);//PoseKey(4), PoseKey(5), PoseKey(3), PoseKey(2), PoseKey(1);
|
2012-02-07 12:58:11 +08:00
|
|
|
|
2012-02-21 08:53:35 +08:00
|
|
|
list<Key> actual = predecessorMap2Keys<Key>(p_map);
|
2010-01-19 03:11:22 +08:00
|
|
|
LONGS_EQUAL(expected.size(), actual.size());
|
|
|
|
|
2012-02-21 08:53:35 +08:00
|
|
|
list<Key>::const_iterator it1 = expected.begin();
|
|
|
|
list<Key>::const_iterator it2 = actual.begin();
|
2010-01-19 03:11:22 +08:00
|
|
|
for(; it1!=expected.end(); it1++, it2++)
|
|
|
|
CHECK(*it1 == *it2)
|
|
|
|
}
|
|
|
|
|
2010-01-13 00:12:25 +08:00
|
|
|
/* ************************************************************************* */
|
2010-01-14 11:21:07 +08:00
|
|
|
TEST( Graph, predecessorMap2Graph )
|
|
|
|
{
|
|
|
|
typedef SGraph<string>::Vertex SVertex;
|
2012-02-21 08:53:35 +08:00
|
|
|
SGraph<Key> graph;
|
2010-01-14 11:21:07 +08:00
|
|
|
SVertex root;
|
2012-02-21 08:53:35 +08:00
|
|
|
map<Key, SVertex> key2vertex;
|
2010-01-14 11:21:07 +08:00
|
|
|
|
2012-02-21 08:53:35 +08:00
|
|
|
PredecessorMap<Key> p_map;
|
|
|
|
p_map.insert(kx(1), kx(2));
|
|
|
|
p_map.insert(kx(2), kx(2));
|
|
|
|
p_map.insert(kx(3), kx(2));
|
|
|
|
tie(graph, root, key2vertex) = predecessorMap2Graph<SGraph<Key>, SVertex, Key>(p_map);
|
2010-01-14 11:21:07 +08:00
|
|
|
|
|
|
|
LONGS_EQUAL(3, boost::num_vertices(graph));
|
2012-02-21 08:53:35 +08:00
|
|
|
CHECK(root == key2vertex[kx(2)]);
|
2010-01-14 11:21:07 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
/* ************************************************************************* */
|
2010-01-13 00:12:25 +08:00
|
|
|
TEST( Graph, composePoses )
|
|
|
|
{
|
2012-01-28 10:31:44 +08:00
|
|
|
pose2SLAM::Graph graph;
|
2010-01-13 00:12:25 +08:00
|
|
|
Matrix cov = eye(3);
|
2010-01-27 10:49:58 +08:00
|
|
|
Pose2 p1(1.0, 2.0, 0.3), p2(4.0, 5.0, 0.6), p3(7.0, 8.0, 0.9), p4(2.0, 2.0, 2.9);
|
2010-08-27 03:55:40 +08:00
|
|
|
Pose2 p12=p1.between(p2), p23=p2.between(p3), p43=p4.between(p3);
|
2012-01-28 10:31:44 +08:00
|
|
|
graph.addOdometry(1,2, p12, cov);
|
|
|
|
graph.addOdometry(2,3, p23, cov);
|
|
|
|
graph.addOdometry(4,3, p43, cov);
|
2010-01-13 00:12:25 +08:00
|
|
|
|
2012-02-21 08:53:35 +08:00
|
|
|
PredecessorMap<Key> tree;
|
2012-02-07 12:02:20 +08:00
|
|
|
tree.insert(pose2SLAM::PoseKey(1),pose2SLAM::PoseKey(2));
|
|
|
|
tree.insert(pose2SLAM::PoseKey(2),pose2SLAM::PoseKey(2));
|
|
|
|
tree.insert(pose2SLAM::PoseKey(3),pose2SLAM::PoseKey(2));
|
|
|
|
tree.insert(pose2SLAM::PoseKey(4),pose2SLAM::PoseKey(3));
|
2010-01-13 00:12:25 +08:00
|
|
|
|
2010-01-27 10:49:58 +08:00
|
|
|
Pose2 rootPose = p2;
|
2010-01-13 00:12:25 +08:00
|
|
|
|
2012-02-04 01:18:32 +08:00
|
|
|
boost::shared_ptr<Values> actual = composePoses<pose2SLAM::Graph, pose2SLAM::Odometry,
|
2012-02-21 08:53:35 +08:00
|
|
|
Pose2, Key> (graph, tree, rootPose);
|
Large gtsam refactoring
To support faster development *and* better performance Richard and I pushed through a large refactoring of NonlinearFactors.
The following are the biggest changes:
1) NonLinearFactor1 and NonLinearFactor2 are now templated on Config, Key type, and X type, where X is the argument to the measurement function.
2) The measurement itself is no longer kept in the nonlinear factor. Instead, a derived class (see testVSLAMFactor, testNonlinearEquality, testPose3Factor etc...) has to implement a function to compute the errors, "evaluateErrors". Instead of (h(x)-z), it needs to return (z-h(x)), so Ax-b is an approximation of the error. IMPORTANT: evaluateErrors needs - if asked - *combine* the calculation of the function value h(x) and the derivatives dh(x)/dx. This was a major performance issue. To do this, boost::optional<Matrix&> arguments are provided, and tin EvaluateErrors you just says something like
if (H) *H = Matrix_(3,6,....);
3) We are no longer using int or strings for nonlinear factors. Instead, the preferred key type is now Symbol, defined in Key.h. This is both fast and cool: you can construct it from an int, and cast it to a strong. It also does type checking: a Symbol<Pose3,'x'> will not match a Symbol<Pose2,'x'>
4) minor: take a look at LieConfig.h: it help you avoid writing a lot of code bu automatically creating configs for a certain type. See e.g. Pose3Config.h. A "double" LieConfig is on the way - Thanks Richard and Manohar !
2010-01-14 06:25:03 +08:00
|
|
|
|
2012-02-03 00:16:46 +08:00
|
|
|
Values expected;
|
2012-02-04 01:18:32 +08:00
|
|
|
expected.insert(pose2SLAM::PoseKey(1), p1);
|
|
|
|
expected.insert(pose2SLAM::PoseKey(2), p2);
|
|
|
|
expected.insert(pose2SLAM::PoseKey(3), p3);
|
|
|
|
expected.insert(pose2SLAM::PoseKey(4), p4);
|
2010-01-13 00:12:25 +08:00
|
|
|
|
2010-01-27 10:49:58 +08:00
|
|
|
LONGS_EQUAL(4, actual->size());
|
2010-01-13 10:09:16 +08:00
|
|
|
CHECK(assert_equal(expected, *actual));
|
2010-01-13 00:12:25 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
/* ************************************************************************* */
|
|
|
|
int main() {
|
|
|
|
TestResult tr;
|
|
|
|
return TestRegistry::runAllTests(tr);
|
|
|
|
}
|
|
|
|
/* ************************************************************************* */
|