From cff977548b25efaa56962a70127bdbe765c685fc Mon Sep 17 00:00:00 2001 From: wolf 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 +inline +double +DerivativeApproximation::SecondDerivative:: +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 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 +inline +double +DerivativeApproximation::SecondDerivative:: +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 inline void -- 2.39.5