1

这是我的函数，它测试两个点 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;
}
4

3 回答 3

0

您的代码中有一些错误。例如，您这样做：

a = Re[j];
b = Im[j];

但是在第一次迭代中j = -1，您将获得数组索引 -1 处的值。那不是你想做的。

另外，为什么是Re和Im数组——你真的需要跟踪计算中的所有中间结果吗？

Wikipedia 包含算法的伪代码，您可能需要检查您自己的代码。

于 2011-04-16T10:04:44.373 回答
0

另一个错误：您的函数需要int参数，因此您的double输入值将被截断（即小数部分将被丢弃）。

于 2011-04-16T10:10:24.943 回答
0

您可能应该在 while 循环中检查转义。也就是说，如果 ((a*a + b*b) > 4)在任何时候，那么该像素已经逃脱，故事结束。通过继续迭代这些像素，以及浪费 CPU 周期，您的值正在无限制地增长，并且似乎超过了可以在 a 中表示的值double- 结果是NaN，因此您的 finaldist 计算正在产生垃圾。

我认为你会从你的主要解决方案中受益。您放在这里的代码没有计算足够的像素来真正看到大部分集合。

于 2011-04-16T12:48:43.510 回答