21 std::vector<PairTimeParameterName>
22 time_parameter_name_mapping)
24 _time_parameter_name_mapping(std::move(time_parameter_name_mapping))
44 "Coordinate system not expected to be set for curve scaled "
65 auto const alpha = (t - t0) / (t1 - t0);
68 std::transform(r0.begin(), r0.end(), r0.begin(),
69 [
alpha](
auto const& v) { return (1 - alpha) * v; });
71 std::transform(r1.begin(), r1.end(), r1.begin(),
72 [
alpha](
auto const& v) { return alpha * v; });
73 std::transform(r0.begin(), r0.end(), r1.begin(), r0.begin(),
74 [](
auto const& v0,
auto const& v1) { return v0 + v1; });
79 std::vector<std::unique_ptr<ParameterBase>>
const& parameters)
81 DBUG(
"TimeDependentHeterogeneousParameter init {:d} time series entries.",
86 &findParameter<double>(time_parameter_map.second, parameters, 0);
97 { return n == p.second->getNumberOfGlobalComponents(); }))
100 "All referenced parameters in time dependent heterogeneous "
101 "parameter '{:s}' have to have the same number of components.",
111 auto const time_series_config =
115 std::vector<TimeDependentHeterogeneousParameter::PairTimeParameterName>
118 for (
auto const p : time_series_config.getConfigSubtreeList(
"pair"))
121 auto time = p.getConfigParameter<
double>(
"time");
122 auto parameter_name =
124 p.getConfigParameter<std::string>(
"parameter_name");
125 time_series.emplace_back(time, parameter_name);
128 if (time_series.empty())
131 "Time dependent heterogeneous parameter '{:s}' doesn't contain "
132 "necessary time series data.",
137 time_series.begin(), time_series.end(),
141 p1) { return p0.first < p1.first; }))
144 "The points in time in the time series '{:s}' aren't in ascending "
149 return std::make_unique<TimeDependentHeterogeneousParameter>(
150 name, std::move(time_series));
void DBUG(char const *fmt, Args const &... args)
void checkConfigParameter(std::string const ¶m, T const &value) const
ConfigTree getConfigSubtree(std::string const &root) const
std::vector< double > operator()(double const t, SpatialPosition const &pos) const override
Returns the parameter value at the given time and position.
TimeDependentHeterogeneousParameter(std::string name, std::vector< PairTimeParameterName > time_parameter_name_mapping)
std::vector< PairTimeParameter > _time_parameter_mapping
void initialize(std::vector< std::unique_ptr< ParameterBase >> const ¶meters) override
bool isTimeDependent() const override
std::pair< double, std::string > PairTimeParameterName
int getNumberOfGlobalComponents() const override
std::vector< PairTimeParameterName > _time_parameter_name_mapping
std::unique_ptr< ParameterBase > createTimeDependentHeterogeneousParameter(std::string const &name, BaseLib::ConfigTree const &config)
std::optional< CoordinateSystem > _coordinate_system