From 6963657b3f3ee4321394953a2cc67cd7386cce2d Mon Sep 17 00:00:00 2001
From: Tim Blechmann <timblech@users.sourceforge.net>
Date: Mon, 27 Dec 2004 22:55:41 +0000
Subject: additions and better code reuse

svn path=/trunk/externals/tb/; revision=2433
---
 chaos/src/coupled_logistic.hpp | 82 ++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 82 insertions(+)
 create mode 100644 chaos/src/coupled_logistic.hpp

(limited to 'chaos/src/coupled_logistic.hpp')

diff --git a/chaos/src/coupled_logistic.hpp b/chaos/src/coupled_logistic.hpp
new file mode 100644
index 0000000..eb5b326
--- /dev/null
+++ b/chaos/src/coupled_logistic.hpp
@@ -0,0 +1,82 @@
+// 
+//  
+//  chaos~
+//  Copyright (C) 2004  Tim Blechmann
+//  
+//  This program is free software; you can redistribute it and/or modify
+//  it under the terms of the GNU General Public License as published by
+//  the Free Software Foundation; either version 2 of the License, or
+//  (at your option) any later version.
+//  
+//  This program is distributed in the hope that it will be useful,
+//  but WITHOUT ANY WARRANTY; without even the implied warranty of
+//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+//  GNU General Public License for more details.
+//  
+//  You should have received a copy of the GNU General Public License
+//  along with this program; see the file COPYING.  If not, write to
+//  the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+//  Boston, MA 02111-1307, USA.
+
+#include "map_base.hpp"
+
+//  coupled_logistic map: x[n+1] = r * x[n] * (1 - x[n]) + e * (y[n] - x[n])
+//                        y[n+1] = r * y[n] * (1 - y[n]) + e * (x[t] - y[t])
+//                        1 <= r <= 4
+//  taken from Willi-Hans Steeb: Chaos and Fractals
+
+class coupled_logistic:
+	public map_base
+{
+public:
+	coupled_logistic()
+	{
+		m_num_eq = 2;
+		m_data = new data_t[2];
+		CHAOS_SYS_INIT(e, 0.06);
+		CHAOS_SYS_INIT(r, 3.7);
+		CHAOS_SYS_INIT(x, 0.1);
+		CHAOS_SYS_INIT(x, 0.2);
+	}
+
+	~coupled_logistic()
+	{
+		delete m_data;
+	}
+
+	virtual void m_step()
+	{
+		data_t x = m_data[0];
+		data_t y = m_data[1];
+		data_t e = CHAOS_PARAMETER(e);
+		data_t r = CHAOS_PARAMETER(r);
+		m_data[0] = r * x * (1.f - x) + e * (y - x);
+		m_data[1] = r * y * (1.f - y) + e * (x - y);
+	}
+
+	CHAOS_SYSPAR_FUNCS(e);
+	CHAOS_SYSPAR_FUNCS_PRED(r, m_pred_r);
+	bool m_pred_r(t_float f)
+	{
+		return (f > 0) && (f < 4);
+	}
+
+	CHAOS_SYSVAR_FUNCS(x, 0);
+	CHAOS_SYSVAR_FUNCS(y, 0);
+	
+};
+
+#define COUPLED_LOGISTIC_CALLBACKS				\
+MAP_CALLBACKS;									\
+CHAOS_SYS_CALLBACKS(e);							\
+CHAOS_SYS_CALLBACKS(r);							\
+CHAOS_SYS_CALLBACKS(x);							\
+CHAOS_SYS_CALLBACKS(y);
+
+#define COUPLED_LOGISTIC_ATTRIBUTES				\
+MAP_ATTRIBUTES;									\
+CHAOS_SYS_ATTRIBUTE(e);							\
+CHAOS_SYS_ATTRIBUTE(r);							\
+CHAOS_SYS_ATTRIBUTE(x);							\
+CHAOS_SYS_ATTRIBUTE(y);
+
-- 
cgit v1.2.1