diff --git a/src/rodsolver.cc b/src/rodsolver.cc
index f2ea5ee58dfe48af2fa8185e93c186e1480cf795..fa7a6b2e9d75c279ce3436da74959bd8e5e17f20 100644
--- a/src/rodsolver.cc
+++ b/src/rodsolver.cc
@@ -230,8 +230,7 @@ void RodSolver<GridType>::solve()
 
         if (this->verbosity_ == FULL)
             std::cout << "Rod energy: " <<rodAssembler_->computeEnergy(x_) << std::endl;
-        //rodAssembler_->assembleGradient(x_, rhs);
-        rodAssembler_->assembleGradientFD(x_, rhs);
+        rodAssembler_->assembleGradient(x_, rhs);
         //rodAssembler_->assembleMatrix(x_, *hessianMatrix_);
         rodAssembler_->assembleMatrixFD(x_, *hessianMatrix_);