josch/mandelbrot
1
0
Fork 0
Code Issues Pull requests Projects Releases Wiki Activity

fix bailout value

Browse source
This commit is contained in:
josch 2013-06-28 15:54:33 +02:00
parent fd57bd5120
commit bcf75e0e2e
2 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
mandel_dd.c
View file

@ -22,12 +22,12 @@ int main(int argc, char **argv) {
DoubleDouble centerx, centery; DoubleDouble centerx, centery;
centerx = dd_new(-0.7436438870371587, -3.628952515063387E-17); centerx = dd_new(-0.7436438870371587, -3.628952515063387E-17);
centery = dd_new(0.13182590420531198, -1.2892807754956678E-17); centery = dd_new(0.13182590420531198, -1.2892807754956678E-17);
double bailout = 4; // the distance must not be greater than 2 (4 = 2*2) double bailout = 128; // with a smaller value there are lines on magn=1
double logLogBailout = log(log(bailout)); double logLogBailout = log(log(bailout));
DoubleDouble magn = dd_new(strtod(argv[5], NULL), 0); DoubleDouble magn = dd_new(strtod(argv[5], NULL), 0);
// maxiter = width * sqrt(magn); /*// maxiter = width * sqrt(magn);
temp1 = dd_sqrt(magn); temp1 = dd_sqrt(magn);
unsigned long maxiter = width * dd_get_ui(temp1); unsigned long maxiter = width * dd_get_ui(temp1);*/
DoubleDouble x2, y2, x0d, y1d; DoubleDouble x2, y2, x0d, y1d;
// x0d = 4 / magn / width; // x0d = 4 / magn / width;
x0d = dd_ui_div(4, magn); x0d = dd_ui_div(4, magn);

2
mandel_mpfr.c
View file

@ -25,7 +25,7 @@ int main(int argc, char **argv) {
mpfr_init2(centery, 300); mpfr_init2(centery, 300);
mpfr_set_str(centerx, argv[3], 10, MPFR_RNDN); mpfr_set_str(centerx, argv[3], 10, MPFR_RNDN);
mpfr_set_str(centery, argv[4], 10, MPFR_RNDN); mpfr_set_str(centery, argv[4], 10, MPFR_RNDN);
double bailout = 4; // the distance must not be greater than 2 (4 = 2*2) double bailout = 128; // with a smaller value there are lines on magn=1
double logLogBailout = log(log(bailout)); double logLogBailout = log(log(bailout));
mpfr_t magn; mpfr_t magn;
mpfr_init2(magn, 300); mpfr_init2(magn, 300);