diff --git a/gtsam/linear/JacobianFactor.cpp b/gtsam/linear/JacobianFactor.cpp
index dbf2757b8..809f54b99 100644
--- a/gtsam/linear/JacobianFactor.cpp
+++ b/gtsam/linear/JacobianFactor.cpp
@@ -131,7 +131,7 @@ namespace gtsam {
         Matrix::Zero(maxrank, Ab_.matrix().cols());
     // FIXME: replace with triangular system
     Ab_.rowEnd() = maxrank;
-    model_ = noiseModel::Unit::Create(maxrank);
+    model_ = SharedDiagonal(); // should be same as Unit::Create(maxrank);
   }
 
   /* ************************************************************************* */
diff --git a/gtsam/linear/tests/testJacobianFactorUnordered.cpp b/gtsam/linear/tests/testJacobianFactorUnordered.cpp
index d280a5827..a417e9229 100644
--- a/gtsam/linear/tests/testJacobianFactorUnordered.cpp
+++ b/gtsam/linear/tests/testJacobianFactorUnordered.cpp
@@ -138,8 +138,7 @@ TEST(JabobianFactor, Hessian_conversion) {
   JacobianFactor expected(0, (Matrix(2,4) <<
       1.2530,   2.1508,   -0.8779,  -1.8755,
            0,   2.5858,    0.4789,  -2.3943).finished(),
-      (Vector(2) << -6.2929, -5.7941).finished(),
-      noiseModel::Unit::Create(2));
+      (Vector(2) << -6.2929, -5.7941).finished());
 
   JacobianFactor actual(hessian);