diff --git a/src/rodsolver.cc b/src/rodsolver.cc
index 5654066e31d6b0cc18724683b6eec66d496691e2..76b0c8681bbd4044ae57603ab35eb80b5fd3b55b 100644
--- a/src/rodsolver.cc
+++ b/src/rodsolver.cc
@@ -227,7 +227,8 @@ void RodSolver<GridType>::solve()
         corr = 0;
 
         std::cout << "Rod energy: " <<rodAssembler_->computeEnergy(x_) << std::endl;
-        rodAssembler_->assembleGradient(x_, rhs);
+        //rodAssembler_->assembleGradient(x_, rhs);
+        rodAssembler_->assembleGradientFD(x_, rhs);
         //rodAssembler_->assembleMatrix(x_, *hessianMatrix_);
         rodAssembler_->assembleMatrixFD(x_, *hessianMatrix_);