home *** CD-ROM | disk | FTP | other *** search
/ Frostbyte's 1980s DOS Shareware Collection / floppyshareware.zip / floppyshareware / FORTH / FRASRC11.ZIP / TESTPT.C < prev    next >
C/C++ Source or Header  |  1989-07-18  |  2KB  |  56 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. The sample code below is a straightforward Mandelbrot routine.
  17.  
  18. */
  19.  
  20. extern int xdots;        /* the screen is this many dots across */
  21. extern int ydots;        /* the screen is this many dots down */
  22. extern int colors;        /* the screen has this many colors */
  23.  
  24. teststart()    /* this routine is called just before the fractal starts */
  25. {
  26. }
  27.  
  28. 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. testpt(initreal,initimag,parm1,parm2,maxit,inside)
  36. double initreal,initimag,parm1,parm2;
  37. int maxit,inside;
  38. {
  39. double oldreal, oldimag, newreal, newimag, magnitude;
  40. register color;
  41.    oldreal=parm1;
  42.    oldimag=parm2;
  43.    magnitude = 0.0;
  44.    color = 0;
  45.    while ((magnitude < 4.0) && (color < maxit)) {
  46.       newreal = oldreal * oldreal - oldimag * oldimag + initreal;
  47.       newimag = 2 * oldreal * oldimag + initimag;
  48.       color++;
  49.       oldreal = newreal;
  50.       oldimag = newimag;
  51.       magnitude = newreal * newreal + newimag * newimag;
  52.       }
  53. if (color >= maxit) color = inside;
  54. return(color);
  55. }
  56.