From 510973ca3a8f650017a1cb339f85705786272fde Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Mon, 24 Jul 2000 08:16:05 +0000
Subject: [PATCH] Add function for general dimension.

git-svn-id: https://svn.dealii.org/trunk@3197 0785d39b-7218-0410-832d-ea1e28bc413d
---
 .../numerics/derivative_approximation.cc      | 23 +++++++++++++++++++
 .../source/numerics/gradient_estimator.cc     | 23 +++++++++++++++++++
 2 files changed, 46 insertions(+)

diff --git a/deal.II/deal.II/source/numerics/derivative_approximation.cc b/deal.II/deal.II/source/numerics/derivative_approximation.cc
index 5b6db2d9eb..7efe887e06 100644
--- a/deal.II/deal.II/source/numerics/derivative_approximation.cc
+++ b/deal.II/deal.II/source/numerics/derivative_approximation.cc
@@ -121,6 +121,29 @@ derivative_norm (const Derivative &d)
 
 
 
+template <int dim>
+inline
+double
+DerivativeApproximation::SecondDerivative<dim>::
+derivative_norm (const Derivative &d)
+{
+				   // computing the spectral norm is
+				   // not so simple in general. it is
+				   // feasible for dim==3, since then
+				   // there are still closed form
+				   // expressions of the roots of the
+				   // third order characteristic
+				   // polynomial, and they can easily
+				   // be computed using
+				   // maple. however, for higher
+				   // dimensions, some other method
+				   // needs to be employed.
+  Assert (false, ExcNotImplemented());
+  return 0;
+};
+
+
+
 template <int dim>
 inline
 void
diff --git a/deal.II/deal.II/source/numerics/gradient_estimator.cc b/deal.II/deal.II/source/numerics/gradient_estimator.cc
index 5b6db2d9eb..7efe887e06 100644
--- a/deal.II/deal.II/source/numerics/gradient_estimator.cc
+++ b/deal.II/deal.II/source/numerics/gradient_estimator.cc
@@ -121,6 +121,29 @@ derivative_norm (const Derivative &d)
 
 
 
+template <int dim>
+inline
+double
+DerivativeApproximation::SecondDerivative<dim>::
+derivative_norm (const Derivative &d)
+{
+				   // computing the spectral norm is
+				   // not so simple in general. it is
+				   // feasible for dim==3, since then
+				   // there are still closed form
+				   // expressions of the roots of the
+				   // third order characteristic
+				   // polynomial, and they can easily
+				   // be computed using
+				   // maple. however, for higher
+				   // dimensions, some other method
+				   // needs to be employed.
+  Assert (false, ExcNotImplemented());
+  return 0;
+};
+
+
+
 template <int dim>
 inline
 void
-- 
2.39.5