set.seed(1)
data <- simulate_data_FRTM(n_obs=100)

dom <- c(0,1)
basis_x <- fda::create.bspline.basis(c(0,1),nbasis =30)
x_fd <- fda::smooth.basis(data$grid_i[[1]]/max(data$grid_i[[1]]),data$x_err[[1]],basis_x)$fd
template_fd <- fda::smooth.basis(data$grid_template,data$template,basis_x)$fd
der_x_fd <- fda::deriv.fd(x_fd,1)
der_template_fd <- fda::deriv.fd(template_fd,1)

mod <- OEBFDTW(x_fd,template_fd,der_x_fd,der_template_fd,get_fd = "x_reg")