mirror of
https://github.com/stevenrobertson/cuburn.git
synced 2025-02-05 11:40:04 -05:00
92 lines
2.7 KiB
Python
92 lines
2.7 KiB
Python
#!/usr/bin/python
|
|
#
|
|
# flam3cuda, one of a surprisingly large number of ports of the fractal flame
|
|
# algorithm to NVIDIA GPUs.
|
|
#
|
|
# This one is copyright 2010 Steven Robertson <steven@strobe.cc>
|
|
#
|
|
# This program is free software; you can redistribute it and/or modify
|
|
# it under the terms of the GNU General Public License version 2 or later
|
|
# as published by the Free Software Foundation.
|
|
|
|
import os
|
|
import sys
|
|
from pprint import pprint
|
|
from ctypes import *
|
|
|
|
import numpy as np
|
|
|
|
from cuburn.device_code import *
|
|
from cuburn.cuda import LaunchContext
|
|
from fr0stlib.pyflam3 import *
|
|
from fr0stlib.pyflam3._flam3 import *
|
|
from cuburn.render import *
|
|
|
|
import pyglet
|
|
|
|
def dump_3d(nda):
|
|
with open('/tmp/data.txt', 'w') as f:
|
|
for row in nda:
|
|
f.write(' | '.join([' '.join(
|
|
['%4.1g\t' % x for x in pt]) for pt in row]) + '\n')
|
|
|
|
def main(args):
|
|
verbose = 1
|
|
if '-d' in args:
|
|
verbose = 3
|
|
|
|
with open(args[-1]) as fp:
|
|
genomes = Genome.from_string(fp.read())
|
|
anim = Animation(genomes)
|
|
anim.compile()
|
|
bins = anim.render_frame()
|
|
|
|
# Embarrassingly poor (and s-l-o-w) approximation of log filtering, because
|
|
# I want to keep working on performance right now
|
|
alpha = np.array([[bins[y][x][3] for x in range(anim.features.hist_width)]
|
|
for y in range(anim.features.hist_height)])
|
|
k2ish = (256./(np.mean(alpha)+1e-9))
|
|
lses = 20 * np.log2(1.0 + alpha * k2ish) / (alpha+1e-6)
|
|
for x in range(anim.features.hist_width):
|
|
for y in range(anim.features.hist_height):
|
|
bins[y][x] *= lses[y][x]
|
|
bins = np.minimum(bins, 255)
|
|
bins = bins.astype(np.uint8)
|
|
|
|
if '-g' not in args:
|
|
return
|
|
|
|
print anim.features.hist_width
|
|
print anim.features.hist_height
|
|
print anim.features.hist_stride
|
|
window = pyglet.window.Window(1600, 900)
|
|
image = pyglet.image.ImageData(anim.features.hist_width,
|
|
anim.features.hist_height,
|
|
'RGBA',
|
|
bins.tostring(),
|
|
-anim.features.hist_stride*4)
|
|
tex = image.texture
|
|
|
|
pal = (anim.ctx.ptx.instances[PaletteLookup].pal * 255.).astype(np.uint8)
|
|
image2 = pyglet.image.ImageData(256, 16, 'RGBA', pal.tostring())
|
|
|
|
@window.event
|
|
def on_draw():
|
|
window.clear()
|
|
tex.blit(0, 0)
|
|
image2.blit(0, 0)
|
|
|
|
@window.event
|
|
def on_key_press(sym, mod):
|
|
if sym == pyglet.window.key.Q:
|
|
pyglet.app.exit()
|
|
|
|
pyglet.app.run()
|
|
|
|
if __name__ == "__main__":
|
|
if len(sys.argv) < 2 or not os.path.isfile(sys.argv[-1]):
|
|
print "Last argument must be a path to a genome file"
|
|
sys.exit(1)
|
|
main(sys.argv)
|
|
|