"""
Ausgehöhltes Dodekaeder
05.05.2021
www.3d-meier.de
"""

import c4d
import math

# Variablen und Konstanten
Titel = 'Ausgehöhltes Dodekaeder' # Name

NP = 32    # Anzahl Punkte
N3 = 60    # Anzahl Dreiecke

Radius = 2 # Radius der Eckpunktkugeln       

a = (3 + math.sqrt(5))/4
b = (1 + math.sqrt(5))/4
c = (math.sqrt(5) - 1)/4

Punkte = [(   b,   -b,    b),   # Dodekaeder
          (   b,   -b,   -b),   # Dodekaeder
          (  -b,   -b,    b),   # Dodekaeder
          (  -b,   -b,   -b),   # Dodekaeder
          (  -b,    b,    b),   # Dodekaeder
          (  -b,    b,   -b),   # Dodekaeder
          (   b,    b,    b),   # Dodekaeder
          (   b,    b,   -b),   # Dodekaeder
          ( 0.0,    a,  0.5),   # Dodekaeder
          ( 0.0,    a, -0.5),   # Dodekaeder
          ( 0.0,   -a,  0.5),   # Dodekaeder
          ( 0.0,   -a, -0.5),   # Dodekaeder
          ( 0.5,  0.0,    a),   # Dodekaeder
          ( 0.5,  0.0,   -a),   # Dodekaeder
          (-0.5,  0.0,    a),   # Dodekaeder
          (-0.5,  0.0,   -a),   # Dodekaeder
          (   a,  0.5,  0.0),   # Dodekaeder
          (   a, -0.5,  0.0),   # Dodekaeder
          (  -a,  0.5,  0.0),   # Dodekaeder
          (  -a, -0.5,  0.0),   # Dodekaeder
          ( 0.0,  c,    0.5),   # Ikosaeder
          ( 0.0,  c,   -0.5),   # Ikosaeder
          ( 0.0, -c,    0.5),   # Ikosaeder
          ( 0.0, -c,   -0.5),   # Ikosaeder
          ( c,    0.5,  0.0),   # Ikosaeder
          ( c,   -0.5,  0.0),   # Ikosaeder
          (-c,    0.5,  0.0),   # Ikosaeder
          (-c,   -0.5,  0.0),   # Ikosaeder
          (   0.5,  0.0,  c),   # Ikosaeder
          (  -0.5,  0.0,  c),   # Ikosaeder
          (   0.5,  0.0, -c),   # Ikosaeder
          (  -0.5,  0.0, -c)]   # Ikosaeder
  
Dreiecke = [(31, 18,  5,  5),
            (31,  5, 15, 15),
            (31, 15,  3,  3),
            (31,  3, 19, 19),
            (31, 19, 18, 18),
            (21,  5,  9,  9),
            (21,  9,  7,  7),
            (21,  7, 13, 13),
            (21, 13, 15, 15),
            (21, 15,  5,  5),
            (23, 15, 13, 13),
            (23, 13,  1,  1),
            (23,  1, 11, 11),
            (23, 11,  3,  3),
            (23,  3, 15, 15),
            (27, 11, 10, 10),
            (27, 10,  2,  2),
            (27,  2, 19, 19),
            (27, 19,  3,  3),
            (27,  3, 11, 11),
            (25,  1, 17, 17),
            (25, 17,  0,  0),
            (25,  0, 10, 10),
            (25, 10, 11, 11),
            (25, 11,  1,  1),
            (29, 18, 19, 19),
            (29, 19,  2,  2),
            (29,  2, 14, 14),
            (29, 14,  4,  4),
            (29,  4, 18, 18),
            (26,  9,  5,  5),
            (26,  5, 18, 18),
            (26, 18,  4,  4),
            (26,  4,  8,  8),
            (26,  8,  9,  9),
            ( 9,  8, 24, 24),
            ( 7,  9, 24, 24),
            (16,  7, 24, 24),
            ( 6, 16, 24, 24),
            ( 8,  6, 24, 24),
            (22,  2, 10, 10),
            (22, 10,  0,  0),
            (22,  0, 12, 12),
            (22, 12, 14, 14),
            (22, 14,  2,  2),
            (30, 17,  1,  1),
            (30,  1, 13, 13),
            (30, 13,  7,  7),
            (30,  7, 16, 16),
            (30, 16, 17, 17),
            (20, 14, 12, 12),
            (20, 12,  6,  6),
            (20,  6,  8,  8),
            (20,  8,  4,  4),
            (20,  4, 14, 14),
            (28, 12,  0,  0),
            (28,  0, 17, 17),
            (28, 17, 16, 16),
            (28, 16,  6,  6),
            (28,  6, 12, 12)]
            
Faktor = 100         # Skalierungsfaktor
k = 1                # Korrekturfaktor

#************************************************************************
def CreateNullobjekt1():
    obj = c4d.BaseObject(c4d.Onull)
    obj.SetName(Titel)
    obj.Message(c4d.MSG_UPDATE)
    return obj

#************************************************************************
def CreateNullobjekt2():
    obj = c4d.BaseObject(c4d.Onull)
    obj.SetName('Ecken')
    obj.Message(c4d.MSG_UPDATE)
    return obj

#************************************************************************
def CreateNullobjekt3():
    obj = c4d.BaseObject(c4d.Onull)
    obj.SetName('Kanten')
    obj.Message(c4d.MSG_UPDATE)
    return obj

#************************************************************************
def CreateNullobjekt4():
    obj = c4d.BaseObject(c4d.Onull)
    obj.SetName('Polygone')
    obj.Message(c4d.MSG_UPDATE)
    return obj

#************************************************************************
def CreateDreiecke():
    obj = c4d.BaseObject(c4d.Opolygon)
    obj.ResizeObject(NP, N3)
    obj.SetName('Dreiecke')
  # Punkte uebergeben  
    for i in xrange(NP):
        obj.SetPoint(i, c4d.Vector(Punkte[i][0]*Faktor*k, Punkte[i][1]*Faktor*k, Punkte[i][2]*Faktor*k))
  # Dreiecke setzen
    for i in xrange(N3):
        obj.SetPolygon( i, c4d.CPolygon(Dreiecke[i][0], Dreiecke[i][1], Dreiecke[i][2], Dreiecke[i][3]))       

    obj.Message(c4d.MSG_UPDATE)
    return obj

#************************************************************************
def CreateSplineDreiecke():
    obj = c4d.BaseObject(c4d.Ospline)
    obj.SetName("Spline-Dreiecke")
    obj.ResizeObject(N3*3)
    zz = 0
    for i in xrange(N3):
        for j in xrange(3):
            x = Punkte[Dreiecke[i][j]][0]
            y = Punkte[Dreiecke[i][j]][1]
            z = Punkte[Dreiecke[i][j]][2]
            obj.SetPoint(zz, c4d.Vector(x*Faktor*k, y*Faktor*k, z*Faktor*k))
            zz = zz + 1
  # Segmente erzeugen
    obj.MakeVariableTag(c4d.Tsegment, N3)
    for i in range(0, N3):
        obj.SetSegment(i, 3, True)
  # Spline schliessen
    obj[c4d.SPLINEOBJECT_CLOSED] = True

    obj.Message(c4d.MSG_UPDATE)
    return obj

#************************************************************************
def main():
    nullobj1 = CreateNullobjekt1()
    nullobj2 = CreateNullobjekt2()   # Ecken
    nullobj3 = CreateNullobjekt3()   # Kanten
    nullobj4 = CreateNullobjekt4()   # Polygone    
        
    plyobj1 = CreateDreiecke()   
    splobj1 = CreateSplineDreiecke()   
 
    doc.InsertObject(nullobj1, None, None, True)  
    doc.InsertObject(nullobj4, nullobj1, None, True)
    doc.InsertObject(nullobj3, nullobj1, None, True)
    doc.InsertObject(nullobj2, nullobj1, None, True) 

    doc.InsertObject(splobj1, nullobj3, None, True)    
    
    doc.InsertObject(plyobj1, nullobj4, None, True)   

 
  # Kugeln auf Eckpunkte setzen
    for i in range(0, NP):
        obj = c4d.BaseObject(c4d.Osphere)
        obj[c4d.PRIM_SPHERE_RAD] = Radius
        obj.SetName(str(i))            
        x = Punkte[i][0]
        y = Punkte[i][1]
        z = Punkte[i][2]
        obj.SetAbsPos(c4d.Vector(x*Faktor*k, y*Faktor*k, z*Faktor*k))
        doc.InsertObject(obj, nullobj2, None, True)
                 
    c4d.EventAdd()

if __name__=='__main__':
    main()