diff --git a/include/teqp/algorithms/VLE.hpp b/include/teqp/algorithms/VLE.hpp
index 9a2ada0f7bbf99eccfcf7af58e7a4507abc050aa..38b36c105840da90acd96b53e271543cdd3b5c42 100644
--- a/include/teqp/algorithms/VLE.hpp
+++ b/include/teqp/algorithms/VLE.hpp
@@ -679,7 +679,7 @@ auto trace_VLE_isotherm_binary(const Model &model, Scalar T, VecType rhovecL0, V
             store_point();
         }
 
-        double dtold = dt;
+        //double dtold = dt;
         auto x0_previous = x0;
 
         if (opt.integration_order == 5) {
diff --git a/include/teqp/algorithms/VLLE.hpp b/include/teqp/algorithms/VLLE.hpp
index 01df6601830d551bf8279b059d6cf95979b27f98..b546789cfac6246c6564218bbadef786508b4975 100644
--- a/include/teqp/algorithms/VLLE.hpp
+++ b/include/teqp/algorithms/VLLE.hpp
@@ -174,7 +174,7 @@ namespace teqp {
                 y.push_back(p);
             }
             auto intersections = get_self_intersections(x, y);
-            auto& trace = traces[0];
+            //auto& trace = traces[0];
 
             auto process_intersection = [&](auto& trace, auto& i) {
                 auto rhoL1_j = traces[0][i.j].at("rhoL / mol/m^3").template get<std::valarray<double>>();