TimeSeries.h

#ifndef _TIMESERIES_H_
#define _TIMESERIES_H_
#include "mspass/seismic/CoreTimeSeries.h"
#include "mspass/utility/ErrorLogger.h"
#include "mspass/utility/ProcessingHistory.h"
namespace mspass::seismic {
class TimeSeries : public mspass::seismic::CoreTimeSeries,
                   public mspass::utility::ProcessingHistory {
public:
  TimeSeries()
      : mspass::seismic::CoreTimeSeries(),
        mspass::utility::ProcessingHistory() {};
  TimeSeries(const size_t nsamples)
      : mspass::seismic::CoreTimeSeries(nsamples),
        mspass::utility::ProcessingHistory() {};
  TimeSeries(const BasicTimeSeries &bts, const Metadata &md)
      : mspass::seismic::CoreTimeSeries(bts, md),
        mspass::utility::ProcessingHistory() {};
  TimeSeries(const Metadata &md);
 
  TimeSeries(const mspass::seismic::CoreTimeSeries &d)
      : mspass::seismic::CoreTimeSeries(d),
        mspass::utility::ProcessingHistory() {};
  TimeSeries(const mspass::seismic::CoreTimeSeries &d, const std::string alg);
  TimeSeries(const mspass::seismic::BasicTimeSeries &b,
             const mspass::utility::Metadata &m,
             const mspass::utility::ProcessingHistory &mcts,
             const std::vector<double> &d);
  TimeSeries(const TimeSeries &parent)
      : mspass::seismic::CoreTimeSeries(parent),
        mspass::utility::ProcessingHistory(parent) {};
  TimeSeries &operator=(const TimeSeries &parent);
  TimeSeries &operator+=(const TimeSeries &d) {
    dynamic_cast<CoreTimeSeries &>(*this) +=
        dynamic_cast<const CoreTimeSeries &>(d);
    this->elog += d.elog;
    return (*this);
  };
  TimeSeries &operator+=(const CoreTimeSeries &d) {
    dynamic_cast<CoreTimeSeries &>(*this) += d;
    return (*this);
  };
  const TimeSeries operator+(const TimeSeries &other) const {
    TimeSeries result(*this);
    result += other;
    return result;
  };
  const CoreTimeSeries operator+(const CoreTimeSeries &other) const {
    CoreTimeSeries result(dynamic_cast<const CoreTimeSeries &>(*this));
    result += other;
    return result;
  };
  TimeSeries &operator*=(const double scale) {
    dynamic_cast<CoreTimeSeries &>(*this) *= scale;
    return *this;
  };
  TimeSeries &operator-=(const TimeSeries &d) {
    dynamic_cast<CoreTimeSeries &>(*this) -=
        dynamic_cast<const CoreTimeSeries &>(d);
    this->elog += d.elog;
    return (*this);
  };
  TimeSeries &operator-=(const CoreTimeSeries &d) {
    dynamic_cast<CoreTimeSeries &>(*this) -= d;
    return (*this);
  };
  const TimeSeries operator-(const TimeSeries &other) const {
    TimeSeries result(*this);
    result -= other;
    return result;
  };
  const CoreTimeSeries operator-(const CoreTimeSeries &other) const {
    CoreTimeSeries result(dynamic_cast<const CoreTimeSeries &>(*this));
    result -= other;
    return result;
  };
  void load_history(const mspass::utility::ProcessingHistory &h);
  size_t memory_use() const;
};
} // namespace mspass::seismic
#endif