View
266 src/fractal.c
@@ -1,5 +1,108 @@
#include "fractal.h"
#include "render_threads.h"
#include <math.h>
int frac_mandel(int depth,
long double wim, long double wre,
long double c_im, long double c_re,
long double wim2, long double wre2 )
{
int wz;
for (wz = 1; wz <= depth; ++wz)
{
wim = 2.0 * wre * wim + c_im;
wre = wre2 - wim2 + c_re;
wim2 = wim * wim;
wre2 = wre * wre;
if (wim2 + wre2 > 4.0F)
return wz;
}
return 0;
}
int frac_mandel_mpfr(int depth,
mpfr_t t1, /* t1 == temporary */
mpfr_t frs_bail,
mpfr_t wim, mpfr_t wre,
mpfr_t c_im, mpfr_t c_re,
mpfr_t wim2, mpfr_t wre2 )
{
int wz;
for (wz = 1; wz <= depth; wz++)
{
/* wim = 2.0 * wre * wim + c_im; */
mpfr_mul( t1, wre, wim, GMP_RNDN);
mpfr_mul_si(t1, t1, 2, GMP_RNDN);
mpfr_add( wim, t1, c_im, GMP_RNDN);
/* wre = wre2 - wim2 + c_re; */
mpfr_sub( t1, wre2, wim2, GMP_RNDN);
mpfr_add( wre, t1, c_re, GMP_RNDN);
/* wim2 = wim * wim; */
mpfr_mul( wim2, wim, wim, GMP_RNDN);
/* wre2 = wre * wre; */
mpfr_mul( wre2, wre, wre, GMP_RNDN);
/* if ((wim2 + wre2) > frs_bail) */
mpfr_add( t1, wim2, wre2, GMP_RNDN);
if (mpfr_greater_p(t1, frs_bail))
return wz;
}
return 0;
}
#ifdef WITH_GMP
int frac_mandel_gmp(int depth, mpf_t frs_bail,
mpf_t t1, /* t1 == temporary */
mpf_t wim, mpf_t wre,
mpf_t c_im, mpf_t c_re,
mpf_t wim2, mpf_t wre2 )
{
int wz;
for (wz = 1; wz <= depth; wz++)
{
/* wim = 2.0 * wre * wim + c_im; */
mpf_mul( t1, wre, wim);
mpf_mul_ui(t1, t1, 2);
mpf_add( wim, t1, c_im);
/* wre = wre2 - wim2 + c_re; */
mpf_sub( t1, wre2, wim2);
mpf_add( wre, t1, c_re);
/* wim2 = wim * wim; */
mpf_mul( wim2, wim, wim);
/* wre2 = wre * wre; */
mpf_mul( wre2, wre, wre);
/* if ((wim2 + wre2) > frs_bail) */
mpf_add( t1, wim2, wre2);
if (mpf_cmp(t1, frs_bail) > 0)
return wz;
}
return 0;
}
#endif
int frac_celtic(int depth,
long double wim, long double wre,
long double c_im, long double c_re,
long double wim2, long double wre2)
{
int wz;
for (wz = 1; wz <= depth; ++wz)
{
wim = 2.0 * wre * wim + c_im;
wre = fabsl(wre2 - wim2) + c_re;
wim2 = wim * wim;
wre2 = wre * wre;
if (wim2 + wre2 > 4.0F)
return wz;
}
return 0;
}
int fractal_calculate_line(image_info* img, int line)
{
@@ -62,25 +165,18 @@ int fractal_calculate_line(image_info* img, int line)
c_im = j_c_im;
break;
}
for (wz = 1; wz <= depth; wz++)
{
wim = 2.0 * wre * wim+ c_im;
wre = wre2 - wim2 + c_re;
wim2 = wim * wim;
wre2 = wre * wre;
if (wim2 + wre2 > 4.0F)
break;
}
if (wz == depth + 1)
wz = 0;
*raw_data = wz;
*raw_data = frac_mandel(depth, wim, wre,
c_im, c_re,
wim2, wre2);
}
if (rth_render_should_stop((rthdata*)img->rth_ptr))
return 0;
}
return 1;
}
int fractal_mpfr_calculate_line(image_info* img, int line)
{
int ret = 1;
@@ -124,7 +220,7 @@ int fractal_mpfr_calculate_line(image_info* img, int line)
mpfr_set_si(frs_bail, 4, GMP_RNDN);
mpfr_set_si(img_rw, img_width, GMP_RNDN);
mpfr_set( img_xmin, img->xmin, GMP_RNDN);
mpfr_set( width, img->pcoords->width, GMP_RNDN);
mpfr_set( width, img->width, GMP_RNDN);
/* y = img->ymax - ((img->xmax - img->xmin)
/ (long double)img->real_width)
@@ -168,27 +264,10 @@ int fractal_mpfr_calculate_line(image_info* img, int line)
break;
}
for (wz = 1; wz <= depth; wz++)
{
/* wim = 2.0 * wre * wim + c_im; */
mpfr_mul( t1, wre, wim, GMP_RNDN);
mpfr_mul_si(t1, t1, 2, GMP_RNDN);
mpfr_add( wim, t1, c_im, GMP_RNDN);
/* wre = wre2 - wim2 + c_re; */
mpfr_sub( t1, wre2, wim2, GMP_RNDN);
mpfr_add( wre, t1, c_re, GMP_RNDN);
/* wim2 = wim * wim; */
mpfr_mul( wim2, wim, wim, GMP_RNDN);
/* wre2 = wre * wre; */
mpfr_mul( wre2, wre, wre, GMP_RNDN);
/* if ((wim2 + wre2) > frs_bail) */
mpfr_add( t1, wim2, wre2, GMP_RNDN);
if (mpfr_greater_p(t1, frs_bail))
break;
}
if (wz == depth + 1)
wz = 0;
*raw_data = wz;
*raw_data = frac_mandel_mpfr(depth, t1, frs_bail,
wim, wre,
c_im, c_re,
wim2, wre2 );
}
if (rth_render_should_stop((rthdata*)img->rth_ptr))
@@ -213,3 +292,122 @@ int fractal_mpfr_calculate_line(image_info* img, int line)
mpfr_clear(t1);
return ret;
}
#ifdef WITH_GMP
int fractal_gmp_calculate_line(image_info* img, int line)
{
int ret = 1;
int ix = 0;
int mx = 0;
int chk_px = ((rthdata*)img->rth_ptr)->check_stop_px;
int img_width = img->real_width;
int* raw_data = &img->raw_data[line * img_width];
depth_t wz;
depth_t depth = img->depth;
mpf_t x, y;
mpf_t x2, y2;
mpf_t c_re, c_im;
/* working variables: */
mpf_t wre, wim;
mpf_t wre2, wim2;
mpf_t frs_bail;
mpf_t width, img_rw, img_xmin;
mpf_t t1;
mpf_init2(x, img->precision);
mpf_init2(y, img->precision);
mpf_init2(x2, img->precision);
mpf_init2(y2, img->precision);
mpf_init2(c_re, img->precision);
mpf_init2(c_im, img->precision);
mpf_init2(wre, img->precision);
mpf_init2(wim, img->precision);
mpf_init2(wre2, img->precision);
mpf_init2(wim2, img->precision);
mpf_init2(frs_bail,img->precision);
mpf_init2(width, img->precision);
mpf_init2(img_rw, img->precision);
mpf_init2(img_xmin,img->precision);
mpf_init2(t1, img->precision);
mpf_set_si(frs_bail, 4);
mpf_set_si(img_rw, img_width);
mpf_set( img_xmin, img->gxmin);
mpf_set( width, img->gwidth);
/* y = img->ymax - ((img->xmax - img->xmin)
/ (long double)img->real_width)
* (long double)img->lines_done; */
mpf_div( t1, width, img_rw);
mpf_mul_ui( t1, t1, line);
mpf_sub( y, img->gymax, t1);
mpf_mul( y2, y, y);
while (ix < img_width)
{
mx += chk_px;
if (mx > img_width)
mx = img_width;
for (; ix < mx; ++ix, ++raw_data)
{
/* x = ((long double)ix / (long double)img->real_width)
* (img->xmax - img->xmin) + img->xmin; */
mpf_ui_div(t1, ix, img_rw);
mpf_mul(x, t1, width);
mpf_add(x, x, img_xmin);
mpf_mul( x2, x, x);
mpf_set( wre, x);
mpf_set( wim, y);
mpf_set( wre2, x2);
mpf_set( wim2, y2);
switch (img->fr_type)
{
case MANDELBROT:
mpf_set(c_re, x);
mpf_set(c_im, y);
break;
case JULIA:
// mpf_set(c_re, img->u.julia.c_re);
// mpf_set(c_im, img->u.julia.c_im);
break;
}
*raw_data = frac_mandel_gmp(depth, frs_bail, t1,
wim, wre,
c_im, c_re,
wim2, wre2 );
}
if (rth_render_should_stop((rthdata*)img->rth_ptr))
{
ret = 0;
break;
}
}
mpf_clear(x);
mpf_clear(y);
mpf_clear(x2);
mpf_clear(y2);
mpf_clear(c_re);
mpf_clear(c_im);
mpf_clear(wre);
mpf_clear(wim);
mpf_clear(wre2);
mpf_clear(wim2);
mpf_clear(frs_bail);
mpf_clear(width);
mpf_clear(img_rw);
mpf_clear(t1);
return ret;
}
#endif /* WITH_GMP */
View
5 src/fractal.h
@@ -17,4 +17,9 @@
int fractal_calculate_line(image_info* img, int line);
int fractal_mpfr_calculate_line(image_info* img, int line);
#ifdef WITH_GMP
int fractal_gmp_calculate_line(image_info* img, int line);
#endif
#endif
View
26 src/image_info.c
@@ -56,6 +56,13 @@ image_info * image_info_create(fractal_type fr)
mpfr_init2( img->u.julia.c_re, DEFAULT_PRECISION);
mpfr_init2( img->u.julia.c_im, DEFAULT_PRECISION);
#ifdef WITH_GMP
mpf_init2( img->gxmin, DEFAULT_PRECISION);
mpf_init2( img->gxmax, DEFAULT_PRECISION);
mpf_init2( img->gymax, DEFAULT_PRECISION);
mpf_init2( img->gwidth, DEFAULT_PRECISION);
#endif
img->pcoords = coords_new(img->real_width, img->real_height,
cx, cy, size);
image_info_reset_view(img);
@@ -71,6 +78,13 @@ void image_info_destroy(image_info* img)
mpfr_clear( img->ymax );
mpfr_clear( img->width );
#ifdef WITH_GMP
mpf_clear( img->gxmin );
mpf_clear( img->gxmax );
mpf_clear( img->gymax );
mpf_clear( img->gwidth );
#endif
mpfr_clear( img->old_cx );
mpfr_clear( img->old_cy );
mpfr_clear( img->old_size );
@@ -207,6 +221,13 @@ void image_info_mpfr_change(image_info * img, mp_prec_t precision)
precision_change( img->u.julia.c_re, precision );
precision_change( img->u.julia.c_re, precision );
#ifdef WITH_GMP
precision_change_gmp( img->gxmin, precision );
precision_change_gmp( img->gxmax, precision );
precision_change_gmp( img->gymax, precision );
precision_change_gmp( img->gwidth, precision );
#endif
coords_set_precision( img->pcoords, precision );
}
@@ -234,7 +255,12 @@ void image_info_set_mpfr(image_info* img, gboolean use_mpfr)
if (use_mpfr)
{
img->using_mpfr = TRUE;
#ifdef WITH_GMP
rth_set_next_line_cb(rth, fractal_gmp_calculate_line);
#else
rth_set_next_line_cb(rth, fractal_mpfr_calculate_line);
#endif
}
else
{
View
4 src/image_info.h
@@ -67,6 +67,10 @@ typedef struct image_info
{
mpfr_t xmin, xmax, ymax, width;
#ifdef WITH_GMP
mpf_t gxmin, gxmax, gymax, gwidth;
#endif
/* saved mandel coords. */
/* we need these when switching back from */
/* julia mode */
View
6 src/main_gui.c
@@ -593,6 +593,12 @@ void gui_start_rendering(image_info* img)
coords_get_rect(img->pcoords, img->xmin, img->xmax,
img->ymax, img->width);
#ifdef WITH_GMP
coords_get_rect_gmp(img->pcoords, img->gxmin, img->gxmax,
img->gymax, img->gwidth);
#endif
if (!img->j_pre)
{
gtk_spin_button_update( GTK_SPIN_BUTTON(depth_spin));
View
4 src/render.c
@@ -33,6 +33,10 @@ int render_to_file(image_info* img)
coords_get_rect(img->pcoords, img->xmin, img->xmax,
img->ymax, img->width);
#ifdef WITH_GMP
coords_get_rect_gmp(img->pcoords, img->gxmin, img->gxmax,
img->gymax, img->gwidth);
#endif
rth_ui_init(rth);