/********************************************************************
Multi Chua II Attraktor
Plugin für Cinema Version 8.0
(c) Jürgen Meier (Juergen-Meier@Hamburg.de)
Version : 1.01
02.09.2003
*********************************************************************/

// Globale Variablen und Konstanten
   // Plugin ID Nummer
      const var cPluginID = 1012251;
   // Pluginname
      const var cPlugName = "Multi Chua II Attraktor";
   // Pluginhilfstext
      const var cPlugHelp = "Erzeugt den Multi Chua II Attraktor";
   // Objektname
      const var cObjName = "Multi Chua II Attraktor";
   // Index
      var i,j,k;
   // Koordinaten der Punkte
      var x,y,z;
   // Differenzen
      var dx,dy,dz;
   // Konstanten
      var alpha=9.0;
      var beta=14.286;
      var a=1.3;
      var b=0.11;
   // Hilfsvariable
      var h;
   // Delta
      var Delta=0.001;
   // Anzahl Iterationen
      var N=500000;
   // Startwerte
      var x0=0.1;
      var y0=-0.2;
      var z0=0.3;
   // Alter Wert
      var xa,ya,za;
   // Abstand zwischen 2 Punkten
      var PP;
      var Abstand=0.2;
   // Skalierungsfaktor
      var FaktorX=10;
      var FaktorY=50;
      var FaktorZ=10;
   // Breite des Menütextes
      var Breite=150;
   // Maximalwerte
      var cMaxFloat=1000000;
      var cMaxInt=100000000;
   // Objektvariable
      var qObj;
   // Dialogvariable
      var gDial;
   // Array für die Konstanten
      var m;
      var c;




// ******************************************************************
// Funktion zum Objekt aktivieren (c) h.g.seib
// ******************************************************************
fSetActiveObject(vDoc,vObj)
{
   // Aktives Objekt ermitteln
      var vAObj = GetActiveObject(vDoc);
   // Active object Bit bzw. Flag setzen
      if (vAObj) vAObj->DelBit(BIT_AOBJ);
          vObj->SetBit(BIT_AOBJ);
   // Änderung im Dokument bekanngeben (neues Objekt selektiert)
      GeEventAdd(NEW_ACTIVE_OBJECT);
}


// ******************************************************************
// Funktion f(x)
// ******************************************************************
   f(x)
   {
    var Summe;
    Summe=0;
    for (k=1;k<6;k++)
        {
        Summe=Summe+((m[k-1] - m[k])*(abs(x + c[k])-abs(x - c[k])));
        }
    h= m[5]*x+0.5*Summe;
    return h;
   }


// ******************************************************************
// Funktion zur Berechnung des Attraktors
   fAttraktor(vDoc)
   {
   // Variablen deklarieren
      // Felder für die Punkte des Attraktors
         var vPointAry;
      // Variable für VariableChanged und BackupTag deklarieren
         var vVarChanged,vBackupTags;
   // Neues Splineobjekt erzeugen und Namen zuweisen
      qObj=new(SplineObject);
      if (!qObj) return FALSE;
      qObj->SetName(cObjName);
   // Undo Funktion aktivieren
      vDoc->StartUndo();
      vDoc->AddUndo(UNDO_OBJECT_NEW,qObj);
      vDoc->EndUndo();
   // Variable für VariableChanged und BackupTag erzeugen
      vVarChanged = new(VariableChanged);
      vBackupTags = new(BackupTags);
   // Feld für die Punkte erzeugen
      vPointAry = new(array,N);
   // Feld für die Punkte dem Splineobjekt zuordnen
      qObj->SetPoints(vPointAry);
   // ???
      vBackupTags->Init(qObj);
      vVarChanged->Init(0,N);
   // Message der Klasse BaseList2D
      if (!qObj->Message(MSG_POINTS_CHANGED, vVarChanged))
         {
             vBackupTags->Restore();
             return FALSE;
         }
   // Objekt in der Dokumenthierarchie positionieren
      vDoc->InsertObject(qObj,NULL,NULL);
      fSetActiveObject(vDoc,qObj);
      GeEventAdd(DOCUMENT_CHANGED);
   // Variable für den Container (Klasse BaseContainer)
      var vContainer = qObj->GetContainer();
   // Spline Typ festlegen
      vContainer->SetData(SPLINEOBJECT_TYPE,3);
   // Containerdaten dem Spline Objekt zurückgeben
      qObj->SetContainer(vContainer);
   // Punkte des Attraktors berechnen
      // Startwerte übergeben
         x=x0;
         y=y0;
         z=z0;
         xa=x0;
         ya=y0;
         za=z0;
      // Schleife durchlaufen
         j=0;
         for (i=0;i<N;i++)
             {
              // Delta Werte berechnen
                 dx= alpha * (y - f(x));
                 dy= x - y + z;
                 dz= -beta * y;
              // Neue Koordinaten berechnen
                 x=x + Delta * dx;
                 y=y + Delta * dy;
                 z=z + Delta * dz;
              // Abstand zum letzten Punkt berechnen
                 PP=sqrt((x-xa)*(x-xa)+(y-ya)*(y-ya)+(z-za)*(z-za));
              // Punkte setzen
                 if (PP > Abstand) {
                     qObj->SetPoint(j,vector(x*FaktorX,z*FaktorZ,y*FaktorY));
                     j=j+1;
                     xa=x;
                     ya=y;
                     za=z;
                 }
             }
     return TRUE;
   }

// ******************************************************************
// Ableitung der Klasse MyDialog von GeDialog
// ******************************************************************
class MyDialog : GeDialog
{
   public:
      MyDialog();
      CreateLayout();
      Init();
      CoreMessage(vID,vMSG);
      Command(vID,vMSG);
}
// ******************************************************************
// Deklaration des Konstruktors
   MyDialog::MyDialog()
   {
   // Konstruktor der parent class aufrufen
      super(cPluginID);
   }

// ******************************************************************
// Deklaration der Methode CreateLayout
   MyDialog::CreateLayout()
   {     
   // Array für die Konstante m
      m = new(array,6);
      m[0]=0.9/7;
      m[1]=-3.0/7;
      m[2]=3.5/7;
      m[3]=-2.7/7;
      m[4]=4.0/7;
      m[5]=-2.4/7;
   // Array für die Konstante c
      c = new(array,6);
      c[0]=0.0;
      c[1]=1.0;
      c[2]=2.15;
      c[3]=3.6;
      c[4]=6.2;
      c[5]=9.0;
   // Plugintitel setzen
      SetTitle(cPlugName);
   // Gruppe für die Konstanten
      AddGroupBeginV(1000,BFH_CENTER,2,"",0);
      {
         AddStaticText(0,BFH_LEFT,Breite,0," Konstante alpha",0);
         AddEditNumber(2014,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," Konstante beta",0);
         AddEditNumber(2015,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," Konstante a",0);
         AddEditNumber(2000,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," Konstante b",0);
         AddEditNumber(2001,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," Konstante c",0);
         AddEditNumber(2002,BFH_LEFT,100,0);
      }
      // Vertikal orientierte Gruppe schliessen
         AddGroupEnd();
   // Waagerechte Trennlinie
      AddSeparatorH(0);
   // Gruppe für die Startwerte
      AddGroupBeginV(1001,BFH_CENTER,2,"",0);
      {
         AddStaticText(0,BFH_LEFT,Breite,0," Start x",0);
         AddEditNumber(2005,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," Start y",0);
         AddEditNumber(2006,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," Start z",0);
         AddEditNumber(2007,BFH_LEFT,100,0);
      }
      // Vertikal orientierte Gruppe schliessen
         AddGroupEnd();
   // Waagerechte Trennlinie
      AddSeparatorH(0);
   // Gruppe für die Parameter der Iteration
      AddGroupBeginV(1002,BFH_CENTER,2,"",0);
      {
         AddStaticText(0,BFH_LEFT,Breite,0," Delta",0);
         AddEditNumber(2008,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," Anzahl N",0);
         AddEditNumber(2009,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," Abstand",0);
         AddEditNumber(2011,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," FaktorX",0);
         AddEditNumber(2010,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," FaktorY",0);
         AddEditNumber(2012,BFH_LEFT,100,0);
         AddStaticText(0,BFH_LEFT,Breite,0," FaktorZ",0);
         AddEditNumber(2013,BFH_LEFT,100,0);
      }
      // Vertikal orientierte Gruppe schliessen
         AddGroupEnd();
   // Waagerechte Trennlinie
      AddSeparatorH(0);
   // Button für den Start der Berechnung
      AddButton (3000,BFH_CENTER,200,0,"Start");
      return TRUE;
   }

// ******************************************************************
// Deklaration der Methode Init
   MyDialog::Init()
   {
   // Startwerte setzen
      SetFloat(2000,a,-cMaxFloat,cMaxFloat,0);
      SetFloat(2001,b,-cMaxFloat,cMaxFloat,0);
      SetFloat(2005,x0,-cMaxFloat,cMaxFloat,0);
      SetFloat(2006,y0,-cMaxFloat,cMaxFloat,0);
      SetFloat(2007,z0,-cMaxFloat,cMaxFloat,0);
      SetFloat(2008,Delta,-cMaxFloat,cMaxFloat,0);
      SetInt(2009,N,0,cMaxInt,0);
      SetFloat(2010,FaktorX,0,cMaxFloat,0);
      SetFloat(2011,Abstand,0,cMaxFloat,0);
      SetFloat(2012,FaktorY,0,cMaxFloat,0);
      SetFloat(2013,FaktorZ,0,cMaxFloat,0);
      SetFloat(2014,alpha,-cMaxFloat,cMaxFloat,0);
      SetFloat(2015,beta,-cMaxFloat,cMaxFloat,0);
      return TRUE;
   }

// ******************************************************************
// Deklaration der Methode CoreMessage
   MyDialog::CoreMessage(vID,vMSG)
   {
      switch (vID)
         {
            case NEW_DOCUMENT:
            case DOCUMENT_CHANGED:
            case NEW_ACTIVE_OBJECT:
            case ACTIVE_OBJECT_CHANGED:
            return TRUE;
            break;
         }
      return FALSE;
   }

// ******************************************************************
// Deklaration der Methode Command
   MyDialog::Command(vID,vMSG)
   {
      var vDoc;
   //
      StopAllThreads();
      switch (vID)
      {
      case 3000 : // Start Button gedrückt
                    {
                       vDoc = GetActiveDocument();
                       if (!vDoc) return FALSE;
                       fAttraktor(vDoc);
                    }
      case 2000 : a=GetFloat(2000);
      case 2001 : b=GetFloat(2001);
//      case 2002 : c=GetFloat(2002);
      case 2005 : x0=GetFloat(2005);
      case 2006 : y0=GetFloat(2006);
      case 2007 : z0=GetFloat(2007);
      case 2008 : Delta=GetFloat(2008);
      case 2009 : N=GetInt(2009);
      case 2010 : FaktorX=GetFloat(2010);
      case 2011 : Abstand=GetFloat(2011);
      case 2012 : FaktorY=GetFloat(2012);
      case 2013 : FaktorZ=GetFloat(2013);
      case 2014 : alpha=GetFloat(2014);
      case 2015 : beta=GetFloat(2015);
      }
      return TRUE;
   }


// ******************************************************************
// Ableitung der Klasse MyMenuPlugin von MenuPlugin
// ******************************************************************
class MyMenuPlugin : MenuPlugin
{     public:
          MyMenuPlugin();
          GetID();
          GetName();
          GetHelp();
          Execute(doc);
          RestoreLayout(secret);
}
// ******************************************************************
// Deklaration des Konstruktors
   MyMenuPlugin::MyMenuPlugin()
   { 
      super();
   }

// ******************************************************************
// Deklaration der Methode GetID
   MyMenuPlugin::GetID()
   {
      return cPluginID;
   }

// ******************************************************************
// Deklaration der Methode GetName
   MyMenuPlugin::GetName()
   {
       return cPlugName;
   }

// ******************************************************************
// Deklaration der Methode GetHelp
   MyMenuPlugin::GetHelp()
   {
       return cPlugHelp;
   }

// ******************************************************************
// Deklaration der Methode Execute
   MyMenuPlugin::Execute(doc)
   {
      if (!gDial)
         gDial=new(MyDialog);
     gDial->Open(TRUE,-1,-1);
   }
// ******************************************************************
// Deklaration der Methode RestoreLayout
   MyMenuPlugin::RestoreLayout(secret)
   {
      if (!gDial)
         gDial=new(MyDialog);
      gDial->RestoreLayout(secret);
   }

// ******************************************************************
main()
{
    Register(MyMenuPlugin);
}