"""
Multi Layer Chua I Attraktor
(c) Juergen Meier (www.3d-meier.de)
24.12.2021
"""

import c4d
import math

a = 1.3             # Konstante Multi Layer Chua I Attraktor
b = 0.11            # Konstante Multi Layer Chua I Attraktor
c = 2.0             # Konstante Multi Layer Chua I Attraktor
alpha = 10.814      # Konstante Multi Layer Chua I Attraktor
beta = 14.0         # Konstante Multi Layer Chua I Attraktor

def Fx(x):
    d=0;
    if (c % 2 ==0):
        d = math.pi
    if (x >= 2*a*c):
        h=((b*math.pi)/(2*a))*(x - 2*a*c)

    if (x <= -2*a*c):
        h=((b*math.pi)/(2*a))*(x + 2*a*c)

    if (x < 2*a*c)  and (x > -2*a*c):
        h=-b*math.sin((math.pi*x)/(2*a)+d)
    return h

def CreateSplineObject():
  # Variablen und Konstanten
    Name = "Multi Layer Chua I Attractor"

    x0 = 1             # Startwerte
    y0 = 1             # Startwerte
    z0 = 1             # Startwerte
    Delta = 0.01       # Delta Wert fuer die Berechnung
    N = 30000          # Anzahl Punkte
    N0 = 900           # Anzahl Punke am Anfang verwerfen
    Abstand = 0.05     # Mindestabstand zwischen zwei Punkten
    FaktorX = 100      # Skalierungsfaktor x
    FaktorY = 500      # Skalierungsfaktor y
    FaktorZ = 100      # Skalierungsfaktor z
    x = x0             # Startwerte zuweisen
    y = y0             # Startwerte zuweisen
    z = z0             # Startwerte zuweisen

    obj = c4d.BaseObject(c4d.Ospline)       # Splineobjekt erzeugen
    obj.ResizeObject(N)                     # Groesse des Spline festlegen
    obj[c4d.SPLINEOBJECT_CLOSED] = False
    #obj[c4d.SPLINEOBJECT_TYPE] = c4d.SPLINETYPE_LINEAR
    obj[c4d.SPLINEOBJECT_TYPE] = c4d.SPLINETYPE_BSPLINE
    #obj[c4d.SPLINEOBJECT_TYPE] = c4d.SPLINETYPE_BEZIER
    obj.SetName(Name)

    xa=x0
    ya=y0
    za=z0
    zz = 0
    for i in range(0,N):
       dx= alpha * (y - Fx(x))
       dy= x - y + z
       dz= -beta * y
       x=x + Delta * dx
       y=y + Delta * dy
       z=z + Delta * dz
       PP=math.sqrt((x-xa)*(x-xa)+(y-ya)*(y-ya)+(z-za)*(z-za))
       if (i > N0):
           if (PP > Abstand):
               obj.SetPoint(zz, c4d.Vector(x*FaktorX, y*FaktorY, z*FaktorZ))
               zz = zz + 1
               xa=x
               ya=y
               za=z

    obj.ResizeObject(zz)

    obj.Message(c4d.MSG_UPDATE)
    return obj

def main():
    slinobj = CreateSplineObject()

    doc.InsertObject(slinobj, None, None, True)

    c4d.EventAdd()

if __name__=='__main__':
    main()