Júlia halmaz - WebCL

Julia fraktál rajzolása.

Készítette Ujj László - 2012

Alaphelyzetben Julia halmaz rajzolása Keszthelyi Balázs példájából kiindulva, de a kernel szabadon átírható kisérletezés céljából!

Iterációk:
a:
b:
Kernel: __kernel void webclSample(__global uchar rgb, uint MaxIterations, float aa, float bb) { const uint x = get_global_id(0), y = get_global_id(1), _width = get_global_size(0), _height = get_global_size(1), i = y_width + x; const float a = convert_float(aa), b = convert_float(bb), MinRe = -1.5f, MaxRe = 1.5f, MinIm = -1.0f, MaxIm = MinIm+(MaxRe-MinRe)convert_float(_height)/convert_float(_width), Re_factor = (MaxRe-MinRe)/convert_float(_width-1), Im_factor = (MaxIm-MinIm)/convert_float(_height-1); const float2 c = (float2)(a,b); float2 Z = (float2)(MinRe + convert_float(x)Re_factor, MaxIm - convert_float(y)Im_factor); uint n; for(n = 0; n < MaxIterations; ++n) { float2 Z2 = ZZ; if(Z2.s0 + Z2.s1 > 4.0f) { break; } Z = (float2)(Z2.s0 - Z2.s1, 2Z.s0Z.s1) + c; } uchar color; if(n == MaxIterations) { rgb[4i + 0] = 255; rgb[4i + 1] = 255; rgb[4i + 2] = 255; rgb[4i + 3] = 255; } else if(convert_float(n)<(convert_float(MaxIterations) / 3.0f)) { color = convert_uchar(convert_float(0xFF) * convert_float(n)3 / convert_float(MaxIterations)); rgb[4i + 0] = 0; rgb[4i + 1] = color; rgb[4i + 2] = 0; rgb[4i + 3] = 255; } else if(convert_float(n)<(convert_float(MaxIterations) 2.0f / 3.0f)) { color = convert_uchar(convert_float(0xFF) * convert_float(n)3.0f /2.0f / convert_float(MaxIterations)); rgb[4i + 0] = 0; rgb[4i + 1] = 255; rgb[4i + 2] = color; rgb[4i + 3] = 255; } else { color = convert_uchar(convert_float(0xFF) convert_float(n) / convert_float(MaxIterations)); rgb[4i + 0] = color; rgb[4i + 1] = 255; rgb[4i + 2] = 255; rgb[4i + 3] = 255; } }
Log: