]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Make the same kind of changes also to FESystem. 1845/head
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 8 Nov 2015 01:49:53 +0000 (19:49 -0600)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 8 Nov 2015 01:49:53 +0000 (19:49 -0600)
source/fe/fe_system.cc

index 0dd710f6629fd10b33b1fc919ddea476ab0d0b48..51697cfda7d189b1c9f868983fb2d5579f10bf6a 100644 (file)
@@ -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];
+
               }
       }
 }

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.