代码帮助确定点是否在Mandelbrot集中(检查我的解决方案)
||
这是我的函数,用于测试两个点x和y在MAX_ITERATION 255之后是否位于mandelbrot集中。如果不存在,则返回0,否则返回1。
int isMandelbrot (int x, int y) {
int i;
int j;
double Re[255];
double Im[255];
double a;
double b;
double dist;
double finaldist;
int check;
i=0;
Re[0]=0;
Im[0]=0;
j=-1;
a=0;
b=0;
while (i < MAX_ITERATION) {
a = Re[j];
b = Im[j];
Re[i]=((a*a)-(b*b))+x;
Im[i]=(2 * a * b) + y;
i++;
j++;
}
finaldist = sqrt(pow(Re[MAX_ITERATION],2)+pow(Im[MAX_ITERATION],2));
if (dist > 2) { //not in mandelbrot
check = 0;
} else if (dist <= 2) { //in mandelbrot set
check = 1;
}
return check;
}
鉴于它是正确的(有人可以验证...还是编写更有效的代码?)。
这是我的打印代码,但是不起作用! (它一直在给出所有点都在集合中)。我在这里做错了什么?
int main(void) {
double col;
double row;
int checkSet;
row = -4;
col = -1;
while (row < 1.0 ) {
while (col < 1.0) {
checkSet = isMandelbrot(row, col);
if (checkSet == 1) {
printf(\"-\");
} else if (checkSet == 0) {
printf(\"*\");
}
col=col+0.5;
}
col=-1;
row=row+0.5;
printf(\"\\n\");
}
return 0;
}
没有找到相关结果
已邀请:
3 个回复
厦惫
但是在第一次迭代中,
,因此您将获得数组索引-1处的值。那不是你想做的。 另外,为什么是
和
数组-您真的需要在计算中跟踪所有中间结果吗? Wikipedia包含该算法的伪代码,您可能需要对照该代码检查自己的代码。
玩翁文醚碱
参数,因此
输入的值将被截断(即小数部分将被丢弃)。
坝硷操
可以表示的结果-结果为
,因此您的finaldist计算会产生垃圾。 我认为您将从主要解决方案中受益。您在此处输入的代码计算的像素不足,无法真正看到其中的大部分内容。