算法之美---由计算机生成的图像

/*
unsigned short RD(int i,int j){
double r=i/256.-2,s=j/256.-2,q=r*r+s*s,n=hypot(r+(.866-r/2)/q,s+(r*.866+s/2)/q),
d=.5/log(n);if(d<0||d>1)d=1;return d*(sin(n*10)*511+512);
}
unsigned short GR(int i,int j){
return 0;
}
unsigned short BL(int i,int j){
double r=i/256.-2,s=j/256.-2,q=r*r+s*s;return RD(i,j)*sqrt(q/40);
}
*/

// --------------------------------------------------------------------------------------

unsigned int    CPixelJoukowsky::CalculatePixel(unsigned int i, unsigned int j)
{
double r=i/256.0-2;
double s=j/256.0-2;
double q=r*r+s*s;
double n=hypot(r+(0.866-r/2)/q,s+(r*.866+s/2)/q);
double d=0.5/log(n);

if(d<0||d>1)
{
d=1;
}

float fr = (float)(d*(sin(n*10)*511+512));

unsigned int br = FLOAT_255_TO_BYTE(fr);
unsigned int bg = 0;

float fb = fr*sqrtf((float)q/40);
unsigned int bb = FLOAT_255_TO_BYTE(fb);

return MAKE_RGB(br,bg,bb);
}