void Generate(image2d_t outputImage)
{
// get id of element in array
int x = get_global_id(0);
int y = get_global_id(1);
int w = get_global_size(0);
int h = get_global_size(1);
float4 result = (float4)(0.0f, 0.0f, 0.0f, 1.0f);
float MinRe = -2.0f;
float MaxRe = 1.0f;
float MinIm = -1.5f;
float MaxIm = MinIm + (MaxRe - MinRe) * h / w;
float Re_factor = (MaxRe - MinRe) / (w - 1);
float Im_factor = (MaxIm - MinIm) / (h - 1);
float MaxIterations = 50;
//C imaginary
float c_im = MaxIm - y * Im_factor;
//C real
float c_re = MinRe + x * Re_factor;
//Z real
float Z_re = c_re, Z_im = c_im;
bool isInside = true;
bool col2 = false;
bool col3 = false;
int iteration = 0;
for (int n = 0; n < MaxIterations; n++)
{
// Z - real and imaginary
float Z_re2 = Z_re * Z_re, Z_im2 = Z_im * Z_im;
//if Z real squared plus Z imaginary squared is greater than c squared
if (Z_re2 + Z_im2 > 4)
{
if (n >= 0 && n <= (MaxIterations / 2 - 1))
{
col2 = true;
isInside = false;
break;
}
else if (n >= MaxIterations / 2 && n <= MaxIterations - 1)
{
col3 = true;
isInside = false;
break;
}
}
Z_im = 2 * Z_re * Z_im + c_im;
Z_re = Z_re2 - Z_im2 + c_re;
iteration++;
}
if (col2)
{
result = (float4)(iteration * 0.05f, 0.0f, 0.0f, 1.0f);
}
else if (col3)
{
result = (float4)(255, iteration * 0.05f, iteration * 0.05f, 1.0f);
}
else if (isInside)
{
result = (float4)(0.0f, 0.0f, 0.0f, 1.0f);
}
write_imagef(outputImage, (int2)(x, y), result);
}
public void write_imagef(image2d_t output, int2 mask, float4 @float4)
{
}