From da609ba8c21e607196daca60b7d994eb2cb6fe4c Mon Sep 17 00:00:00 2001
From: Michael Bien <mbien@fh-landshut.de>
Date: Tue, 26 Jan 2010 23:42:24 +0100
Subject: initial import of GL-CL Mandelbrot Set example.

---
 src/com/mbien/opencl/demos/fractal/Mandelbrot.cl | 44 ++++++++++++++++++++++++
 1 file changed, 44 insertions(+)
 create mode 100644 src/com/mbien/opencl/demos/fractal/Mandelbrot.cl

(limited to 'src/com/mbien/opencl/demos/fractal/Mandelbrot.cl')

diff --git a/src/com/mbien/opencl/demos/fractal/Mandelbrot.cl b/src/com/mbien/opencl/demos/fractal/Mandelbrot.cl
new file mode 100644
index 0000000..e6583f4
--- /dev/null
+++ b/src/com/mbien/opencl/demos/fractal/Mandelbrot.cl
@@ -0,0 +1,44 @@
+/**
+ * For a description of this algorithm please refer to
+ * http://en.wikipedia.org/wiki/Mandelbrot_set
+ * @author Michael Bien
+ */
+kernel void mandelbrot(
+        global uint *output,
+        const int width, const int height,
+        const float x0, const float y0,
+        const float x1, const float y1,
+        global uint *colorMap, const int colorMapSize, const int maxIterations) {
+
+    unsigned int ix = get_global_id(0);
+    unsigned int iy = get_global_id(1);
+
+    float r = x0 + ix * (x1-x0) / width;
+    float i = y0 + iy * (y1-y0) / height;
+
+    float x = 0;
+    float y = 0;
+
+    float magnitudeSquared = 0;
+    int iteration = 0;
+
+    while (magnitudeSquared < 4 && iteration < maxIterations) {
+        float x2 = x*x;
+        float y2 = y*y;
+        y = 2 * x * y + i;
+        x = x2 - y2 + r;
+        magnitudeSquared = x2+y2;
+        iteration++;
+    }
+
+    if (iteration == maxIterations)  {
+        output[iy * width + ix] = 0;
+    }else {
+        float alpha = (float)iteration / maxIterations;
+        int colorIndex = (int)(alpha * colorMapSize);
+        output[iy * width + ix] = colorMap[colorIndex];
+      // monochrom
+      //  output[iy * width + ix] = 255*iteration/maxIterations;
+    }
+
+}
\ No newline at end of file
-- 
cgit v1.2.3