head	1.1;
access;
symbols;
locks; strict;
comment	@# @;


1.1
date	97.12.01.03.59.30;	author hans;	state Exp;
branches;
next	;


desc
@animation of chaotic behavior
the number crunching thread
@


1.1
log
@Initial revision
@
text
@// the number crunching class - separate thread that calculates
// position and velocity of pendulum and returns it

// last changed 05/08/97 
// status: constructor takes too many arguments


import java.util.*;
import java.lang.*;

class animation extends Thread
{

// variables used in the static rk4 routine have to be static   
   static double y1[]=new double[2];
   static double y2[]=new double[2];
   
   static double damp, ext, phi, w, x;
   
// I need this to reset x0, v0 to initial values     
   double oldy1[]=new double[2];
   double oldy2[]=new double[2];

   double duration, xstart=0;
   int delay;
   
   Control_Win control;
   graphic out1, out2;
   P_Graph phase1, phase2;
   
   animation(Control_Win incontrol, double iny1[], double iny2[], double inphi, 
             double indamp, 
             double inw, double inext, double induration, int indelay, 
             graphic inout1, graphic inout2, P_Graph inphase1, P_Graph inphase2)
   {
      
      control = incontrol;
      y1=iny1;
      oldy1=iny1;
      y2=iny2;
      oldy2=iny2;
      
      damp=indamp;
      delay=indelay;
      duration=induration;
      phi=inphi;
      ext=inext;
      w=inw;
      out1=inout1;
      out2=inout2;
      phase1=inphase1;
      phase2=inphase2;
   }
   
   public void cleanup()
   {
      y1=oldy1;
      y2=oldy2;
   } 
   
   public void run()
   {
      int f_data=0;
      double k;
      for(x=0; x<duration; x+=0.05)
      {     
         for (k=x; k<=x+0.05; k+=0.005)
         {
            y1=rk4(x, y1, 0.005);
            y2=rk4(x, y2, 0.005);
         }

         out1.draw(y1[0], y1[1]);
         out2.draw(y2[0], y2[1]);
         
         phase1.draw(y1[0], y1[1]);
         phase2.draw(y2[0], y2[1]);
        
         try
         { 
            if (control.kill_thread == true)
            {
               System.out.println("kill thread");
               this.stop();
               
            }
            else
            {
               this.sleep(delay);
            }
         }catch (InterruptedException e)
         {}
      }

      this.cleanup();
   }    

// the Runge-Kutta routine
   
   static double[] rk4(double x, double y[], double step)
   {
       double h=step/2;
       double temp[]= new double[2];
       double k1[]= new double[2];
       double k2[]= new double[2];
       double k3[]= new double[2];
       double k4[]= new double[2];
       int i;         
             
       for (i=0;i<2;i++)   k1[i]=step*f(x, y, i);
       for (i=0;i<2;i++)   temp[i]=y[i]+0.5*k1[i];
       for (i=0;i<2;i++)  k2[i]=step*f(x+h, temp, i);
       for (i=0;i<2;i++)  temp[i]=y[i]+0.5*k2[i];
       for (i=0;i<2;i++)  k3[i]=step*f(x+h, temp, i);
       for (i=0;i<2;i++)  temp[i]=y[i]+k3[i];
       for (i=0;i<2;i++)  k4[i]=step*f(x+step, temp, i);
       for (i=0;i<2;i++) y[i]+=(k1[i]+2*k2[i]+2*k3[i]+k4[i])/6.0;
       return y;
   }

// the routine which returns the derivatives/definition of the system	
      
   static double f(double x, double y[], int i)
   {
      double tmp=0;
      if (i==0) tmp=(y[1]);
      if (i==1) tmp=(-Math.sin(y[0])-damp*y[1]+ext*Math.cos(w*x+phi));
      
      return tmp;
   }
}       
@