series.hh
No OneTemporary
Actions

Subscribers

None

File Metadata

	#ifndef SERIES_HPP
	#define SERIES_HPP

	#include <cmath>

	class Series {
	public:
	unsigned int current_index = 0;
	double current_value = 0.;
	virtual double compute(unsigned int N) {
	if (current_index <= N) {
	N = N - current_index;
	}
	else { current_value = 0.;
	current_index = 0;
	}
	for (int k = 0;k<N; ++k){
	addIndex();
	}
	return current_value;
	}

	virtual double getAnalyticPrediction() {
	return std::nan("");
	}

	void addIndex() {
	current_index = current_index + 1;
	current_value = current_value + computeValue(current_index);
	}

	virtual double computeValue(int k) = 0;
	};

	#endif