Basic Mandelbrot Loop

loop.cpp

unsigned char iter;
for(unsigned int y=0; y < h; y++)
{
    for(unsigned int x=0; x < w; x++)
    {
        iter  = 0;
        z_re  = 0;
        z_im  = 0;
        z_re2 = 0;
        z_im2 = 0;

        while(z_re2 + z_im2 < THRESHOLD && iter++ < MAX_ITERS)
        {
            z_re2 = z_re * z_re;
            z_im2 = z_im * z_im;
            z_im  = (2.0 * z_re * z_im) + c_im;
            z_re  = z_re2 - z_im2 + c_re;
        }
        c_re += inc_re;
        ofs << iter << iter << iter;
    }
    c_re = init_re;
    c_im += inc_im;
}