##############################
###  DM2_w24  hu_03_setUp  ###
###     Lea Koraimann      ###
##############################


import rhinoscriptsyntax as rs
import random

def create_cube_with_points(center, size, num_points, color=None):
    
    points = []
    step = size / (num_points ** (1/3))
    
    for x in range(int(num_points ** (1/3))):
        for y in range(int(num_points ** (1/3))):
            for z in range(int(num_points ** (1/3))):
                px = center[0] - size / 2 + x * step
                py = center[1] - size / 2 + y * step
                pz = center[2] - size / 2 + z * step
                points.append((px, py, pz))

    for point in points:
        rs.AddPoint(point)

    if color:
        for point in points:
            obj_id = rs.coerceguid(rs.AddPoint(point))
            rs.ObjectColor(obj_id, color)

total_points = 30000
points_per_cube = total_points // 3
size = 15

centers = [
    (0, 0, 0),
    (9, 8, 12),
    (18, 16, 2)
]

create_cube_with_points(centers[0], size, points_per_cube)
create_cube_with_points(centers[1], size, points_per_cube, color=(128, 0, 128))
create_cube_with_points(centers[2], size, points_per_cube)

################
rs.ZoomExtents()