From 52b64a86e71c10a377981e9d251d35d42453bf20 Mon Sep 17 00:00:00 2001
From: guido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 13 Nov 2000 11:19:48 +0000
Subject: [PATCH] stupid error removed

git-svn-id: https://svn.dealii.org/trunk@3481 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/base/source/function_derivative.cc | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/deal.II/base/source/function_derivative.cc b/deal.II/base/source/function_derivative.cc
index 5866eb9747..ad76305d02 100644
--- a/deal.II/base/source/function_derivative.cc
+++ b/deal.II/base/source/function_derivative.cc
@@ -58,7 +58,7 @@ FunctionDerivative<dim>::value (const Point<dim>   &p,
   switch (formula)
     {
     case Euler:
-      return (f.value(p+incr, component)-f.value(p-incr, component))/2*h;
+      return (f.value(p+incr, component)-f.value(p-incr, component))/(2*h);
     case UpwindEuler:
       return (f.value(p, component)-f.value(p-incr, component))/h;
     default:
@@ -91,8 +91,11 @@ FunctionDerivative<dim>::value_list (const vector<Point<dim> > &points,
 	}
       f.value_list(p1, values, component);
       f.value_list(p2, e2, component);
+      
       for (unsigned int i=0;i<n;++i)
-	values[i] = (values[i]-e2[i])/2*h;
+	{
+	  values[i] = (values[i]-e2[i])/(2*h);
+	}
       break;
     }    
     case UpwindEuler:
-- 
2.39.5