]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fixed bug that prevented mapping q eulerian to work properly in codim one.
authorheltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 23 Sep 2010 15:12:02 +0000 (15:12 +0000)
committerheltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 23 Sep 2010 15:12:02 +0000 (15:12 +0000)
git-svn-id: https://svn.dealii.org/trunk@22133 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/mapping_q_eulerian.cc

index d3d936a2ea0b22c34e286c5c7b0aa872ca52ca1f..314270e7a61cdcd681a1ab61b8453a70802250b1 100644 (file)
@@ -153,7 +153,7 @@ compute_mapping_support_points
   const unsigned int n_support_pts = support_quadrature.n_quadrature_points;
   const unsigned int n_components  = euler_dof_handler->get_fe().n_components();
 
-  Assert (n_components >= dim, ExcWrongNoOfComponents() );
+  Assert (n_components >= spacedim, ExcWrongNoOfComponents() );
 
   std::vector<Vector<double> > shift_vector(n_support_pts,Vector<double>(n_components));
 
@@ -175,7 +175,7 @@ compute_mapping_support_points
   for(unsigned int q=0; q<n_support_pts; ++q)
     {
       a[q] = fe_values.quadrature_point(q);
-      for(unsigned int d=0; d<dim; ++d)
+      for(unsigned int d=0; d<spacedim; ++d)
        a[q](d) += shift_vector[q](d);
     }
 }

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.