diff --git a/gtsam_unstable/linear/tests/testLPSolver.cpp b/gtsam_unstable/linear/tests/testLPSolver.cpp
index 7045040fe..6d64c6314 100644
--- a/gtsam_unstable/linear/tests/testLPSolver.cpp
+++ b/gtsam_unstable/linear/tests/testLPSolver.cpp
@@ -341,7 +341,7 @@ public:
     if (Aterms.size() > 0) {
       Vector b = zero(delta.at(key).size());
       Factor::const_iterator it = lp_.cost.find(key);
-      if (it != lp_.cost.end()) b = lp_.cost.getA(it);
+      if (it != lp_.cost.end()) b = lp_.cost.getA(it).transpose();
       return boost::make_shared<JacobianFactor>(Aterms, b); // compute the least-square approximation of dual variables
     }
     else {