Cyclic Quine

2017-08-12, post № 177

Programming, Python, Quine, #self-referential

A classic quine is a program which outputs its own source code.
At first, such a program’s existence seems weird if not impossible, as it has to be so self-referential that it knows about itself everything, including how to know about itself. However, writing quines is possible, if not [1] trivial.

A cyclic quine therefore is a program which outputs source code which differs from its own source code, yet outputs the original source code when run (the cycle length could be greater than one). So when running source codes 𝛹 and 𝛷, they output source codes 𝛷 and 𝛹.

Therefore, when one saves the first program as q0.py and the second as q1.py, one can create both source codes from one another (the following bash commands [2] will not change the files’ contents).

$ python q0.py > q1.py
$ python q1.py > q0.py
$ python q0.py | python > q0.py
$ python q0.py | python | python | python | python > q1.py
Source code: cyclic-quine_q0.py
Source code: cyclic-quine_q1.py

Mandelbrot Set III

2017-07-29, post № 176

Java, mathematics, programming, #dragging, #fractal, #fractals, #Mandelweed, #mouse control, #multithreading, #scrolling, #thread, #viewer

I wrote my first ever Mandelbrot Set renderer back in 2015 and used Python to slowly create fractal images. Over a year later, I revisited the project with a Java version which — due to its code being actually compiled — ran much faster, yet had the same clunky interface; a rectangle the user had to draw and a key they had to press to the view change to the selected region.
In this post, over half a year later, I present my newest Mandelbrot Set fractal renderer (download the .jar), written in Java, which both runs fast and allows a much more intuitive and immersive walk through the complex plane by utilizing mouse dragging and scrolling.
The still time demanding task of rendering fractals — even in compiled languages — is split up into a low quality preview rendering, a normal quality display rendering and a high quality 4K (UHD-1 at 𝟥𝟪𝟦𝟢 ⨉ 𝟤𝟣𝟨𝟢 pixels to keep a 𝟣𝟨 : 𝟫 image ratio) rendering, all running in seperate threads.

Rainbow spiral

The color schemes where also updated, apart from the usual black-and-white look there are multiple rainbow color schemes which rely on the HSB color space, zebra color schemes which use the iterations taken modulo some constant to define the color and a prime color scheme which tests if the number of iterations taken is prime.

Zebra spiral

Apart from the mouse and keyboard control, there is also a menu bar (implemented using Java’s JMenuBar) which allows for more conventional user input through a proper GUI.


  • Left mouse dragging: pan view,
  • Left mouse double click: set cursor’s complex number to image center,
  • Mouse scrolling: zoom view,
  • Mouse scrolling +CTRL: pan view,
  • ‘p’: render high definition fractal,
  • ‘r’: reset view to default,
  • ‘w’, ‘s’: zoom frame,
  • Arrow keys: pan view,
  • Arrow keys +CTRL: zoom view,
  • Menu bar

    • “Fractal”: extra info about current fractal rendering,
    • “Color Scheme”: change color scheme and maximum iteration depth,
    • “HD”: controls for high definition rendering,
    • “Extra”: help and about.
Blue spiral

A bit more on how the three threads are implemented.
Whenever the user changes the current view, the main program thread renders a low quality preview and immediately draws it to the screen. In the background, the normal quality thread (its pixel dimensions match the frame’s pixel dimensions) is told to start working. Once this medium quality rendering is finished, it is preferred to the low quality rendering and gets drawn on the screen.
If the user likes a particular frame, they can initiate a high quality rendering (4K UHD-1, 𝟥𝟪𝟦𝟢 ⨉ 𝟤𝟣𝟨𝟢 pixels( either by pressing ‘q’ or selecting “HD ❯ Render current frame”. This high quality rendering obviously takes some time and a lot of processing power, so this thread is throttled by default to allow the user to further explore the fractal. Throttling can be disabled through the menu option “HD ❯ Fast rendering”. There is also the option to tell the program to exit upon having finished the last queued high definition rendering (“HD ❯ Quit when done”).
The high definition renderings are saved as .png files and named with their four defining constants. Zim and Zre define the image’s complex center, Zom defines the complex length above the image’s center. Clr defines the number of maximum iterations.

Another blue spiral

Just to illustrate how resource intensive fractal rendering really is.
A 4K fractal at 𝟥𝟪𝟦𝟢 ⨉ 𝟤𝟣𝟨𝟢 pixels with an iteration depth of 𝟤𝟧𝟨 would in the worst case scenario (no complex numbers actually escape) require 3840\cdot 2160\cdot 256\cdot 4=8493465600 double multiplications. If you had a super-optimized CPU which could do one double multiplication a clock tick (which current CPUs definitely cannot) and ran at 𝟦.𝟢𝟢 GHz, it would still take that massively overpowered machine \frac{8493465600}{4\cdot 10^9}=2.123 seconds. [1] Larger images and higher maximum iterations would only increase the generated overhead.
The program’s source code is listed below and can also be downloaded (.java), though the compiled .jar can also be downloaded.

Green self-similarity

Unrelated to algorithmically generating fractal renderings, I recently found a weed which seemed to be related to the Mandelbrot Set and makes nature’s intertwined relationship with fractals blatently obvious. I call it the Mandel Weed.

Mandel Weed
Source code: MandelbrotSetExplorer.java
Extra assets: mandelbrot-set-iii_mk.sh, mandelbrot-set-iii_3840x2160_Zre-0.7451119128483528_Zim0.0980857735166215_Zom0.0015327239951060212_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.3091727307741978_Zim0.0798169059094002_Zom4.0442023600132944E-8_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.3203805004961733_Zim0.08763686694940123_Zom6.070130782934575E-9_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.3716081094521764_Zim0.0758961984414791_Zom8.091648167718248E-4_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.99181475237973_Zim-3.982534247685533E-7_Zom3.0401148534882E-6_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.9963766117665036_Zim-8.944472584209038E-7_Zom7.062361041362908E-6_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.996376645843924_Zim3.7516943208736117E-6_Zom2.859109719673806E-10_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.1907127228021151_Zim-0.5505294608801341_Zom5.308930362839995E-4_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.2428560033896521_Zim0.5073262577216794_Zom5.8089779367154003E-5_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.3411570576688668_Zim0.5750384378500207_Zom1.100339471773629E-4_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.35167270726867744_Zim-0.06184250438694921_Zom9.837555143550957E-5_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.35940345191970086_Zim-0.06373194819231794_Zom8.586182289196699E-7_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.4443477244180509_Zim-0.38358230098736174_Zom1.4540775100673537E-6_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.44503254319557867_Zim-0.38339448301161105_Zom2.4624930313254245E-6_Clr256.png


2017-07-15, post № 175

art, ascii, programming, Python PIL, #ANSI, #ASCII, #characters, #image processing, #low resolution, #retro look

Most images nowadays are represented using pixels. They are square, often relatively small and numerous, come in (2^8)^3 different colors and thereby do a good job being the fundamental building block of images. But one can imagine more coarse-grained and differently shaped pixels.
An interesting fact is, that in most monotype fonts two characters placed right next to each other (for example ‘$$’) occupy roughly a square area. So simple ASCII characters can indeed be used to approximately describe any ordinary image.
Asciify does exactly this; it takes in an image and some optional parameters and maps the pixels’ intensity onto a character set. Both the large and small default character sets are taken from a post by Paul Bourke.

In conjunction with asciify.py, I wrote index.py, which asciifies a bunch of images and results in their html form; it also creates an index. All images asciified for this post can be viewed through this index [1].

Converting an image to its asciified form works best when there is a lot of contrast in the image. Because of this, some pre-processing of the image may be required for best results (all images shown where only cropped or rotated). The built-in color functionality also only knows of 𝟪 colors, so bright and different colors look the best, as they interestingly differentiate from one another. The asciified image’s size also plays a role, the larger it is, the better the characters blend into one and appear to be one image.

Asciify is operated on a command prompt; python asciify.py img.png. To parse arguments, the built-in Python module argparse is used. The images are opened and read using the Python Imaging Library module PIL, which needs to be installed for this program to work.
Optional arguments include --size N, where the maximum size can be specified, --invert and --smallcharset, which can sometimes increase the asciified image’s visual appeal and --html, which will output an html file to be viewed in a browser. To see the program’s full potential, simply run python asciify.py --help [2].
Source code for both asciify.py and index.py can be downloaded, the first is also listed below.

The two examples above use the color mode, though certain images also work in default black and white mode, such as this spider I photographed.

Then again, the colored text also has its charm, especially when the source image has bright colors and a lot of contrast.

Source code: asciify.py


2017-07-01, post № 174

mathematics, #17, #2017, #integer, #integer sequences, #July, #sequences

Today it is the first day of July in the year 2017. On this day there is a point in time which can be represented as 1.7.2017, 17:17:17.
To celebrate this symbolically speaking 17-heavy day, I created a list of 17 integer sequences which all contain the number 17.
All sequences were generated using a Python program; the source code can be viewed below or downloaded. Because the following list is formatted using LaTex, the program’s plaintext output can also be downloaded.

  1. Prime numbers 𝑛.

    \{2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, \dots\}
  2. Odd positive integers 𝑛 whose number of goldbach sums (all possible sums of two primes) of 𝑛 + 𝟣 and 𝑛 - 𝟣 are equal to one another.

    \{5, 7, 15, 17, 19, 23, 25, 35, 75, 117, 177, 207, 225, 237, 321, 393, 453, 495, 555, 567, \dots\}
  3. Positive integers n who are part of a Pythagorean triple excluding 𝟢: n^2=a^2+b^2 with integers a,b>0.

    \{5, 10, 13, 15, 17, 20, 25, 26, 29, 30, 34, 35, 37, 39, 40, 41, 45, 50, 51, 52, \dots\}
  4. Positive integers 𝑛 where \lfloor(n!)^{\frac{1}{n}}\rfloor is prime

    \{4, 5, 6, 7, 8, 12, 13, 17, 18, 19, 28, 29, 33, 34, 35, 44, 45, 46, 49, 50, \dots\}
  5. Positive integers 𝑛 with distance 𝟣 to a perfect square.

    \{1, 2, 3, 5, 8, 10, 15, 17, 24, 26, 35, 37, 48, 50, 63, 65, 80, 82, 99, 101, \dots\}
  6. Positive integers 𝑛 where the number of perfect squares including 𝟢 less than 𝑛 is prime.

    \{2, 3, 4, 5, 6, 7, 8, 9, 17, 18, 19, 20, 21, 22, 23, 24, 25, 37, 38, 39, \dots\}
  7. Prime numbers 𝑛 where either 𝑛 - 𝟤 or 𝑛 + 𝟤 (exclusive) are prime.

    \{3, 7, 11, 13, 17, 19, 29, 31, 41, 43, 59, 61, 71, 73, 101, 103, 107, 109, 137, 139, \dots\}
  8. Positive integers 𝑛 whose three-dimensional vector’s (n,n,n) floored length is prime, \lfloor\sqrt{3\cdot n^2}\rfloor is prime.

    \{2, 3, 8, 10, 11, 17, 18, 24, 25, 31, 39, 41, 46, 48, 60, 62, 63, 76, 91, 100, \dots\}
  9. Positive integers 𝑛 who are the sum of a perfect square and a perfect cube (excluding 𝟢).

    \{2, 5, 9, 10, 12, 17, 24, 26, 28, 31, 33, 36, 37, 43, 44, 50, 52, 57, 63, 65, \dots\}
  10. Positive integers 𝑛 whose decimal digit sum is the cube of a prime.

    \{8, 17, 26, 35, 44, 53, 62, 71, 80, 107, 116, 125, 134, 143, 152, 161, 170, 206, 215, 224, \dots\}
  11. Positive integers 𝑛 for which \text{decimal\_digitsum}(n)+n is a perfect square.

    \{2, 8, 17, 27, 38, 72, 86, 135, 161, 179, 216, 245, 275, 315, 347, 432, 467, 521, 558, 614, \dots\}
  12. Prime numbers 𝑛 for which \text{decimal\_digitsum}(n^4) is prime.

    \{2, 5, 7, 17, 23, 41, 47, 53, 67, 73, 97, 103, 113, 151, 157, 163, 173, 179, 197, 199, \dots\}
  13. Positive integers 𝑛 where \text{decimal\_digitsum}(2 \cdot n) is a substring of 𝑛.

    \{9, 17, 25, 52, 58, 66, 71, 85, 90, 104, 107, 115, 118, 123, 137, 142, 151, 156, 165, 170, \dots\}
  14. Positive integers 𝑛 whose decimal reverse is prime.

    \{2, 3, 5, 7, 11, 13, 14, 16, 17, 20, 30, 31, 32, 34, 35, 37, 38, 50, 70, 71, \dots\}
  15. Positive integers 𝑛 who are a decimal substring of n^n.