Что не так с моим кодом набора Мандельброта? - PullRequest
Новая
       5

Что не так с моим кодом набора Мандельброта?

6 голосов
/ 20 октября 2011

Я пытаюсь реализовать множество Мандельброта в C, но у меня странная проблема.Мой код выглядит следующим образом: 

#include <stdio.h>
#include <math.h>
#include <complex.h>

int iterate_pt(complex c);

int main() {
FILE *fp;
fp = fopen("mand.ppm", "w+");


double crmin = -.75;
double crmax = -.74;
double cimin = -.138;
double cimax = -.75; //Changing this value to -.127 fixed my problem.

int ncols = 256;
int nrows = 256;
int mand[ncols][nrows];
int x, y, color;
double complex c;

double dx = (crmax-crmin)/ncols;
double dy = (cimax-cimin)/nrows;

for (x = 0; x < ncols; x++){
    for (y = 0; y < nrows; y++){
        double complex imaginary = 0+1.0i;
        c = crmin+(x*dx) + (cimin+(y*dy)) * imaginary;
        mand[x][y] = iterate_pt(c);
    }
}

printf("Printing ppm header.");
fprintf(fp, "P3\n");
fprintf(fp, "%d %d\n255\n\n", ncols, nrows);

for (x = 0; x < ncols; x++) {
    for (y = 0; y < nrows; y++){
        color = mand[x][y];
        fprintf(fp, "%d\n", color);
        fprintf(fp, "%d\n", color);
        fprintf(fp, "%d\n\n", color); //Extra new line added, telling the ppm to go to next pixel.
    }
}
fclose(fp);

return 0;
}

int iterate_pt(double complex c){
double complex z = 0+0.0i;
int iterations = 0;
int k;
for (k = 1; k <= 255; k++) {
    z = z*z + c;
    if (sqrt( z*conj(z) ) > 50){
        break;
    }
    else
        ++iterations;
}
return iterations;
}

Тем не менее, вывод этой программы, которая хранится в виде ppm-файла выглядит так:за вашу помощь!

Ответы [ 3 ]

3 голосов
/ 20 октября 2011

Попробуйте установить cimax на -0.127, я тоже работаю над этим проектом, и, похоже, он добился цели;)

2 голосов
/ 20 октября 2011

Код выглядит хорошо. Но ваш начальный прямоугольник выглядит неправильно!

вы используете 

Real ranage [  -.75 ,  -.74 ]
Imag range  [ -.138 ,  -.75 ]

Вы уверены, что это то, что вы хотели? Мне кажется, что это ужасно растянутая шкала y.

Кроме того, стандартные алгоритмы Мандельброта имеют тенденцию использовать 

magnitude > 2

, а не 50. в качестве спасательного чека. Хотя это не должно влиять на фактическую форму набора.

0 голосов
/ 28 октября 2011

Между прочим, нет смысла вычислять площадь z * con (z).Просто возведите в квадрат выражения по обе стороны неравенства, получив if (z*conj(z) > 2500), и вы повысите производительность.

...