/*
*  Coupled Map Lattices
*  Por Rogério Neves
*  LSI- USP
*/

import java.awt.*;
import java.awt.Graphics;
import java.awt.Color;
import java.applet.Applet;
import java.util.Random;

public class cml extends Applet implements Runnable {
    Random rnd = new Random();
    Thread tread;
    int i,j,cont;
    int Xx=100,Xy=100;
    double[][] X;
    double E=.4;
    double mu=3.85;
    int step=8;
    int delay=10;
    char map='g';

    public synchronized void init() {
        String par;
        //  Get Parameters
        par = getParameter("xsize");
        if (par!=null) Xx = Integer.parseInt(par);
        par = getParameter("ysize");
        if (par!=null) Xy = Integer.parseInt(par);
        par =  getParameter("cc");
        if (par!=null) E = Double.valueOf(par).doubleValue();
        par = getParameter("mu");
        if (par!=null) mu = Double.valueOf(par).doubleValue();
        par = getParameter("step");
        if (par!=null) step = Integer.parseInt(par);
        par = getParameter("delay");
        if (par!=null) delay = Integer.parseInt(par);
        par = getParameter("palette");
        if (par!=null) map = par.toLowerCase().charAt(0);
        X = new double [Xx][Xy];
        //  Initializes X randomly
        for(i=0;i<Xx;i++) for(j=0;j<Xy;j++)
            X[i][j] = rnd.nextDouble();
        //  Insulate borders
        for(i=0;i<Xx;i++) {
            X[i][0] = 0;
            X[i][Xy-1] = 0;
        }
        for(j=0;j<Xy;j++) {
            X[0][j] = 0;
            X[Xx-1][j] = 0;
        }
    }
    public synchronized void paint(Graphics g) {
        setBackground(Color.white);
        int w = (getSize().width)/Xx;
        int h = (getSize().height)/Xy;
        for(i=0;i<Xx;i++) for(j=0;j<Xy;j++) {
            if (map=='g')
                g.setColor(Color.getHSBColor((float)1,(float)0,(float)X[i][j]));
            if (map=='h')
                g.setColor(Color.getHSBColor((float)(.65+.59*X[i][j]),(float)1,(float)(.5+.5*X[i][j])));
            if (map=='c')
                g.setColor(Color.getHSBColor((float)(X[i][j]),(float)1,(float)1));
            if (map=='b') {
                float c=0;
                if (X[i][j]>.5) c=1;
                g.setColor(Color.getHSBColor((float)1,(float)0,c));
            }
            g.fillRect(1+i*w,1+j*h,w,h);
        }
    }
    public void start() {
	setEnabled(true);
        tread = new Thread(this);
        tread.start();
    }
    public void stop() {
	setEnabled(false);
        tread = null;
    }
    public void dispose() {
    }
    public void run() {
        Thread me = Thread.currentThread();
        while (tread == me) {
            try {
                Thread.currentThread().sleep(delay);
            } catch (InterruptedException e) {   }
            cont++;
            //  Main Interactor
            for(i=0;i<Xx;i++) for(j=0;j<Xy;j++)
                X[i][j]=mu*X[i][j]*(1-X[i][j]);
            for(i=0;i<Xx;i++) for(j=0;j<Xy;j++)
                if ((i>0)&(i<Xx-1)&(j>0)&(j<Xy-1))
                    X[i][j]=(1-E)*X[i][j]+
                               E*(X[(i-1)][(j)]+
                                  X[(i+1)][(j)]+
                                  X[(i)][(j-1)]+
                                  X[(i)][(j+1)]+
                                  X[(i-1)][(j-1)]+
                                  X[(i-1)][(j+1)]+
                                  X[(i+1)][(j-1)]+
                                  X[(i+1)][(j+1)])/8;
            if (cont>step) {
                repaint();
                cont=1;
            }
        }
    }

    public void update(Graphics g) {
        paint(g);
    }
    public String getAppletInfo() {
        return "Logistic Map for\nTime-series\nBy Rogério Neves.";
    }
    public String[][] getParameterInfo() {
        String[][] info = {
            {"xsize", "integer number", "The cell matrix x size in pixels."},
            {"ysize", "integer number", "The cell matrix y size in pixels."},
            {"cc", "double between 0 and 1", "The 8-neighbour coupling constant."},
            {"mu", "double between 0 and 1", "The cell time-series parameter."},
            {"step", "integer greater than 0", "The Nth. interaction to show, 1=show every frame"},
            {"delay", "integer greater than 0", "The frame sleep time in microseconds"},
            {"palette", "string", "The color oscheme to use (binary, gray, color, hot."}
        };
        return info;
    }
}