find_harmonics(y, m)

initial_parameters(k_vector, use_trend, use_damped_trend, ar_coeffs, ma_coeffs)

makeXMatrix(b, s_vector, d_vector, epsilon_vector)

findPQ(ar_coeffs, ma_coeffs)

makeTBATSWMatrix(phi, k_vector, ar_coeffs, ma_coeffs, tau)

makeTBATSGMatrix(
    k_vector,
    alpha,
    adj_beta,
    beta,
    gamma_one_v,
    gamma_two_v,
    p,
    q,
    tau
)

extract_params(
    params,
    use_boxcox,
    use_trend,
    use_damped_trend,
    use_arma_errors,
    seasonal_periods,
    p,
    q
)

updateTBATSWMatrix(w_transpose, phi, tau, ar_coeffs, ma_coeffs, p, q)

updateTBATSGMatrix(
    g,
    gamma_bold,
    alpha,
    beta,
    k_vector,
    gamma_one_v,
    gamma_two_v
)

updateTBATSFMatrix(
    F,
    phi,
    alpha,
    beta,
    gamma_bold,
    ar_coeffs,
    ma_coeffs,
    p,
    q,
    tau
)

checkAdmissibility(
    BoxCox_lambda,
    bc_lower_bound,
    bc_upper_bound,
    phi,
    ar_coeffs,
    ma_coeffs,
    D
)

calcLikelihoodTBATS(
    params,
    use_boxcox,
    use_trend,
    use_damped_trend,
    use_arma_errors,
    y,
    y_trans,
    seasonal_periods,
    k_vector,
    tau,
    w_transpose,
    F,
    g,
    gamma_bold,
    x_nought,
    x_nought_untransformed,
    bc_lower_bound,
    bc_upper_bound,
    p,
    q,
    scale
)

tbats_model_generator(
    y,
    seasonal_periods,
    k_vector,
    use_boxcox,
    bc_lower_bound,
    bc_upper_bound,
    use_trend,
    use_damped_trend,
    use_arma_errors,
    ar_coeffs,
    ma_coeffs
)

tbats_model(
    y,
    seasonal_periods,
    k_vector,
    use_boxcox,
    bc_lower_bound,
    bc_upper_bound,
    use_trend,
    use_damped_trend,
    use_arma_errors
)

tbats_selection(
    y,
    seasonal_periods,
    use_boxcox,
    bc_lower_bound,
    bc_upper_bound,
    use_trend,
    use_damped_trend,
    use_arma_errors
)

tbats_forecast(mod, h)