[dt] Reduced the tolerance in the time step size calculation under the...

[dt] Reduced the tolerance in the time step size calculation under the constraint of specified times

NumLib/TimeStepping/Algorithms/EvolutionaryPIDcontroller.cpp

@@ -147,9 +147,8 @@ double EvolutionaryPIDcontroller::checkSpecificTimeReached(const double h_new)
         return h_new;
 
     const double specific_time = _specified_times.back();
-    const double zero_threshold = std::numeric_limits<double>::epsilon();
     if ((specific_time > _ts_current.current()) &&
-        (_ts_current.current() + h_new - specific_time > zero_threshold))
+        (_ts_current.current() + h_new - specific_time > 0.0))
     {
         _specified_times.pop_back();
         return specific_time - _ts_current.current();

ProcessLib/Output/Output.cpp

@@ -71,7 +71,7 @@ bool Output::shallDoOutput(unsigned timestep, double const t)
         return make_output;
 
     const double specific_time = _specified_times.back();
-    const double zero_threshold = std::numeric_limits<double>::epsilon();
+    const double zero_threshold = std::numeric_limits<double>::min();
     if (std::fabs(specific_time - t) < zero_threshold)
     {
         _specified_times.pop_back();
@@ -85,7 +85,7 @@ Output::Output(std::string output_directory, std::string prefix,
                bool const compress_output, std::string const& data_mode,
                bool const output_nonlinear_iteration_results,
                std::vector<PairRepeatEachSteps> repeats_each_steps,
-               const std::vector<double>&& specified_times)
+               std::vector<double>&& specified_times)
     : _output_directory(std::move(output_directory)),
       _output_file_prefix(std::move(prefix)),
       _output_file_compression(compress_output),

ProcessLib/Output/Output.h

@@ -43,7 +43,7 @@ public:
            bool const compress_output, std::string const& data_mode,
            bool const output_nonlinear_iteration_results,
            std::vector<PairRepeatEachSteps> repeats_each_steps,
-           const std::vector<double>&& specified_times);
+           std::vector<double>&& specified_times);
 
     //! TODO doc. Opens a PVD file for each process.
     void addProcess(ProcessLib::Process const& process, const int process_id);