/*
 *
 * (C) 2013-25 - ntop.org
 *
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software Foundation,
 * Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
 *
 */

#ifndef _DES_COUNTER_H_
#define _DES_COUNTER_H_

#include "ntop_includes.h"

/* ************************ */

/* Counter based on Double Exponential smoothing algorithm */
class DESCounter : public BehaviouralCounter {
 private:
  struct ndpi_des_struct des;

 public:
  DESCounter(double alpha = 0.9, double beta = 0.035, float significance = 0.05)
      : BehaviouralCounter() {
    if (ndpi_des_init(&des, alpha, beta, significance) != 0)
      throw "Error while creating DES";
  }

  DESCounter(struct ndpi_des_struct *_des) {
    memcpy(&des, _des, sizeof(struct ndpi_des_struct));
  }

  bool addObservation(u_int64_t value) {
    double forecast, confidence_band;
    bool rc =
        (ndpi_des_add_value(&des, value, &forecast, &confidence_band) == 1)
            ? true
            : false;
    double l_forecast = forecast - confidence_band;
    double h_forecast = forecast + confidence_band;

    last_value = value;
    last_lower = (u_int64_t)floor(((l_forecast < 0) ? 0 : l_forecast));
    last_upper = (u_int64_t)round(h_forecast + 0.5);

    if (rc) {
      is_anomaly =
          ((value < last_lower) || (value > last_upper)) ? true : false;

      if (is_anomaly) tot_num_anomalies++;

      return (is_anomaly);
    }

    return (rc);
  }

  inline void resetStats() { ndpi_des_reset(&des); }
  inline DESCounter *clone() { return (new (std::nothrow) DESCounter(&des)); }
  inline void set(DESCounter *c) {
    memcpy(&des, &c->des, sizeof(struct ndpi_des_struct));
  }
};

#endif /* _DES_COUNTER_H_ */
