From b66e7577770f6d4730a330d903e53e46eadd6083 Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Mon, 24 Feb 2003 15:53:02 +0000
Subject: [PATCH] One more place where we have to do something special for
 non-primitive elements.

git-svn-id: https://svn.dealii.org/trunk@7231 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/deal.II/source/dofs/dof_tools.cc | 27 ++++++++++++++++++++++--
 1 file changed, 25 insertions(+), 2 deletions(-)

diff --git a/deal.II/deal.II/source/dofs/dof_tools.cc b/deal.II/deal.II/source/dofs/dof_tools.cc
index 8b7daae6c5..a36065742a 100644
--- a/deal.II/deal.II/source/dofs/dof_tools.cc
+++ b/deal.II/deal.II/source/dofs/dof_tools.cc
@@ -509,10 +509,33 @@ DoFTools::make_flux_sparsity_pattern (const DoFHandler<dim> &dof,
 
   for (unsigned int i=0; i<total_dofs; ++i)
     {
+      const unsigned int ii
+        = (dof.get_fe().is_primitive(i) ?
+           dof.get_fe().system_to_component_index(i).first
+           :
+           (std::find (dof.get_fe().get_nonzero_components(i).begin(),
+                       dof.get_fe().get_nonzero_components(i).end(),
+                       true)
+            -
+           dof.get_fe().get_nonzero_components(i).begin())
+           );
+      Assert (ii < dof.get_fe().n_components(),
+              ExcInternalError());
+
       for (unsigned int j=0; j<total_dofs; ++j)
 	{
-	  unsigned int ii = fe.system_to_component_index(i).first;
-	  unsigned int jj = fe.system_to_component_index(j).first;
+          const unsigned int jj
+            = (dof.get_fe().is_primitive(j) ?
+               dof.get_fe().system_to_component_index(j).first
+               :
+               (std::find (dof.get_fe().get_nonzero_components(j).begin(),
+                           dof.get_fe().get_nonzero_components(j).end(),
+                           true)
+                -
+                dof.get_fe().get_nonzero_components(j).begin())
+               );
+          Assert (jj < dof.get_fe().n_components(),
+                  ExcInternalError());          
 	  
 	  if (int_mask(ii,jj) != 0)
 	    int_dof_mask(i,j) = 'f';
-- 
2.39.5