18 std::set<std::string>&& output_variables_,
19 std::vector<double>&& fixed_output_times_,
20 std::vector<PairRepeatEachSteps>&& repeats_each_steps_,
21 bool const output_residuals_)
22 : output_variables(std::move(output_variables_)),
23 fixed_output_times(std::move(fixed_output_times_)),
24 repeats_each_steps(std::move(repeats_each_steps_)),
25 output_residuals(output_residuals_)
30 "Vector of fixed output time steps passed to the "
31 "OutputDataSpecification constructor must be sorted");
36 if (pair.each_steps == 0)
39 "Step in pair of <repeats><steps> is zero but has to be "
40 "greater than zero.");
50 double const time)
const
52 auto isFixedOutputStep = [
this](
double const time) ->
bool
55 [&](
auto fixed_output_time)
57 return (std::abs(fixed_output_time - time) <
62 auto isPairRepeatsEachTimeStepOutput = [
this](
int timestep) ->
bool
68 each_steps = pair.each_steps;
70 if (timestep > pair.repeat * each_steps)
72 timestep -= pair.repeat * each_steps;
80 return timestep % each_steps == 0;
83 return isFixedOutputStep(time) || isPairRepeatsEachTimeStepOutput(timestep);
95 os <<
"OuputDataSpecification" << std::endl;
96 os <<
"\toutput_variables: ";
98 std::ostream_iterator<std::string>(os,
" "));
100 os <<
"\tfixed_output_times: ";
102 std::ostream_iterator<double>(os,
" "));
104 os <<
"\trepeats_each_steps: ";
106 std::ostream_iterator<PairRepeatEachSteps>(os,
" "));
static constexpr double minimalTimeStepSize
std::ostream & operator<<(std::ostream &os, Output const &output)
Holds information about which variables to write to output files.
std::vector< PairRepeatEachSteps > repeats_each_steps
Describes after which timesteps to write output.
std::vector< double > fixed_output_times
Given times that steps have to reach.
bool output_residuals
Tells if also to output extrapolation residuals.
std::set< std::string > output_variables
All variables that shall be output.
bool isOutputStep(int timestep, double const time) const
OutputDataSpecification(std::set< std::string > &&output_variables_, std::vector< double > &&fixed_output_times_, std::vector< PairRepeatEachSteps > &&repeats_each_steps_, bool const output_residuals_)
const int repeat
Apply each_steps repeat times.
const int each_steps
Do output every each_steps timestep.