"""
Triakisikosaeder
16.01.2021
www.3d-meier.de
http://dmccooey.com/polyhedra/TriakisIcosahedron.html
"""

import c4d
import math

# Variablen und Konstanten
Titel = 'Triakisikosaeder' # Name

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

Radius = 2.0 # Radius der Eckpunktkugeln    

C0 = 5 * (7 + math.sqrt(5)) / 44
C1 = 5 * (3 + 2 * math.sqrt(5)) / 22
C2 = (5 + math.sqrt(5)) / 4
C3 = 5 * (13 + 5 * math.sqrt(5)) / 44
C4 = (5 + 3 * math.sqrt(5)) / 4

Punkte = [( C2, 0.0,  C4),
          ( C2, 0.0, -C4),
          (-C2, 0.0,  C4),
          (-C2, 0.0, -C4),
          ( C4,  C2, 0.0),
          ( C4, -C2, 0.0),
          (-C4,  C2, 0.0),
          (-C4, -C2, 0.0),
          (0.0,  C4,  C2),
          (0.0,  C4, -C2),
          (0.0, -C4,  C2),
          (0.0, -C4, -C2),
          (0.0,  C0,  C3),
          (0.0,  C0, -C3),
          (0.0, -C0,  C3),
          (0.0, -C0, -C3),
          ( C3, 0.0,  C0),
          ( C3, 0.0, -C0),
          (-C3, 0.0,  C0),
          (-C3, 0.0, -C0),
          ( C0,  C3, 0.0),
          ( C0, -C3, 0.0),
          (-C0,  C3, 0.0),
          (-C0, -C3, 0.0),
          ( C1,  C1,  C1),
          ( C1,  C1, -C1),
          ( C1, -C1,  C1),
          ( C1, -C1, -C1),
          (-C1,  C1,  C1),
          (-C1,  C1, -C1),
          (-C1, -C1,  C1),
          (-C1, -C1, -C1)]


Dreiecke = [( 12,  0,  8,  8 ),
            ( 12,  8,  2,  2 ),
            ( 12,  2,  0,  0 ),
            ( 13,  1,  3,  3 ),
            ( 13,  3,  9,  9 ),
            ( 13,  9,  1,  1 ),
            ( 14,  0,  2,  2 ),
            ( 14,  2, 10, 10 ),
            ( 14, 10,  0,  0 ),
            ( 15,  1, 11, 11 ),
            ( 15, 11,  3,  3 ),
            ( 15,  3,  1,  1 ),
            ( 16,  0,  5,  5 ),
            ( 16,  5,  4,  4 ),
            ( 16,  4,  0,  0 ),
            ( 17,  1,  4,  4 ),
            ( 17,  4,  5,  5 ),
            ( 17,  5,  1,  1 ),
            ( 18,  2,  6,  6 ),
            ( 18,  6,  7,  7 ),
            ( 18,  7,  2,  2 ),
            ( 19,  3,  7,  7 ),
            ( 19,  7,  6,  6 ),
            ( 19,  6,  3,  3 ),
            ( 20,  4,  9,  9 ),
            ( 20,  9,  8,  8 ),
            ( 20,  8,  4,  4 ),
            ( 21,  5, 10, 10 ),
            ( 21, 10, 11, 11 ),
            ( 21, 11,  5,  5 ),
            ( 22,  6,  8,  8 ),
            ( 22,  8,  9,  9 ),
            ( 22,  9,  6,  6 ),
            ( 23,  7, 11, 11 ),
            ( 23, 11, 10, 10 ),
            ( 23, 10,  7,  7 ),
            ( 24,  0,  4,  4 ),
            ( 24,  4,  8,  8 ),
            ( 24,  8,  0,  0 ),
            ( 25,  1,  9,  9 ),
            ( 25,  9,  4,  4 ),
            ( 25,  4,  1,  1 ),
            ( 26,  0, 10, 10 ),
            ( 26, 10,  5,  5 ),
            ( 26,  5,  0,  0 ),
            ( 27,  1,  5,  5 ),
            ( 27,  5, 11, 11 ),
            ( 27, 11,  1,  1 ),
            ( 28,  2,  8,  8 ),
            ( 28,  8,  6,  6 ),
            ( 28,  6,  2,  2 ),
            ( 29,  3,  6,  6 ),
            ( 29,  6,  9,  9 ),
            ( 29,  9,  3,  3 ),
            ( 30,  2,  7,  7 ),
            ( 30,  7, 10, 10 ),
            ( 30, 10,  2,  2 ),
            ( 31,  3, 11, 11 ),
            ( 31, 11,  7,  7 ),
            ( 31,  7,  3,  3 )]
                        

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()