From: Denis Davydov <davydden@gmail.com>
Date: Thu, 19 May 2016 12:38:27 +0000 (+0200)
Subject: fix component_to_base_table for non-tensor product case
X-Git-Tag: v8.5.0-rc1~1025^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F2614%2Fhead;p=dealii.git

fix component_to_base_table for non-tensor product case
---

diff --git a/source/fe/fe_tools.cc b/source/fe/fe_tools.cc
index b0a6e7e2a6..7f87499c2b 100644
--- a/source/fe/fe_tools.cc
+++ b/source/fe/fe_tools.cc
@@ -568,15 +568,26 @@ namespace FETools
   {
     unsigned int total_index = 0;
 
-    for (unsigned int base=0; base < fe.n_base_elements(); ++base)
-      for (unsigned int m = 0; m < fe.element_multiplicity(base); ++m)
-        {
-          for (unsigned int k=0; k<fe.base_element(base).n_components(); ++k)
-            component_to_base_table[total_index++]
-              = std::make_pair(std::make_pair(base,k), m);
-        }
-    Assert (total_index == component_to_base_table.size(),
-            ExcInternalError());
+    if (do_tensor_product)
+      {
+        for (unsigned int base=0; base < fe.n_base_elements(); ++base)
+          for (unsigned int m = 0; m < fe.element_multiplicity(base); ++m)
+            {
+              for (unsigned int k=0; k<fe.base_element(base).n_components(); ++k)
+                component_to_base_table[total_index++]
+                  = std::make_pair(std::make_pair(base,k), m);
+            }
+        Assert (total_index == component_to_base_table.size(),
+                ExcInternalError());
+      }
+    else
+      {
+        // The base element establishing a component does not make sense in this case.
+        // Set up to something meaningless:
+        for (unsigned int i = 0; i < component_to_base_table.size(); i++)
+          component_to_base_table[i] = std::make_pair(std::make_pair(numbers::invalid_unsigned_int,numbers::invalid_unsigned_int), numbers::invalid_unsigned_int);
+
+      }
 
 
     // Initialize index tables.  Multi-component base elements have to be