source: util/src/AnimateStereo.py.in@ a3ffb44

#!@PYTHON@
#
# Creates a series of stereo header (left and right) files to be included in a raster3d with a rotational coordinate transformation

import sys, random, math, re
wrerr=sys.stderr.write
wrout=sys.stdout.write

# check arguments
if len(sys.argv) < 2:
        print "Usage: "+sys.argv[0]+" <prefix> [AtomCount]"
        sys.exit(1)

# get atomcount
if len(sys.argv) < 2:
        Atoms=10
else:
        Atoms=int(sys.argv[2])

# deduce r-factor
factor=15.0 + (40.0-15.0)/2000*Atoms

for i in range(36):
        for type in range(2):
                c = math.cos(2.*math.pi/36.*i);
                s = math.sin(2.*math.pi/36.*i);
                if type == 0: # left
                        output = open("%s-%03d-l.r3d" % (sys.argv[1],i), "w")
                else:
                        output = open("%s-%03d-r.r3d" % (sys.argv[1],i), "w")
                output.write("Example of material properties and file indirection\n");
                output.write("80 64     tiles in x,y\n");
                output.write(" 8  8     pixels (x,y) per tile\n");
                output.write("4         3x3 virtual pixels -> 2x2 pixels\n");
                output.write("0 0 0   background colour\n");
                output.write("T         cast shadows\n");
                output.write("25        Phong power\n");
                output.write("0.15      secondary light contribution\n");
                output.write("0.05      ambient light contribution\n");
                output.write("0.25      specular reflection component\n");
                output.write("0.0       eye position\n");
                output.write("1 1 1     main light source position\n");
                output.write("%lg %lg 0 0   input coordinate, radius transformation\n" % (c, s));
                output.write("%lg %lg 0 0\n" % (-s, c));
                if type == 0:   #left
                        output.write("0.03 0 1 0\n");
                else: # right
                        output.write("-0.03 0 1 0\n");
                output.write("0 0 0 %5.2f\n" % (factor));
                output.write("3         mixed objects\n");
                output.write("*\n");
                output.write("*\n");
                output.write("*\n");
                output.close()

for i in range(36):
        for type in range(2):
                c = math.cos(2.*math.pi/36.*i);
                s = math.sin(2.*math.pi/36.*i);
                if type == 0: # left
                        output = open("%s-%03d-l.r3d" % (sys.argv[1],i+36), "w")
                else:
                        output = open("%s-%03d-r.r3d" % (sys.argv[1],i+36), "w")
                output.write("Example of material properties and file indirection\n");
                output.write("80 64     tiles in x,y\n");
                output.write(" 8  8     pixels (x,y) per tile\n");
                output.write("4         3x3 virtual pixels -> 2x2 pixels\n");
                output.write("0 0 0   background colour\n");
                output.write("T         cast shadows\n");
                output.write("25        Phong power\n");
                output.write("0.15      secondary light contribution\n");
                output.write("0.05      ambient light contribution\n");
                output.write("0.25      specular reflection component\n");
                output.write("0.0       eye position\n");
                output.write("1 1 1     main light source position\n");
                output.write("%lg 0 %lg 0   input coordinate, radius transformation\n" % (c, s));
                output.write("0 1 0 0\n");
                if type == 0:   #left
                        output.write("%lg 0 %lg 0\n" % (-s+0.03, c));
                else: # right
                        output.write("%lg 0 %lg 0\n" % (-s-0.03, c));
                output.write("0 0 0 %5.2f\n" % (factor));
                output.write("3         mixed objects\n");
                output.write("*\n");
                output.write("*\n");
                output.write("*\n");
                output.close()

for i in range(36):
        for type in range(2):
                c = math.cos(2.*math.pi/36.*i);
                s = math.sin(2.*math.pi/36.*i);
                if type == 0: # left
                        output = open("%s-%03d-l.r3d" % (sys.argv[1],i+72), "w")
                else:
                        output = open("%s-%03d-r.r3d" % (sys.argv[1],i+72), "w")
                output.write("Example of material properties and file indirection\n");
                output.write("80 64     tiles in x,y\n");
                output.write(" 8  8     pixels (x,y) per tile\n");
                output.write("4         3x3 virtual pixels -> 2x2 pixels\n");
                output.write("0 0 0   background colour\n");
                output.write("T         cast shadows\n");
                output.write("25        Phong power\n");
                output.write("0.15      secondary light contribution\n");
                output.write("0.05      ambient light contribution\n");
                output.write("0.25      specular reflection component\n");
                output.write("0.0       eye position\n");
                output.write("1 1 1     main light source position\n");
                output.write("1 0 0 0   input coordinate, radius transformation\n");
                output.write("0 %lg %lg 0\n" % (c, s));
                if type == 0:   #left
                        output.write("0.03 %lg %lg 0\n" % (-s, c));
                else: # right
                        output.write("-0.03 %lg %lg 0\n" % (-s, c));
                output.write("0 0 0 %5.2f\n" % (factor));
                output.write("3         mixed objects\n");
                output.write("*\n");
                output.write("*\n");
                output.write("*\n");
                output.close()