12345qiupeng
/
gtsam
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Comment out the test for raw memory access

release/4.3a0
Sungtae An 2014-11-12 04:19:44 -05:00
parent a143815e79
commit 102974588a
1 changed files with 17 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

33
gtsam/linear/tests/testGaussianFactorGraph.cpp
View File

@ -316,27 +316,28 @@ TEST( GaussianFactorGraph, multiplyHessianAdd2 )
}
/* ************************************************************************* */
TEST( GaussianFactorGraph, multiplyHessianAdd3 )
{
GaussianFactorGraph gfg = createGaussianFactorGraphWithHessianFactor();
// Raw memory access
//TEST( GaussianFactorGraph, multiplyHessianAdd3 )
//{
// GaussianFactorGraph gfg = createGaussianFactorGraphWithHessianFactor();
// brute force
Matrix AtA; Vector eta; boost::tie(AtA,eta) = gfg.hessian();
Vector X(6); X<<1,2,3,4,5,6;
Vector Y(6); Y<<-450, -450, 300, 400, 2950, 3450;
EXPECT(assert_equal(Y,AtA*X));
// // brute force
// Matrix AtA; Vector eta; boost::tie(AtA,eta) = gfg.hessian();
// Vector X(6); X<<1,2,3,4,5,6;
// Vector Y(6); Y<<-450, -450, 300, 400, 2950, 3450;
// EXPECT(assert_equal(Y,AtA*X));
double* x = &X[0];
// double* x = &X[0];
Vector fast_y = gtsam::zero(6);
gfg.multiplyHessianAdd(1.0, x, fast_y.data());
EXPECT(assert_equal(Y, fast_y));
// Vector fast_y = gtsam::zero(6);
// gfg.multiplyHessianAdd(1.0, x, fast_y.data());
// EXPECT(assert_equal(Y, fast_y));
// now, do it with non-zero y
gfg.multiplyHessianAdd(1.0, x, fast_y.data());
EXPECT(assert_equal(2*Y, fast_y));
// // now, do it with non-zero y
// gfg.multiplyHessianAdd(1.0, x, fast_y.data());
// EXPECT(assert_equal(2*Y, fast_y));
}
//}
/* ************************************************************************* */