"""
Kachelung 46
12.10.2015
"""

import c4d
import math

# Variablen und Konstanten
N = 5      # Anzahl Punkte x
r = 1.0    # Umkreisradius Sechseck

ris = 0.5 * math.sqrt(3.0) * r             # Inkreisradius Sechseck
a = 60                                     # Winkel Raute
hr = r*math.sin(a*math.pi/180.0)           # Hoehe Raute
b = r*math.cos(a*math.pi/180.0)
c = r*math.sin((60-a)*math.pi/180.0)
d = r*math.cos((60-a)*math.pi/180.0)

Abstand_x = 0.5*r + d + r
Abstand_y = 2*ris + hr
Versatz_x = b
Versatz_y = ris - c

NS = 5       # Anzahl Sechsecke von der Mitte in einer Richtung
if NS%2 ==0: # nur ungerade NS zulassen
    NS = NS +1

Faktor = 1000        # Skalierungsfaktor

Sexeck = []
Dreieck = []
Raute = []

Inversion = True
#Inversion = False

def Invers(z):
    return 1/z

def CreateNullobjekt():
    obj = c4d.BaseObject(c4d.Onull)   # Nullobjekt erzeugen
    obj.SetName("Kachelung 46")
    obj.Message(c4d.MSG_UPDATE)
    return obj

def Drehung(z, Grad):
    Winkel = Grad * math.pi/180
    x = z.real
    y = z.imag
    x1 = math.cos(Winkel) * x + math.sin(Winkel) * y
    y1 = -math.sin(Winkel) * x + math.cos(Winkel) * y
    return complex(x1, y1)


# Rauten  ***************************************************************************************
def CreateSplineObjectR():
    obj = c4d.BaseObject(c4d.Ospline)   # Splineobjekt erzeugen
    obj.ResizeObject(N*4*(2*NS+1)*(2*NS+1)*3)
    dt= r/N
    zz = 0   # Zaehler Punkte

    # Rautendaten erzeugen
    for i in range(0,N):
        z = complex(i * dt, 0)
        Raute.append(z)
    for i in range(0,N):
        z = complex(i * dt, 0)
        z = Drehung(z, -a) + complex(r, 0)
        Raute.append(z)
    for i in range(0,N):
        z = complex(i * dt, 0)
        z = Drehung(z, -180) + complex(r + b, hr)
        Raute.append(z)
    for i in range(0,N):
        z = complex(i * dt, 0)
        z = Drehung(z, 180 - a) + complex(b, hr)
        Raute.append(z)
    # Raute verschieben
    for i in range(0,4*N):
        Raute[i] = Raute[i] + + complex(-0.5*r, ris)
    # Rauten duplizieren und drehen
    for k in range(1,3):
        for i in range(0,4*N):
            z = Drehung(Raute[i], -60*k)
            Raute.append(z)

    for k in range(0,2*NS+1):
        x1 = -Versatz_x*NS + k*Versatz_x
        y1 = -Abstand_y*NS + k*Abstand_y
        for j in range(0,2*NS+1):
            x0 = -Abstand_x*NS + j*Abstand_x
            y0 = -Versatz_y*NS + j*Versatz_y
            for i in range(0,4*N*3):
                z = Raute[i] + complex(x0, y0) + complex(x1, y1)
                if Inversion:
                    z = Invers(z)
                obj.SetPoint(zz, c4d.Vector(z.real*Faktor, z.imag*Faktor, 0.0))
                zz = zz + 1

    obj.MakeVariableTag(c4d.Tsegment, (2*NS+1)*(2*NS+1)*3)
    for i in range(0, (2*NS+1)*(2*NS+1)*3):
        obj.SetSegment(i, 4*N, True)

    obj[c4d.SPLINEOBJECT_CLOSED] = True
    obj.SetName("Rauten")

    obj.Message(c4d.MSG_UPDATE)
    return obj


# Sechsecke *************************************************************
def CreateSplineObjectS():
    # Splineobjekte erzeugen
    obj = c4d.BaseObject(c4d.Ospline)
    obj.ResizeObject(N*6*(2*NS+1)*(2*NS+1))
    dt= r/N  # Unterteilung Kante
    zz = 0   # Zaehler Punkte

    # Sechseckdaten erzeugen
    for i in range(0,N):
        z = complex(-0.5 * r + i * dt, -ris)
        Sexeck.append(z)
    for k in range(1,6):
        for j in range(0,N):
            z = Drehung(Sexeck[j], -60*k)
            Sexeck.append(z)

    for k in range(0,2*NS+1):
        x1 = -Versatz_x*NS + k*Versatz_x
        y1 = -Abstand_y*NS + k*Abstand_y
        for j in range(0,2*NS+1):
            x0 = -Abstand_x*NS + j*Abstand_x
            y0 = -Versatz_y*NS + j*Versatz_y
            for i in range(0,6*N):
                z = Sexeck[i] + complex(x0, y0) + complex(x1, y1)
                if Inversion:
                    z = Invers(z)
                obj.SetPoint(zz, c4d.Vector(z.real*Faktor, z.imag*Faktor, 0.0))
                zz = zz + 1

    obj.MakeVariableTag(c4d.Tsegment, (2*NS+1)*(2*NS+1))
    for i in range(0, (2*NS+1)*(2*NS+1)):
        obj.SetSegment(i, 6*N, True)

    obj[c4d.SPLINEOBJECT_CLOSED] = True
    obj.SetName("Sechsecke")

    obj.Message(c4d.MSG_UPDATE)
    return obj


# Dreiecke ***************************************************************
def CreateSplineObjectD():
    obj = c4d.BaseObject(c4d.Ospline)   # Splineobjekt erzeugen
    obj.ResizeObject(N*3*(2*NS+1)*(2*NS+1)*2)
    dt= r/N
    zz = 0   # Zaehler Punkte

    # Dreiecksdaten erzeugen
    for i in range(0,N):
        z = complex(i * dt,)
        Dreieck.append(z)
    for i in range(0,N):
        z = complex(i * dt, 0)
        z = Drehung(z, -120) + complex(r,  0)
        Dreieck.append(z)
    for i in range(0,N):
        z = complex(i * dt, 0)
        z = Drehung(z, 120) + complex(0.5*r, ris)
        Dreieck.append(z)
    # Dreieck drehen und verschieben
    for i in range(0,3*N):
        Dreieck[i] = Drehung(Dreieck[i], 60 - a) + complex(0.5*r, ris)
    # Dreieck duplizieren und drehen
    for i in range(0,3*N):
        z = Drehung(Dreieck[i], 60)
        Dreieck.append(z)

    for k in range(0,2*NS+1):
        x1 = -Versatz_x*NS + k*Versatz_x
        y1 = -Abstand_y*NS + k*Abstand_y
        for j in range(0,2*NS+1):
            x0 = -Abstand_x*NS + j*Abstand_x
            y0 = -Versatz_y*NS + j*Versatz_y
            for i in range(0,3*N*2):
                z = Dreieck[i] + complex(x0, y0) + complex(x1, y1)
                if Inversion:
                    z = Invers(z)
                obj.SetPoint(zz, c4d.Vector(z.real*Faktor, z.imag*Faktor, 0.0))
                zz = zz + 1

    obj.MakeVariableTag(c4d.Tsegment, (2*NS+1)*(2*NS+1)*2)
    for i in range(0, (2*NS+1)*(2*NS+1)*2):
        obj.SetSegment(i, 3*N, True)

    obj[c4d.SPLINEOBJECT_CLOSED] = True
    obj.SetName("Dreiecke")

    obj.Message(c4d.MSG_UPDATE)
    return obj


def main():
    nullobj = CreateNullobjekt()
    slinobj1 = CreateSplineObjectS()
    slinobj2 = CreateSplineObjectD()
    slinobj3 = CreateSplineObjectR()

    doc.InsertObject(nullobj, None, None, True)
    doc.InsertObject(slinobj1, nullobj, None, True)
    doc.InsertObject(slinobj2, nullobj, None, True)
    doc.InsertObject(slinobj3, nullobj, None, True)

    c4d.EventAdd()

if __name__=='__main__':
    main()