###################################
### DM2_w24 # AGruber@tugraz.at ###
### ue_05 UN_headquaters  NYC   ###
###################################


##############################
import rhinoscriptsyntax as rs
import random, time, sys   ###
sys.path.append("P:/")
sys.path.append("C:/Users/annas/OneDrive/Dokumente/uni-dokumente/dm2")    ###
import DM_lib as dm        ### 
reload(dm)
##############################      

rs.UnitSystem(4)                                        # km = 5, meters = 4, cm = 3 etc
rs.ShowGrid(None, 0)                                    # grid > 0 = off
rs.ShowGridAxes(None, 1)                                # y/y/z axen display > 0/1 = off/on
rs.ViewDisplayMode(rs.CurrentView(), "wireframe")
rs.Command("cplane w t enter", 0)						# cPlane World Top
dm.PointRadius(displayModeX=0, rad=3, styl=3)
dm.printDisplay(state=1) ### nomen est omen / printDisplay(state=1, scale=1000, thickness=1, color="Display")
rs.EnableRedraw(0)       ### MAC >> problem !



###_________________________________________#
### basic settings for grid to fit UN_slab  #
###                                         # ! no need 2 change !
floors = H  = dm.H = 40                     # default=40 / incl roof
slabs  = L  = dm.L = 11                     # default=11
depth  = D  = dm.D =  4                     # default= 4 / division in building_depth
floorHeight = fH = dm.fH = 4.0              # default= 4.0 / 4.0*(H-1) = 156 meters
                                            #
#############                               # get gridCoords L*D*H = 1760
UnoGridCoords = dm.getUnoGridCoords()       # get from DM_lib as dm
#################___________________________#


dm.newEmptyLayer("UNO::setUp::points", [150,150,200]  )
dm.newEmptyLayer("UNO::setUp::geo",    [150, 20,150]  )
dm.newEmptyLayer("UNO::setUp::crv",    [200, 100, 100])
dm.newEmptyLayer("UNO::setUp",         [120,120,140]  )


################
lengthVec = lVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) )) ## rs.AddPoint( dm.getUnoCoord(1, 0, 0) )
depthVec  = dVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 1, 0), dm.getUnoCoord(0, 0, 0) )) ## rs.AddPoint( dm.getUnoCoord(0, 1, 0) )
################

if 1: ### DEMOs / SETUP >> dont' exec @ homework !
    if 0:
        rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(10, 0,  0) )
        rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord( 0, 3,  0) )
        rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord( 0, 0, 39) )
        
        rs.CurveArrows(rs.AllObjects()[0], arrow_style=2)
        rs.CurveArrows(rs.AllObjects()[1], arrow_style=2)
        rs.CurveArrows(rs.AllObjects()[2], arrow_style=2)
        
        rs.ObjectColor(rs.AllObjects()[2], [0, 222, 0] )
        rs.ObjectColor(rs.AllObjects()[1], [222, 0, 0] )
        rs.ObjectColor(rs.AllObjects()[0], [0, 0, 222] )
        
        rs.ObjectPrintWidth( rs.AllObjects()[0:3], 1.0 )

        p0 = dm.getUnoCoord(  0, 0, 0)
        p1 = dm.getUnoCoord( 10, 0, 0)
        p2 = dm.getUnoCoord( 10, 3, 0)
        p3 = dm.getUnoCoord(  0, 3, 0)
        unoBaseCoords = [p0, p1, p2, p3, p0]
        unoBaseCrv = rs.AddCurve( unoBaseCoords, 1)
        unoTopCoords = [dm.getUnoCoord(  0, 0, 39),
                        dm.getUnoCoord( 10, 0, 39),
                        dm.getUnoCoord( 10, 3, 39),
                        dm.getUnoCoord(  0, 3, 39),
                        dm.getUnoCoord(  0, 0, 39),
                       ]
        unoTopCrv = rs.AddCurve( unoTopCoords, 1)
        rs.ZoomExtents()
    if 0:
        ##############################
        rs.CurrentLayer("UNO::setUp::points")
        rs.AddPoints( UnoGridCoords )
        for i in range(L): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(i, 0, 0) ), [0, 222, 0] )
        for i in range(D): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(0, i, 0) ), [222, 0, 0] )
        for i in range(H): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(0, 0, i) ), [0, 0, 222] )
        
        ##############################
    
    
    if 0:
        rs.CurrentLayer("UNO::setUp::crv")
        dm.UN_slab(showAll=1)
        for i in range(8):
            p0 = dm.getUnoCoord( random.randint(0,L-1), random.randint(0,D-1), random.randint(0,H-1) )
            p1 = dm.getUnoCoord( random.randint(0,L-1), random.randint(0,D-1), random.randint(0,H-1) )
            rs.AddCurve( [p0, p1] )
            dm.textDots(  [p0, p1] )
        rs.ZoomExtents()
    if 0:
        coords = UnoGridCoords[:] ### make COPY - not equal !
        random.shuffle( coords )
        rs.AddCurve( coords[0:], random.randint(1, 3) )

#################
### abbreveation:
def gC (l=1, d=2, h=3):
    cor = dm.getUnoCoord(l, d, h)
    return cor
## rs.AddPoint( gC(1,0,1) )
#################



#_____________________________here you go:
dm.newEmptyLayer("myPROJ", [100,100,100])

rs.CurrentLayer("UNO::setUp::crv")

#wiederholung der uebung (hab schon wieder vergessen was wir gemacht haben)

if(0):
    p0 = gC(l = 10, d = 0, h = 39) ###holen von verschiedenen Punkten
    rs.AddPoint(p0)
    p1 = rs.VectorAdd(p0, rs.VectorScale( dVec, 100)) ###holen von Punkten ausserhalb des Gebaeudes
    p2 = gC( 5,3,39 )
    p3 = gC( 8,0,20 )
    p4 = rs.VectorAdd(p3, rs.VectorScale( dVec, -100))
    p5 = gC( 0,0,0 )
    
    coords = [p0, p1, p2,p3,p4,p5] ###Punkte werden zu Liste hinzugefuegt 
    #rs.AddCurve( coords,1)
    crvDiv = 44                    ### Curve segments werden defined 
    crv = rs.AddCurve( coords, 2)
    coords = rs.DivideCurve(crv, crvDiv, create_points = 0) ### Curve wird in 44 teile geteilt
    rs.AddPoints( coords ) 
    
    for cor in coords:
        pass 
        #rs.AddCircle( cor, 5.5 ) ###erzeugt an jedem punkt einen kreis mit radius 5.5
    
    dom = rs.CurveDomain( crv ) ####Returns the domain (Bereich) of a curve object
    print "dom =", dom   ### um curve domain anzeigen zu lassen
    dom = rs.CurveDomain( crv )[1] ### Returns the domain maximum of a curve 
    print "dom =", dom 
    #pX = rs.EvaluateCurve( crv, dom*0.5 ) 
    #planeX = rs.CurvePerpFrame( crv, dom*0.5) ####returns the perpendicular plane at a parameter of a curve
    #rs.AddCircle( planeX, 5.5) 
    
    #paraX = rs.CurveClosestPoint( crv, coords[10] ) ###Punkt am naechsten zum zehnten Punkt in der Liste
    #planeX = rs.CurvePerpFrame( crv, paraX) ### plane definieren 
    #rs.AddCircle(planeX, 5.5)
    
    circDiv = 16
    allCoordsLists = [] ###leere Liste definieren
    for cor in coords:
        paraX = rs.CurveClosestPoint( crv, cor ) 
        planeX = rs.CurvePerpFrame( crv, paraX ) 
        circ = rs.AddCircle( planeX, 5.5)
        circCoords = rs.DivideCurve( circ, circDiv, 0) ###teilt die kreise in segmente
        rs.DeleteObject( circ) ####kreise werden nicht mehr gebraucht
        allCoordsLists.append( circCoords ) 
    
    i = 0 
    for i in range(circDiv):
        coords = []
        for list in allCoordsLists:
            pX = list[i]
            coords.append( pX )
        rs.AddCurve( coords, 2 ) 
            

#######ab jetzt hu

####Holen der Punkte fuer Rollercoaster 
p0 = gC(l = 8, d = 0, h = 0)
p1 = gC(8,0,35)
p2 = gC(8,3,39)
p3 = rs.VectorAdd(gC(8,3,30), rs.VectorScale(dVec, 35))
p4 = rs.VectorAdd(gC(5,3,28), rs.VectorScale(dVec, 60))
p5 = gC(0,3,25)
p6 = rs.VectorAdd(gC(5,1,22), rs.VectorScale(dVec, -60))
p7 = rs.VectorAdd(gC(10,3,19), rs.VectorScale(lVec, -35))
p8 = rs.VectorAdd(gC(5,3,16), rs.VectorScale(dVec, 40))
p9 = gC(3,3,13)
p10 = rs.VectorAdd(gC(5,1,10), rs.VectorScale(dVec, -40))
p11 = rs.VectorAdd(gC(10,3,7), rs.VectorScale(lVec, -35))
p12 = rs.VectorAdd(gC(5,3,4), rs.VectorScale(dVec, 20))
p13 = gC(5,3,1)
p14 = rs.VectorAdd(gC(5,1,0), rs.VectorScale(dVec, -35))
p15 = rs.VectorAdd(gC(10,3,0), rs.VectorScale(lVec, -10))

crv = rs.AddCurve( [p0,p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15],3 )

########Unterteilung der Curve in Segments
crvDiv = 100
coords = rs.DivideCurve(crv, crvDiv, create_points = 0)
#rs.AddPoints(coords)

#####Returning domain of curve
dom = rs.CurveDomain(crv)

#####Erstellen von Kreisen an Curve Segments und Unterteilen in segments
circDiv = 18
allCoordsLists = []
for cor in coords:
    paraX = rs.CurveClosestPoint(crv, cor)
    planeX = rs.CurvePerpFrame(crv, paraX)
    circ = rs.AddCircle(planeX, 5.5)
    circCoords = rs.DivideCurve( circ, circDiv, 0)
    #rs.AddPoints(circCoords)
    rs.DeleteObject(circ)
    allCoordsLists.append(circCoords)
    
#######Stuetzen
#Kreise an bestimmten Segments erzeugen und diese in segments unterteilen
paraA = rs.CurveClosestPoint(crv,coords[26]) 
planeA = rs.CurvePerpFrame(crv, paraA)
circA = rs.AddCircle(planeA, 5)

paraB = rs.CurveClosestPoint(crv,coords[36]) 
planeB = rs.CurvePerpFrame(crv, paraB)
circB = rs.AddCircle(planeB, 5)

paraC = rs.CurveClosestPoint(crv,coords[44]) 
planeC = rs.CurvePerpFrame(crv, paraC)
circC = rs.AddCircle(planeC, 5)

paraD = rs.CurveClosestPoint(crv,coords[59]) 
planeD = rs.CurvePerpFrame(crv, paraD)
circD = rs.AddCircle(planeD, 5)

paraE = rs.CurveClosestPoint(crv,coords[66]) 
planeE = rs.CurvePerpFrame(crv, paraE)
circE = rs.AddCircle(planeE, 5)

paraF = rs.CurveClosestPoint(crv,coords[72]) 
planeF = rs.CurvePerpFrame(crv, paraF)
circF = rs.AddCircle(planeF, 5)

paraG = rs.CurveClosestPoint(crv,coords[83]) 
planeG = rs.CurvePerpFrame(crv, paraG)
circG = rs.AddCircle(planeG, 5)

allCirc = [circA, circB, circC, circD, circE, circF, circG]

allcircDiv = 4
allCirccoords = []
circCoordsA = rs.DivideCurve(circA, allcircDiv, 0)
circCoordsB = rs.DivideCurve(circB, allcircDiv, 0)
circCoordsC = rs.DivideCurve(circC, allcircDiv, 0)
circCoordsD = rs.DivideCurve(circD, allcircDiv, 0)
circCoordsE = rs.DivideCurve(circE, allcircDiv, 0)
circCoordsF = rs.DivideCurve(circF, allcircDiv, 0)
circCoordsG = rs.DivideCurve(circG, allcircDiv, 0)


rs.ObjectColor(allCirc, [200,200,200])
rs.ObjectPrintWidth(allCirc,3)

#####Stuetzen schaffen an den definierten Kreisen
lineA = rs.AddLine(circCoordsA[3], gC(8,0,0))
lineB = rs.AddLine(circCoordsB[3], gC(8,0,0))
lineC = rs.AddLine(circCoordsC[3], gC(8,0,0))
lineD = rs.AddLine(circCoordsD[3], gC(8,0,0))
lineE = rs.AddLine(circCoordsE[3], gC(8,0,0))
lineF = rs.AddLine(circCoordsF[3], gC(8,0,0))
lineG = rs.AddLine(circCoordsG[3], gC(8,0,0))

allLines = [lineA, lineB, lineC, lineD, lineE, lineF, lineG]
rs.ObjectColor(allLines, [200,200,200])
rs.ObjectPrintWidth(allLines, 1)

for i in range(circDiv):
    coords =[]
    for list in allCoordsLists:
        pX = list[i]
        coords.append( pX )
    rs.AddCurve(coords,2)
    

####Umrahmung
pkt0 = gC(0,0,0)
pkt1 = gC(10,0,0)
pkt3 = gC(0,3,0)
pkt4 = gC(10,3,0)
pkt5 = gC(0,0,39)
pkt6 = gC(10,0,39)
pkt7 = gC(0,3,39)
pkt8 = gC(10,3,39)

unten = [pkt0, pkt1, pkt4, pkt3, pkt0]
CrvUnten = rs.AddCurve(unten,1)
oben = [pkt5, pkt6, pkt8, pkt7, pkt5]
CrvOben = rs.AddCurve(oben, 1)
links = [pkt0, pkt3, pkt7, pkt5, pkt0]
CrvLinks = rs.AddCurve(links, 1)
rechts = [pkt1, pkt4, pkt8, pkt6, pkt1]
CrvRechts = rs.AddCurve(rechts,1)

Umrahmung = [CrvUnten, CrvOben, CrvLinks, CrvRechts]
rs.ObjectColor(Umrahmung, [225,197,192])



















#____________________________:here you end

###################################
dm.zA( proz=0.95 )
rs.EnableRedraw(1) ### 4 the MACs
#dm.printDisplay(1)
dm.newEmptyLayer("Default")