import rhinoscriptsyntax as rs
import random as ran
import flipped_classroom_lib as fc
reload(fc)

rs.DeleteObjects(rs.AllObjects())

###############################################################
# Your Panel functions...
###############################################################
def my_panel(points, scale):
    #rs.AddSrfPt(points)
    (P0,P1,P2,P3) = points
    rs.AddLine(P0,P2)
    d_cen1 = rs.VectorDivide(rs.VectorAdd(P0,P2),2)
    rs.AddPoint(d_cen1)
    rs.AddLine(P1,P3)
    d_cen2 = rs.VectorDivide(rs.VectorAdd(P1,P3),2)
    rs.AddPoint(d_cen2)
    center_pt = rs.VectorDivide(rs.VectorAdd(d_cen1,d_cen2),2)
    rs.AddPoint(center_pt)
    if not(rs.PointCompare(d_cen1,d_cen2)):
        rs.AddLine(d_cen1,d_cen2)
    vec1 = rs.VectorSubtract(P0, center_pt)
    vec2 = rs.VectorSubtract(center_pt, P3)
    normal = rs.VectorCrossProduct(vec1, vec2)
    u_normal = rs.VectorUnitize(normal)
    su_normal = rs.VectorScale(u_normal, scale)
    normal_pt = rs.VectorAdd(center_pt, su_normal)
    rs.AddPoint(normal_pt)
    rs.AddLine(center_pt, normal_pt)
    rs.AddSrfPt([P0, normal_pt, P1])
    rs.AddSrfPt([P1, normal_pt, P2])
    rs.AddSrfPt([P2, normal_pt, P3])
    rs.AddSrfPt([P3, normal_pt, P0])

###############################################################

rs.EnableRedraw(False)

#quad = [[0,0,0], [10,0,0], [10,0,12], [0,0,12]]
#quad = [[0,0,0], [10,4,0], [10,4,12], [0,0,12]]
#quad = [[0,0,0], [10,0,0], [10,2,12], [0,2,12]]
#quad = [[0,0,0], [10,0,0], [8,2,12], [2,2,12]]
#quad = [[0,0,0], [10,2,0], [10,-2,12], [0,2,10]]
#     
#result_list = fc.PEF_single_face(quad)
#result_list = fc.PEF_face(7, 10, 8.0, 5.0, 30)
#result_list = fc.PEF_face_w(7, 10, 8.0, 5.0, 30, 2.0)
#result_list = fc.PEF_pantheon()
quad = [[0,0,0], [10,0,0], [10,0,12], [0,0,12]]
result_list = fc.PEF_pantheon()


p_list = result_list[0]  # Pointlist, contains single points
e_list = result_list[1]  # Horizontal Edgelist, contains lists with two points
ve_list = result_list[2] # Vertical Edgelist, contains lists with two points
f_list = result_list[3]  # Facelist, contans lists with four points (quads)
zcol = result_list[4]    # number of levels of faces
xcol = result_list[5]    # number of faces in one level
print "there are " + str(len(p_list)) + " points and " + str(len(f_list)) + " faces on " + str(zcol) + " levels, "+str(xcol)+" per level!"




if 1:
    for i in range(zcol):
        for j in range(xcol):
            points = f_list[i*xcol+j]
            scale = ran.uniform(0.5, 8.5)
            myp= my_panel(points, scale)
            #rs.AddSrfPt(f_list[i*xcol+j])


#######oben########
if 1:
    for i in range(zcol):
        for j in range(xcol):
            points = f_list[i*xcol+j]
            (P0,P1,P2,P3) = points
            cpoints=[P0,P1,P2,P3,P0]
            plinea = rs.AddCurve(cpoints,2)
            plineb = rs.AddCurve(cpoints,3)
            #loftsurf = rs.AddLoftSrf([plinea, plineb])
            rs.Command("-_Loft selid {} selid {} _Enter _Enter _Enter".format(plinea, plineb), False)
            loftsurf = rs.FirstObject()
            tkt = rs.OffsetSurface(loftsurf, 2, both_sides=True, create_solid=True)
            tdt = rs.CopyObject(tkt,[0,0,0])
            two = rs.RotateObject(tkt, [0,0,0],180,[1,0,0])
            rs.MoveObject(two, [0,0,-18])
            rs.MoveObject(tdt, [0,0,-19])