Merge pull request #121 from svkucheryavski/revert-119-unimod

Revert "fix bugs in constraintUnimod"

R/constraints.R

@@ -132,18 +132,13 @@ constraintUnimod <- function(x, d, tol = 0) {
    # process each component separately
    for (a in seq_len(ncol(x))) {
 
-      if (peak.ind[a] != 1) {
-          # flatten peaks to the left of maximum
-          left_part <- (peak.ind[a] - 1):1
-          x[, a] <- f(x[, a], max = x[peak.ind[a], a], indseq = left_part, step = +1)
-      }
-
-      if (peak.ind[a] != nvar) {
-          # flatten peaks to the right of maximum
-          right_part <- (peak.ind[a] + 1):nvar
-          x[, a] <- f(x[, a], max = x[peak.ind[a], a], indseq = right_part, step = -1)
-      }
+      # flatten peaks to the left of maximum
+      left_part <- (peak.ind[a] - 1):1
+      x[, a] <- f(x[, a], max = x[peak.ind[a], a], indseq = left_part, step = +1)
 
+      # flatten peaks to the right of maximum
+      right_part <- (peak.ind[a] + 1):nvar
+      x[, a] <- f(x[, a], max = x[peak.ind[a], a], indseq = right_part, step = -1)
    }
 
    return(x)