diff --git a/include/evalhyd/evald.hpp b/include/evalhyd/evald.hpp
index 671d1a1bb292da5a15b4e5e45b93648b69ada5fd..ae72beab38001e33ed971acedb4e233812214192 100644
--- a/include/evalhyd/evald.hpp
+++ b/include/evalhyd/evald.hpp
@@ -3,7 +3,6 @@
 
 #include <unordered_map>
 #include <vector>
-#include <optional>
 
 #include <xtensor/xexpression.hpp>
 #include <xtensor/xtensor.hpp>
@@ -160,6 +159,7 @@ namespace evalhyd
             const std::vector<std::string>& dts = {}
     )
     {
+        // retrieve real types of the expressions
         const D2& q_obs_ = q_obs.derived_cast();
         const D2& q_prd_ = q_prd.derived_cast();