diff --git a/heat/src/wSolver.cpp b/heat/src/wSolver.cpp
index 22a68553c23fe93fb27a419cea5dfc0975062d77..cc3f53dd6a95dd7fbe26545863a127895831d097 100644
--- a/heat/src/wSolver.cpp
+++ b/heat/src/wSolver.cpp
@@ -256,7 +256,7 @@ void Solver::start(MshExport *mshWriter)
     Eigen::Vector2d qM = Eigen::Vector2d::Zero();
     for (auto mat : pbl->media)
     {
-        tbb::parallel_for_each(mat->tag->elems.begin(), mat->tag->elems.end(),[&](Element *e) {
+        tbb::parallel_for_each(mat->tag->elems.begin(), mat->tag->elems.end(), [&](Element *e) {
             if (e->type() != ELTYPE::TRI3)
                 return;
             //std::cout << "processing element #" << e->no << "\n";
@@ -372,23 +372,23 @@ void Solver::buildK(Eigen::SparseMatrix<double, Eigen::RowMajor> &K2, Eigen::Map
             if (verbose > 1)
                 std::cout << "\tprocessing " << *mat << '\n';
             tbb::parallel_for_each(mat->tag->elems.begin(), mat->tag->elems.end(), [&](Element *e) {
-                    if (e->type() != ELTYPE::TRI3)
-                        return;
-                    //std::cout << "processing element #" << e->no << "\n";
+                if (e->type() != ELTYPE::TRI3)
+                    return;
+                //std::cout << "processing element #" << e->no << "\n";
 
-                    Eigen::MatrixXd Ke = HeatTerm::build(*e, T1, *mat->k, false);
+                Eigen::MatrixXd Ke = HeatTerm::build(*e, T1, *mat->k, false);
 
-                    // assembly
-                    tbb::spin_mutex::scoped_lock lock(mutex);
-                    for (size_t ii = 0; ii < e->nodes.size(); ++ii)
+                // assembly
+                tbb::spin_mutex::scoped_lock lock(mutex);
+                for (size_t ii = 0; ii < e->nodes.size(); ++ii)
+                {
+                    Node *nodi = e->nodes[ii];
+                    for (size_t jj = 0; jj < e->nodes.size(); ++jj)
                     {
-                        Node *nodi = e->nodes[ii];
-                        for (size_t jj = 0; jj < e->nodes.size(); ++jj)
-                        {
-                            Node *nodj = e->nodes[jj];
-                            T.push_back(Eigen::Triplet<double>(rows[nodi->row], nodj->row, Ke(ii, jj)));
-                        }
+                        Node *nodj = e->nodes[jj];
+                        T.push_back(Eigen::Triplet<double>(rows[nodi->row], nodj->row, Ke(ii, jj)));
                     }
+                }
             });
         }
     }