diff --git a/src/NumericalInitialization.cc b/src/NumericalInitialization.cc
index e004adc402caa63b31525903f9ad7eace9fdb4df..f9b496a67f5e3909d77f2bfb65effc07475f6afe 100644
--- a/src/NumericalInitialization.cc
+++ b/src/NumericalInitialization.cc
@@ -619,8 +619,8 @@ void
 LoadParamsAndSteadyStateStatement::fillEvalContext(eval_context_t &eval_context) const
 {
   for (const auto & it : content)
-    /* We use strtod() instead of stod() because we want overflows and
-       underflows to respectively yield 0 and ±Inf. See also the comment in
+    /* We use strtod() instead of stod() because we want underflows and
+       overflows to respectively yield 0 and ±Inf. See also the comment in
        NumericalConstants.cc */
     eval_context[it.first] = strtod(it.second.c_str(), nullptr);
 }