TimeSeriesData.cxx

Go to the documentation of this file.

// C++ includes
#include <algorithm>
#include <ctime>
#include <iomanip>
#include <iostream>
#include <string>
#include <map>
#include <vector>

// NOvASoft includes
#include "TimeSeriesData.h"

namespace nl {
    
///////////////////////////////////////////////////////////////////////////
TimeSeriesData::TimeSeriesData(std::vector<std::time_t> t, std::vector<double> x)
{
    // Check that the two inputs are the same size
    if (t.size() != x.size()) {
        std::cerr << "ERROR. Vectors provided are different sizes. Initializing this object's properties to 0." << std::endl;
        this->T = {0};
        this->X = {0};
    } else {
        this->T = t;
        this->X = x;
    }
}

///////////////////////////////////////////////////////////////////////////
TimeSeriesData::~TimeSeriesData() {}
    

///////////////////////////////////////////////////////////////////////////
unsigned int TimeSeriesData::N() const
{
    return this->T.size();
}

///////////////////////////////////////////////////////////////////////////
std::time_t TimeSeriesData::MinTime() const
{
    return *std::min_element(this->T.begin(), this->T.end());
}

///////////////////////////////////////////////////////////////////////////
std::time_t TimeSeriesData::MaxTime() const
{
    return *std::max_element(this->T.begin(), this->T.end());
}

///////////////////////////////////////////////////////////////////////////
int TimeSeriesData::RangeTime() const
{
    time_t t1 = this->MinTime();
    time_t t2 = this->MaxTime();
    
    return t2-t1;
}

///////////////////////////////////////////////////////////////////////////
double TimeSeriesData::MinValue() const
{
    return *std::min_element(this->X.begin(), this->X.end());
}

///////////////////////////////////////////////////////////////////////////
double TimeSeriesData::MaxValue() const
{
    return *std::max_element(this->X.begin(), this->X.end());
}

///////////////////////////////////////////////////////////////////////////
double TimeSeriesData::MeanValue() const
{
    auto meanLambda = [&](double a, double b){return a + b / this->X.size();};
    return std::accumulate(this->X.begin(), this->X.end(), 0.0, meanLambda);
}

///////////////////////////////////////////////////////////////////////////
double TimeSeriesData::RangeValue() const
{
    double x1 = this->MinValue();
    double x2 = this->MaxValue();
    
    return x2-x1;
}

///////////////////////////////////////////////////////////////////////////
std::ostream& operator<< (std::ostream& out, const TimeSeriesData& ts)
{
    out << std::setiosflags(std::ios::fixed) << std::setprecision(2);
    out << "Time series data at " << ts << " contains " << ts.N() << " entries.";
    out << "Time span: " << ts.MinTime() << " -- " << ts.MaxTime();
    out << "Values: MIN=" << ts.MinValue() << ", MAX=" << ts.MaxValue() << ", MEAN=" << ts.MeanValue();

    return out;
}

} //end namespace nl