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

import c4d
import math

# Variablen und Konstanten
Titel = 'Pentagonikositetraeder' # Name

NP = 38    # Anzahl Punkte
N5 = 24    # Anzahl Fuenfecke

Radius = 1 # Radius der Eckpunktkugeln       


C0 = 0.218796643000480441021525033712
C1 = 0.740183741369857222808508165943
C2 = 1.02365617811269018236347687527
C3 = 1.36141015192644253450100912043

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


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

                        

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 CreateFuenfecke():
    obj = c4d.BaseObject(c4d.Opolygon)
    obj.ResizeObject(NP, N5*2)
    obj.SetName('Fünfecke')
  # 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))
    zz = 0 # Zaehler fuer Polygone zuruecksetzen
    for i in xrange(N5):
        obj.SetPolygon(zz, c4d.CPolygon(Fuenfecke[i][0], Fuenfecke[i][1], Fuenfecke[i][2], Fuenfecke[i][4]))
        zz = zz + 1
        obj.SetPolygon(zz, c4d.CPolygon(Fuenfecke[i][2], Fuenfecke[i][3], Fuenfecke[i][4], Fuenfecke[i][4]))        
        zz = zz + 1
   
    obj.Message(c4d.MSG_UPDATE)
    return obj


#************************************************************************
def CreateSplineFuenfecke():
    obj = c4d.BaseObject(c4d.Ospline)
    obj.SetName("Spline-Fuenfecke")
    obj.ResizeObject(N5*5)
    zz = 0
    for i in xrange(N5):
        for j in xrange(5):
            x = Punkte[Fuenfecke[i][j]][0]
            y = Punkte[Fuenfecke[i][j]][1]
            z = Punkte[Fuenfecke[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, N5)
    for i in range(0, N5):
        obj.SetSegment(i, 5, 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 = CreateFuenfecke()   
    splobj1 = CreateSplineFuenfecke() 
 
    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()