void add()

in fatal/math/statistical_moments.h [79:99]


  void add(value_type const &sample) {
    auto const n = samples_++;

    auto const delta = sample - moment_1_;
    auto const normalized_delta = delta / samples_;
    auto const normalized_delta_squared = normalized_delta * normalized_delta;

    auto const term1 = delta * normalized_delta * n;

    moment_4_ += term1 * normalized_delta_squared
      * (samples_ * samples_ - 3 * samples_ + 3)
      + 6 * normalized_delta_squared * moment_2_
      - 4 * normalized_delta * moment_3_;

    moment_3_ += term1 * normalized_delta * (n - 1)
      - 3 * normalized_delta * moment_2_;

    moment_2_ += term1;

    moment_1_ += normalized_delta;
  }