pub type TimeFunction<F> = Arc<dyn Fn(F) -> Array2<F> + Send + Sync>;
Type alias for time-dependent matrix function
pub struct TimeFunction<F> { /* private fields */ }