Indent

source/material_model/rheology/elasticity.cc

@@ -443,7 +443,7 @@ namespace aspect
                 // for the fields, which for the first nonlinear iteration means the extrapolated solution from
                 // the last and previous to last timesteps. In later iterations, it holds the current solution.
                 // This can lead to alternately computing a reaction term that is correct and one that is zero.
-                // Therefore we subtract the old stress stress_t. 
+                // Therefore we subtract the old stress stress_t.
                 for (unsigned int j = 0; j < SymmetricTensor<2, dim>::n_independent_components; ++j)
                   {
                     //out.reaction_terms[i][j] = -in.composition[i][j] + stress_0[SymmetricTensor<2, dim>::unrolled_to_component_indices(j)];

source/material_model/viscoelastic.cc

@@ -38,7 +38,7 @@ namespace aspect
 
       // Store which components to exclude during volume fraction computation.
       ComponentMask composition_mask(this->n_compositional_fields(), true);
-      // Assign compositional fields associated with the viscoelastic stress 
+      // Assign compositional fields associated with the viscoelastic stress
       // or the old viscoelastic stress a value of 0.
       // Assume these fields are listed first.
       for (unsigned int i=0; i < 2*SymmetricTensor<2,dim>::n_independent_components; ++i)

source/particle/property/elastic_stress.cc

@@ -119,7 +119,7 @@ namespace aspect
           // plugins, the data position is moved forward each time. However, here we want to apply
           // an update to properties that have already been updated by another plugin. I.e., we
           // want to update the 'composition' properties, which carry the stress tensor components.
-          // This plugin is necessary because it updates the composition properties with the 
+          // This plugin is necessary because it updates the composition properties with the
           // newest solution, i.e. the one including the operator splitting update.
           // The current plugin then applies the reaction term (the rotation). Note that the reaction_terms
           // will be zero for ve_stress_**_old, so we only loop over n_independent_components instead of
@@ -151,8 +151,8 @@ namespace aspect
       std::vector<std::pair<std::string, unsigned int>>
       ElasticStress<dim>::get_property_information() const
       {
-         // No need to return anything.
-         return std::vector<std::pair<std::string, unsigned int>>();
+        // No need to return anything.
+        return std::vector<std::pair<std::string, unsigned int>>();
       }
     }
   }

source/simulator/solver_schemes.cc

@@ -344,14 +344,14 @@ namespace aspect
     stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_yy"));
     if (dim == 2)
       {
-       stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_xy"));
+        stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_xy"));
       }
     else if (dim == 3)
       {
-       stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_zz"));
-       stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_xy"));
-       stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_xz"));
-       stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_yz"));
+        stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_zz"));
+        stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_xy"));
+        stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_xz"));
+        stress_indices.push_back(introspection.compositional_index_for_name("ve_stress_yz"));
       }