From: guido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Fri, 9 Jul 2004 14:58:26 +0000 (+0000)
Subject: attempt to W1infty-norm
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=0fdb91244641f4aba8462c5d989676815e072fba;p=dealii-svn.git

attempt to W1infty-norm


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

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 <i>p</i>th power is
+					    * integrated and the pth
+					    * root is computed on each
+					    * cell. The exponent
+					    * <i>p</i> 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 <i>p</i>th power is
-					    * integrated and the pth
-					    * root is computed on each
-					    * cell. The exponent
-					    * <i>p</i> 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
 					    * <i>L<sup>p</sup></i>.
 					    */
-	  W1p_norm
+	  W1p_norm,
+					   /**
+					    * #Linfty_norm of the gradient.
+					    */
+	  W1infty_seminorm,
+					   /**
+					    * same as #H1_norm for
+					    * <i>L<sup>infty</sup></i>.
+					    */
+	  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<dim>    &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<dim>    &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<n_components; ++k)
 	      for (unsigned int q=0; q<n_q_points; ++q)
@@ -1454,6 +1457,24 @@ VectorTools::integrate_difference (const Mapping<dim>    &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<n_components; ++k)
+	      for (unsigned int q=0; q<n_q_points; ++q)
+		{
+		  double t = 0.;
+		  for (unsigned int d=0;d<dim;++d)
+		    t = std::max(t,std::fabs(psi_grads[q][k][d])
+				 * weight_vectors[q](k));
+		  
+		  psi_scalar[q] = std::max(psi_scalar[q],t);
+		}
+
+	    for (unsigned int i=0;i<psi_scalar.size();++i)
+	      diff = std::max (diff, psi_scalar[i]);
+	    break;
 	  default:
 	    break;
  	}