From 0fdb91244641f4aba8462c5d989676815e072fba Mon Sep 17 00:00:00 2001 From: guido Date: Fri, 9 Jul 2004 14:58:26 +0000 Subject: [PATCH] attempt to W1infty-norm git-svn-id: https://svn.dealii.org/trunk@9496 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/include/numerics/vectors.h | 34 ++++++++++++++-------- deal.II/deal.II/source/numerics/vectors.cc | 21 +++++++++++++ 2 files changed, 43 insertions(+), 12 deletions(-) diff --git a/deal.II/deal.II/include/numerics/vectors.h b/deal.II/deal.II/include/numerics/vectors.h index 2a0584c162..71673f92cb 100644 --- a/deal.II/deal.II/include/numerics/vectors.h +++ b/deal.II/deal.II/include/numerics/vectors.h @@ -325,6 +325,17 @@ class VectorTools * computed on each cell. */ L2_norm, + /** + * The absolute value to + * the pth power is + * integrated and the pth + * root is computed on each + * cell. The exponent + * p is the last + * parameter of the + * function. + */ + Lp_norm, /** * The maximum absolute * value of the function. @@ -341,17 +352,6 @@ class VectorTools * of the #H1_seminorm. */ H1_norm, - /** - * The absolute value to - * the pth power is - * integrated and the pth - * root is computed on each - * cell. The exponent - * p is the last - * parameter of the - * function. - */ - Lp_norm, /** * #Lp_norm of the gradient. */ @@ -360,7 +360,17 @@ class VectorTools * same as #H1_norm for * Lp. */ - W1p_norm + W1p_norm, + /** + * #Linfty_norm of the gradient. + */ + W1infty_seminorm, + /** + * same as #H1_norm for + * Linfty. + */ + W1infty_norm + }; /** diff --git a/deal.II/deal.II/source/numerics/vectors.cc b/deal.II/deal.II/source/numerics/vectors.cc index cededd77a9..73ae7cb54c 100644 --- a/deal.II/deal.II/source/numerics/vectors.cc +++ b/deal.II/deal.II/source/numerics/vectors.cc @@ -1229,10 +1229,12 @@ VectorTools::integrate_difference (const Mapping &mapping, { case H1_seminorm: case W1p_seminorm: + case W1infty_seminorm: update_flags |= UpdateFlags (update_gradients); break; case H1_norm: case W1p_norm: + case W1infty_norm: update_flags |= UpdateFlags (update_gradients); // no break! default: @@ -1402,6 +1404,7 @@ VectorTools::integrate_difference (const Mapping &mapping, diff=std::sqrt(diff); break; case Linfty_norm: + case W1infty_norm: std::fill_n (psi_scalar.begin(), n_q_points, 0.0); for (unsigned int k=0; k &mapping, 0.0); diff = std::sqrt(diff); break; + case W1infty_seminorm: + case W1infty_norm: + Assert(false, ExcNotImplemented()); + std::fill_n (psi_scalar.begin(), n_q_points, 0.0); + for (unsigned int k=0; k