]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix a problem where we weren't prepared to deal with hp cells in 1d.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 4 May 2010 18:24:24 +0000 (18:24 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 4 May 2010 18:24:24 +0000 (18:24 +0000)
git-svn-id: https://svn.dealii.org/trunk@21064 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/include/numerics/vectors.templates.h

index 4adc1b6ccc7f4f0dedab2edb9a894000c3f9bcee..132b0fd96fc0ab042a83ad2c7a385affa3f5bb77 100644 (file)
@@ -1301,7 +1301,9 @@ VectorTools::interpolate_boundary_values (const Mapping<DH::dimension, DH::space
 
   for (unsigned int i=0; i<fe.dofs_per_vertex; ++i)
     if (component_mask[fe.face_system_to_component_index(i).first])
-      boundary_values[outermost_cell->vertex_dof_index(direction,i)]
+      boundary_values[outermost_cell->
+                     vertex_dof_index(direction,i,
+                                      outermost_cell->active_fe_index())]
        = function_values(fe.face_system_to_component_index(i).first);
 }
 
@@ -1818,7 +1820,8 @@ VectorTools::interpolate_boundary_values (const Mapping<DH::dimension, DH::space
       {
                                   // TODO: should we clear the other
                                   // entries in the line here?
-       const unsigned int row = outermost_cell->vertex_dof_index(direction,i);
+       const unsigned int row = outermost_cell->vertex_dof_index(direction,i,
+                                                                 outermost_cell->active_fe_index());
        constraints.add_line (row);
        constraints.set_inhomogeneity (row,
                    function_values(fe.face_system_to_component_index(i).first));

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.