From 534dbceeccdaa65d41b0f704def17fbdfd35a1bc Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 22 Oct 2007 16:24:31 +0000
Subject: [PATCH] compute gradient analytically and not be fd approximation

[[Imported from SVN: r1711]]
---
 src/rodsolver.cc | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/src/rodsolver.cc b/src/rodsolver.cc
index f2ea5ee5..fa7a6b2e 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_);
 
-- 
GitLab