jesli zrezygnowac z testu z kwadratem
to tutaj ten kernel z trzema mnozeniami
mozna zamianic na wersje tylko z dowma
mnozeniami

for(n=0; n<=max_iter; n++)
{

reim2 = (re + re) * im;

re = (re - im) * (re + im) + cRe;
im = reim2 + cIm;

// if( fabs(re) > 2.0 ) break;
// if( fabs(im) > 2.0 ) break;

if( re > 2.0 || re < -2.0
&& im > 2.0 || im < -2.0 ) break;


}

to zauwazylem ze szkoly a*a - b*b = (a+b)*(a-b) u mnie w wersji c to akurat
przyspiesza

nie wiadomo jednak czy wersja z dwoma
mul i czterema cmp (statystycznie szczesliwie mniej ) bylaby szybsza
niz kernel z 3 mul i 1 cmp

nalezaloby sprawdzic ale to musialbym
znowu pogrzebac w asmie