Added to repo the julia kernel too
This commit is contained in:
parent
43ecf7464f
commit
f2f4a1c9cd
1 changed files with 48 additions and 0 deletions
48
julia_kernel.cl
Normal file
48
julia_kernel.cl
Normal file
|
@ -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;
|
||||
}
|
||||
}
|
Loading…
Reference in a new issue