From f2f4a1c9cd553f82944d889338d9aff7903c519b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jos=C3=A9=20Carlos=20Cuevas?= <reset.reboot@gmail.com>
Date: Sun, 1 Sep 2013 02:58:48 +0200
Subject: [PATCH] Added to repo the julia kernel too

---
 julia_kernel.cl | 48 ++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 48 insertions(+)
 create mode 100644 julia_kernel.cl

diff --git a/julia_kernel.cl b/julia_kernel.cl
new file mode 100644
index 0000000..71636ea
--- /dev/null
+++ b/julia_kernel.cl
@@ -0,0 +1,48 @@
+float map_x(int x, int width, float zoom)
+{
+    return (((float)x / (float)width) * (3.5 * zoom)) - (1.75 - (1.0 - zoom));
+}
+
+float map_y(int y, int height, float zoom)
+{
+    return (((float)y / (float)height) * (2.0 * zoom)) - (1.00001 - (1.0 - zoom));
+}
+
+__kernel void fractal_point(__global const int *res_x, 
+                               __global const int *res_y, 
+                               __global const int *line, 
+                               __global const float *zoom, 
+                               __global int *graph_line) 
+{
+    // Get the index of the current element
+    int image_x = get_global_id(0);
+    int image_y = *line;
+    float x = map_x(image_x, *res_x, 1.0);
+    float y = map_y(image_y, *res_y, 1.0);
+
+    int iteration = 0;
+    int max_iteration = 256;
+    float xtemp, xx, yy;
+
+    while (iteration < max_iteration)
+    {
+       xx = x * x;
+       yy = y * y;
+       if ((xx) + (yy) > (4.0)) break;
+
+       xtemp = xx - yy + 0.353 + *zoom;
+       y = 2.0 * x * y + 0.288;
+
+       x = xtemp;
+       iteration++;
+    }
+
+    if (iteration > max_iteration)
+    {
+       graph_line[image_x] = 0;
+    }
+    else
+    {
+       graph_line[image_x] = iteration;
+    }
+}