diff --git a/include/micm/solver/rosenbrock.inl b/include/micm/solver/rosenbrock.inl index b50314bee..bb2226819 100644 --- a/include/micm/solver/rosenbrock.inl +++ b/include/micm/solver/rosenbrock.inl @@ -21,7 +21,7 @@ namespace micm auto& K = derived_class_temporary_variables->K_; auto& Yerror = derived_class_temporary_variables->Yerror_; const double h_max = parameters_.h_max_ == 0.0 ? time_step : std::min(time_step, parameters_.h_max_); - const double h_start = + const double h_start = parameters_.h_start_ == 0.0 ? std::max(parameters_.h_min_, DELTA_MIN) : std::min(h_max, parameters_.h_start_); SolverStats stats;