from pythonscad import *
from pylibfive import *

levels = 4
diameter = 7
c = lv_coord()

ang = (lv.atan2(c[0], c[1])) % (3.1415 / 20.0)
mag = lv.sqrt(c[0] * c[0] + c[1] * c[1])
cir_c = [mag * lv.cos(ang), mag * lv.sin(ang), (lv_clamp(c[2] / 1.5, 0, 3)) % 1.0]
cir_c1 = [cir_c[0], cir_c[1], lv.min(cir_c[2], 1 - cir_c[2])]
cir_c2 = [mag * lv.cos(ang), mag * lv.sin(ang), (lv_clamp(c[2] / 1.5, 3, 4)) % 1.0]
out = lv_cylinder(
    c,
    levels * 1.5,
    diameter + 0.2 * lv.pow(0.5 * lv.cos(2 * 3.1415 * c[2] / 1.5) + 0.5, 10.0) - 0.2,
)
inn = lv_cylinder(lv_trans(c, [0, 0, -0.2]), levels * 1.5 + 0.4, diameter - 0.5)

# Saeulen
cyl1 = lv_cylinder(lv_trans(cir_c1, [diameter + 0.10, 0.6, 0.1]), 0.98, 0.05)
cyl2 = lv_box(lv_trans(cir_c1, [diameter + 0.10, 0.6, 0.1]), [0.1, 0.1, 0.05])
cyl = lv_union(cyl1, cyl2)

# Doors
door1 = lv_box(lv_trans(cir_c, [diameter - 1, +0.15, 0.4]), [2, 0.3, 0.3])
door2 = lv_cylinder([cir_c[2] - 0.5, cir_c[1] - 0.15, cir_c[0]], diameter + 2, 0.3)
door = lv_union(door1, door2)

# Windows
window = lv_box(lv_trans(cir_c2, [diameter - 1, +0.15, 0.5]), [2, 0.3, 0.2])

# Compile
notbuilt = lv_box(lv_trans(c, [0, -3.9, 1.5 * 5]), [9, 4 + 0.2 * c[2], 4.5])
negatives = lv_union(inn, lv_union(lv_union(door, window), notbuilt))
collosseum = lv_union(lv_difference(out, negatives) - 0.1, lv_difference(cyl, notbuilt))
fobj = frep(collosseum, [-15, -15, -1], [15, 15, 8], 40)

output(fobj)