/* gcc  -o bs.x baksneppen.c -L/usr/X11R6/lib -lplot -lXaw -lXmu -lXt -lSM  -lICE -lXext -lX11 -lm */
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<plot.h>
#include<sys/types.h>
#include<unistd.h>
#include<termios.h>
#include<term.h>
#include<ncurses.h>
#define MAX 500
long int rng1(long int nrand);

long int ia=16807;
long int RMAX=2147483647;
long int iq=127773;
long int ir=2836;
long int irandom;


main(int argc, char *argv[]){
	long int vmin,v[MAX];
	char sbuff[7];
	int N,i,j,n,nmin,m,pl_handle;
	long seed=1003;
    	int ch=0;
    	long int vmax=0;
    	int time=0;
    	int avatime=0;
    	float x,y,dx;
    	if(argc != 3){
		printf("./bs.x  seed N(<500)\n");
        	exit(0);
    	}
    	i=atoi(argv[1]);
	rng1(i);
    	N=atoi(argv[2]);     /* number of particles */
    	if(N >= 500){
		printf("need N<500 \n");
		exit(0);
    	}
    	/* randomly initialize the fitnesses*/
    	vmin=RMAX;
    	for(n=0;n<N;n++){
		v[n]=rng1(0);
		if(v[n] < vmin){ 
			vmin=v[n];
			nmin=n;
		}
    	}

  	/* set up a plotter window            */
  	/* an animation with double buffering */
  	pl_parampl ("BITMAPSIZE", "700x700");
  	pl_parampl ("USE_DOUBLE_BUFFERING", "yes"); 
  	pl_handle = pl_newpl ("X", stdin, stdout, stderr);
  	pl_selectpl (pl_handle);
  	pl_openpl ();
  	pl_fspace (-52.0, -60.0, 52.0, 52.0);
  	pl_flinewidth (0.01);
  	dx=45.0/N; /* Set symbolsize */
	
  while(time<10000000) {  /* the time loop */
	if(vmin>vmax){
		vmax=vmin;
		avatime=0;
	}
   	++time;
   	++avatime;
    for(i=0;i<00000;i++); /* waiting loop */
      	pl_erase();
      	pl_fmove(-50.0,-58.0);
      	pl_label("TIME= ");
      	pl_fmove(-40.0,-58.0);
      	pl_label(gcvt(time,6,sbuff));
      	pl_fmove(-20.0,-58.0);
      	pl_label("AVA-TIME= ");
      	pl_fmove(0.0,-58.0);
      	pl_label(gcvt(avatime,6,sbuff));
      	pl_fmove(15.0,-58.0);
      	pl_label("Thresh= ");
      	pl_fmove(30.0,-58.0);
	y=vmax*4.6566128e-10;
      	pl_label(gcvt(y,6,sbuff));
      	pl_filltype(0);
      	pl_fbox(-51.0,-51.0,51.0,51.0); 
	y=(vmax*4.6566128e-8)-52;
	pl_fline(-52,y,52,y);
      	x=-51+dx;
      	for(n=0;n<N;n++){
	  	if(((N+n-1)%N) == nmin)
	  		pl_fillcolor(65535,0,65535);
	  	else if(((N+n+1)%N) == nmin)
	  		pl_fillcolor(65535,0,65535);
	  	else if(n == nmin)
	  		pl_fillcolor(65535,0,0);
	  	else
	  		pl_fillcolor(0,0,65535);
      	  	pl_filltype(1);
		y= (v[n]*4.6566128e-8)-52;
	  	pl_fcircle(x,y,dx);
	  	x=x+2*dx;
      	}
    	for(n=nmin-1;n<nmin+2;n++){
		v[(N+n)%N]=rng1(0);
    	}
    	vmin=RMAX;
    	for(n=0;n<N;n++){
		if(v[n] < vmin){ 
			vmin=v[n];
			nmin=n;
		}
    	}
	
  }
}



long int rng1(long int nrand){
/* This one does the multiplicative LCG of Lewis, Goodman and Miller,
    using the Schrage modulus trick 
*/
long int l;
if(nrand>0) irandom=nrand;
l = irandom/iq;
irandom = ia*(irandom-iq*l) - ir*l;
if (irandom<0) irandom=irandom+RMAX;
return(irandom);
}