/*
 * @(#)Graph.java	1.9 99/08/04
 * MODIFS PYL juillet Aout 2005!!!!!!!!!!
 *
 * This software is not designed or intended for use in on-line control of
 * aircraft, air traffic, aircraft navigation or aircraft communications; or in
 * the design, construction, operation or maintenance of any nuclear
 * facility. Licensee represents and warrants that it will not use or
 * redistribute the Software for such purposes.
 */
//javac guiChal.java Kul.java ;java guiChal  
//javac guiChal.java Kul.java; open index.html       


import java.util.*;
import java.awt.*;
import java.applet.Applet;
import java.awt.event.*;
import javax.swing.*;
import java.io.*;

 
 
 
class Kul extends Thread    
 {
 public double t;
 private double dt=0;
 public double dx,dy,xb=.2;
 public double T[][];
 public double To[][];
 public double Too[][];
 public double Ti[][];
 public double T0=0,T1;
 public double Tx0,Tx1,T0y,T1y;
 public boolean slopex,slopey,mixtex,mixtey,valx,valy,symx,symy;
 public int ib;
 public boolean doublet=true,bouge=false;
 public double Ts=0;
 public double q=0;
 public double Bi=0;
 public double omega=6,vitesse=0,amp=0;
 public int n; 
// public tridiag tridiagC;
 boolean pleaseWait = true;
 
private int property;
public int getProperty  (         ){ return property; }
public void setProperty (int value){ property = value; }
public void setdt (double value){ dt = value;System.out.println("------------>>>dt"+dt); }
public double getdt  (         ){/*System.out.println("-------------<<<dt"+dt);*/return dt; }

public void Kul( )
 {    
 }
 
 public void cl0000( )
 {   
 	doublet =true; 
 	bouge=false;

   }

 public void cl1000( )
 {   
 	doublet=false; 
 	bouge=false;
 }
 
public void cl0100( )
 {    
  ib=n/3;	
  //xb=.2;
  amp=1;
  
  System.out.println("   PYL !!! " + xb + " t" + t);
  bouge=true;
 } 

public void cl0010( )
 {    
 Tx0=0;
 T1y=0;
 Tx1=1;
 T0y=0;
 Ts=0;
 q=0;
 
 } 
  
 
 public void clt0( )
 {    
 Tx0=0;
 T1y=0;
 Tx1=0;
 T0y=0;
 Ts=0;
 Ts=1;
 q=0;
 
 } 
 
 
 
 public void clB1( )
 {    
 Tx0=0;
 T1y=1;
 Tx1=1;
 T0y=1;
 Ts=0;
 q=0;
 
 } 

 
  
 public void prems()
 { int i,j;
    //pleaseWait = true;
    t=0;
	n=300;//370;
	dx=1./n;
	dy=1./n;
	omega=0;
	xb=0.2;
   System.out.println(" initialize "+dx); 
   To = new double[n+1][n+1];
   Too= new double[n+1][n+1];
   T  = new double[n+1][n+1];
   Ti = new double[n+1][n+1];
   for(i=0;i<=n;i++){
      for(j=0;j<=n;j++)
		{To[i][j]=0;T[i][j]=To[i][j];Ti[i][j]=To[i][j];Too[i][j]=To[i][j];}}
	  System.out.println("toutazero");
	   for(i=1;i<n;i++){
      for(j=1;j<n;j++)
		{To[i][j]=0;T[i][j]=To[i][j];Ti[i][j]=To[i][j];Too[i][j]=To[i][j];}}
	
	  
	
	
	
 
    }
 
   public void run() {
   int i,j,it=0;
   double phase=0;
   //System.out.println("params " + dt + " " +dx +" dt/dx^2="+(dt/dx/dx));
	
	while (!interrupted()) {
	   
	    
		if(dt>0){
		phase=phase+omega*dt;
		t=t+dt;
		xb=xb+(dt*vitesse);
		ib=2+(int)(xb*n);
		// System.out.println("params " + ib + " t" + t);
	   
		if(ib>n-2){ib=2;xb=0;}
			To[ib][n/2]=amp*0.3*Math.cos(phase);
		
		for(i=1;i<=n-1;i++){
		for(j=1;j<=n-1;j++){
        T[i][j]= 2*To[i][j] -Too[i][j] +
                          dt*dt*((To[i-1][j]-2*To[i][j]+To[i+1][j])+
                                 (To[i][j-1]-2*To[i][j]+To[i][j+1])) /dx/dx;
    
		}}

					
		 
	//	tridiag.u0=Tx0+Ts*Math.sin(phase);
		 for(i=1;i<=n-1;i++){
		for(j=1;j<=n-1;j++){	 
		Too[i][j]=To[i][j];
		To[i][j]=T[i][j];}}
				}
		
		try {
		Thread.sleep(10);
	    } catch (InterruptedException e) {}
				// Do work
                // Check if should wait
                synchronized (this) {
                    while (pleaseWait) {
                        try {
                            wait();
                        } catch (Exception e) {}
                    }
                }
    }
	
}
		
 }