| 
									
										
										
										
											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
										 |  |  |  */ | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2012-06-03 00:18:40 +08:00
										 |  |  | #include <gtsam/slam/pose2SLAM.h>
 | 
					
						
							|  |  |  | #include <gtsam/inference/graph-inl.h>
 | 
					
						
							| 
									
										
										
										
											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
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2012-06-03 00:18:40 +08:00
										 |  |  | #include <boost/shared_ptr.hpp>
 | 
					
						
							|  |  |  | #include <boost/assign/std/list.hpp> // for operator +=
 | 
					
						
							|  |  |  | using namespace boost::assign; | 
					
						
							| 
									
										
										
										
											2010-01-18 13:38:53 +08:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2012-06-03 00:18:40 +08:00
										 |  |  | #include <iostream>
 | 
					
						
							| 
									
										
										
										
											2010-01-13 00:12:25 +08:00
										 |  |  | 
 | 
					
						
							|  |  |  | using namespace std; | 
					
						
							|  |  |  | using namespace gtsam; | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											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; | 
					
						
							| 
									
										
										
										
											2012-06-03 00:18:40 +08:00
										 |  |  | 	p_map.insert(1,1); | 
					
						
							|  |  |  | 	p_map.insert(2,1); | 
					
						
							|  |  |  | 	p_map.insert(3,1); | 
					
						
							|  |  |  | 	p_map.insert(4,3); | 
					
						
							|  |  |  | 	p_map.insert(5,1); | 
					
						
							| 
									
										
										
										
											2012-02-07 12:58:11 +08:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2012-02-21 08:53:35 +08:00
										 |  |  | 	list<Key> expected; | 
					
						
							| 
									
										
										
										
											2012-06-03 00:18:40 +08:00
										 |  |  | 	expected += 4,5,3,2,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; | 
					
						
							| 
									
										
										
										
											2012-06-03 00:18:40 +08:00
										 |  |  | 	p_map.insert(1, 2); | 
					
						
							|  |  |  | 	p_map.insert(2, 2); | 
					
						
							|  |  |  | 	p_map.insert(3, 2); | 
					
						
							| 
									
										
										
										
											2012-06-05 02:46:05 +08:00
										 |  |  | 	boost::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-06-03 00:18:40 +08:00
										 |  |  | 	CHECK(root == key2vertex[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; | 
					
						
							| 
									
										
										
										
											2012-06-23 03:36:49 +08:00
										 |  |  | 	SharedNoiseModel cov = noiseModel::Unit::Create(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-06-24 10:48:12 +08:00
										 |  |  | 	graph.addRelativePose(1,2, p12, cov); | 
					
						
							|  |  |  | 	graph.addRelativePose(2,3, p23, cov); | 
					
						
							|  |  |  | 	graph.addRelativePose(4,3, p43, cov); | 
					
						
							| 
									
										
										
										
											2010-01-13 00:12:25 +08:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2012-02-21 08:53:35 +08:00
										 |  |  | 	PredecessorMap<Key> tree; | 
					
						
							| 
									
										
										
										
											2012-06-03 00:18:40 +08:00
										 |  |  | 	tree.insert(1,2); | 
					
						
							|  |  |  | 	tree.insert(2,2); | 
					
						
							|  |  |  | 	tree.insert(3,2); | 
					
						
							|  |  |  | 	tree.insert(4,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-06-03 00:18:40 +08:00
										 |  |  | 	expected.insert(1, p1); | 
					
						
							|  |  |  | 	expected.insert(2, p2); | 
					
						
							|  |  |  | 	expected.insert(3, p3); | 
					
						
							|  |  |  | 	expected.insert(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); | 
					
						
							|  |  |  | } | 
					
						
							|  |  |  | /* ************************************************************************* */ |