Iterointi luuppi vaikka Mandelbrotin joukon laskemiseen. Koodi on optimoitu nopeimpaan kuosiin, vai onko sinulla ideaa, miten luuppiin saisi enemmän tehoa?
/*
Jouni Aro
22.2.2000
*/
#ifndef __LONG_DOUBLE
#define __LONG_DOUBLE
typedef long double lreal;
#define lr(ex) (lreal)(ex)
#endif
#define MAX_ITER 256
int iteroi(lreal a, lreal b)
{
int loop=0;
lreal rr=lr(0);
lreal x=0, y=0, r, i;
while (rr < lr(4096))
{
r=x*x;
i=y*y;
rr=r+i;
y=2*x*y+b;
x=r-i+a;
if (++loop==MAX_ITER) return 0;
}
return loop;
}