rl_algo.hpp

template <typename OperatorType,
typename SignalType>
unsigned int rl(OperatorType &op, SignalType &gn, SignalType &bg, SignalType &x, unsigned int maxit) {
  unsigned int it = 0;

  SignalType temp = x;
  SignalType temp2 = x;

  while (it < maxit) {
    op.blur(x, temp);
    temp += bg;
    SignalType::transform(gn.begin(), gn.end(), temp.begin(), temp2.begin(), typename SignalType::div_op());
    op.blurT(temp2, temp);
    x *= temp;

    ++it;
  }
  return it;
}