##############################
###  DM2_w24  hu_03_setUp  ###
###  _diag  /  2024 10 17  ###
##############################
import rhinoscriptsyntax as rs
import random, time, sys   ###      
sys.path.append("P:/WWW/dreisi23")     ### add path where "DM_lib.py" can be found !!!                                              
import DM_lib as dm        ### reload(dm)    
############################## 



####Cube

length = 40                                         

if 1:                                           
    for i in range(length):                        
        for j in range(length):
            for k in range(length):
                pnt = rs.CreatePoint(i, j, k)
                points = rs.AddPoint(pnt)
                rs.ObjectColor(points, [148, 0, 211])  




####Sphere1

anz = 5000
rad = 1
center = [ 20, 20, 20]
deltaAngle = 360/anz
vec = [6,6,6]
for i in range(anz):
    angX = deltaAngle*i
    #print angX
    vecX = rs.VectorRotate( vec, angX , [random.uniform(-1,1),random.uniform(-1,1),0]  ) #vector, angle_degrees, axis)
    vecX = rs.VectorScale( vecX, random.uniform( 4,rad) )
    vecX = rs.VectorAdd( vecX, center)
   
    rs.AddPoint( vecX )
    
####Sphere2

anz = 400
rad = 1
center = [ 20, 20, 20]
deltaAngle = 360/anz
vec = [20, 20, 20]
for i in range(anz):
    angX = deltaAngle*i
    #print angX
    vecX = rs.VectorRotate( vec, angX , [random.uniform(-1,1),random.uniform(-1,1),0]  ) #vector, angle_degrees, axis)
    vecX = rs.VectorScale( vecX, random.uniform( 4,rad) )
    vecX = rs.VectorAdd( vecX, center)
   
    rs.AddPoint( vecX )