From b5876c653a7f4be18ccd8a74e19697da6f3807ca Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 15 Aug 2007 09:46:30 +0000
Subject: [PATCH] use fd approximation for the gradient as well

[[Imported from SVN: r1525]]
---
 src/rodsolver.cc | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/rodsolver.cc b/src/rodsolver.cc
index 5654066e..76b0c868 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_);
 
-- 
GitLab