M"Gc@sjdZdZdZyeWnd@\ZZnXyddkZddkZddkZddkZddk Z ei i ei i e i eiddd eeei Zdd klZlZWn1yeiid eWn eGHnXnXei ieiddZei iedZeid ieiZd ZdZ dZhdd<dd<dd<dd<dd<dd<eid<dd<dd<dd<d d!e=o3e?d<e=doeidZ<eid=Z=nWq:q:Xne?d>e<ee@e=fpd?GHndS(Asbsimple program to create Mandelbrot set this takes command-line arguments to create a PPM files- mandelgen -- generate image of Mandelbrot set Copyright (C) 2005 John Comeau This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. s7Not all needed libraries found, upgrade or check path: iiiNitlibtpython(tgplt jclicenses%s s [0-9A-Za-z]+cGstS(N(tFalse(twhatever((s./mandelgen.pyt DebugPrint2scGs|\}}t|tijot|}nttdo|i|Sn^d}xPtdt|D]9}|||}|t|djo||}qoqoW|SdS(sfor pythons without str.jointjointiiN( ttypettypest StringTypetevalthasattrtstrRtrangetlen(targststringtarraytjoinedtindex((s./mandelgen.pyR5s &cGsg|\}}t|p d}ntddo|i|Sn titi|i|SdS(sfor pythons without str.splitRtsplitN(RtNoneR Rtretcompiletescape(RRtstring_to_split((s./mandelgen.pyRCs  gtxgtyi@twidthitheighttratiog?tzoomtoutputtcolorst maxiterationsitmaxvalueit ppm_maxvalc Gs]td|t|ttdtdtdtdtdtdtdtd S( smain routine creates PPM image of Mandelbrot set at a particular location and zoom level each pixel is z -> z^2 + z0 RtcurrentRRRRR R!N(Rtinitt ppmheaderR't gen_by_pixel(R((s./mandelgen.pyt mandelgen]s   cGstd}}tiidtdodtd|i d|i dtd jo|d 7}n t d WqPqXqWt ||o)|tddjotiid n|d 7}q;W|otiid ndS(Nis#Showing one dot per line processed sgen_by_pixel() startingiR$RRiR%itescapedt.s (tTruetsyststderrtwriteRR'tinttcomplext widthvaluet heightvaluetrealtimagt Exceptiontppmpixel(Rt showprogressRRRtctzt iterations((s./mandelgen.pyR*js*  #2#   cCs;||\}}\}}||||||||gS(N((titjtatbR;td((s./mandelgen.pytcomplex_multiplyscCs+||\}}\}}||||gS(N((R>R?R@RAR;RB((s./mandelgen.pyt complex_addscGs4tdtd}|d}td||tdS(NttoptbottomiR(R'(Rt heightspant pixelheight((s./mandelgen.pyR5s cGs4tdtd}|d}td||tdS(NtrighttleftiR(R'(Rt widthspant pixelwidth((s./mandelgen.pyR4s cCs7ttttd}t|}td||S(NR#tmaxval(R(tmaptmaxR'R(tmaxrgbRM((s./mandelgen.pyt find_maxvals   cGstddjoot|djo|d}ntd|t|djo-t|dt|dtd! --    "cCs5d}x(|D] }|dK}||djO}q W|S(Niit1((t binary_stringtnumbert character((s./mandelgen.pytbinarys  cGs dddgS(s6simple 3-pass coloring scheme to match WikiHow articleii(iii(iii(iii((R((s./mandelgen.pytrgbscGsvg}g}tdddD]1}||ddfd|dfdd|fgq~}x|D]}||7}q^W|S(s.alternate r, g, b darker and darker to 0, 0, 0iii(R(Rtvaluest_[1]tcolorR#ttriplet((s./mandelgen.pyt rgb2blacksJcGsbtd}td}g}t|ddD]&}|d||d>||d>fq/~}|S(spdefine start color as cyan: r=0 g>=b decrement for each iteration such that 0 is black to match the colorForth program, g must be 6 bits and b 5, with the low bits of each cleared, and two subtracted each time also, found out that the g and b values must be shifted into the high-order bits of their respective bytes, to match colorForth results t 11111011110t100000iiii(R]R(RtcyantdivisorR`RaR_((s./mandelgen.pyRTs  ?cGs:g}tdddD]}|dd|fq~}|S(suse blue values from 1023 to 0iii(R(RR`RaR_((s./mandelgen.pyt bluevaluess6c Gsg}tdddD]}||ddfq~g}tdddD]}|d|dfqJ~g}tdddD]}|dd|fq~~}|S(s,progress from white to cyan to blue to blackiiii(R(RR`tredt_[2]tgreent_[3]tblueR_((s./mandelgen.pyt shadesofbluesg7c Gsng}tdddD]}|d|dfq~g}tdddD]}|dd|fqJ~}|S(s#progress from cyan to blue to blackiiii(R(RR`RkRjRmR_((s./mandelgen.pyt shadesofblue2s37cGs>g}tddD]}|d||fq~dg}|S(s<go from dark to light until the black of the Mandelbrot lakeiiii(iii(R(RR`RaR_((s./mandelgen.pyt inverse_reds:cGs>g}tddD]}|||dfq~dg}|S(s<go from dark to light until the black of the Mandelbrot lakeiii(iii(R(RR`RaR_((s./mandelgen.pyt inverse_blues:cGs^g}xQtdddD]=}x4t|ddD] }|||||fg7}q2WqW|S(sJFn's alternating color scheme e.g. yellow rg0 fading into magenta r0biii(R(RR_tstepRa((s./mandelgen.pytfnbrots"c Gsng}tdddD]}|d||fq~g}tdddD]}|dd|fqJ~}|S(scyan to black, blue to blackiii(R(RR`RaRjRmR_((s./mandelgen.pyt twotoneblue s37c Gsg}tdddD]a}||||f||df|d|fd||fdd|f|ddfd|dfgq~}g}x|D]}|i|qW|dg7}|dS(sgswitch between different color patterns from white to black it looks Peter Max-ish, hence the nameiiii(iii(Rtextend(RR`RaR_R#Rb((s./mandelgen.pyt maxelbrotsz csttdodtd0sis len(colormap)(R(RR'Rtfiltertappend(RRW((Rgs./mandelgen.pyRW&s      cCs1tdiddtdtdtdfdS(NR"s %s %d %d %d tP6RRR&(R'R1(((s./mandelgen.pyR)9sc Cstd|\}}}tddjo"tdid|||fn@tdid|d?|d|d?|d|d?|dfdS(NRRR&iR"s%c%c%cs %c%c%c%c%c%ci(R'R1(R=trtgRA((s./mandelgen.pyR9=s "cGs0tddtds<5                      #                 "" ,