small optimizations

• replace pow(x, 2) by Optim::pow2(x)
  • implement running average on a vector (start_pos and len given)