From 43dd4efedf1ecfe721cde5830bdcee67ffa48907 Mon Sep 17 00:00:00 2001
From: Tim Blechmann <timblech@users.sourceforge.net>
Date: Tue, 28 Dec 2004 12:38:45 +0000
Subject: a few new attractors ...

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

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

diff --git a/chaos/src/latoomutalpha.hpp b/chaos/src/latoomutalpha.hpp
new file mode 100644
index 0000000..6926d1b
--- /dev/null
+++ b/chaos/src/latoomutalpha.hpp
@@ -0,0 +1,88 @@
+// 
+//  
+//  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"
+#include <cmath>
+
+//  latoocarfian model, mutation alpha:
+//           x1[n+1] = sin(x2[n]*b) + pow(sin(x1[n]*b),2) + pow(sin(x1[n]*b),3)
+//           x2[n+1] = sin(x1[n]*a) + pow(sin(x2[n]*a),2) + pow(sin(x2[n]*c),3)
+//  taken from Pickover: Chaos In Wonderland
+
+class latoomutalpha
+	: public map_base
+{
+public:
+	latoomutalpha()
+	{
+		m_num_eq = 2;
+		m_data = new data_t[m_num_eq];
+
+		CHAOS_SYS_INIT(x1,0.5);
+		CHAOS_SYS_INIT(x2,0.2);
+		CHAOS_SYS_INIT(a,-0.966918);
+		CHAOS_SYS_INIT(b,2.879879);
+		CHAOS_SYS_INIT(c,0.765145);
+	}
+	
+	~latoomutalpha()
+	{
+		delete m_data;
+	}
+
+	virtual void m_step()
+	{
+		data_t x1 = m_data[0], x2 = m_data[1];
+		data_t a = CHAOS_PARAMETER(a), b = CHAOS_PARAMETER(b),
+			c = CHAOS_PARAMETER(c);
+		data_t tmp1, tmp2;
+		
+		tmp1 = sin(x1*b);
+		m_data[0] = sin(x2*b) + tmp1*tmp1 + tmp1*tmp1*tmp1;
+		tmp1 = sin(x2*a);
+		tmp2 = sin(x2*c);
+		m_data[1] = sin(x1*a) + tmp1*tmp1 + tmp2*tmp2*tmp2;
+	}
+
+	CHAOS_SYSVAR_FUNCS(x1, 0);
+	CHAOS_SYSVAR_FUNCS(x2, 1);
+
+	CHAOS_SYSPAR_FUNCS(a);
+	CHAOS_SYSPAR_FUNCS(b);
+	CHAOS_SYSPAR_FUNCS(c);
+};
+
+
+#define LATOOMUTALPHA_CALLBACKS					\
+MAP_CALLBACKS;									\
+CHAOS_SYS_CALLBACKS(x1);						\
+CHAOS_SYS_CALLBACKS(x2);						\
+CHAOS_SYS_CALLBACKS(a);							\
+CHAOS_SYS_CALLBACKS(b);							\
+CHAOS_SYS_CALLBACKS(c);
+
+
+#define LATOOMUTALPHA_ATTRIBUTES				\
+MAP_ATTRIBUTES;									\
+CHAOS_SYS_ATTRIBUTE(x1);						\
+CHAOS_SYS_ATTRIBUTE(x2);						\
+CHAOS_SYS_ATTRIBUTE(a);							\
+CHAOS_SYS_ATTRIBUTE(b);							\
+CHAOS_SYS_ATTRIBUTE(c);
-- 
cgit v1.2.1