diff --git a/ChemistryLib/CreateChemicalSolverInterface.cpp b/ChemistryLib/CreateChemicalSolverInterface.cpp
index 9da6a484ab594fbaa9cf3190f68d20c4285edea9..8ad2d0acc01bc295d681d6c247ec3336577daba9 100644
--- a/ChemistryLib/CreateChemicalSolverInterface.cpp
+++ b/ChemistryLib/CreateChemicalSolverInterface.cpp
@@ -91,7 +91,7 @@ createChemicalSolverInterface<ChemicalSolver::Phreeqc>(
     auto const ls_name =
         //! \ogs_file_param{prj__chemical_system__linear_solver}
         config.getConfigParameter<std::string>("linear_solver");
-    auto& linear_solver = BaseLib::getOrError(
+    auto const& linear_solver = BaseLib::getOrError(
         linear_solvers, ls_name,
         "A linear solver with the given name does not exist.");
 
diff --git a/ChemistryLib/PhreeqcKernelData/ReactionRate.cpp b/ChemistryLib/PhreeqcKernelData/ReactionRate.cpp
index 9375a7709e38764ebeae24ec8a3e3d9cc96085bc..c3f7a83329a33cacd9553c10bfac3b36ab7ab8a8 100644
--- a/ChemistryLib/PhreeqcKernelData/ReactionRate.cpp
+++ b/ChemistryLib/PhreeqcKernelData/ReactionRate.cpp
@@ -14,8 +14,8 @@ namespace ChemistryLib
 {
 namespace PhreeqcKernelData
 {
-ReactionRate::ReactionRate(
-    std::string kinetic_reactant_, std::vector<std::string> statements)
+ReactionRate::ReactionRate(std::string kinetic_reactant_,
+                           std::vector<std::string> const& statements)
     : kinetic_reactant(std::move(kinetic_reactant_))
 {
     int line_number = 1;
diff --git a/ChemistryLib/PhreeqcKernelData/ReactionRate.h b/ChemistryLib/PhreeqcKernelData/ReactionRate.h
index 0050de792267f7766e4c9f91f61b266c831a6d0b..796dc6ec6db5bef9c0d76389b7cac3a2e02a8c6d 100644
--- a/ChemistryLib/PhreeqcKernelData/ReactionRate.h
+++ b/ChemistryLib/PhreeqcKernelData/ReactionRate.h
@@ -21,7 +21,7 @@ class ReactionRate
 {
 public:
     ReactionRate(std::string kinetic_reactant_,
-                 std::vector<std::string> statements);
+                 std::vector<std::string> const& statements);
 
     std::string const& commands() const { return _commands; }
 
diff --git a/GeoLib/GEOObjects.cpp b/GeoLib/GEOObjects.cpp
index ad16d68e75dbe9a67021f72b8872c089ee32ef5c..a062996c3f90e886ab9d718d9ba246eb58568118 100644
--- a/GeoLib/GEOObjects.cpp
+++ b/GeoLib/GEOObjects.cpp
@@ -330,7 +330,7 @@ bool GEOObjects::isPntVecUsed(const std::string& name) const
 
 void GEOObjects::getStationVectorNames(std::vector<std::string>& names) const
 {
-    for (auto point : _pnt_vecs)
+    for (auto const* point : _pnt_vecs)
     {
         if (point->getType() == PointVec::PointType::STATION)
         {
@@ -342,7 +342,7 @@ void GEOObjects::getStationVectorNames(std::vector<std::string>& names) const
 std::vector<std::string> GEOObjects::getGeometryNames() const
 {
     std::vector<std::string> names;
-    for (auto const point : _pnt_vecs)
+    for (auto const* point : _pnt_vecs)
     {
         if (point->getType() == PointVec::PointType::POINT)
         {
diff --git a/MeshLib/IO/XDMF/HdfWriter.cpp b/MeshLib/IO/XDMF/HdfWriter.cpp
index 2728c1d1637d825653397c019464a583045e5a8f..f6e2290b988f4536380c469b860011ecc6769c55 100644
--- a/MeshLib/IO/XDMF/HdfWriter.cpp
+++ b/MeshLib/IO/XDMF/HdfWriter.cpp
@@ -206,7 +206,7 @@ struct HdfWriter::HdfMesh final
     std::vector<HdfData> const variable_attributes;
 };
 
-HdfWriter::HdfWriter(std::vector<MeshHdfData> meshes,
+HdfWriter::HdfWriter(std::vector<MeshHdfData> const& meshes,
                      unsigned long long const initial_step,
                      std::filesystem::path const& filepath,
                      bool const use_compression,
diff --git a/MeshLib/IO/XDMF/HdfWriter.h b/MeshLib/IO/XDMF/HdfWriter.h
index 3e400b8e4036786c5435a07e5f48e0c339f70c18..48a7d8396edec7cbc9c0f2964ce62bd434b6867c 100644
--- a/MeshLib/IO/XDMF/HdfWriter.h
+++ b/MeshLib/IO/XDMF/HdfWriter.h
@@ -45,7 +45,7 @@ public:
      * @param is_file_manager True if process (in parallel execution) is
      * @param n_files Number of output files
      */
-    HdfWriter(std::vector<MeshHdfData> meshes,
+    HdfWriter(std::vector<MeshHdfData> const& meshes,
               unsigned long long initial_step,
               std::filesystem::path const& filepath,
               bool use_compression,
diff --git a/MeshLib/Vtk/VtkMappedMeshSource.cpp b/MeshLib/Vtk/VtkMappedMeshSource.cpp
index 435c24d90cbe8246a195265ffa2134f5b71cb255..8d10c9da60901bd47049252596086c43dbdfbf1a 100644
--- a/MeshLib/Vtk/VtkMappedMeshSource.cpp
+++ b/MeshLib/Vtk/VtkMappedMeshSource.cpp
@@ -123,54 +123,59 @@ int VtkMappedMeshSource::RequestData(vtkInformation* /*request*/,
         {
             continue;
         }
-        if (auto p = dynamic_cast<PropertyVector<double>*>(property))
+        if (auto const* p = dynamic_cast<PropertyVector<double>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<float>*>(property))
+        else if (auto const* p = dynamic_cast<PropertyVector<float>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<int>*>(property))
+        else if (auto const* p = dynamic_cast<PropertyVector<int>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<unsigned>*>(property))
+        else if (auto const* p =
+                     dynamic_cast<PropertyVector<unsigned>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<long>*>(property))
+        else if (auto const* p = dynamic_cast<PropertyVector<long>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<long long>*>(property))
+        else if (auto const* p =
+                     dynamic_cast<PropertyVector<long long>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<PropertyVector<unsigned long>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<unsigned long long>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<PropertyVector<unsigned long long>*>(
+                         property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<std::size_t>*>(property))
+        else if (auto const* p =
+                     dynamic_cast<PropertyVector<std::size_t>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<char>*>(property))
+        else if (auto const* p = dynamic_cast<PropertyVector<char>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<PropertyVector<unsigned char>*>(property))
         {
             addProperty(*p);
         }
-        else if (auto p = dynamic_cast<PropertyVector<uint8_t>*>(property))
+        else if (auto const* p =
+                     dynamic_cast<PropertyVector<uint8_t>*>(property))
         {
             addProperty(*p);
         }
diff --git a/MeshToolsLib/IntegrationPointDataTools.cpp b/MeshToolsLib/IntegrationPointDataTools.cpp
index 3320a43d821e3a2efc0bf5ce11d27843ddbf690a..1dfcc6f946f09904e41aa699b3460548c9e2110b 100644
--- a/MeshToolsLib/IntegrationPointDataTools.cpp
+++ b/MeshToolsLib/IntegrationPointDataTools.cpp
@@ -115,7 +115,7 @@ std::vector<std::size_t> getIntegrationPointDataOffsetsOfMeshElements(
         MeshLib::getIntegrationPointMetaData(properties, pv.getPropertyName());
     for (std::size_t i = 0; i < mesh_elements.size(); i++)
     {
-        auto const element = mesh_elements[i];
+        auto const* const element = mesh_elements[i];
 
         // Assuming that the order of elements in mesh_elements is not touched.
         element_ip_data_offsets[i] = counter;
diff --git a/MeshToolsLib/MeshEditing/MeshRevision.cpp b/MeshToolsLib/MeshEditing/MeshRevision.cpp
index a72e0d6e60b95c6980a00b5ad1e6c61ccbf9b551..e2474e67119789763d00594774b1c8d0353a9825 100644
--- a/MeshToolsLib/MeshEditing/MeshRevision.cpp
+++ b/MeshToolsLib/MeshEditing/MeshRevision.cpp
@@ -937,7 +937,7 @@ MeshLib::Properties copyProperties(MeshLib::Properties const& props,
         if (props.existsPropertyVector<int>(name, MeshLib::MeshItemType::Node,
                                             1))
         {
-            auto p = props.getPropertyVector<int>(
+            auto const* p = props.getPropertyVector<int>(
                 name, MeshLib::MeshItemType::Node, 1);
             auto new_node_vec = new_properties.createNewPropertyVector<int>(
                 name, MeshLib::MeshItemType::Node, 1);
@@ -947,7 +947,7 @@ MeshLib::Properties copyProperties(MeshLib::Properties const& props,
         if (props.existsPropertyVector<float>(name, MeshLib::MeshItemType::Node,
                                               1))
         {
-            auto p = props.getPropertyVector<float>(
+            auto const* p = props.getPropertyVector<float>(
                 name, MeshLib::MeshItemType::Node, 1);
             auto new_node_vec = new_properties.createNewPropertyVector<float>(
                 name, MeshLib::MeshItemType::Node, 1);
@@ -957,7 +957,7 @@ MeshLib::Properties copyProperties(MeshLib::Properties const& props,
         if (props.existsPropertyVector<double>(name,
                                                MeshLib::MeshItemType::Node, 1))
         {
-            auto p = props.getPropertyVector<double>(
+            auto const* p = props.getPropertyVector<double>(
                 name, MeshLib::MeshItemType::Node, 1);
             auto new_node_vec = new_properties.createNewPropertyVector<double>(
                 name, MeshLib::MeshItemType::Node, 1);
@@ -967,7 +967,7 @@ MeshLib::Properties copyProperties(MeshLib::Properties const& props,
         if (props.existsPropertyVector<int>(name, MeshLib::MeshItemType::Cell,
                                             1))
         {
-            auto p = props.getPropertyVector<int>(
+            auto const* p = props.getPropertyVector<int>(
                 name, MeshLib::MeshItemType::Cell, 1);
             auto new_cell_vec = new_properties.createNewPropertyVector<int>(
                 name, MeshLib::MeshItemType::Cell, 1);
@@ -977,7 +977,7 @@ MeshLib::Properties copyProperties(MeshLib::Properties const& props,
         if (props.existsPropertyVector<float>(name, MeshLib::MeshItemType::Cell,
                                               1))
         {
-            auto p = props.getPropertyVector<float>(
+            auto const* p = props.getPropertyVector<float>(
                 name, MeshLib::MeshItemType::Cell, 1);
             auto new_cell_vec = new_properties.createNewPropertyVector<float>(
                 name, MeshLib::MeshItemType::Cell, 1);
@@ -987,7 +987,7 @@ MeshLib::Properties copyProperties(MeshLib::Properties const& props,
         if (props.existsPropertyVector<double>(name,
                                                MeshLib::MeshItemType::Cell, 1))
         {
-            auto p = props.getPropertyVector<double>(
+            auto const* p = props.getPropertyVector<double>(
                 name, MeshLib::MeshItemType::Cell, 1);
             auto new_cell_vec = new_properties.createNewPropertyVector<double>(
                 name, MeshLib::MeshItemType::Cell, 1);
diff --git a/MeshToolsLib/MeshInformation.cpp b/MeshToolsLib/MeshInformation.cpp
index fe9b43eb97b0c9e34ff56e998af17f3b4b043e9e..015f140d765e133a6d65a26341c3da107eba3e83 100644
--- a/MeshToolsLib/MeshInformation.cpp
+++ b/MeshToolsLib/MeshInformation.cpp
@@ -77,51 +77,56 @@ void MeshInformation::writePropertyVectorInformation(const MeshLib::Mesh& mesh)
 
     for (auto [name, property] : properties)
     {
-        if (auto p = dynamic_cast<MeshLib::PropertyVector<double>*>(property))
+        if (auto const* p =
+                dynamic_cast<MeshLib::PropertyVector<double>*>(property))
         {
             printBounds(*p);
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<float>*>(property))
         {
             printBounds(*p);
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<int>*>(property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<int>*>(property))
         {
             printBounds(*p);
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<unsigned>*>(property))
         {
             printBounds(*p);
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<long>*>(property))
         {
             printBounds(*p);
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<long long>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<long long>*>(
+                         property))
         {
             printBounds(*p);
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<unsigned long>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<unsigned long>*>(
+                         property))
         {
             printBounds(*p);
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<unsigned long long>*>(
                          property))
         {
             printBounds(*p);
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<std::size_t>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<std::size_t>*>(
+                         property))
         {
             printBounds(*p);
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<char>*>(property))
         {
             printBounds(*p);
diff --git a/MeshToolsLib/MeshQuality/MeshValidation.cpp b/MeshToolsLib/MeshQuality/MeshValidation.cpp
index 3db81708983cb1fccda818685034d15e8ff46936..76895d4cf4ce295e2095233683e8efb991b2960a 100644
--- a/MeshToolsLib/MeshQuality/MeshValidation.cpp
+++ b/MeshToolsLib/MeshQuality/MeshValidation.cpp
@@ -145,7 +145,7 @@ std::vector<ElementErrorCode> MeshValidation::testElementGeometry(
         std::accumulate(error_count, error_count + nErrorCodes, 0.0)));
     if (error_sum != 0)
     {
-        ElementErrorFlag flags[nErrorCodes] = {
+        ElementErrorFlag const flags[nErrorCodes] = {
             ElementErrorFlag::ZeroVolume, ElementErrorFlag::NonCoplanar,
             ElementErrorFlag::NonConvex, ElementErrorFlag::NodeOrder};
         for (std::size_t i = 0; i < nErrorCodes; ++i)
diff --git a/MeshToolsLib/MeshSurfaceExtraction.cpp b/MeshToolsLib/MeshSurfaceExtraction.cpp
index a24971642502be38fe7f87841522f5d9c58d672d..7b426d62c203f98040d10af119942ea23706fcca 100644
--- a/MeshToolsLib/MeshSurfaceExtraction.cpp
+++ b/MeshToolsLib/MeshSurfaceExtraction.cpp
@@ -93,60 +93,65 @@ bool createSfcMeshProperties(MeshLib::Mesh& sfc_mesh,
         }
 
         auto const& id_map = *id_maps.at(property->getMeshItemType());
-        if (auto p = dynamic_cast<MeshLib::PropertyVector<double>*>(property))
+        if (auto const* p =
+                dynamic_cast<MeshLib::PropertyVector<double>*>(property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<float>*>(property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<int>*>(property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<int>*>(property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<unsigned>*>(property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<long>*>(property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<long long>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<long long>*>(
+                         property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<unsigned long>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<unsigned long>*>(
+                         property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<unsigned long long>*>(
                          property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<std::size_t>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<std::size_t>*>(
+                         property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
             vectors_copied++;
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<char>*>(property))
         {
             processPropertyVector(*p, id_map, sfc_mesh);
diff --git a/ProcessLib/LIE/Common/PostUtils.cpp b/ProcessLib/LIE/Common/PostUtils.cpp
index 2970b6a472062914bad30c5c7f8de60dbaac8f87..d9f550e263fcff26be585e029084cec02bc11972 100644
--- a/ProcessLib/LIE/Common/PostUtils.cpp
+++ b/ProcessLib/LIE/Common/PostUtils.cpp
@@ -212,51 +212,56 @@ PostProcessTool::PostProcessTool(
 
     for (auto [name, property] : _org_mesh.getProperties())
     {
-        if (auto p = dynamic_cast<MeshLib::PropertyVector<double>*>(property))
+        if (auto const* p =
+                dynamic_cast<MeshLib::PropertyVector<double>*>(property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<float>*>(property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<int>*>(property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<int>*>(property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<unsigned>*>(property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<long>*>(property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<long long>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<long long>*>(
+                         property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<unsigned long>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<unsigned long>*>(
+                         property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<unsigned long long>*>(
                          property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p = dynamic_cast<MeshLib::PropertyVector<std::size_t>*>(
-                     property))
+        else if (auto const* p =
+                     dynamic_cast<MeshLib::PropertyVector<std::size_t>*>(
+                         property))
         {
             copyPropertyValues(*p, createProperty(*p));
         }
-        else if (auto p =
+        else if (auto const* p =
                      dynamic_cast<MeshLib::PropertyVector<char>*>(property))
         {
             copyPropertyValues(*p, createProperty(*p));
diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture-impl.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture-impl.h
index b214c225e2d98c7a046e669b85fef8356fef6f82..84ad8e6353c1c43c0a086a403d9999979cd48dab 100644
--- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture-impl.h
+++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture-impl.h
@@ -184,7 +184,7 @@ void SmallDeformationLocalAssemblerFracture<ShapeFunction, DisplacementDim>::
         auto const& w_prev = ip_data.w_prev;
         auto& C = ip_data.C;
         auto& state = *ip_data.material_state_variables;
-        auto& N = _secondary_data.N[ip];
+        auto const& N = _secondary_data.N[ip];
 
         auto const ip_physical_coords(computePhysicalCoordinates(_element, N));
         std::vector<double> const levelsets(duGlobalEnrichments(
diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h
index a026e2d42d33a3056c8ec9c2aae9fdc19dd1020b..895e6f9b98db8ab23e3f1990ac968ec984b88c2b 100644
--- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h
+++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h
@@ -183,7 +183,7 @@ void SmallDeformationLocalAssemblerMatrix<ShapeFunction, DisplacementDim>::
         _integration_method.getNumberOfPoints();
     for (unsigned ip = 0; ip < n_integration_points; ip++)
     {
-        auto& ip_data = _ip_data[ip];
+        auto const& ip_data = _ip_data[ip];
 
         ele_stress += ip_data._sigma;
         ele_strain += ip_data._eps;
diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrixNearFracture-impl.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrixNearFracture-impl.h
index ad101a3be4d4fb38d683f4cfa313c564ed143da5..aeccadf93cd01e176a6c47a0335b470b58563a4e 100644
--- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrixNearFracture-impl.h
+++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrixNearFracture-impl.h
@@ -315,7 +315,7 @@ void SmallDeformationLocalAssemblerMatrixNearFracture<ShapeFunction,
         _integration_method.getNumberOfPoints();
     for (unsigned ip = 0; ip < n_integration_points; ip++)
     {
-        auto& ip_data = _ip_data[ip];
+        auto const& ip_data = _ip_data[ip];
 
         ele_stress += ip_data._sigma;
         ele_strain += ip_data._eps;
diff --git a/ProcessLib/TimeLoop.cpp b/ProcessLib/TimeLoop.cpp
index 28633f000deafddb0ff415bb4bfdbf13f11d00f5..ddd7026da8160826542f102c9343f8d1b775b844 100644
--- a/ProcessLib/TimeLoop.cpp
+++ b/ProcessLib/TimeLoop.cpp
@@ -151,7 +151,7 @@ setInitialConditions(
     std::vector<GlobalVector*> process_solutions;
     std::vector<GlobalVector*> process_solutions_prev;
 
-    for (auto& process_data : per_process_data)
+    for (auto const& process_data : per_process_data)
     {
         auto const process_id = process_data->process_id;
         auto& ode_sys = *process_data->tdisc_ode_sys;
@@ -165,7 +165,7 @@ setInitialConditions(
                 ode_sys.getMatrixSpecifications(process_id)));
     }
 
-    for (auto& process_data : per_process_data)
+    for (auto const& process_data : per_process_data)
     {
         auto& pcs = process_data->process;
         auto const process_id = process_data->process_id;
@@ -184,7 +184,7 @@ void calculateNonEquilibriumInitialResiduum(
     std::vector<GlobalVector*> const& process_solutions,
     std::vector<GlobalVector*> const& process_solutions_prev)
 {
-    for (auto& process_data : per_process_data)
+    for (auto const& process_data : per_process_data)
     {
         auto& nonlinear_solver = process_data->nonlinear_solver;
 
@@ -255,7 +255,7 @@ TimeLoop::TimeLoop(std::vector<Output>&& outputs,
 
 void TimeLoop::setCoupledSolutions()
 {
-    for (auto& process_data : _per_process_data)
+    for (auto const& process_data : _per_process_data)
     {
         auto const& x = *_process_solutions[process_data->process_id];
 
@@ -465,7 +465,7 @@ TimeLoop::generateOutputTimeStepConstraints(
 /// initialize output, convergence criterion, etc.
 void TimeLoop::initialize()
 {
-    for (auto& process_data : _per_process_data)
+    for (auto const& process_data : _per_process_data)
     {
         auto& pcs = process_data->process;
         for (auto& output : _outputs)
@@ -656,7 +656,7 @@ NumLib::NonlinearSolverStatus TimeLoop::solveUncoupledEquationSystems(
 {
     NumLib::NonlinearSolverStatus nonlinear_solver_status;
 
-    for (auto& process_data : _per_process_data)
+    for (auto const& process_data : _per_process_data)
     {
         auto const process_id = process_data->process_id;
         nonlinear_solver_status = solveMonolithicProcess(
@@ -796,7 +796,7 @@ TimeLoop::solveCoupledEquationSystemsByStaggeredScheme(
     }
 
     {
-        for (auto& process_data : _per_process_data)
+        for (auto const& process_data : _per_process_data)
         {
             auto& pcs = process_data->process;
             int const process_id = process_data->process_id;
@@ -826,7 +826,7 @@ void TimeLoop::outputSolutions(bool const output_initial_condition,
     bool const is_staggered_coupling =
         !isMonolithicProcess(*_per_process_data[0]);
 
-    for (auto& process_data : _per_process_data)
+    for (auto const& process_data : _per_process_data)
     {
         // If nonlinear solver diverged, the solution has already been
         // saved.