- start of external Imager API access:
[imager.git] / Mandelbrot / mandel.c
1 #include "imext.h"
2 #include <stdlib.h>
3
4 char evalstr[]="Mandlebrot renderer";
5
6 /* Example Mandlebrot generator */
7
8 /* input parameters
9    image is the image object.
10 */
11
12
13 static
14 int
15 mandel(double x, double y, int max_iter) {
16   double xn, yn;
17   double xo, yo;
18   double dist;
19   int iter = 1;
20   /*    Z(n+1) = Z(n) ^2 + c */
21
22   /* printf("(%.2f, %.2f) -> \n", x,y);   */
23
24   xo = x;
25   yo = y;
26
27   while( xo*xo+yo*yo <= 10 && iter < max_iter) {
28     xn = xo*xo-yo*yo + x;
29     yn = 2*xo*yo     + y;
30     xo=xn;
31     yo=yn;
32     iter++;
33   }
34   return (iter == max_iter)?0:iter;
35 }
36
37 void 
38 mandelbrot(i_img *im, double minx, double miny, double maxx, double maxy, int max_iter) {
39
40   i_color vl;
41   int i,bytes,x,y;
42   int idx;
43   double divx, divy;
44
45   i_color icl[256];
46   srand(12235);
47   for(i=1;i<256; i++) {
48     icl[i].rgb.r = 100+(int) (155.0*rand()/(RAND_MAX+1.0));
49     icl[i].rgb.g = 100+(int) (155.0*rand()/(RAND_MAX+1.0));
50     icl[i].rgb.g = 100+(int) (155.0*rand()/(RAND_MAX+1.0));
51   }
52
53   icl[0].rgb.r = 0;
54   icl[0].rgb.g = 0;
55   icl[0].rgb.g = 0;
56     
57   if (maxx <= minx)
58     maxx = minx + 1.0;
59   if (maxy <= miny)
60     maxy = miny + 1.0;
61
62   divx = (maxx - minx) / im->xsize;
63   divy = (maxy - miny) / im->ysize;
64
65   for(y = 0; y < im->ysize; y ++) {
66     for(x = 0; x < im->xsize; x ++ ) {
67       idx = mandel(minx + x*divx , miny + y*divy, max_iter);
68       idx = idx % 256;
69       i_ppix(im,x,y,&icl[idx]); 
70     }
71   }
72 }