From: Wolfgang Bangerth Date: Mon, 24 Jul 2000 08:16:05 +0000 (+0000) Subject: Add function for general dimension. X-Git-Tag: v8.0.0~20229 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=510973ca3a8f650017a1cb339f85705786272fde;p=dealii.git 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 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