"""
Bouali Attraktor Typ III
(c) Juergen Meier (www.3d-meier.de)
26.12.2021
"""

import c4d
import math

def CreateSplineObject():
  # Variablen und Konstanten
    Name = "Bouali Attraktor Typ III"
    a = 3.0            # Konstante Bouali Attraktor Typ III
    b = 2.2            # Konstante Bouali Attraktor Typ III
    c = 1.0            # Konstante Bouali Attraktor Typ III
    d = 1.51           # Konstante Bouali Attraktor Typ III
    x0 = 1.0           # Startwerte
    y0 = 1.0           # Startwerte
    z0 = 0.0           # Startwerte
    Delta = 0.0001     # Delta Wert fuer die Berechnung
    N = 3600000        # Anzahl Punkte
    N0 = 1             # Anzahl Punke am Anfang verwerfen
    Abstand = 0.01     # Abstand zwischen zwei Punkten
    FaktorX = 10       # Skalierungsfaktor x
    FaktorY = 10       # Skalierungsfaktor y
    FaktorZ = 10       # 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 = a*x*(1 - y) - b*z
       dy = -c*y*(1 - x*x)
       dz = d*x
       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()