def create_data():
global offset
DB.connect()
make_tables()
save_coords()
offset = 0
o = Tetrahedron()
poly_faces(o, "TETRA")
o = InvTetrahedron()
poly_faces(o, "INVTET")
o = Cube()
poly_faces(o, "CUBE")
o = Octahedron()
poly_faces(o, "OCTA")
o = RD()
poly_faces(o, "RD")
o = Cuboctahedron()
poly_faces(o, "VE")
#save_faces()
save_polys()
DB.disconnect()
def animation3(targdir="anim3"):
targdir = os.path.join(".",targdir)
if not os.path.isdir(targdir):
os.mkdir(targdir)
radii = np.linspace(0, 0.5, 8).tolist()
radii += reversed(radii[1:-1])
for frame_id, radius in enumerate(radii, start=1):
filename = f"balls{frame_id:03}.pov"
with open(os.path.join(targdir,filename), "w") as target:
target.write(pov_header)
cubocta = Cuboctahedron()
cubocta.vert_radius = radius
cubocta.vert_color = "rgb <1,0,0>"
if frame_id <= 9:
draw_poly(cubocta, target, v=True, e=False)
draw_vert(Qvector((0,0,0,0)), c=cubocta.vert_color,
r=cubocta.vert_radius, t=target)
if frame_id > 9:
cube = Cube()
ico = Icosahedron() * (1/PHI * math.sqrt(2)/2) * 1.2 * radius
st = Struts(cube, ico)
draw_poly(ico, target)
draw_poly(st, target)
# draw_poly(cube, target)
for ik in twelve_around_one(ico):
draw_poly(ik, target)
# for c in twelve_around_one(cube):
# draw_poly(c, target)
for s in twelve_around_one(st):
draw_poly(s, target)
def cubocta_in_cube_lattice():
with open("cubocta.pov", "w") as target:
target.write(pov_header)
cubocta = Cuboctahedron()
cube = Cube()
draw_poly(cubocta, target)
for c in twelve_around_one(cube):
draw_poly(c, target)
def ch():
with open("concentric_hierarchy.pov", "w") as target:
target.write(pov_header)
ico = Icosahedron() * 2
cubo = Cuboctahedron() * 2
cub = Cube() * 2
rd = RD() * 2
tet = Tetrahedron() * 2
octa = Octahedron() * 2
draw_poly(ico, target)
draw_poly(cubo, target)
draw_poly(rd, target)
draw_poly(cub, target)
draw_poly(tet, target)
draw_poly(octa, target)
def animation4(targdir="anim4"):
targdir = os.path.join(".",targdir)
if not os.path.isdir(targdir):
os.mkdir(targdir)
radii = np.linspace(0.6, 1.3, 8).tolist()
radii += reversed(radii[1:-1])
for frame_id, radius in enumerate(radii, start=1):
filename = f"balls{frame_id:03}.pov"
with open(os.path.join(targdir,filename), "w") as target:
target.write(pov_header)
ico = Icosahedron() * radius
cube= Cube() * PHI * 1.5
s = Struts(cube, ico)
draw_poly(ico, target)
draw_poly(cube, target)
draw_poly(s, target)
def spheres_bulging():
with open("sphere_pack2.pov", "w") as target:
target.write(pov_header)
cube = Cube()
rd = RD()
octa = Octahedron()
cubocta = Cuboctahedron()
draw_poly(cube, target)
draw_poly(octa, target)
draw_poly(rd, target)
for p in twelve_around_one(cube):
draw_poly(p, target)
for p in twelve_around_one(rd):
draw_poly(p, target)
for p in twelve_around_one(octa):
draw_poly(p, target)
cubocta.vert_radius = Qvector((2,1,1,0)).length()/2 # grow the balls
cubocta.vert_color = "rgb <1,0,0>"
draw_poly(cubocta, target, v=True, e=False)
draw_vert(Qvector((0,0,0,0)), c=cubocta.vert_color,
r=cubocta.vert_radius, t=target)