diff --git a/NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.cpp b/NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.cpp
index b9ec34a2b6622d3b5fe9d6524464dbd989cb6316..1f8a692ac896a51a3b3cd64149557bf2345ef28c 100644
--- a/NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.cpp
+++ b/NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.cpp
@@ -56,7 +56,8 @@ bool EvolutionaryPIDcontroller::next(double const solution_error,
: 0.5 * _ts_current.dt();
h_new = limitStepSize(h_new, is_previous_step_accepted);
- h_new = possiblyClampToNextFixedTime(h_new);
+ h_new = possiblyClampDtToNextFixedTime(_ts_current.current(), h_new,
+ _fixed_output_times);
_ts_current = _ts_prev;
_ts_current += h_new;
@@ -116,7 +117,8 @@ bool EvolutionaryPIDcontroller::next(double const solution_error,
}
h_new = limitStepSize(h_new, is_previous_step_accepted);
- h_new = possiblyClampToNextFixedTime(h_new);
+ h_new = possiblyClampDtToNextFixedTime(_ts_current.current(), h_new,
+ _fixed_output_times);
_dt_vector.push_back(h_new);
_ts_prev = _ts_current;
@@ -165,27 +167,6 @@ double EvolutionaryPIDcontroller::limitStepSize(
return limited_h;
}
-double EvolutionaryPIDcontroller::possiblyClampToNextFixedTime(
- const double h_new) const
-{
- auto const specific_time =
- std::upper_bound(std::cbegin(_fixed_output_times),
- std::cend(_fixed_output_times), _ts_current.current());
-
- if (specific_time == std::cend(_fixed_output_times))
- {
- return h_new;
- }
-
- if ((*specific_time > _ts_current.current()) &&
- (_ts_current.current() + h_new - *specific_time > 0.0))
- {
- return *specific_time - _ts_current.current();
- }
-
- return h_new;
-}
-
void EvolutionaryPIDcontroller::addFixedOutputTimes(
std::vector<double> const& extra_fixed_output_times)
{
diff --git a/NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.h b/NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.h
index 069daca9aab88a8c5416ad041d521bd8d02f4af7..af5a598e642cd72eb6fa23e7c7456a10fdf0af4f 100644
--- a/NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.h
+++ b/NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.h
@@ -111,10 +111,6 @@ private:
*/
double limitStepSize(const double h_new,
const bool previous_step_accepted) const;
-
- /// If any time will be reached with given time increment, it will be
- /// reduced, otherwise the input will be returned.
- double possiblyClampToNextFixedTime(const double h_new) const;
};
/// Create an EvolutionaryPIDcontroller time stepper from the given
diff --git a/NumLib/TimeStepping/Algorithms/IterationNumberBasedTimeStepping.cpp b/NumLib/TimeStepping/Algorithms/IterationNumberBasedTimeStepping.cpp
index b240bca0a67fb2db8a7050de71670b90533a6d07..a204a053b02c3621d184a1adb338d8d7c1e66a7b 100644
--- a/NumLib/TimeStepping/Algorithms/IterationNumberBasedTimeStepping.cpp
+++ b/NumLib/TimeStepping/Algorithms/IterationNumberBasedTimeStepping.cpp
@@ -79,7 +79,8 @@ bool IterationNumberBasedTimeStepping::next(double const /*solution_error*/,
// prepare the next time step info
_ts_current = _ts_prev;
- _ts_current += possiblyClampToNextFixedTime(getNextTimeStepSize());
+ _ts_current += possiblyClampDtToNextFixedTime(
+ _ts_current.current(), getNextTimeStepSize(), _fixed_output_times);
return true;
}
@@ -133,27 +134,6 @@ double IterationNumberBasedTimeStepping::getNextTimeStepSize() const
return dt;
}
-double IterationNumberBasedTimeStepping::possiblyClampToNextFixedTime(
- const double h_new) const
-{
- auto const specific_time =
- std::upper_bound(std::cbegin(_fixed_output_times),
- std::cend(_fixed_output_times), _ts_current.current());
-
- if (specific_time == std::cend(_fixed_output_times))
- {
- return h_new;
- }
-
- if ((*specific_time > _ts_current.current()) &&
- (_ts_current.current() + h_new - *specific_time > 0.0))
- {
- return *specific_time - _ts_current.current();
- }
-
- return h_new;
-}
-
void IterationNumberBasedTimeStepping::addFixedOutputTimes(
std::vector<double> const& extra_fixed_output_times)
{
diff --git a/NumLib/TimeStepping/Algorithms/IterationNumberBasedTimeStepping.h b/NumLib/TimeStepping/Algorithms/IterationNumberBasedTimeStepping.h
index f911e773d12a6579f2b7156696c1597bf1d12060..9dcf45eacbc2e783bd0ec650f69dab6cd246e144 100644
--- a/NumLib/TimeStepping/Algorithms/IterationNumberBasedTimeStepping.h
+++ b/NumLib/TimeStepping/Algorithms/IterationNumberBasedTimeStepping.h
@@ -116,10 +116,6 @@ private:
/// Find a multiplier for the given number of iterations.
double findMultiplier(int number_iterations) const;
- /// If any time will be reached with given time increment, it will be
- /// reduced, otherwise the input will be returned.
- double possiblyClampToNextFixedTime(const double h_new) const;
-
/// This vector stores the number of iterations to which the respective
/// multiplier coefficient will be applied.
const std::vector<int> _iter_times_vector;
diff --git a/NumLib/TimeStepping/Algorithms/TimeStepAlgorithm.cpp b/NumLib/TimeStepping/Algorithms/TimeStepAlgorithm.cpp
new file mode 100644
index 0000000000000000000000000000000000000000..0f627b3bfc035153f41036d6279366a187e2cf60
--- /dev/null
+++ b/NumLib/TimeStepping/Algorithms/TimeStepAlgorithm.cpp
@@ -0,0 +1,35 @@
+/**
+ * \file
+ * \copyright
+ * Copyright (c) 2012-2020, OpenGeoSys Community (http://www.opengeosys.org)
+ * Distributed under a Modified BSD License.
+ * See accompanying file LICENSE.txt or
+ * http://www.opengeosys.org/project/license
+ */
+
+#include "TimeStepAlgorithm.h"
+
+#include <algorithm>
+
+namespace NumLib
+{
+double possiblyClampDtToNextFixedTime(
+ double const t, double const dt,
+ std::vector<double> const& fixed_output_times)
+{
+ auto const specific_time = std::upper_bound(
+ std::cbegin(fixed_output_times), std::cend(fixed_output_times), t);
+
+ if (specific_time == std::cend(fixed_output_times))
+ {
+ return dt;
+ }
+
+ if ((*specific_time > t) && (t + dt - *specific_time > 0.0))
+ {
+ return *specific_time - t;
+ }
+
+ return dt;
+}
+} // namespace NumLib
diff --git a/NumLib/TimeStepping/Algorithms/TimeStepAlgorithm.h b/NumLib/TimeStepping/Algorithms/TimeStepAlgorithm.h
index d0285855b379bf1410dc7448b6f3df3886b119cd..faf0a559f57379d1a7729ebd6300b5442d0b3b8f 100644
--- a/NumLib/TimeStepping/Algorithms/TimeStepAlgorithm.h
+++ b/NumLib/TimeStepping/Algorithms/TimeStepAlgorithm.h
@@ -141,4 +141,14 @@ protected:
std::vector<double> _dt_vector;
};
+/// If any of the fixed times will be reached with given time increment, it will
+/// be reduced, otherwise the input will be returned.
+/// \pre The input vector of fixed times must be sorted.
+/// \param t Current time.
+/// \param dt Suggested time increment.
+/// \param fixed_output_times Sorted list of times which are to be reached.
+double possiblyClampDtToNextFixedTime(
+ double const t, double const dt,
+ std::vector<double> const& fixed_output_times);
+
} // namespace NumLib