From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Sun, 8 Nov 2015 01:49:53 +0000 (-0600)
Subject: Make the same kind of changes also to FESystem.
X-Git-Tag: v8.4.0-rc2~242^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F1845%2Fhead;p=dealii.git

Make the same kind of changes also to FESystem.
---

diff --git a/source/fe/fe_system.cc b/source/fe/fe_system.cc
index 0dd710f662..51697cfda7 100644
--- a/source/fe/fe_system.cc
+++ b/source/fe/fe_system.cc
@@ -1197,29 +1197,31 @@ compute_fill (const Mapping<dim,spacedim>                      &mapping,
                 Assert (this->n_nonzero_components(system_index) ==
                         base_fe.n_nonzero_components(base_index),
                         ExcInternalError());
-                for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s)
-                  {
-                    if (base_flags & update_values)
-                      for (unsigned int q=0; q<n_q_points; ++q)
-                        output_data.shape_values[out_index+s][q] =
-                          base_data.shape_values(in_index+s,q);
-
-                    if (base_flags & update_gradients)
-                      for (unsigned int q=0; q<n_q_points; ++q)
-                        output_data.shape_gradients[out_index+s][q] =
-                          base_data.shape_gradients[in_index+s][q];
-
-                    if (base_flags & update_hessians)
-                      for (unsigned int q=0; q<n_q_points; ++q)
-                        output_data.shape_hessians[out_index+s][q] =
-                          base_data.shape_hessians[in_index+s][q];
-
-                    if (base_flags & update_3rd_derivatives)
-                      for (unsigned int q=0; q<n_q_points; ++q)
-                        output_data.shape_3rd_derivatives[out_index+s][q] =
-                          base_data.shape_3rd_derivatives[in_index+s][q];
 
-                  }
+                if (base_flags & update_values)
+                  for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s)
+                    for (unsigned int q=0; q<n_q_points; ++q)
+                      output_data.shape_values[out_index+s][q] =
+                        base_data.shape_values(in_index+s,q);
+
+                if (base_flags & update_gradients)
+                  for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s)
+                    for (unsigned int q=0; q<n_q_points; ++q)
+                      output_data.shape_gradients[out_index+s][q] =
+                        base_data.shape_gradients[in_index+s][q];
+
+                if (base_flags & update_hessians)
+                  for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s)
+                    for (unsigned int q=0; q<n_q_points; ++q)
+                      output_data.shape_hessians[out_index+s][q] =
+                        base_data.shape_hessians[in_index+s][q];
+
+                if (base_flags & update_3rd_derivatives)
+                  for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s)
+                    for (unsigned int q=0; q<n_q_points; ++q)
+                      output_data.shape_3rd_derivatives[out_index+s][q] =
+                        base_data.shape_3rd_derivatives[in_index+s][q];
+
               }
       }
 }