import rhinoscriptsyntax as rs
import random as ran

allobjs = rs.AllObjects()
rs.DeleteObjects(allobjs)

rs.EnableRedraw(False)

A=5
B=A/3
thick=0.2
hgt=2.7
xcol=3
ycol=4
levels=3
f_height=0.5   
f_size=0.8    

center_pt = [A*(xcol-1)/2, A*(ycol-1)/2, f_height]                 #Mittelpunkt der Bodenplatte
p_width = A*(xcol-1)+2*B                                           #Breite der Bodenplatte
p_lenght = A*(ycol-1) + f_size                                  #Laenge der Bodenplatte

#create box at centre point
def make_box(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,(-xsize/2,-ysize/2,0))
    rs.MoveObject(box,insertion)
    return(box)


#create box at corner point
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)



#make_box(12,16,35)
#make_box([10,3,4],ysize=9)
#make_box([10,3,4],12,16,35)

#Funktion um Fundament zu kreieren

def make_foundations(A=A, f_size=f_size, f_height=f_height, xcol=xcol, ycol=ycol):
    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)

#Funktion um Saeulen zu kreieren

def make_columns(A=A, level=0.7, thick=thick, hgt=hgt, xcol=xcol, ycol=ycol):
    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)

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)

def make_terrace(A=A, level=0.7, thick=thick, hgt=hgt, xcol=xcol, ycol=ycol):
    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)


#stairs
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[-1] [1],pointlist[0] [2]-thick])
    pointlist.append([pointlist[0] [0],pointlist[-1] [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)
    return(hull)


#Dom-Ino wird gebaut

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 = []                                                        #Liste der Fundamente
    c_list = []                                                        #Liste der Saeulen
    p_list = []                                                        #Liste der Platten
    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))
            c_list.extend(make_terrace(A, level, thick, hgt, xcol, ycol))
        p_list.append(make_box(center_pt, p_width, p_lenght, thick))
    level = f_height + thick +(levels-1)*(hgt+thick)
    c_list.extend(make_terrace(A, level, thick, hgt, xcol, ycol))
    ##################################################################
    #calculate stair values#
    ########################################################
    steps = int((hgt+thick)/.17)
    if steps%2:
        steps = steps-1
    th = (hgt+thick)/steps
    if(th>.19):
        steps = steps+2
        th = (hgt+thick)/steps
    ###################################################################
    #step parameters
    #########################################################
    tt = .3                         #step size
    s_width = 1.2                   #stair width
    pod_w = B                       #depht of landing
    start = [pod_w,-(s_width*2+f_size/2), f_height+thick]        #startpoint of stair
    
    #make_stair(start, th, tt, steps, thick, s_width):
    stair_1 = []
    for i in range(levels):
        start[2] = f_height+thick + i*(thick+hgt)                #z-wert bei jeder Iteration neu gesetzt
        if i==levels-1:         #letztes Podest
            stair_1.append(make_podest([start[0]-pod_w, start[1]+s_width, start[2]-thick], pod_w, s_width, thick))
        else:
            stair_1.append(make_podest([start[0]-pod_w, start[1],start[2]-thick], pod_w, s_width*2, thick))
            stair_1.append(make_stair(start, th, tt, int(steps/2), thick, s_width))
            stair_1.append(make_podest([start[0]+(steps/2)*tt, start[1], start[2]+(steps/2)*th-thick], pod_w, s_width*2, thick))
            stair_1.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_1)

(f_list, c_list, p_list, stair_1) = make_domino()




rs.AddLayer("Fundament")
rs.LayerColor("Fundament", (220,60,60))
rs.ObjectLayer(f_list, "Fundament")

rs.AddLayer("Saeule")
rs.LayerColor("Saeule", (60,220,6))
rs.ObjectLayer(c_list, "Saeule")

rs.AddLayer("Platten")
rs.LayerColor("Platten", (60,60,220))
rs.ObjectLayer(p_list, "Platten")

# Function to create lamella facade while skipping the area on the staircase side where the first podium is located
def create_lamella_facade(A=A, xcol=xcol, ycol=ycol, thick=0.1, lamella_height=hgt * levels + thick * (levels - 1), gap=0.3):
    lamellas = []
    lamella_width = 0.2
    x_start = -B
    y_start = -f_size / 2
    z_start = 0

    # Define the specific area to skip (first podium position)
    staircase_x_min = B  # Adjusted based on the podium's approximate x-position
    staircase_x_max = staircase_x_min + 1.2  # Width of the podium's area to skip

    num_lamellas_x = int((A * (xcol - 1) + 2 * B) / (lamella_width + gap))
    num_lamellas_y = int((A * (ycol - 1) + f_size) / (lamella_width + gap))

    # Front Side (Staircase Side) - Skip lamellae overlapping the first podium only
    for i in range(num_lamellas_x):
        x_pos = x_start + i * (lamella_width + gap)
        if not (staircase_x_min <= x_pos <= staircase_x_max):  # Exclude lamellae in the range of the first podium
            lamellas.append(make_box([x_pos, y_start, z_start], lamella_width, thick, lamella_height))

    # Back Side
    for i in range(num_lamellas_x):
        x_pos = x_start + i * (lamella_width + gap)
        lamellas.append(make_box([x_pos, y_start + (A * (ycol - 1) + f_size), z_start], lamella_width, thick, lamella_height))

    # Left and Right Sides
    for j in range(num_lamellas_y):
        y_pos = y_start + j * (lamella_width + gap)
        lamellas.append(make_box([x_start, y_pos, z_start], thick, lamella_width, lamella_height))
        lamellas.append(make_box([x_start + (A * (xcol - 1) + 2 * B), y_pos, z_start], thick, lamella_width, lamella_height))
        
    return lamellas

lamella_facade = create_lamella_facade()

# Add layer for the lamella facade
rs.AddLayer("Fassade", (150, 150, 150))
rs.ObjectLayer(lamella_facade, "Fassade")

import math

# Function to duplicate building in random positions within a larger area and ensure they do not overlap
def duplicate_building_randomly(f_list, c_list, p_list, stair_1, lamella_facade, num_copies=20, min_distance=15, area_range=200):
    building_parts = f_list + c_list + p_list + stair_1 + lamella_facade
    positions = []  # Track positions to avoid overlaps

    for _ in range(num_copies):
        valid_position = False

        while not valid_position:
            # Generate random x, y coordinates within the area range
            x_offset = ran.uniform(-area_range, area_range)
            y_offset = ran.uniform(-area_range, area_range)
            offset = [x_offset, y_offset, 0]  # z = 0 for ground level placement

            # Check if this position is far enough from all previous positions
            valid_position = all(
                math.sqrt((x_offset - pos[0])**2 + (y_offset - pos[1])**2) >= min_distance for pos in positions
            )

        # If valid, add to positions list and copy the building parts
        positions.append((x_offset, y_offset))
        rs.CopyObjects(building_parts, offset)

# Duplicate the building 5 times at random positions within a large area, ensuring no overlap
duplicate_building_randomly(f_list, c_list, p_list, stair_1, lamella_facade, num_copies=20, min_distance=15, area_range=200)

# Enable redraw to update view
rs.EnableRedraw(True)