Welcome, guest | Sign In | My Account | Store | Cart
```# True-color plasma fractal using non-recursive algorithm
# FB - 20160122
import random
from PIL import Image
imgx = 512; imgy = 512
image = Image.new("RGB", (imgx, imgy))
mx = imgx - 1; my = imgy - 1
f = 2.0 # roughness

def putpixel(x, y, r, g, b):
global pixels
pixels[int(round(x)), int(round(y))] = \
(int(round(r)), int(round(g)), int(round(b)))

def getpixel(x, y):
return pixels[int(round(x)), int(round(y))]

putpixel(0, 0, random.randint(0, 255), \
random.randint(0, 255), random.randint(0, 255))
putpixel(mx, 0, random.randint(0, 255), \
random.randint(0, 255), random.randint(0, 255))
putpixel(mx, my, random.randint(0, 255), \
random.randint(0, 255), random.randint(0, 255))
putpixel(0, my, random.randint(0, 255), \
random.randint(0, 255), random.randint(0, 255))

j = -1
while True:
j += 1; j2 = 2 ** j
jx = float(mx) / j2; jy = float(my) / j2
if jx < 1 and jy < 1: break
for i in range(j2):
y0 = i * jy; y1 = y0 + jy; y = y0 + jy / 2.0
for k in range(j2):
x0 = k * jx; x1 = x0 + jx; x = x0 + jx / 2.0

a = getpixel(x0, y0); b = getpixel(x1, y0)
c = getpixel(x0, y1); d = getpixel(x1, y1)

putpixel(x, y, (a[0] + b[0] + c[0] + d[0]) / 4.0, \
(a[1] + b[1] + c[1] + d[1]) / 4.0, (a[2] + b[2] + c[2] + d[2]) / 4.0)
putpixel(x, y0, (a[0] + b[0]) / 2.0 + jx * (random.random() - .5) * f, \
(a[1] + b[1]) / 2.0 + jx * (random.random() - .5) * f, \
(a[2] + b[2]) / 2.0 + jx * (random.random() - .5) * f)
putpixel(x0, y, (a[0] + c[0]) / 2.0 + jy * (random.random() - .5) * f, \
(a[1] + c[1]) / 2.0 + jy * (random.random() - .5) * f, \
(a[2] + c[2]) / 2.0 + jy * (random.random() - .5) * f)
putpixel(x1, y, (b[0] + d[0]) / 2.0 + jy * (random.random() - .5) * f, \
(b[1] + d[1]) / 2.0 + jy * (random.random() - .5) * f, \
(b[2] + d[2]) / 2.0 + jy * (random.random() - .5) * f)
putpixel(x, y1, (c[0] + d[0]) / 2.0 + jx * (random.random() - .5) * f, \
(c[1] + d[1]) / 2.0 + jx * (random.random() - .5) * f, \
(c[2] + d[2]) / 2.0 + jx * (random.random() - .5) * f)

image.save("plasma.png", "PNG")
```

#### Diff to Previous Revision

```--- revision 3 2010-03-29 02:28:35
+++ revision 4 2016-01-22 22:44:33
@@ -1,61 +1,56 @@
# True-color plasma fractal using non-recursive algorithm
-# FB - 201003287
+# FB - 20160122
import random
from PIL import Image
-imgx = 512
-imgy = 512
+imgx = 512; imgy = 512
image = Image.new("RGB", (imgx, imgy))
-mx = imgx - 1
-my = imgy - 1
+mx = imgx - 1; my = imgy - 1
f = 2.0 # roughness

-image.putpixel((0, 0), (random.randint(0, 255),
-                        random.randint(0, 255),
-                        random.randint(0, 255)))
-image.putpixel((mx, 0), (random.randint(0, 255),
-                         random.randint(0, 255),
-                         random.randint(0, 255)))
-image.putpixel((mx, my), (random.randint(0, 255),
-                          random.randint(0, 255),
-                          random.randint(0, 255)))
-image.putpixel((0, my), (random.randint(0, 255),
-                         random.randint(0, 255),
-                         random.randint(0, 255)))
+def putpixel(x, y, r, g, b):
+    global pixels
+    pixels[int(round(x)), int(round(y))] = \
+        (int(round(r)), int(round(g)), int(round(b)))
+
+def getpixel(x, y):
+    return pixels[int(round(x)), int(round(y))]
+
+putpixel(0, 0, random.randint(0, 255), \
+    random.randint(0, 255), random.randint(0, 255))
+putpixel(mx, 0, random.randint(0, 255), \
+    random.randint(0, 255), random.randint(0, 255))
+putpixel(mx, my, random.randint(0, 255), \
+    random.randint(0, 255), random.randint(0, 255))
+putpixel(0, my, random.randint(0, 255), \
+    random.randint(0, 255), random.randint(0, 255))
+
j = -1
while True:
-    j += 1
-    j2 = 2 ** j
-    jx = float(mx) / j2
-    jy = float(my) / j2
+    j += 1; j2 = 2 ** j
+    jx = float(mx) / j2; jy = float(my) / j2
if jx < 1 and jy < 1: break
for i in range(j2):
-        y0 = i * jy
-        y1 = y0 + jy
-        y = y0 + jy / 2.0
+        y0 = i * jy; y1 = y0 + jy; y = y0 + jy / 2.0
for k in range(j2):
-            x0 = k * jx
-            x1 = x0 + jx
-            x = x0 + jx / 2.0
+            x0 = k * jx; x1 = x0 + jx; x = x0 + jx / 2.0

-            a = image.getpixel((x0, y0))
-            b = image.getpixel((x1, y0))
-            c = image.getpixel((x0, y1))
-            d = image.getpixel((x1, y1))
+            a = getpixel(x0, y0); b = getpixel(x1, y0)
+            c = getpixel(x0, y1); d = getpixel(x1, y1)

-            image.putpixel((x, y), ((a[0] + b[0] + c[0] + d[0]) / 4.0,
-                                    (a[1] + b[1] + c[1] + d[1]) / 4.0,
-                                    (a[2] + b[2] + c[2] + d[2]) / 4.0))
-            image.putpixel((x, y0), ((a[0] + b[0]) / 2.0 + jx * (random.random() - .5) * f,
-                                     (a[1] + b[1]) / 2.0 + jx * (random.random() - .5) * f,
-                                     (a[2] + b[2]) / 2.0 + jx * (random.random() - .5) * f))
-            image.putpixel((x0, y), ((a[0] + c[0]) / 2.0 + jy * (random.random() - .5) * f,
-                                     (a[1] + c[1]) / 2.0 + jy * (random.random() - .5) * f,
-                                     (a[2] + c[2]) / 2.0 + jy * (random.random() - .5) * f))
-            image.putpixel((x1, y), ((b[0] + d[0]) / 2.0 + jy * (random.random() - .5) * f,
-                                     (b[1] + d[1]) / 2.0 + jy * (random.random() - .5) * f,
-                                     (b[2] + d[2]) / 2.0 + jy * (random.random() - .5) * f))
-            image.putpixel((x, y1), ((c[0] + d[0]) / 2.0 + jx * (random.random() - .5) * f,
-                                     (c[1] + d[1]) / 2.0 + jx * (random.random() - .5) * f,
-                                     (c[2] + d[2]) / 2.0 + jx * (random.random() - .5) * f))
+            putpixel(x, y, (a[0] + b[0] + c[0] + d[0]) / 4.0, \
+                (a[1] + b[1] + c[1] + d[1]) / 4.0, (a[2] + b[2] + c[2] + d[2]) / 4.0)
+            putpixel(x, y0, (a[0] + b[0]) / 2.0 + jx * (random.random() - .5) * f, \
+                (a[1] + b[1]) / 2.0 + jx * (random.random() - .5) * f, \
+                (a[2] + b[2]) / 2.0 + jx * (random.random() - .5) * f)
+            putpixel(x0, y, (a[0] + c[0]) / 2.0 + jy * (random.random() - .5) * f, \
+                (a[1] + c[1]) / 2.0 + jy * (random.random() - .5) * f, \
+                (a[2] + c[2]) / 2.0 + jy * (random.random() - .5) * f)
+            putpixel(x1, y, (b[0] + d[0]) / 2.0 + jy * (random.random() - .5) * f, \
+                (b[1] + d[1]) / 2.0 + jy * (random.random() - .5) * f, \
+                (b[2] + d[2]) / 2.0 + jy * (random.random() - .5) * f)
+            putpixel(x, y1, (c[0] + d[0]) / 2.0 + jx * (random.random() - .5) * f, \
+                (c[1] + d[1]) / 2.0 + jx * (random.random() - .5) * f, \
+                (c[2] + d[2]) / 2.0 + jx * (random.random() - .5) * f)

-image.save("plasma3.png", "PNG")
+image.save("plasma.png", "PNG")
```