home *** CD-ROM | disk | FTP | other *** search
/ ftp.xmission.com / 2014.06.ftp.xmission.com.tar / ftp.xmission.com / pub / lists / fractint / archive / v01.n311 < prev    next >
Internet Message Format  |  1998-10-08  |  42KB
  1. From: owner-fractint-digest@lists.xmission.com (fractint-digest)
  2. To: fractint-digest@lists.xmission.com
  3. Subject: fractint-digest V1 #311
  4. Reply-To: fractint-digest
  5. Sender: owner-fractint-digest@lists.xmission.com
  6. Errors-To: owner-fractint-digest@lists.xmission.com
  7. Precedence: bulk
  8.  
  9.  
  10. fractint-digest        Friday, October 9 1998        Volume 01 : Number 311
  11.  
  12.  
  13.  
  14.  
  15. ----------------------------------------------------------------------
  16.  
  17. Date: Fri, 9 Oct 1998 10:46:01 -0700 (PDT)
  18. From: Ken Childress <kchildre@uccs.jpl.nasa.gov>
  19. Subject: Re: (fractint) Re: [fractal-art] cost of images
  20.  
  21.  
  22. Kerry,
  23.  
  24. > Ken,
  25. > 
  26. > I agree with you in principle, but I think that you can still make good
  27. > enlargements from monitor shots, depending on your equipment, image, and
  28. > standards.  Current 17" monitors deliver 1600 x 1200 resolutions, and a
  29. > good video card and shoot that in true color, elminating any jpeg
  30. > artifacts.
  31.  
  32. I don't mean to say that this isn't a method worth consideration.  It
  33. just isn't worth the effort for me because I can get better results for
  34. less effort.  If you like the results, by all means use it.
  35.  
  36. > Also, actual dpi numbers lose relevance when you're talking
  37. > about putting a picture on the wall where people generally don't get their
  38. > noses into it.  Of course, you're right about making big and/or detailed
  39. > enlargements--you just can't do that from a monitor shot.
  40.  
  41. This is true.  Which is why religious arguments about Epson vs. HP vs.
  42. Canon printers are a waste of time.  Bigger enlargements can suffer less
  43. DPI since you are viewing them from farther away.  However, you can't
  44. take a 1600x1200 pixel image and get a very good 30" by 40" print from
  45. it, IMO.
  46.  
  47. > To print from transparencies like this, you should find a lab that does
  48. > Ilfochrome processing (used to be called, "Cibachrome").  It's a process
  49. > that prints directly from the trans, without an internegative.  There is a
  50. > bunch of other reasons why it's good (like color saturation and archival
  51. > life), but the main one is the "wow" factor.  Expensive, but definitely
  52. > worth it.
  53.  
  54. I agree completely here.  If anyone hasn't seen an Ilfochrome print, you
  55. should make an effort to see the results.  They can be absolutely
  56. stunning.  With the brilliant colors in so many of the fractal images,
  57. an Ilfochrome print would be dazzling.  Definitely worth it, plus the
  58. image will last for a long, long time.
  59.  
  60.  
  61. Ken...
  62.  
  63.  
  64.  
  65.  
  66.  
  67. - --------------------------------------------------------------
  68. Thanks for using Fractint, The Fractals and Fractint Discussion List
  69. Post Message:   fractint@lists.xmission.com
  70. Get Commands:   majordomo@lists.xmission.com "help"
  71. Administrator:  twegner@phoenix.net
  72. Unsubscribe:    majordomo@lists.xmission.com "unsubscribe fractint"
  73.  
  74. ------------------------------
  75.  
  76. Date: Fri, 9 Oct 1998 18:49:42 +0100
  77. From: "Les St Clair" <les_stclair@crosstrees.prestel.co.uk>
  78. Subject: (fractint) September par collection
  79.  
  80. Hi everyone,
  81.  
  82. A little later than usual, but the September par collection is now available.
  83.  
  84. As usual there are two formats:
  85. 1. Just the pars - http://ourworld.compuserve.com/homepages/Les_StClair/
  86. 2. The pars with the text of the original messages -
  87. http://www.homeusers.prestel.co.uk/crosstrees/fractasi.htm
  88.  
  89. Now for a shameless plug...
  90. For some really nice fractal animations (Java applets) check out my new page at
  91. Prestel
  92. http://www.homeusers.prestel.co.uk/crosstrees/cube.htm
  93. (I first saw these on Sharon Webb's site and fell in love with them - thanks
  94. Sharon!).
  95.  
  96. cheers, Les
  97.  
  98.  
  99. - --------------------------------------------------------------
  100. Thanks for using Fractint, The Fractals and Fractint Discussion List
  101. Post Message:   fractint@lists.xmission.com
  102. Get Commands:   majordomo@lists.xmission.com "help"
  103. Administrator:  twegner@phoenix.net
  104. Unsubscribe:    majordomo@lists.xmission.com "unsubscribe fractint"
  105.  
  106. ------------------------------
  107.  
  108. Date: Fri, 9 Oct 1998 11:52:11 -0700 (PDT)
  109. From: Ken Childress <kchildre@uccs.jpl.nasa.gov>
  110. Subject: Re: (fractint) September par collection
  111.  
  112. > 
  113. > Hi everyone,
  114. > 
  115. > A little later than usual, but the September par collection is now available.
  116. > 
  117. > As usual there are two formats:
  118. > 1. Just the pars - http://ourworld.compuserve.com/homepages/Les_StClair/
  119.  
  120. WOW!!!  I have arrived.  :-)
  121.  
  122. I made one of the thumbnail images on Les's page.
  123.  
  124. Thanks Les, I'm truly honored.
  125.  
  126.  
  127.  
  128. Ken...
  129. (I'm not worthy. :-) )
  130.  
  131.  
  132.  
  133. - --------------------------------------------------------------
  134. Thanks for using Fractint, The Fractals and Fractint Discussion List
  135. Post Message:   fractint@lists.xmission.com
  136. Get Commands:   majordomo@lists.xmission.com "help"
  137. Administrator:  twegner@phoenix.net
  138. Unsubscribe:    majordomo@lists.xmission.com "unsubscribe fractint"
  139.  
  140. ------------------------------
  141.  
  142. Date: Fri, 9 Oct 1998 21:00:52 +0100
  143. From: "Kim Bach Petersen" <kimb@post8.tele.dk>
  144. Subject: Sv: (fractint) Re: [fractal-art] cost of images
  145.  
  146. Hi Steve - and all of you other great people outthere,
  147.  
  148. >I would LOVE....I say again LOVE...to hear a detailes description of the
  149. whole
  150. >process....starting with "This is the size I render to print a ??X??
  151. photo."
  152.  
  153.  
  154. I've been doing some 15"x10" digital prints (A3 paper). Many a photo shop
  155. provides such prints, usually made on a digital color photocopying machine,
  156. that can print as well. The machine my local printshop uses, prints at 400
  157. dpi allowing lots of detail. The colors are vivid and clear, I'm very
  158. satisfied with the results. Still, intense colors seems to do best, while
  159. more pale or dusty colors have a lesser "wow-factor".
  160.  
  161. A print cost some 5 $ (30 danish kroner), but if you want multible prints of
  162. the same picture, the price can be reduced greatly, even halfed. The method
  163. is easy, all you need is fractint, a disk and maybe a photoeditor, and you
  164. can make impressive posters.
  165.  
  166. In need of a portable harddisk or the like, I usually stick to 100 dpi
  167. prints, since it allows me to have the file on a standard floppy disk. I
  168. render a picture on 4500x3000 or 6000x4000 and then reduces it to 1500x1000
  169. to anti-alias it. Then I use the photoeditor to specify the wanted size
  170. (15"x10") and resolution (100 dpi), and the saves the ready-to-print file as
  171. jpeg (with lowest compression, unlike gif jpeg allows the specification of
  172. parameters such as size and resolution). Some photoeditors allows all this
  173. to be done in one process. Notice that I use a resolution that is dividable
  174. with the printers resolution (100 = 400 / 4). If I did not, it would
  175. generate aliasing while printing - and it's nicer to avoid.
  176.  
  177. Kind regards, Kim.
  178.  
  179.  
  180.  
  181. - --------------------------------------------------------------
  182. Thanks for using Fractint, The Fractals and Fractint Discussion List
  183. Post Message:   fractint@lists.xmission.com
  184. Get Commands:   majordomo@lists.xmission.com "help"
  185. Administrator:  twegner@phoenix.net
  186. Unsubscribe:    majordomo@lists.xmission.com "unsubscribe fractint"
  187.  
  188. ------------------------------
  189.  
  190. Date: Fri, 9 Oct 1998 13:19:48 -0700 (MST)
  191. From: Kerry Mitchell <lkmitch@primenet.com>
  192. Subject: (fractint) 5 point star coloring methods
  193.  
  194. Fresh from the, "Just because you can doesn't mean that you should" files,
  195. :-) here are some formulas built around 5 point stars
  196.  
  197. - -------------------------------------------------------------------------------
  198. Kerry Mitchell
  199. lkmitch@primenet.com    www.primenet.com/~lkmitch/
  200. - -------------------------------------------------------------------------------
  201.  
  202. comment { ; narrative copyright Kerry Mitchell 09oct98
  203.  
  204. Seeing Stars
  205.  
  206. Most fractal images involve circles in some respect:  either stopping
  207. the iteration when the orbit moves outside of a given circle, or
  208. coloring by how close the orbit comes to a certain circle, or some
  209. variation.  These formulas use a 5 point star instead of a circle.
  210. The "starnear" formulas color by how close the orbit comes to a star.
  211. The "starbail" formulas are set up to stop iterating either when the
  212. orbit comes into a star or when it leaves a star.
  213.  
  214. The star is represented by the 5 outer points.  These are equally spaced
  215. on a circle.  The center and size of the circle are user-specified, as
  216. is the rotation angle of the star.
  217.  
  218. How z at any iteration relates to the star (inside, outside, how close)
  219. is determined by looking at each of the 10 sides, one at a time.  Each
  220. side can be represented by a line A*x + B*y + C = 0, where A, B, and C
  221. come from the coordinates of the 2 outer points that are joined to make
  222. the side.  Given the numbers A, B, and C, the quantity
  223.  
  224.         q = A*real(z) + B*imag(z) + C
  225.  
  226. is computed.  If q is positive, then z is on one side of the line, and
  227. if q is negative, then z is on the other side of the line.  If |q| is
  228. very small, that means that z is very close to that side of the star.
  229. Taking the signs of q for all 10 sides will determine if z is inside or
  230. outside of the star.  Using the smallest value of |q| gives the distance
  231. from z to the star.
  232.  
  233. To see this in action, use the formula "star".  It will draw the outline
  234. of a star, with the position, orientation and size that you choose.
  235. This formula, and the "starnear" formulas that are built on it, use a
  236. logarithmic transfer function from distance to angle (for decomposition
  237. coloring in Fractint).  This helps highlight the star without needing to
  238. find out if the point is exactly on the star or not.  The "star-inout"
  239. formula will draw a star in 2 colors.  Everything inside the star is one
  240. color, and everything outside of the star is the other color.  These 2
  241. formulas are simply to give an idea of how the method works.
  242.  
  243. The fractal formulas need a bit of extra explanation.  In the
  244. "starnear_jul" formula, you have the option of whether or not to check
  245. the pixel (initial value of z) for its distance from the star.  If the
  246. star is centered at the origin (p2=0/0), this can make quite a difference.
  247. For example, use the "starnear2" parameter set.  As written, the real
  248. part of p3 is -0.28. The magnitude of real(p3) is the radius of the circle
  249. containing the 5 outer points.  If real(p3) is negative, that tells the
  250. formula to check the pixel for its distance before iterating.  In the
  251. image, this shows up as a star in the center.  Now, make real(p3) +0.28.
  252. You'll see that the outer portion of the image is unchanged, but the star
  253. in the middle has gone away.  In the "starbail" formulas, the flag to
  254. check the pixel before iterating (pixelflag) has been set in the formula,
  255. and is not interactively changeable (too many parameters).
  256.  
  257. The "starnear" formulas use a standard bailout condition--stop when the
  258. iterate gets too big, or leaves a very large circle.  Conversely, in the
  259. "starbail" formulas, the bailout condition is built around a star.  If 
  260. the star is centered at the origin and is big enough, then you can
  261. essentially use the star instead of the traditional bailout circle.  This
  262. is what has been done in parameter set "starbail1".  The magnification
  263. was decreased to show the bailout star, and the fractal can clearly be
  264. seen in the middle.  It is also possible to bailout when the orbit enters
  265. a star.  This is similar to the orbit trap methods that others have
  266. written, and works particularly well for Julia sets that are dusts.  An
  267. example of this is in parameter set "starbail2".  Some care must be taken
  268. if the star is not centered on the origin.  The "starbail" formulas set
  269. z to a default (1,0) if the orbit was not caught by a star.  If the star
  270. is too far away or too small to catch many orbits, the majority of the
  271. screen will be one color (color 0 when using decomposition coloring).  If
  272. you use stars off of the origin, you may want to increase the size and
  273. set "pixelflag=0" in the formula, to prevent the image from being taken
  274. over by one huge rogue star.
  275.  
  276. One final note:  This method uses angles a lot to figure out which side
  277. belongs to whom.  Consequently, the ability to rotate the star is a bit
  278. limited.  If the rotations are kept to the range of 0 to 36 degrees, there
  279. shouldn't be any problems.  Rotations outside this range (try it with the
  280. "star" formula) will result in branch cuts as the formula gets its angles
  281. confused.  Fortunately, since the star is symmetric, there shouldn't be
  282. a need to rotate it beyond those bounds.
  283.  
  284. }
  285.  
  286. starnear1 { ; copyright Kerry Mitchell 09oct98
  287.   ;
  288.   ; sample of "near to 5 point star" coloring
  289.   ;
  290.   reset=1960 type=formula formulafile=fractint.frm
  291.   formulaname=starnear_jul center-mag=-0.358/0.651/8
  292.   params=0/1/0/0/0.25/18 float=y maxiter=256 inside=0
  293.   decomp=256 periodicity=0 colors=zrM<40>zzzzzzyyz<60>QQzPPzPP\
  294.   zPPy<61>000000110<60>yn0zo0zo0zo1<19>zrL cyclerange=0/255
  295.   }
  296.  
  297. starnear2 { ; copyright Kerry Mitchell 09oct98
  298.   ;
  299.   ; sample of "near to 5 point star" coloring
  300.   ;
  301.   reset=1960 type=formula formulafile=fractint.frm
  302.   formulaname=starnear_jul center-mag=0/0/1.1/1/-90
  303.   params=0.28/0/0/0/-0.28/36 float=y maxiter=256 inside=0
  304.   decomp=256 periodicity=0 colors=vlv<23>i1ih0hf0f<28>303\
  305.   101100<29>x99zAAzBB<30>zzz<30>i2ih0hg0g<29>202000000<30>\
  306.   y99zAAzCC<30>zzz<5>wnw cyclerange=0/255
  307.   }
  308.  
  309. starnear3 { ; copyright Kerry Mitchell 09oct98
  310.   ;
  311.   ;  Mandelbrot sample of coloring by nearest approach to 5 point star
  312.   ;
  313.   reset=1960 type=formula formulafile=fractint.frm
  314.   formulaname=starnear_man center-mag=-1/0/0.56/1/-90
  315.   params=-2/0/0.25/36 float=y maxiter=100 inside=0 decomp=256
  316.   periodicity=0 colors=000<35>x00z00z10<34>zx0zz0zz1<36>zzz<51>\
  317.   zz2zz0zy0<51>y20x00w00<30>200000000000000 cyclerange=0/255
  318.   }
  319.  
  320. starbail1 { ; copyright Kerry Mitchell 09oct98
  321.   ;
  322.   ; bails out when orbit leaves a 5 point star
  323.   ;
  324.   reset=1960 type=formula formulafile=fractint.frm
  325.   formulaname=starbail_jul center-mag=0/0/0.25 cyclerange=0/255
  326.   params=0.3795135923319061/0.3349323012073617/0/0/5/18 float=y
  327.   maxiter=256 inside=0 decomp=256 periodicity=0 colors=000<40>x00z\
  328.   00z00<40>zy0zz0zz1<39>zzxzzzzzz<40>1zz0zz0yz<39>02z00z00z<41>000
  329.   }
  330.  
  331. starbail2 { ; copyright Kerry Mitchell 09oct98
  332.   ;
  333.   ; bails out when orbit enters a 5 point star
  334.   ;
  335.   reset=1960 type=formula formulafile=fractint.frm
  336.   formulaname=starbail_jul passes=1 center-mag=0/0/1.1/1/-90
  337.   params=0.3/0/0/0/-0.5/36 float=y maxiter=256 inside=0
  338.   decomp=256 periodicity=0 colors=000<15>zo0<14>zyvzzzxxz<13>\
  339.   SSzPPzOOw<14>000<15>zo0<14>zywzzzxxz<13>RRzPPzOOv<11>55B337\
  340.   223110<13>tj0xm0zo1<13>zxtzyxyyz<12>VVzTTzRRzPPyNNu<12>3361\
  341.   12220<14>yn0zo2<14>zyyyyz<13>SSzQQzOOx<14>111 cyclerange=0/255
  342.   }
  343.  
  344. starbail3 { ; copyright Kerry Mitchell 09oct98
  345.   ;
  346.   ;  Mandelbrot sample of bailout out when moving inside a 5 point star
  347.   ;
  348.   reset=1960 type=formula formulafile=fractint.frm
  349.   formulaname=starbail_man center-mag=-0.407/0/0.813
  350.   params=1/0/-1/36 float=y maxiter=100 inside=0 decomp=256
  351.   periodicity=0 colors=000<35>x00z00z10<34>zx0zz0zz1<36>zzz<51>\
  352.   zz2zz0zy0<51>y20x00w00<30>200000000000000 cyclerange=0/255
  353.   }
  354.  
  355. frm:starbail_man { ; Kerry Mitchell 09oct98
  356.         ;
  357.         ; Mandelbrot, bails out when orbit enters/leaves 5 point star
  358.         ; p1 = center of star
  359.         ; cabs(real(p2)) = star size of star (try 1)
  360.         ; sign(real(p2)) = in/out flag:
  361.         ;   + = bailout when orbit leaves star
  362.         ;   - = bailout when orbit enters star
  363.         ; imag(p2) = star rotation angle, degrees
  364.         ;   only use angles from 0 to 36 degrees
  365.         ;
  366.         ;  initialize iteration parameters
  367.         ;
  368.         zc=0, c=pixel, done=0, iter=1, bailout=1e12
  369.         pixelflag=0  ; set by hand
  370.         ;
  371.         ;  star parameters
  372.         ;
  373.         center=p1, xcen=real(center), ycen=imag(center)
  374.         r=cabs(real(p2)), inout=1
  375.         if(real(p2)<0)
  376.           inout=0
  377.         endif
  378.         phi=imag(p2)/180*pi, twopi=2*pi, temp=twopi/10
  379.         t0=phi, t1=t0+temp
  380.         t2=t1+temp, t3=t2+temp, t4=t3+temp, t5=t4+temp
  381.         t6=t5+temp, t7=t6+temp, t8=t7+temp, t9=t8+temp
  382.         ;
  383.         ;  set up control points
  384.         ;
  385.         x0=r*cos(t0)+xcen, y0=r*sin(t0)+ycen
  386.         x1=r*cos(t2)+xcen, y1=r*sin(t2)+ycen
  387.         x2=r*cos(t4)+xcen, y2=r*sin(t4)+ycen
  388.         x3=r*cos(t6)+xcen, y3=r*sin(t6)+ycen
  389.         x4=r*cos(t8)+xcen, y4=r*sin(t8)+ycen
  390.         ;
  391.         ;  check pixel to see if it bailed
  392.         ;
  393.         if(pixelflag!=0)
  394.         x=real(zc), y=imag(zc), t=imag(log(zc-center))
  395.         if(t<0)
  396.           t=t+twopi
  397.         endif
  398.         flag=0
  399.         if((t>t0)&&(t<=t1))
  400.           f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  401.           if(f<0)
  402.             flag=1
  403.           endif
  404.         elseif((t>t1)&&(t<=t2))
  405.           f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  406.           if(f>0)
  407.             flag=1
  408.           endif
  409.         elseif((t>t2)&&(t<=t3))
  410.           f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  411.           if(f<0)
  412.             flag=1
  413.           endif
  414.         elseif((t>t3)&&(t<=t4))
  415.           f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  416.           if(f>0)
  417.             flag=1
  418.           endif
  419.         elseif((t>t4)&&(t<=t5))
  420.           f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  421.           if(f<0)
  422.             flag=1
  423.           endif
  424.         elseif((t>t5)&&(t<=t6))
  425.           f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  426.           if(f>0)
  427.             flag=1
  428.           endif
  429.         elseif((t>t6)&&(t<=t7))
  430.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  431.           if(f<0)
  432.             flag=1
  433.           endif
  434.         elseif((t>t7)&&(t<=t8))
  435.           f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  436.           if(f>0)
  437.             flag=1
  438.           endif
  439.         elseif((t>t8)&&(t<=t9))
  440.           f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  441.           if(f<0)
  442.             flag=1
  443.           endif
  444.         else
  445.           f=x*(y0-y3)+y*(x3-x0)-x3*y0+x0*y3
  446.           if(f>0)
  447.             flag=1
  448.           endif
  449.         endif
  450.         if(flag==inout)
  451.           done=1
  452.           z=zc-center
  453.         endif
  454.         endif
  455.         :
  456.         ;  standard iteration if pixel didn't bail
  457.         ;
  458.         if(done==0)
  459.         iter=iter+1, zc=sqr(zc)+c
  460.         x=real(zc), y=imag(zc), t=imag(log(zc-center))
  461.         if(t<0)
  462.           t=t+twopi
  463.         endif
  464.         ;
  465.         ;  see if orbit is inside or outside of star
  466.         ;    by checking each side one at a time
  467.         ;
  468.         flag=0
  469.         if((t>t0)&&(t<=t1))
  470.           f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  471.           if(f<0)
  472.             flag=1
  473.           endif
  474.         elseif((t>t1)&&(t<=t2))
  475.           f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  476.           if(f>0)
  477.             flag=1
  478.           endif
  479.         elseif((t>t2)&&(t<=t3))
  480.           f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  481.           if(f<0)
  482.             flag=1
  483.           endif
  484.         elseif((t>t3)&&(t<=t4))
  485.           f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  486.           if(f>0)
  487.             flag=1
  488.           endif
  489.         elseif((t>t4)&&(t<=t5))
  490.           f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  491.           if(f<0)
  492.             flag=1
  493.           endif
  494.         elseif((t>t5)&&(t<=t6))
  495.           f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  496.           if(f>0)
  497.             flag=1
  498.           endif
  499.         elseif((t>t6)&&(t<=t7))
  500.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  501.           if(f<0)
  502.             flag=1
  503.           endif
  504.         elseif((t>t7)&&(t<=t8))
  505.           f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  506.           if(f>0)
  507.             flag=1
  508.           endif
  509.         elseif((t>t8)&&(t<=t9))
  510.           f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  511.           if(f<0)
  512.             flag=1
  513.           endif
  514.         else
  515.           f=x*(y0-y3)+y*(x3-x0)-x3*y0+x0*y3
  516.           if(f>0)
  517.             flag=1
  518.           endif
  519.         endif
  520.         ;
  521.         ;  if the orbit was on the appropriate side,
  522.         ;    set "done" flag
  523.         ;    set z to iteration variable for coloring
  524.         ;
  525.         if(flag==inout)
  526.           done=1
  527.           z=zc-center
  528.         ;
  529.         ;  safety valve in case star doesn't catch orbit
  530.         ;  if orbit goes to infinity or max iterations reached:
  531.         ;     set "done" flag
  532.         ;     set z to 1 to differentiate it from star trap
  533.         ;
  534.         elseif((|zc|>bailout)||(iter==maxit))
  535.           done=1
  536.           z=1
  537.         endif
  538.         endif
  539.         done==0
  540.         }
  541.  
  542. frm:starbail_jul { ; Kerry Mitchell 09oct98
  543.         ;
  544.         ; Julia set, bails out when orbit enters/leaves 5 point star
  545.         ; p1 = Julia parameter
  546.         ; p2 = center of star
  547.         ; cabs(real(p3)) = star size of star (try 1)
  548.         ; sign(real(p3)) = in/out flag:
  549.         ;   + = bailout when orbit leaves star
  550.         ;   - = bailout when orbit enters star
  551.         ; imag(p3) = star rotation angle, degrees
  552.         ;   only use angles from 0 to 36 degrees
  553.         ;
  554.         ;  initialize iteration parameters
  555.         ;
  556.         c=p1, zc=pixel, done=0, iter=1, bailout=1e12
  557.         pixelflag=1  ; set by hand
  558.         ;
  559.         ;  star parameters
  560.         ;
  561.         center=p2, xcen=real(center), ycen=imag(center)
  562.         r=cabs(real(p3)), inout=1
  563.         if(real(p3)<0)
  564.           inout=0
  565.         endif
  566.         phi=imag(p3)/180*pi, twopi=2*pi, temp=twopi/10
  567.         t0=phi, t1=t0+temp
  568.         t2=t1+temp, t3=t2+temp, t4=t3+temp, t5=t4+temp
  569.         t6=t5+temp, t7=t6+temp, t8=t7+temp, t9=t8+temp
  570.         ;
  571.         ;  set up control points
  572.         ;
  573.         x0=r*cos(t0)+xcen, y0=r*sin(t0)+ycen
  574.         x1=r*cos(t2)+xcen, y1=r*sin(t2)+ycen
  575.         x2=r*cos(t4)+xcen, y2=r*sin(t4)+ycen
  576.         x3=r*cos(t6)+xcen, y3=r*sin(t6)+ycen
  577.         x4=r*cos(t8)+xcen, y4=r*sin(t8)+ycen
  578.         ;
  579.         ;  check pixel to see if it bailed
  580.         ;
  581.         if(pixelflag!=0)
  582.         x=real(zc), y=imag(zc), t=imag(log(zc-center))
  583.         if(t<0)
  584.           t=t+twopi
  585.         endif
  586.         flag=0
  587.         if((t>t0)&&(t<=t1))
  588.           f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  589.           if(f<0)
  590.             flag=1
  591.           endif
  592.         elseif((t>t1)&&(t<=t2))
  593.           f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  594.           if(f>0)
  595.             flag=1
  596.           endif
  597.         elseif((t>t2)&&(t<=t3))
  598.           f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  599.           if(f<0)
  600.             flag=1
  601.           endif
  602.         elseif((t>t3)&&(t<=t4))
  603.           f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  604.           if(f>0)
  605.             flag=1
  606.           endif
  607.         elseif((t>t4)&&(t<=t5))
  608.           f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  609.           if(f<0)
  610.             flag=1
  611.           endif
  612.         elseif((t>t5)&&(t<=t6))
  613.           f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  614.           if(f>0)
  615.             flag=1
  616.           endif
  617.         elseif((t>t6)&&(t<=t7))
  618.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  619.           if(f<0)
  620.             flag=1
  621.           endif
  622.         elseif((t>t7)&&(t<=t8))
  623.           f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  624.           if(f>0)
  625.             flag=1
  626.           endif
  627.         elseif((t>t8)&&(t<=t9))
  628.           f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  629.           if(f<0)
  630.             flag=1
  631.           endif
  632.         else
  633.           f=x*(y0-y3)+y*(x3-x0)-x3*y0+x0*y3
  634.           if(f>0)
  635.             flag=1
  636.           endif
  637.         endif
  638.         if(flag==inout)
  639.           done=1
  640.           z=zc-center
  641.         endif
  642.         endif
  643.         :
  644.         ;  standard iteration if pixel didn't bail
  645.         ;
  646.         if(done==0)
  647.         iter=iter+1, zc=sqr(zc)+c
  648.         x=real(zc), y=imag(zc), t=imag(log(zc-center))
  649.         if(t<0)
  650.           t=t+twopi
  651.         endif
  652.         ;
  653.         ;  see if orbit is inside or outside of star
  654.         ;    by checking each side one at a time
  655.         ;
  656.         flag=0
  657.         if((t>t0)&&(t<=t1))
  658.           f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  659.           if(f<0)
  660.             flag=1
  661.           endif
  662.         elseif((t>t1)&&(t<=t2))
  663.           f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  664.           if(f>0)
  665.             flag=1
  666.           endif
  667.         elseif((t>t2)&&(t<=t3))
  668.           f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  669.           if(f<0)
  670.             flag=1
  671.           endif
  672.         elseif((t>t3)&&(t<=t4))
  673.           f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  674.           if(f>0)
  675.             flag=1
  676.           endif
  677.         elseif((t>t4)&&(t<=t5))
  678.           f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  679.           if(f<0)
  680.             flag=1
  681.           endif
  682.         elseif((t>t5)&&(t<=t6))
  683.           f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  684.           if(f>0)
  685.             flag=1
  686.           endif
  687.         elseif((t>t6)&&(t<=t7))
  688.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  689.           if(f<0)
  690.             flag=1
  691.           endif
  692.         elseif((t>t7)&&(t<=t8))
  693.           f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  694.           if(f>0)
  695.             flag=1
  696.           endif
  697.         elseif((t>t8)&&(t<=t9))
  698.           f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  699.           if(f<0)
  700.             flag=1
  701.           endif
  702.         else
  703.           f=x*(y0-y3)+y*(x3-x0)-x3*y0+x0*y3
  704.           if(f>0)
  705.             flag=1
  706.           endif
  707.         endif
  708.         ;
  709.         ;  if the orbit was on the appropriate side,
  710.         ;    set "done" flag
  711.         ;    set z to iteration variable for coloring
  712.         ;
  713.         if(flag==inout)
  714.           done=1
  715.           z=zc-center
  716.         ;
  717.         ;  safety valve in case star doesn't catch orbit
  718.         ;  if orbit goes to infinity or max iterations reached:
  719.         ;     set "done" flag
  720.         ;     set z to 1 to differentiate it from star trap
  721.         ;
  722.         elseif((|zc|>bailout)||(iter==maxit))
  723.           done=1
  724.           z=1
  725.         endif
  726.         endif
  727.         done==0
  728.         }
  729.  
  730. frm:starnear_man { ; Kerry Mitchell 09oct98
  731.         ;
  732.         ; Mandelbrot, colors by nearest approach to 5 point star
  733.         ; p1 = center of star
  734.         ; cabs(real(p2)) = star size of star (try 1)
  735.         ; sign(real(p2)) = pixel flag:
  736.         ;   + = don't consider pixel--start after first iteration
  737.         ;   - = consider pixel
  738.         ; imag(p2) = star rotation angle, degrees
  739.         ;   only use angles from 0 to 36 degrees
  740.         ; use decomp=256, bailout hardcoded to 10^12
  741.         ;
  742.         ;  initialize iteration parameters
  743.         ;
  744.         zc=0, c=pixel, done=0, iter=1
  745.         bailout=1e12, fmin=bailout
  746.         ;
  747.         ;  star parameters
  748.         ;
  749.         center=p1, xcen=real(center), ycen=imag(center)
  750.         r=cabs(real(p2)), pixelflag=0
  751.         if(real(p2)<0)
  752.           pixelflag=1
  753.         endif
  754.         phi=imag(p2)/180*pi, twopi=2*pi, temp=twopi/10
  755.         t0=phi, t1=t0+temp
  756.         t2=t1+temp, t3=t2+temp, t4=t3+temp, t5=t4+temp
  757.         t6=t5+temp, t7=t6+temp, t8=t7+temp, t9=t8+temp
  758.         ;
  759.         ;  set up control points
  760.         ;
  761.         x0=r*cos(t0)+xcen, y0=r*sin(t0)+ycen
  762.         x1=r*cos(t2)+xcen, y1=r*sin(t2)+ycen
  763.         x2=r*cos(t4)+xcen, y2=r*sin(t4)+ycen
  764.         x3=r*cos(t6)+xcen, y3=r*sin(t6)+ycen
  765.         x4=r*cos(t8)+xcen, y4=r*sin(t8)+ycen
  766.         ;
  767.         ;  if pixel flag =/= 0, then check pixel for how close
  768.         ;    it is to star
  769.         ;
  770.         if(pixelflag!=0)
  771.           x=real(zc), y=imag(zc), t=imag(log(zc-center))
  772.           if(t<0)
  773.             t=t+twopi
  774.           endif
  775.           if((t>t0)&&(t<=t1))
  776.             f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  777.           elseif((t>t1)&&(t<=t2))
  778.             f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  779.           elseif((t>t2)&&(t<=t3))
  780.             f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  781.           elseif((t>t3)&&(t<=t4))
  782.             f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  783.           elseif((t>t4)&&(t<=t5))
  784.             f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  785.           elseif((t>t5)&&(t<=t6))
  786.             f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  787.           elseif((t>t6)&&(t<=t7))
  788.             f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  789.           elseif((t>t7)&&(t<=t8))
  790.             f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  791.           elseif((t>t8)&&(t<=t9))
  792.             f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  793.           else
  794.             f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  795.           endif
  796.           fmin=cabs(f)
  797.         endif
  798.         :
  799.         ;  standard iteration, find polar angle of iterate
  800.         ;
  801.         iter=iter+1, zc=sqr(zc)+c
  802.         x=real(zc), y=imag(zc), t=imag(log(zc-center))
  803.         if(t<0)
  804.           t=t+twopi
  805.         endif
  806.         ;
  807.         ;  compute how close iterate is to each side of star
  808.         ;
  809.         if((t>t0)&&(t<=t1))
  810.           f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  811.         elseif((t>t1)&&(t<=t2))
  812.           f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  813.         elseif((t>t2)&&(t<=t3))
  814.           f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  815.         elseif((t>t3)&&(t<=t4))
  816.           f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  817.         elseif((t>t4)&&(t<=t5))
  818.           f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  819.         elseif((t>t5)&&(t<=t6))
  820.           f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  821.         elseif((t>t6)&&(t<=t7))
  822.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  823.         elseif((t>t7)&&(t<=t8))
  824.           f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  825.         elseif((t>t8)&&(t<=t9))
  826.           f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  827.         else
  828.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  829.         endif
  830.         f=cabs(f)
  831.         ;
  832.         ;  update minimum distance
  833.         ;
  834.         if(f<fmin)
  835.           fmin=f
  836.           zmin=zc
  837.         endif
  838.         ;
  839.         ;  upon escape or maximum iterations:
  840.         ;    set "done" flag
  841.         ;    use minimum distance from star as polar angle of z
  842.         ;      for use with decomp coloring
  843.         ;
  844.         if((|zc|>bailout)||(iter==maxit))
  845.           done=1
  846.           t=log(fmin)
  847.           z=cos(t)+flip(sin(t))
  848.         endif
  849.         done==0
  850.         }
  851.  
  852. frm:starnear_jul { ; Kerry Mitchell 09oct98
  853.         ;
  854.         ; Julia set, colors by nearest approach to 5 point star
  855.         ; p1 = Julia parameter
  856.         ; p2 = center of star
  857.         ; cabs(real(p3)) = star size of star (try 1)
  858.         ; sign(real(p3)) = pixel flag:
  859.         ;   + = don't consider pixel--start after first iteration
  860.         ;   - = consider pixel
  861.         ; imag(p3) = star rotation angle, degrees
  862.         ;   only use angles from 0 to 36 degrees
  863.         ; use decomp=256, bailout hardcoded to 10^12
  864.         ;
  865.         ;  initialize iteration parameters
  866.         ;
  867.         c=p1, zc=pixel, done=0, iter=1
  868.         bailout=1e12, fmin=bailout
  869.         ;
  870.         ;  star parameters
  871.         ;
  872.         center=p2, xcen=real(center), ycen=imag(center)
  873.         r=cabs(real(p3)), pixelflag=0
  874.         if(real(p3)<0)
  875.           pixelflag=1
  876.         endif
  877.         phi=imag(p3)/180*pi, twopi=2*pi, temp=twopi/10
  878.         t0=phi, t1=t0+temp
  879.         t2=t1+temp, t3=t2+temp, t4=t3+temp, t5=t4+temp
  880.         t6=t5+temp, t7=t6+temp, t8=t7+temp, t9=t8+temp
  881.         ;
  882.         ;  set up control points
  883.         ;
  884.         x0=r*cos(t0)+xcen, y0=r*sin(t0)+ycen
  885.         x1=r*cos(t2)+xcen, y1=r*sin(t2)+ycen
  886.         x2=r*cos(t4)+xcen, y2=r*sin(t4)+ycen
  887.         x3=r*cos(t6)+xcen, y3=r*sin(t6)+ycen
  888.         x4=r*cos(t8)+xcen, y4=r*sin(t8)+ycen
  889.         ;
  890.         ;  if pixel flag =/= 0, then check pixel for how close
  891.         ;    it is to star
  892.         ;
  893.         if(pixelflag!=0)
  894.           x=real(zc), y=imag(zc), t=imag(log(zc-center))
  895.           if(t<0)
  896.             t=t+twopi
  897.           endif
  898.           if((t>t0)&&(t<=t1))
  899.             f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  900.           elseif((t>t1)&&(t<=t2))
  901.             f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  902.           elseif((t>t2)&&(t<=t3))
  903.             f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  904.           elseif((t>t3)&&(t<=t4))
  905.             f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  906.           elseif((t>t4)&&(t<=t5))
  907.             f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  908.           elseif((t>t5)&&(t<=t6))
  909.             f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  910.           elseif((t>t6)&&(t<=t7))
  911.             f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  912.           elseif((t>t7)&&(t<=t8))
  913.             f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  914.           elseif((t>t8)&&(t<=t9))
  915.             f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  916.           else
  917.             f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  918.           endif
  919.           fmin=cabs(f)
  920.         endif
  921.         :
  922.         ;  standard iteration, find polar angle of iterate
  923.         ;
  924.         iter=iter+1, zc=sqr(zc)+c
  925.         x=real(zc), y=imag(zc), t=imag(log(zc-center))
  926.         if(t<0)
  927.           t=t+twopi
  928.         endif
  929.         ;
  930.         ;  compute how close iterate is to each side of star
  931.         ;
  932.         if((t>t0)&&(t<=t1))
  933.           f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  934.         elseif((t>t1)&&(t<=t2))
  935.           f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  936.         elseif((t>t2)&&(t<=t3))
  937.           f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  938.         elseif((t>t3)&&(t<=t4))
  939.           f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  940.         elseif((t>t4)&&(t<=t5))
  941.           f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  942.         elseif((t>t5)&&(t<=t6))
  943.           f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  944.         elseif((t>t6)&&(t<=t7))
  945.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  946.         elseif((t>t7)&&(t<=t8))
  947.           f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  948.         elseif((t>t8)&&(t<=t9))
  949.           f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  950.         else
  951.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  952.         endif
  953.         f=cabs(f)
  954.         ;
  955.         ;  update minimum distance
  956.         ;
  957.         if(f<fmin)
  958.           fmin=f
  959.           zmin=zc
  960.         endif
  961.         ;
  962.         ;  upon escape or maximum iterations:
  963.         ;    set "done" flag
  964.         ;    use minimum distance from star as polar angle of z
  965.         ;      for use with decomp coloring
  966.         ;
  967.         if((|zc|>bailout)||(iter==maxit))
  968.           done=1
  969.           t=log(fmin)
  970.           z=cos(t)+flip(sin(t))
  971.         endif
  972.         done==0
  973.         }
  974.  
  975. frm:star-inout { ; Kerry Mitchell 09oct98
  976.         ;
  977.         ; draws 5 point star--1 color inside, 1 color outside
  978.         ; p1 = center of star
  979.         ; real(p2) = size of star (try 1)
  980.         ; imag(p2) = rotation angle, degrees (0 for point at top)
  981.         ; use decomp=256
  982.         ;
  983.         zc=pixel, done=1
  984.         center=p1, xcen=real(center), ycen=imag(center), r=real(p2)
  985.         phi=(imag(p2)+18)/180*pi, twopi=2*pi, temp=twopi/10
  986.         if((phi<0)||(phi>=temp))
  987.           phi=0
  988.         endif
  989.         t0=phi, t1=temp+phi, t2=2*temp+phi, t3=3*temp+phi
  990.         t4=4*temp+phi, t5=5*temp+phi, t6=6*temp+phi
  991.         t7=7*temp+phi, t8=8*temp+phi, t9=9*temp+phi
  992.         x0=r*cos(t0)+xcen, y0=r*sin(t0)+ycen
  993.         x1=r*cos(t2)+xcen, y1=r*sin(t2)+ycen
  994.         x2=r*cos(t4)+xcen, y2=r*sin(t4)+ycen
  995.         x3=r*cos(t6)+xcen, y3=r*sin(t6)+ycen
  996.         x4=r*cos(t8)+xcen, y4=r*sin(t8)+ycen
  997.         :
  998.         x=real(zc), y=imag(zc), t=imag(log(zc-center))
  999.         if(t<0)
  1000.           t=t+twopi
  1001.         endif
  1002.         flag=0
  1003.         if((t>t0)&&(t<=t1))
  1004.           f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  1005.           if(f<0)
  1006.             flag=1
  1007.           endif
  1008.         elseif((t>t1)&&(t<=t2))
  1009.           f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  1010.           if(f>0)
  1011.             flag=1
  1012.           endif
  1013.         elseif((t>t2)&&(t<=t3))
  1014.           f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  1015.           if(f<0)
  1016.             flag=1
  1017.           endif
  1018.         elseif((t>t3)&&(t<=t4))
  1019.           f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  1020.           if(f>0)
  1021.             flag=1
  1022.           endif
  1023.         elseif((t>t4)&&(t<=t5))
  1024.           f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  1025.           if(f<0)
  1026.             flag=1
  1027.           endif
  1028.         elseif((t>t5)&&(t<=t6))
  1029.           f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  1030.           if(f>0)
  1031.             flag=1
  1032.           endif
  1033.         elseif((t>t6)&&(t<=t7))
  1034.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  1035.           if(f<0)
  1036.             flag=1
  1037.           endif
  1038.         elseif((t>t7)&&(t<=t8))
  1039.           f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  1040.           if(f>0)
  1041.             flag=1
  1042.           endif
  1043.         elseif((t>t8)&&(t<=t9))
  1044.           f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  1045.           if(f<0)
  1046.             flag=1
  1047.           endif
  1048.         else
  1049.           f=x*(y0-y3)+y*(x3-x0)-x3*y0+x0*y3
  1050.           if(f>0)
  1051.             flag=1
  1052.           endif
  1053.         endif
  1054.         z=0
  1055.         if(flag==1)
  1056.           z=(0.0,1.0)
  1057.         endif
  1058.         done==0
  1059.         }
  1060.  
  1061. frm:star { ; Kerry Mitchell 09oct98
  1062.         ;
  1063.         ; draws 5-point star, not a fractal
  1064.         ;
  1065.         ; p1 = center of star
  1066.         ; real(p2) = size of star (try 1)
  1067.         ; imag(p2) = rotation angle, degrees (0 for point at top)
  1068.         ; use decomp=256
  1069.         ;
  1070.         zc=pixel, done=1
  1071.         center=p1, xcen=real(center), ycen=imag(center), r=real(p2)
  1072.         phi=imag(p2)/180*pi, twopi=2*pi, temp=twopi/10
  1073.         t0=phi, t1=temp+phi, t2=2*temp+phi, t3=3*temp+phi
  1074.         t4=4*temp+phi, t5=5*temp+phi, t6=6*temp+phi
  1075.         t7=7*temp+phi, t8=8*temp+phi, t9=9*temp+phi
  1076.         x0=r*cos(t0)+xcen, y0=r*sin(t0)+ycen
  1077.         x1=r*cos(t2)+xcen, y1=r*sin(t2)+ycen
  1078.         x2=r*cos(t4)+xcen, y2=r*sin(t4)+ycen
  1079.         x3=r*cos(t6)+xcen, y3=r*sin(t6)+ycen
  1080.         x4=r*cos(t8)+xcen, y4=r*sin(t8)+ycen
  1081.         :
  1082.         x=real(zc), y=imag(zc), t=imag(log(zc-center))
  1083.         if(t<0)
  1084.           t=t+twopi
  1085.         endif
  1086.         if((t>t0)&&(t<=t1))
  1087.           f=x*(y0-y2)+y*(x2-x0)-x2*y0+x0*y2
  1088.         elseif((t>t1)&&(t<=t2))
  1089.           f=x*(y1-y4)+y*(x4-x1)-x4*y1+x1*y4
  1090.         elseif((t>t2)&&(t<=t3))
  1091.           f=x*(y1-y3)+y*(x3-x1)-x3*y1+x1*y3
  1092.         elseif((t>t3)&&(t<=t4))
  1093.           f=x*(y2-y0)+y*(x0-x2)-x0*y2+x2*y0
  1094.         elseif((t>t4)&&(t<=t5))
  1095.           f=x*(y2-y4)+y*(x4-x2)-x4*y2+x2*y4
  1096.         elseif((t>t5)&&(t<=t6))
  1097.           f=x*(y3-y1)+y*(x1-x3)-x1*y3+x3*y1
  1098.         elseif((t>t6)&&(t<=t7))
  1099.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  1100.         elseif((t>t7)&&(t<=t8))
  1101.           f=x*(y4-y2)+y*(x2-x4)-x2*y4+x4*y2
  1102.         elseif((t>t8)&&(t<=t9))
  1103.           f=x*(y4-y1)+y*(x1-x4)-x1*y4+x4*y1
  1104.         else
  1105.           f=x*(y3-y0)+y*(x0-x3)-x0*y3+x3*y0
  1106.         endif
  1107.         t=log(cabs(f))
  1108.         z=cos(t)+flip(sin(t))
  1109.         done==0
  1110.         }
  1111.  
  1112.  
  1113.  
  1114. - --------------------------------------------------------------
  1115. Thanks for using Fractint, The Fractals and Fractint Discussion List
  1116. Post Message:   fractint@lists.xmission.com
  1117. Get Commands:   majordomo@lists.xmission.com "help"
  1118. Administrator:  twegner@phoenix.net
  1119. Unsubscribe:    majordomo@lists.xmission.com "unsubscribe fractint"
  1120.  
  1121. ------------------------------
  1122.  
  1123. Date: Fri, 9 Oct 1998 17:41:38 EDT
  1124. From: SKarl52884@aol.com
  1125. Subject: Re: (fractint) Re: [fractal-art] cost of images
  1126.  
  1127. In a message dated 10/9/98 1:44:48 PM Eastern Daylight Time,
  1128. kchildre@uccs.jpl.nasa.gov writes:
  1129.  
  1130. << The best way I know to do
  1131.  that is to use slides for photographs for underwater photographs, and
  1132.  go from digital to a print for digital images. >>
  1133.  
  1134. Hi there...
  1135.  
  1136. Which brings to mind the digital camera.
  1137.  
  1138. My brother is bringing it over tonight.
  1139. It produces jpg and he claims that to get prints on paper all one needs to do
  1140. is deliver them to the right service agency.
  1141. Appearently the files aren't huge as he sent me a few shots of the of the
  1142. mountains last week while on vacation. They seemed to be at
  1143. 1280x768...downloaded fast
  1144. and looked very very good!
  1145. Steve
  1146.  
  1147.  
  1148.  
  1149.  
  1150.  
  1151.  
  1152.  
  1153.  
  1154. - --------------------------------------------------------------
  1155. Thanks for using Fractint, The Fractals and Fractint Discussion List
  1156. Post Message:   fractint@lists.xmission.com
  1157. Get Commands:   majordomo@lists.xmission.com "help"
  1158. Administrator:  twegner@phoenix.net
  1159. Unsubscribe:    majordomo@lists.xmission.com "unsubscribe fractint"
  1160.  
  1161. ------------------------------
  1162.  
  1163. Date: Fri, 9 Oct 1998 17:51:19 EDT
  1164. From: SKarl52884@aol.com
  1165. Subject: Re: Sv: (fractint) Re: [fractal-art] cost of images/videomode?
  1166.  
  1167. In a message dated 10/9/98 3:06:47 PM Eastern Daylight Time,
  1168. kimb@post8.tele.dk writes:
  1169.  
  1170. << I
  1171.  render a picture on 4500x3000 or 6000x4000 and then reduces it to 1500x1000
  1172.  to anti-alias it.  >>
  1173.  
  1174. Hi Kim
  1175.  
  1176. Have i missed a fractint up-date???  The highest videomode I have
  1177. [in 19.6.... disk ram ] is 2048x2048 or 2048x1536.
  1178. How are you getting 4500x3000 and 6000x4000?
  1179. Steve
  1180.  
  1181. - --------------------------------------------------------------
  1182. Thanks for using Fractint, The Fractals and Fractint Discussion List
  1183. Post Message:   fractint@lists.xmission.com
  1184. Get Commands:   majordomo@lists.xmission.com "help"
  1185. Administrator:  twegner@phoenix.net
  1186. Unsubscribe:    majordomo@lists.xmission.com "unsubscribe fractint"
  1187.  
  1188. ------------------------------
  1189.  
  1190. Date: Fri, 09 Oct 1998 17:02:28 -0500
  1191. From: "Paul N. Lee" <Paul.N.Lee@Worldnet.att.net>
  1192. Subject: Re: Sv: (fractint) Re: [fractal-art] cost of images/videomode?
  1193.  
  1194. SKarl52884@aol.com wrote:
  1195. >
  1196. > Have i missed a fractint up-date???  The highest videomode
  1197. > I have [in 19.6.... disk ram ] is 2048x2048 or 2048x1536.
  1198. > How are you getting 4500x3000 and 6000x4000?
  1199.  
  1200. Fractint still generates images at 2048x2048, but if you want higher
  1201. resolution images for printing quality prints, then use Fractint's
  1202. "pieces" divide-and-conquer feature.  You can create multiple PAR
  1203. entries that break an image up into pieces so that you can generate the
  1204. image pieces one by one.  There are two reasons for doing this.  The
  1205. first is in case the fractal is very slow, and you want to generate
  1206. parts of the image at the same time on several computers.  The second is
  1207. that you might want to make an image greater than 2048x2048.  The 
  1208. parameters for this feature are:
  1209.     X Multiples - How many divisions of final image in the x direction
  1210.     Y Multiples - How many divisions of final image in the y direction
  1211.     Video mode  - Fractint video mode for each piece (e.g. "F3")
  1212.  
  1213. The last item defaults to the current video mode.  If either X Multiples
  1214. or Y Multiples are greater than 1, then multiple numbered PAR entries
  1215. for the pieces are added to the PAR file, and a MAKEMIG.BAT file is
  1216. created that builds all of the component pieces and then stitches them
  1217. together into a "multi-image" GIF.  The current limitations of the
  1218. "divide and conquer" algorithm are 36 or fewer X and Y multiples (so you
  1219. are limited to "only" 36x36=1296 component images), and a final
  1220. resolution limit in both the X and Y directions of 65,535 (a limitation
  1221. of "only" four billion pixels or so).
  1222.  
  1223. The final image generated by MAKEMIG is a "multi-image" GIF file called
  1224. FRACTMIG.GIF.  In case you have other software that can't handle
  1225. multi-image GIF files, MAKEMIG includes a final (but commented out) call
  1226. to SIMPLGIF, a companion program that reads a GIF file that may contain
  1227. little tricks like multiple images and creates a simple GIF from it. 
  1228. Fair warning: SIMPLGIF needs room to build a composite image while it
  1229. works, and it does that using a temporary disk file equal to the size of
  1230. the final image - and a 64Kx64K GIF image requires a 4GB temporary disk
  1231. file!
  1232.  
  1233. P.N.L.
  1234. - -------------------------------------------------
  1235. http://www.fractalus.com/cgi-bin/theway?ring=fractals&id=43&go
  1236.  
  1237. - --------------------------------------------------------------
  1238. Thanks for using Fractint, The Fractals and Fractint Discussion List
  1239. Post Message:   fractint@lists.xmission.com
  1240. Get Commands:   majordomo@lists.xmission.com "help"
  1241. Administrator:  twegner@phoenix.net
  1242. Unsubscribe:    majordomo@lists.xmission.com "unsubscribe fractint"
  1243.  
  1244. ------------------------------
  1245.  
  1246. End of fractint-digest V1 #311
  1247. ******************************
  1248.  
  1249.