1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105
| #include <iostream>
#include <string>
#include <vector>
#include <algorithm> // std::transform
using namespace std;
class SpatioTemporalDataSet
{
public:
SpatioTemporalDataSet(const std::string & name, std::size_t nbDemes, std::size_t nbTimeSteps, std::vector<double>& data);
virtual ~SpatioTemporalDataSet();
typedef std::vector<double>::iterator iterator;
iterator begin() { return m_data.begin(); }
iterator end() { return m_data.end(); }
typedef std::vector<double>::const_iterator const_iterator;
const_iterator begin() const { return m_data.begin(); }
const_iterator end() const { return m_data.end(); }
// STL non constant iterator to iterate over a same temporal step
iterator begin(std::size_t time) { return m_data.begin() + (time)*m_nbrOfDemes ; }
iterator end(std::size_t time) { return begin(time + 1) ; }
// STL constant iterator to iterate over a same temporal step
const_iterator begin(std::size_t time) const { return m_data.begin() + (time)*m_nbrOfDemes ; }
const_iterator end(std::size_t time) const { return begin(time + 1) ; }
// Getters
int GetNumberOfDemes(void) const { return m_nbrOfDemes; }
int GetNumberOfTimeSteps(void) const { return m_nbrOfTimeSteps; }
std::string GetNameOfEnvironmentalVariable(void) const { return m_name; }
private:
std::vector<double> m_data;
std::string m_name;
std::size_t m_nbrOfDemes;
std::size_t m_nbrOfTimeSteps;
};
SpatioTemporalDataSet::SpatioTemporalDataSet(const std::string & name, std::size_t nbDemes, std::size_t nbTimeSteps, std::vector<double>& data) :
m_data(data),
m_name(name),
m_nbrOfDemes(nbDemes),
m_nbrOfTimeSteps(nbTimeSteps)
{}
SpatioTemporalDataSet::~SpatioTemporalDataSet(){}
double op_increase (double i) { return i += 1; }
int main()
{
vector<double> v = {1,2,3,4,5,6,7,8,9,10,11,12};
SpatioTemporalDataSet dataSet("pluviometry", 4, 3, v);
std::cout << "This is a spatio-temporal dataSet for "
<< dataSet.GetNameOfEnvironmentalVariable() << " with "
<< dataSet.GetNumberOfDemes() << " demes and "
<< dataSet.GetNumberOfTimeSteps() << " time steps "
<< std::endl;
// Print all values :
for (std::vector<double>::const_iterator i = dataSet.begin(); i != dataSet.end(); ++i)
std::cout << *i << ' ';
std::cout << std::endl;
for (auto i : dataSet) // works too
std::cout << i << ' ';
std::cout << std::endl;
// Access temporal steps one by one:
for (std::vector<double>::const_iterator i = dataSet.begin(0); i != dataSet.end(0); ++i) // and not 0 : need an assert
std::cout << *i << ' ';
std::cout << std::endl;
for (std::vector<double>::const_iterator i = dataSet.begin(1); i != dataSet.end(1); ++i)
std::cout << *i << ' ';
std::cout << std::endl;
for (std::vector<double>::const_iterator i = dataSet.begin(2); i != dataSet.end(2); ++i)
std::cout << *i << ' ';
std::cout << std::endl;
// Add 1 for all data
std::vector<double> bar;
int d = dataSet.GetNumberOfDemes();
int t = dataSet.GetNumberOfTimeSteps();
bar.resize(d*t); // allocate space
transform (dataSet.begin(), dataSet.end(), bar.begin(), op_increase);
for (auto i : bar) // works too
std::cout << i << ' ';
std::cout << std::endl;
cout << "Hello world!" << endl;
return 0;
} |
Partager