From 1b7c9b6e2eb89e2d17fd0b76d3e3f51308dbc37a Mon Sep 17 00:00:00 2001
From: Florian Jung <flo@windfisch.org>
Date: Thu, 8 Jan 2015 00:40:23 +0100
Subject: panorama!

---
 ringbuf.h | 126 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 126 insertions(+)
 create mode 100644 ringbuf.h

(limited to 'ringbuf.h')

diff --git a/ringbuf.h b/ringbuf.h
new file mode 100644
index 0000000..5756153
--- /dev/null
+++ b/ringbuf.h
@@ -0,0 +1,126 @@
+#include <math.h>
+
+class Ringbuffer
+{
+	public:
+
+	Ringbuffer(int size)
+	{
+		this->size = size;
+		idx = 0;
+		buf = new double[size];
+		for (int i=0; i<size; i++)
+			buf[i] = 0;
+		avg = 0.0;
+		avg_valid = false;
+	}
+
+	~Ringbuffer()
+	{
+		delete [] buf;
+	}
+
+	double get()
+	{
+		if (!avg_valid)
+		{
+			avg=0.0;
+			for (int i=0; i<size; i++)
+				avg += buf[i];
+			avg/=size;
+			avg_valid=true;
+		}
+
+		return avg;
+	}
+
+	void put(double val)
+	{
+		buf[idx] = val;
+		idx = (idx+1) % size;
+		avg_valid = false;
+	}
+
+	void set(double val)
+	{
+		for (int i=0; i<size; i++)
+			buf[i]=val;
+		avg = val;
+		avg_valid=true;
+	}
+
+	void add(double val)
+	{
+		for (int i=0; i<size; i++)
+			buf[i]+=val;
+		avg += val;
+	}
+
+	private:
+		double* buf;
+		int idx;
+		int size;
+		double avg;
+		bool avg_valid;
+};
+
+class ModuloRingbuffer
+{
+	private:
+		Ringbuffer* rb;
+		double low,upp,span;
+
+	public:
+	
+	ModuloRingbuffer(int size, double low, double upp)
+	{
+		rb = new Ringbuffer(size);
+		this->low=low;
+		this->upp=upp;
+		this->span=upp-low;
+		rb->set(_fixup_range(0.0));
+	}
+
+	~ModuloRingbuffer()
+	{
+		delete rb;
+	}
+
+	double _fixup_range(double val)
+	{
+		while (val < low) val+=span;
+		while (val>= upp) val-=span;
+		return val;
+	}
+
+	double get()
+	{
+		return rb->get();
+	}
+
+	void put(double val)
+	{
+		val=_fixup_range(val);
+
+		// direct way
+		double dist1 = val - get();
+		// over borders
+		double dist2 = dist1 + span;
+		double dist3 = dist1 - span;
+
+		if (fabs(dist1) <= fabs(dist2) && fabs(dist1) <= fabs(dist3))
+			rb->put(val);
+		else if (fabs(dist2) <= fabs(dist3))
+			rb->put(val+span);
+		else
+			rb->put(val-span);
+
+		while (get() < low) rb->add(span);
+		while (get()>= upp) rb->add(-span);
+	}
+
+	void set(double val)
+	{
+		rb->set(_fixup_range(val));
+	}
+};
-- 
cgit v1.2.3