R/randmut2.R

In adana: Adaptive Nature-Inspired Algorithms for Hybrid Genetic Optimization

# random mutation 2
randmut2 = function(y, lb, ub, mutpm, ...){
  m = length(y)
  if(missing(lb) | missing(ub))
    stop("randmut2 requires lb and ub!")
  if(missing(mutpm)) mutpm=0.05
  for(j in 1:m){
    r = runif(1, 0, 1)
    val = rnorm(1, mean=0, sd=(ub[j]-lb[j])/10)
    val = ifelse(r < mutpm, val, 0)
    y[j] = y[j] + val
  }
  return(list(mutant=y))
}