R/powmut.R

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

# Power mutation
powmut = function(y, lb, ub, mutpow, ...){
  m = length(y)
  if(missing(lb) | missing(ub)) 
    stop("powmut requires lb and ub !")
  if(missing(mutpow)) mutpow = 2
  v = sample(1:m, 1)
  r = runif(1)
  p = runif(1)^mutpow
  t = (y[v]-lb[v])/(ub[v]-y[v])
  y[v] = y[v] + ifelse(r>t, -p*(y[v]-lb[v]), +p*(ub[v]-y[v]))
  return(list(mutant=y, mutgen=v))
}