##################################
### DM2_w22 # AGruber@tugraz.at ###
### ue_07 UN_headquaters  NYC   ###
###################################


##############################
import rhinoscriptsyntax as rs      ###
import random, time, sys            ###
sys.path.append("P:/")              ### add LW P:/ to rhino's default search path AND
sys.path.append("P:/WWW/default114/dm2/")     ### add YOUR LW P:/WWW/user/ to rhino's default search path ...
import DM2_lib as dm                ### ... thus DM2_lib.py can be found !
#import OSM_lib as osm              ### "osm" is alias 4 calling defs @ _osm_lib
##############################      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.printDisplay(1)                                      # nomen est omen
rs.EnableRedraw(0) ### MAC >> uncomment !


###_________________________________________# import OSM objects

###_________________________________________#
### basic settings for grid to fit UN_slab  #
###                                         # ! no need 2 change !
floors = H  = dm.H = 40                     # def=40 / incl roof
slabs  = L  = dm.L = 11                     # def=11
depth  = D  = dm.D =  4                     # def= 4 / division in building_depth
floorHeight = fH = dm.fH = 4.0              # def= 4.0 / 4.0*(H-1) = 156 meters
                                            #
#############                               # get gridCoords L*D*H = 1760
UnoGridCoords = dm.UnoGridCoords            # get from DM2_lib as dm
#################___________________________#


dm.newEmptyLayer("UNO::slabs",  [ 20, 60,100])
dm.newEmptyLayer("UNO::points", [150,150,200])
dm.newEmptyLayer("UNO::geo", [150, 20,150])
dm.newEmptyLayer("UNO::crv", [220, 20, 20])
dm.newEmptyLayer("UNO::tmp", [120,120,140])


##################
depthVec  = dVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 2, 0) ))
lengthVec = lVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) ))
##################

if 1:
    rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) )
    rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 1, 0) )
    rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 0, 1) )
    
    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.5 )
    
    p0 = dm.getUnoCoord( 0,   0,   0)
    p1 = dm.getUnoCoord( L-1, 0,   0)
    p2 = dm.getUnoCoord( L-1, D-1, 0)
    p3 = dm.getUnoCoord( 0,   D-1, 0)
    unoBaseCoords = [p0, p1, p2, p3, p0]
    unoBaseCrv = rs.AddCurve( unoBaseCoords, 1)
if 0:
    ################################
    pXX = dm.getUnoCoord(5, 0, 40-1)
    lenx = 30.1
    rs.AddLine( pXX , rs.VectorAdd(pXX, rs.VectorScale( dVec, lenx ) ) )
    ####################################################################

if 0:
    rs.CurrentLayer("UNO::points")
    rs.AddPoints( UnoGridCoords )

if 0:
    rs.CurrentLayer("UNO::crv")
    for i in range(64):
        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#[:]
    random.shuffle( coords )
    rs.ObjectColor(rs.AddCurve( coords[0:], 2 ), [50,50,50] )


#############################
##### demo karla m. / DM2_w20
def karla(anzahl = 128, makeAll = 0, delCircs = 1):
    rs.CurrentLayer("UNO::crv")  
    p1 = dm.getUnoCoord( L-1, 0,   0)
    if 0:
        pX = dm.getUnoCoord( 4, 2, 20 )
        dirDepthvec = rs.VectorSubtract( dm.getUnoCoord( 0,   0,   0), dm.getUnoCoord( 0,   D-1, 0) )
        lenDepth = rs.VectorLength( dirDepthvec )
        print "lenDepth",lenDepth
        dirDepthvec = rs.VectorScale( rs.VectorUnitize( dirDepthvec ), lenDepth )
        p1 = rs.VectorAdd( p1, [0,0,lenDepth*2.5] ) 
        #crv = rs.AddCurve( [p0, pX, pVorFass, p1, p7, p2, p3] , 2 )
    karlaCoords = []
    p0 = rs.VectorAdd( dm.getUnoCoord( 0, 0, 10 ), rs.VectorScale(lVec, 10) )
    karlaCoords = [p0]
    counter = 0
    for h in range(0, H-1, 2):
        counter = counter+1
        if counter%4 == 1:
            karlaCoor = rs.VectorAdd(dm.getUnoCoord( int((L-1)/2), 0, h ),  rs.VectorScale(dVec, random.uniform(0,50)))
        if counter%4 == 2:
            karlaCoor = dm.getUnoCoord( L-1, int((D-1)/2), h )
        if counter%4 == 3:
            karlaCoor = rs.VectorAdd(dm.getUnoCoord( int((L-1)/2), D-1, h ), rs.VectorScale(dVec, random.uniform(0,50)*-1))
        if counter%4 == 0:
            karlaCoor = dm.getUnoCoord( 0, int((D-1)/2), h )
        karlaCoords.append(karlaCoor)
    
    karlaCoords.append(karlaCoords[0])
    crv = rs.AddCurve(karlaCoords, 2)

    anz = 110
    circDiv = 64
    allCircCoords = []
    dom = rs.CurveDomain(crv)[1]
    for p in range(anzahl+1):
        fac = p/anzahl
        coord = rs.EvaluateCurve( crv, dom*fac )
        #pnt = rs.AddPoint( coord )
        para = rs.CurveClosestPoint(unoBaseCrv, coord)
        coordOnCrv = rs.EvaluateCurve( unoBaseCrv, para )
        vecXY = dm.vecCurvePerpXY( crv, dom*fac )
        if makeAll:
            rs.AddCurve( [coord, coordOnCrv ] )
        planeX = dm.plnCurvePerp( crv, dom*fac )
        #planeX = rs.CurvePerpFrame(crv, dom*fac )
        angX = rs.Angle([0,0,0], vecXY )
        #print "ang", angX
        #rs.AddLine( planeX[0], rs.VectorAdd(planeX[0], planeX[1]) )
        #print planeX[0], "...", planeX[1]
        #planeX= rs.RotatePlane(planeX, angX[2], planeX.ZAxis)
        circ = rs.AddCircle(planeX, 3.4 )
        coords = rs.DivideCurve( circ, circDiv, 0)
        rs.AddLine(coord, coords[0])
        if delCircs:
            rs.DeleteObject(circ)
        allCircCoords.append(coords)
        vecXY = dm.vecCurvePerpXY( crv, dom*fac )
        #rs.AddLine( planeX[0], rs.VectorAdd(planeX[0], rs.VectorScale(vecXY, 3.4)) )
    if makeAll:
        for i in range(24):
            corrdsX = []
            for circList in allCircCoords:
                corrdsX.append(circList[i])
            rs.ObjectColor(rs.AddCurve(corrdsX, 1), [100,160,160] )

    anzPoints = len(rs.ObjectsByType(1))
    print "anzPoints", anzPoints

#karla( anzahl = 2**8, makeAll=1, delCircs=0 )

#_____________________________here you go:
dm.newEmptyLayer("UNO::PROJ", [100,100,100])




#____________________________:here you end

###################################
rs.ZoomExtents()
rs.EnableRedraw(1) ### 4 the MACs
dm.printDisplay(1)
dm.newEmptyLayer("Default")