From: bangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Wed, 15 Dec 2010 15:37:41 +0000 (+0000)
Subject: Simplify code some more.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=25e0233bb42f52a586159c13984dc3ebb87565db;p=dealii-svn.git

Simplify code some more.


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

diff --git a/deal.II/source/fe/mapping_q1.cc b/deal.II/source/fe/mapping_q1.cc
index d272ec8515..376d74f174 100644
--- a/deal.II/source/fe/mapping_q1.cc
+++ b/deal.II/source/fe/mapping_q1.cc
@@ -995,11 +995,22 @@ namespace internal
       	    Assert (JxW_values.size() == n_q_points,
       		    ExcDimensionMismatch(JxW_values.size(), n_q_points));
 
-      	  mapping.transform (data.unit_tangentials[face_no],
-	  		     data.aux[0],
-	  		     data,
-	  		     mapping_contravariant);
-
+					   // map the unit tangentials to the real cell
+	  for (unsigned int d=0; d<dim-1; ++d)
+	    {
+	      Assert (face_no+GeometryInfo<dim>::faces_per_cell*d <
+		      data.unit_tangentials.size(),
+		      ExcInternalError());
+	      Assert (data.aux[d].size() <=
+		      data.unit_tangentials[face_no+GeometryInfo<dim>::faces_per_cell*d].size(),
+		      ExcInternalError());
+	      
+	      mapping.transform (data.unit_tangentials[face_no+GeometryInfo<dim>::faces_per_cell*d],
+				 data.aux[d],
+				 data,
+				 mapping_contravariant);
+	    }
+	  
       	  switch (dim)
       	    {
       	      case 2:
@@ -1011,19 +1022,8 @@ namespace internal
 
       	      case 3:
       	      {
-      	  	Assert (face_no+GeometryInfo<dim>::faces_per_cell <
-      	  		data.unit_tangentials.size(),
-      	  		ExcInternalError());
-      	  	Assert (data.aux[1].size() <=
-      	  		data.unit_tangentials[face_no+GeometryInfo<dim>::faces_per_cell].size(),
-      	  		ExcInternalError());
-      	  	mapping.transform (data.unit_tangentials[face_no+GeometryInfo<dim>::faces_per_cell],
-				   data.aux[1],
-				   data,
-				   mapping_contravariant);
       	  	for (unsigned int i=0; i<n_q_points; ++i)
       	  	  cross_product (boundary_forms[i], data.aux[0][i], data.aux[1][i]);
-
       	  	break;
       	      }