import rhinoscriptsyntax as rs
import random as ran
import math

# Delete everything and start from scratch
allobjs = rs.AllObjects()
rs.DeleteObjects(allobjs)

rs.EnableRedraw(False)

#-------------------------------------------------------------------------------
# dom-ino variables:
#-------------------------------------------------------------------------------
A       = 5                         # A = Module size (distance between columns)
B       = A/3                       # B = Distance of columns to end of plate
thick   = 0.2                       # thickness of all slabs
hgt     = 2.7                       # height of room
xcol    = 2                         # columns in x direction
ycol    = 3                         # columns in y direction
levels  = 5                         # number of floor plates
f_height= 0.5                       # f_height = foundation height
f_size  = 0.8                       # f_size = foundation edge size
w_pts   = 360                       # number of Facade Elements
w_height= hgt*(levels-f_size)       # w_height = hight of the facade
w_size  = 0.02                      # w_size = facade edge size
frequency = 0.5                     # distance of the wave
#insertion = [A*(xcol-1)/2,A*(ycol-1)/2, f_height]

#-------------------------------------------------------------------------------
# derived values:
#-------------------------------------------------------------------------------
center_pt   = [A*(xcol-1)/2,A*(ycol-1)/2, f_height] # insertion point of floor plate
p_width     = A*(xcol-1) + 2*B                      # width of floor plate (x)
p_length    = A*(ycol-1) + f_size                   # length of floor plate (y)

# function to make box
def make_box(insertion=[0,0,0],xsize=10,ysize=10,zsize=10):
    # Define the corners of the box
    corners = [[0, 0, 0], [xsize, 0, 0], [xsize, ysize, 0], [0, ysize, 0],
            [0, 0, zsize], [xsize, 0, zsize], [xsize, ysize, zsize], [0, ysize, zsize]]
    box=rs.AddBox(corners)
    rs.MoveObject(box, (-xsize/2,-ysize/2,0))
    rs.MoveObject(box, insertion)
    return(box)

#function to create podest
def make_podest(insertion=[0,0,0],xsize=10,ysize=10,zsize=10):
    corners = [[0,0,0], [xsize,0,0], [xsize,ysize,0], [0,ysize,0],
               [0,0,zsize],[xsize,0,zsize],[xsize,ysize,zsize],[0,ysize,zsize]]
    box=rs.AddBox(corners)
    rs.MoveObject(box, insertion)
    return(box)

# function to create a field of foundations
def make_foundations(A=5.0, f_size=0.8, f_height=0.5, xcol=2, ycol=3):
    fns = []
    for i in range(xcol):
        for j in range(ycol):
            fns.append(make_box([i*A,j*A,0], f_size, f_size, f_height))
    return(fns)

# function to create a field of columns
def make_columns(A=5.0, level=0.7, thick=0.2, hgt=3.0, xcol=2, ycol=3):
    cls =[]
    for i in range(xcol):
        for j in range(ycol):
            cls.append(make_box([i*A,j*A,level], thick, thick, hgt))
    return(cls)

# Function to create a facade with height adapting to the building height
def make_facade(insertion=[0, 0, 0], p_width=p_width, p_length=p_length, w_size=w_size, building_height=10, frequency=frequency):
    fcd = []
    plane = rs.WorldXYPlane()
    plane = rs.RotatePlane(plane, 0, [0, 0, 0])
    rec = rs.AddRectangle(plane, p_width, p_length)
    move_vector = rs.VectorCreate(
        [center_pt[0] - p_width / 2, center_pt[1] - p_length / 2, center_pt[2] + thick],
        [0, 0, 0]
    )
    rs.MoveObject(rec, move_vector)
    fcd.append(rec)
    pts = rs.DivideCurve(rec, w_pts, False, True)

    # Create the "facade" layer if it doesnt exist
    if not rs.IsLayer("facade"):
        rs.AddLayer("facade")
        rs.LayerColor("facade", (139, 69, 19))  # Brown color for facade

    for i, p in enumerate(pts):
        # Adjust height of each facade element based on sine function and building height
        random_factor = ran.uniform(0.6, 1)
        element_height = max(thick, building_height * math.sin(math.radians(frequency * i)) * random_factor)
        random_size = w_size * ran.uniform(0.8, 1.5)

        # Create facade element with some random rotation for variety
        w_box = make_box(insertion, xsize=random_size, ysize=random_size, zsize=-element_height)
        rs.MoveObject(w_box, p)

        angle = ran.uniform(-10, 10)  # Random rotation angle
        rs.RotateObject(w_box, p, angle)

        # Set alternating colors for the facade elements
        if i % 2 == 0:
            rs.ObjectColor(w_box, (139, 69, 19))
        else:
            rs.ObjectColor(w_box, (205, 133, 63))

        # Assign facade elements to "facade" layer
        rs.ObjectLayer(w_box, "facade")

        fcd.append(w_box)

    rs.DeleteObjects(rec)
    return fcd, pts

# Function to create a Arc
def make_arc(insertion, rad, thick, hgt, orientation):
    segs=[]
    segs.append(rs.AddArc([0,0,0], rad, 180))
    segs.append(rs.AddArc([0,0,0], rad-thick, 180))
    segs.append(rs.AddLine([rad,0,0],[rad-thick,0,0]))
    segs.append(rs.AddLine([-rad,0,0],[-(rad-thick),0,0]))
    crv=rs.JoinCurves(segs, delete_input=True)
    path=rs.AddLine([0,0,0],[0,0,hgt])
    arc=rs.ExtrudeCurve(crv, path)
    rs.CapPlanarHoles(arc)
    rs.DeleteObject(crv)
    rs.DeleteObject(path)
    rs.RotateObject(arc, [0,0,0], orientation)
    rs.MoveObject(arc, insertion)
    return(arc)

# Function to create a terrace
def make_terrace(A=5.0, level=0.7, thick=0.2, hgt=3.0, xcol=2, ycol=3):
    cls =[]
    for i in range(xcol-1):
        for j in range(ycol-1):
            ori = ran.randint(0, 6)
            if ori<4:
                orientation = ori*90
                cls.append(make_arc([i*A+A/2,j*A+A/2,level], A/2, thick, hgt, orientation))
    return(cls)

#function to create stair
def make_stair(start, th, tt, steps, thick, s_width):
    pointlist=[start]
    for i in range(steps):
        pointlist.append([pointlist[-1][0],pointlist[-1][1],pointlist[-1][2]+th])
        pointlist.append([pointlist[-1][0]+tt,pointlist[-1][1],pointlist[-1][2]])
    pointlist.append([pointlist[-1][0],pointlist[-1][1],pointlist[-1][2]-thick])
    pointlist.append([pointlist[0][0],pointlist[0][1],pointlist[0][2]-thick])
    pointlist.append([pointlist[0][0],pointlist[0][1],pointlist[0][2]])
    s_outline=rs.AddPolyline(pointlist)
    path = rs.AddLine(start,[start[0],start[1]+s_width,start[2]]) 
    hull = rs.ExtrudeCurve(s_outline, path)
    rs.CapPlanarHoles(hull)
    rs.DeleteObjects((s_outline,path))
    return(hull)

# Function to create dom-ino
def make_domino(A=A, B=B, thick=thick, hgt=hgt, levels=levels, xcol=xcol, ycol=ycol, f_height=f_height, f_size=f_size):
    f_list = []                          # list of foundations
    c_list = []                          # list of columns
    p_list = []                          # list of plates
    facade_elements = []                 # list of facade elements
    building_height = f_height + (levels-1) * (hgt + thick)  # Calculate building height

    for i in range(levels):
        center_pt[2] = f_height + i * (thick + hgt)
        level = f_height + thick + (i - 1) * (hgt + thick)
        if i == 0:
            f_list = make_foundations(A, f_size, f_height, xcol, ycol)
        else:
            c_list.extend(make_columns(A, level, thick, hgt, xcol, ycol))
            if i % 2:
                c_list.extend(make_terrace(A, level, thick, hgt, xcol, ycol))
        p_list.append(make_box(center_pt, p_width, p_length, thick))

    # Pass building_height to make_facade
    facade_elements, _ = make_facade(building_height=building_height)
    level = f_height + thick + (levels - 1) * (hgt + thick)
    
    #---------------------------------------------------------------------------
    # Staircase parameters
    #---------------------------------------------------------------------------
    steps = int((hgt + thick) / 0.17)
    steps = steps - 1 if steps % 2 else steps
    th = (hgt + thick) / steps
    steps += 2 if th > 0.19 else 0
    th = (hgt + thick) / steps
    tt = 0.28
    s_width = 1.2
    pod_w = B
    start = [0, -(s_width * 2 + f_size / 2), f_height + thick]

    # Create staircase for each level
    stair_l = []
    for i in range(levels):
        start[2] = f_height + thick + i * (thick + hgt)
        if i == levels - 1:
            stair_l.append(make_podest([start[0] - pod_w, start[1] + s_width, start[2] - thick], pod_w, s_width, thick))
        else:
            stair_l.append(make_podest([start[0] - pod_w, start[1], start[2] - thick], pod_w, s_width * 2, thick))
            stair_l.append(make_stair(start, th, tt, int(steps / 2), thick, s_width))
            stair_l.append(make_podest([start[0] + (steps / 2) * tt, start[1], start[2] + (steps / 2) * th - thick], pod_w, s_width * 2, thick))
            stair_l.append(make_stair([start[0] + (steps / 2) * tt, start[1] + s_width, start[2] + (steps / 2) * th], th, -tt, int(steps / 2), thick, s_width))

    return f_list, c_list, p_list, stair_l, facade_elements

# Ensure that the necessary layers exist
def setup_layers():
    layers = {
        "foundation": (100, 40, 19),
        "columns": (160, 90, 90),
        "plates": (139, 110, 60),
        "stairs": (139, 69, 19),
    }
    for layer, color in layers.items():
        if not rs.IsLayer(layer):
            rs.AddLayer(layer)
            rs.LayerColor(layer, color)

# Set up layers
setup_layers()

# Create buildings and assign layers at the end
for i in range(10):
    xran = ran.randint(0, 150)
    yran = ran.randint(0, 150)
    levels = ran.randint(2, 15)
    
    # Create the building
    f_list, c_list, p_list, stair_l, facade_elements = make_domino(levels=levels)
    
    # Move the entire building
    rs.MoveObjects(f_list + c_list + p_list + stair_l + facade_elements, (xran, yran, 0))
    
    # Assign layers at the end
    rs.ObjectLayer(f_list, "foundation")
    rs.ObjectLayer(c_list, "columns")
    rs.ObjectLayer(p_list, "plates")
    rs.ObjectLayer(stair_l, "stairs")