"""
Rikitake Attraktor
(c) Juergen Meier (www.3d-meier.de)
25.12.2021
"""

import c4d
import math

def CreateSplineObject():
  # Variablen und Konstanten
    Name = "Rikitake Attractor"
    a = 5              # Konstante Rikitake Attraktor
    b = 2              # Konstante Rikitake Attraktor
    x0 = 1             # Startwerte
    y0 = 1             # Startwerte
    z0 = 1             # Startwerte
    Delta = 0.001      # Delta Wert fuer die Berechnung
    N = 50000           # Anzahl Punkte
    Faktor = 100       # Skalierungsfaktor
    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)


    for i in range(0,N):
       h=math.sqrt(x*x+y*y)
       dx=-b*x+z*y
       dy=-b*y+(z-a)*x
       dz=1-x*y
       x=x + Delta * dx
       y=y + Delta * dy
       z=z + Delta * dz
       obj.SetPoint(i, c4d.Vector(x*Faktor, y*Faktor, z*Faktor))

    obj.Message(c4d.MSG_UPDATE)
    return obj

def main():
    slinobj = CreateSplineObject()

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

    c4d.EventAdd()

if __name__=='__main__':
    main()