Add function for general dimension.

git-svn-id: https://svn.dealii.org/trunk@3197 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/source/numerics/derivative_approximation.cc b/deal.II/deal.II/source/numerics/derivative_approximation.cc
index 5b6db2d9eba4c675e7b439bc201d7fbb7f8bf478..7efe887e0691ba6a2894f02965fe9daa2dbb86da 100644 (file)
--- 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 5b6db2d9eba4c675e7b439bc201d7fbb7f8bf478..7efe887e0691ba6a2894f02965fe9daa2dbb86da 100644 (file)
--- 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