diff --git a/R/fit_task0.R b/R/fit_task0.R
index 72f1b14..9a15ccd 100644
--- a/R/fit_task0.R
+++ b/R/fit_task0.R
@@ -161,7 +161,7 @@ find_breakpoints = function(d, avg_points_per_window, max_breakpoints, norm, n_t
 
   max_breakpoints = min(max_breakpoints, as.integer(length(x) / avg_points_per_window))
   available_breakpoints <- 1:max_breakpoints
-  message("Initial proposals")
+  message("Intializing breakpoints")
   proposed_breakpoints <- parallel::mclapply(available_breakpoints, FUN = function(n_breakpoints) {
     print(n_breakpoints)
     convergence <<- FALSE
@@ -200,11 +200,12 @@ find_breakpoints = function(d, avg_points_per_window, max_breakpoints, norm, n_t
   }
 
   m <- biPOD:::get_model("fit_breakpoints")
-  message("Proposals' optimization")
+  message("Breakpoints optimization")
   fits <- list()
   plots <- list()
   proposed_breakpoints$idx <- c(1:nrow(proposed_breakpoints)) + 1
 
+
   loos <- lapply(0:nrow(proposed_breakpoints), function(j) {
     #print(j)
     if (j == 0) {
@@ -268,6 +269,8 @@ find_breakpoints = function(d, avg_points_per_window, max_breakpoints, norm, n_t
     return(list(best_bp=NULL, best_fit=NULL))
   }
 
+  message("Choosing optimal breakpoints")
+
   bic_comp <- dplyr::tibble(bic = loos, n_breakpoints = 1:length(loos) - 1)
   best_j <- bic_comp %>%
     dplyr::filter(bic == min(bic)) %>%