diff --git a/R/time-warping.R b/R/time-warping.R
index baa93cd..956c821 100644
--- a/R/time-warping.R
+++ b/R/time-warping.R
@@ -24,7 +24,7 @@
 #' @param parallel A boolean specifying whether to run calculations in parallel.
 #'   Defaults to `FALSE`.
 #' @param optim_method A string specifying the algorithm used for optimization.
-#'   Choices are `"DP"`, `"DP2"` and `"RBFGS"`. Defaults to `"DP"`.
+#'   Choices are `"DP"`, `"DPo"`, `"DP2"`, and `"RBFGS"`. Defaults to `"DP"`.
 #' @param max_iter An integer value specifying the maximum number of iterations.
 #'   Defaults to `20L`.
 #'
@@ -80,7 +80,7 @@ time_warping <- function(f, time,
                          smooth_data = FALSE,
                          sparam = 25L,
                          parallel = FALSE,
-                         optim_method = c("DP", "DP2", "RBFGS"),
+                         optim_method = c("DP", "DPo", "DP2", "RBFGS"),
                          max_iter = 20L) {
   penalty_method <- rlang::arg_match(penalty_method)
   centroid_type <- rlang::arg_match(centroid_type)
diff --git a/man/time_warping.Rd b/man/time_warping.Rd
index 33df8ff..84210bd 100644
--- a/man/time_warping.Rd
+++ b/man/time_warping.Rd
@@ -14,7 +14,7 @@ time_warping(
   smooth_data = FALSE,
   sparam = 25L,
   parallel = FALSE,
-  optim_method = c("DP", "DP2", "RBFGS"),
+  optim_method = c("DP", "DPo", "DP2", "RBFGS"),
   max_iter = 20L
 )
 }
@@ -49,7 +49,7 @@ box filter. Defaults to \code{25L}. This is used only when \code{smooth_data = T
 Defaults to \code{FALSE}.}
 
 \item{optim_method}{A string specifying the algorithm used for optimization.
-Choices are \code{"DP"}, \code{"DP2"} and \code{"RBFGS"}. Defaults to \code{"DP"}.}
+Choices are \code{"DP"}, \code{"DPo"}, \code{"DP2"}, and \code{"RBFGS"}. Defaults to \code{"DP"}.}
 
 \item{max_iter}{An integer value specifying the maximum number of iterations.
 Defaults to \code{20L}.}