double ewma(double avg, double sample, double weight)
{
   return avg==0 ? sample : weight*avg+(1.0-weight)*sample;
}