import rhinoscriptsyntax as rs
import random

alleobjekte = rs.AllObjects()
rs.DeleteObjects(alleobjekte)
rs.EnableRedraw(0)

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

coords = []
anzahl = 20
xVal = 4
yVal = 4
zVal = 4

for i in range(anzahl):
    for j in range(anzahl):
        for k in range(anzahl):
            x = i*xVal + random.uniform(0.01, xVal) -40
            y = j*yVal + random.uniform(0.01, yVal) -40
            z = k*zVal + random.uniform(0.01, zVal) -40
            
            cor = [x, y, z]
            coords.append(cor)

points = rs.AddPoints(coords)

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

coords2 = []
anzahl2 = 2
xVal2 = 120
yVal2 = 120
zVal2 = 120

for a in range(anzahl2):
    for b in range(anzahl2):
        for c in range(anzahl2):
            x2 = a*xVal2 - 60
            y2 = b*yVal2 - 60
            z2 = c*zVal2 -60
            
            cor2 = [x2, y2, z2 ]
            coords2.append(cor2)

points2 = rs.AddPoints(coords2)

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

coordsrod = []
for a1 in range(5):
    for b1 in range(5):
        for c1 in range(40):
            x3 = a1*xVal-8
            y3 = b1*yVal-8
            z3 = c1*zVal-80
            
            corrod  = (x3+80, y3+80, z3 )
            corrod1 = (x3+80, y3-80, z3 )
            corrod2 = (x3-80, y3+80, z3 )
            corrod3 = (x3-80, y3-80, z3 )
            coordsrod.append(corrod)
            coordsrod.append(corrod1)
            coordsrod.append(corrod2)
            coordsrod.append(corrod3)

pointsrod = rs.AddPoints(coordsrod)

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

coordsrod0 = []
for a2 in range(40):
    for b2 in range(5):
        for c2 in range(5):
            x4 = a2*xVal-80
            y4 = b2*yVal-8
            z4 = c2*zVal-8
            
            corrod00 = (x4, y4+80, z4+80 )
            corrod01 = (x4, y4-80, z4+80 )
            corrod02 = (x4, y4+80, z4-80 )
            corrod03 = (x4, y4-80, z4-80 )
            coordsrod0.append(corrod00)
            coordsrod0.append(corrod01)
            coordsrod0.append(corrod02)
            coordsrod0.append(corrod03)

pointsrod0 = rs.AddPoints(coordsrod0)

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

coordsrod00 = []
for a3 in range(5):
    for b3 in range(40):
        for c3 in range(5):
            x5 = a3*xVal-8
            y5 = b3*yVal-80
            z5 = c3*zVal-8
            
            corrod000 = (x5+80, y5, z5+80 )
            corrod001 = (x5-80, y5, z5+80 )
            corrod002 = (x5+80, y5, z5-80 )
            corrod003 = (x5-80, y5, z5-80 )
            coordsrod00.append(corrod000)
            coordsrod00.append(corrod001)
            coordsrod00.append(corrod002)
            coordsrod00.append(corrod003)

pointsrod00 = rs.AddPoints(coordsrod00)

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

curve1 = rs.AddCurve(coordsrod, 1)
curve2 = rs.AddCurve(coordsrod0, 1)
curve3 = rs.AddCurve(coordsrod00, 1)

rs.ObjectColor(points, [150, 0, 0] )
rs.ObjectColor(pointsrod, [150, 150, 0])
rs.ObjectColor(pointsrod0, [0, 150, 150])
rs.ObjectColor(pointsrod00, [150, 0 ,150])

rs.ObjectColor(curve1, [150, 150, 0])
rs.ObjectColor(curve2, [0, 150, 150])
rs.ObjectColor(curve3, [150, 0 ,150])