From 64bb04964733535b9566ef883a5eaf9ac689f7a3 Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Fri, 1 Jul 2016 22:37:19 -0400 Subject: [PATCH] [Numerics] Handle zero-size timestep in CVodesIntegrator::advance This permits simpler integration loops starting from t=0 --- src/numerics/CVodesIntegrator.cpp | 3 +++ 1 file changed, 3 insertions(+) diff --git a/src/numerics/CVodesIntegrator.cpp b/src/numerics/CVodesIntegrator.cpp index 66cf616f8..dcc6b9a5f 100644 --- a/src/numerics/CVodesIntegrator.cpp +++ b/src/numerics/CVodesIntegrator.cpp @@ -386,6 +386,9 @@ void CVodesIntegrator::applyOptions() void CVodesIntegrator::integrate(double tout) { + if (tout == m_time) { + return; + } int flag = CVode(m_cvode_mem, tout, m_y, &m_time, CV_NORMAL); if (flag != CV_SUCCESS) { throw CanteraError("CVodesIntegrator::integrate",