home *** CD-ROM | disk | FTP | other *** search
/ Computerworld 1996 March / Computerworld_1996-03_cd.bin / idg_cd3 / grafika / fraktaly / frasr192 / testpt.c < prev    next >
Text File  |  1994-08-18  |  2KB  |  55 lines
  1. /*
  2.  
  3. Write your fractal program here. initreal and initimag are the values in
  4. the complex plane; parm1, and parm2 are paramaters to be entered with the
  5. "params=" option (if needed). The function should return the color associated
  6. with initreal and initimag.  FRACTINT will repeatedly call your function with
  7. the values of initreal and initimag ranging over the rectangle defined by the
  8. "corners=" option. Assuming your formula is iterative, "maxit" is the maximum
  9. iteration. If "maxit" is hit, color "inside" should be returned.
  10.  
  11. Note that this routine could be sped up using external variables/arrays
  12. rather than the current parameter-passing scheme.  The goal, however was
  13. to make it as easy as possible to add fractal types, and this looked like
  14. the easiest way.
  15.  
  16. This module is part of an overlay, with calcfrac.c.  The routines in it
  17. must not be called by any part of Fractint other than calcfrac.
  18.  
  19. The sample code below is a straightforward Mandelbrot routine.
  20.  
  21. */
  22.  
  23. int teststart()     /* this routine is called just before the fractal starts */
  24. {
  25.     return( 0 );
  26. }
  27.  
  28. void testend()         /* this routine is called just after the fractal ends */
  29. {
  30. }
  31.  
  32.         /* this routine is called once for every pixel */
  33.     /* (note: possibly using the dual-pass / solif-guessing options */
  34.  
  35. int testpt(double initreal,double initimag,double parm1,double parm2,
  36.            long maxit,int inside)
  37. {
  38. double oldreal, oldimag, newreal, newimag, magnitude;
  39. long color;
  40.    oldreal=parm1;
  41.    oldimag=parm2;
  42.    magnitude = 0.0;
  43.    color = 0;
  44.    while ((magnitude < 4.0) && (color < maxit)) {
  45.       newreal = oldreal * oldreal - oldimag * oldimag + initreal;
  46.       newimag = 2 * oldreal * oldimag + initimag;
  47.       color++;
  48.       oldreal = newreal;
  49.       oldimag = newimag;
  50.       magnitude = newreal * newreal + newimag * newimag;
  51.       }
  52. if (color >= maxit) color = inside;
  53. return((int)color);
  54. }
  55.