import java.applet.Applet;
import java.awt.*;
//import OsmosisCanvas;
//import gjt.*;

/*
 *
 * Osmosis
 *
 */
public class OsmPressure4 extends Nernst
{
        int Particles[] = {0,2,5,10,20,30};
        OsmPressureReportPanel4 pReport;
        OsmPressureControlPanel4 pControl;
        //OsmCanvas CDiff;
        Panel p;
        GraphCanvas4 c;
        int PLOT_AV=0, PLOT_REAL=0, PLT=0, SX, SY, SXY, SX2, Sum, i, j, N, ERR=0, n, is=0, is1=0, X[]= new int[10], Y[]= new int[10], xs[]=new int[600], ys[]= new int[600], zs[]=new int[600];
        float X1, Y1, a, b;
        public void init(){

                Panel pnl=new Panel();
           CDiff=new OsmPCanvas(this);
           //setLayout(new BorderLayout()); //
           //add("Center",CDiff);         //
                Panel pnl2=new Panel();
           tdbDisplay=new ThreeDBorder(pnl2,5);
           pnl2.setLayout(new BorderLayout());
           pnl2.add("Center",CDiff);
           pnl2.add("East", OutBar);

           //tdbDisplay=new ThreeDBorder(CDiff,5);
           tdbDisplay.inset();
           pReport=new OsmPressureReportPanel4(this);
           pControl=new OsmPressureControlPanel4(this);
           setLayout(new BorderLayout());
           add("Center",new ThreeDBorder(pnl));
           pnl.setLayout(new BorderLayout());
           Panel pnl3=new Panel();
           pnl3.setLayout(new BorderLayout());
           pnl3.add("North", pReport);
           pnl3.add("Center",tdbDisplay);
           pnl3.add("South", pControl);
           c = new GraphCanvas4(this);
           c.resize(300,300);
           pnl.add("Center", pnl3);
           pnl.add("East", c);
    }

        public boolean action(Event event, Object what){
                if(event.target==Start_Stop){
                        if (running){
                                CDiff.stop();
                                OutBar.stop();
                                running=false;
                                Start_Stop.setLabel("Start");
                        }else{
                                CDiff.start();
                                OutBar.start();
                                running=true;
                                Start_Stop.setLabel("Stop");
                        }
                }else if (event.target==Num_Particles){
                        int NewNum=Particles[Num_Particles.getSelectedIndex()];
                        if (running) CDiff.stop();
                        ((OsmPCanvas)CDiff).OutIon=NewNum;
                        ((OsmPCanvas)CDiff).OutWater=150-(((OsmPCanvas)CDiff).Hydration+1)*NewNum;
                        CDiff.init();
                        if (running) CDiff.start();
                }else if (event.target==Plot){
                        PLOT_AV=1;
                        c.repaint();
                }else if (event.target==Plot2){
                        PLOT_REAL=1;
                        c.repaint();
                }else if (event.target==Fit){
                        N=0;
                        for(i=0;i<6;i++){
                                n=0; Sum=0;
                                for(j=0;j<is;j++){
                                        if((Particles[i])==xs[j]){
                                                if(zs[j]==2){
                                                        Sum=Sum+ys[j];
                                                        n++;
                                                }
                                        }
                                }
                                if(n > 0){
                                        N++;
                                        X[N]=Particles[i];
                                        Y[N]=Sum/n;
                                }
                        }
                        if (N < 2){
                                ERR=1;
                        }
                        else{
                                SX=0;
                                SY=0;
                                SXY=0;
                                SX2=0;
                                for (i=1;i<=N;i++){
                                        SX=SX+X[i];
                                        SY=SY+Y[i];
                                        SXY=SXY+X[i]*Y[i];
                                        SX2=SX2+X[i]*X[i];
                                   //     SY2=SY2+Y[i]*Y[i];
                                }
                                X1=(float)SX/(float)N;
                                Y1=(float)SY/(float)N;
                                b=((float)SXY-(float)N*X1*Y1)/((float)SX2-(float)N*X1*X1);
                                a=Y1-b*X1;
                                PLT=1;
                        }        
                        c.repaint();
                }
                                
                return true;                    

        }

}


class OsmPressureReportPanel4 extends Panel{
        OsmPressure4 app;

        public OsmPressureReportPanel4(Applet applet){
                app=(OsmPressure4)applet;
                Label l_Vol_i=new Label("vol.",Label.CENTER);
                Label l_Conc_i=new Label("conc.",Label.CENTER);
                Label l_Vol_o=new Label("vol.",Label.CENTER);
                Label l_Conc_o=new Label("conc.",Label.CENTER);
                Label l_Voltage=new Label("real",Label.CENTER);
                Label l_Average=new Label("aver.",Label.CENTER);

                Panel pIn=new Panel();
                Label inside=new Label("Inside",Label.CENTER);
                app.bInside= new RBox(pIn,inside);
                pIn.setLayout(new GridLayout(2,2,2,0));
                pIn.add(l_Vol_i);
                pIn.add(l_Conc_i);
                pIn.add(app.CDiff.R_in);
                pIn.add(app.CDiff.B_in);

                Panel pOut=new Panel();
                Label outside=new Label("Outside",Label.CENTER);
                app.bOutside= new RBox(pOut,outside);
                pOut.setLayout(new GridLayout(2,2,2,0));
                pOut.add(l_Vol_o);
                pOut.add(l_Conc_o);
                pOut.add(app.CDiff.R_out);
                pOut.add(app.CDiff.B_out); 

                Panel pVoltage=new Panel();
                Label Volt=new Label("Pressure",Label.CENTER);
                app.bVoltage= new RBox(pVoltage,Volt);
                pVoltage.setLayout(new GridLayout(2,2,2,0));
                pVoltage.add(l_Voltage);
                pVoltage.add(l_Average);
                pVoltage.add(app.CDiff.Voltage);
                pVoltage.add(app.CDiff.Average);

             //   app.bInside.resize(90,5);
             //   app.bVoltage.resize(90,5);
             //   app.bOutside.resize(90,5);
                
              //  app.bVoltage= new RBox(app.CDiff.Voltage,"Pressure");
              //  app.bAverage= new RBox(app.CDiff.Average,"Average");
               setLayout(new GridLayout(1,3,0,0));

                add(new ThreeDBorder(app.bInside));
                add(new ThreeDBorder(app.bVoltage));
                add(new ThreeDBorder(app.bOutside));
                super.resize(300,5);
        }

        public Insets insets(){
                return new Insets(5,5,5,5);
        }

}

class OsmPressureControlPanel4 extends Panel{

       OsmPressure4 app;

        public OsmPressureControlPanel4(Applet applet){
              
                app=(OsmPressure4)applet;
                app.p=new Panel();
                app.Num_Particles= new Choice();
                app.Num_Particles.addItem("0 Particles");
                app.Num_Particles.addItem("2 Particles");
                app.Num_Particles.addItem("5 Particles");
                app.Num_Particles.addItem("10 Particles");
                app.Num_Particles.addItem("20 Particles");
                app.Num_Particles.addItem("30 Particles");
                app.Num_Particles.select(3);

                app.Start_Stop= new Button("Stop");
                app.Plot= new Button("Average");
                app.Plot2 = new Button("Real");
                app.Fit= new Button("CurveFit");
                GridBagLayout gbl= new GridBagLayout();
                app.p.setLayout(gbl);
                GridBagConstraints c= new GridBagConstraints();
                c.gridwidth=1;
                c.gridheight=1;
                c.fill= GridBagConstraints.BOTH;
                c.weightx=1;
                c.weighty=1;
                ThreeDBorder b1= new ThreeDBorder(app.Start_Stop);
                gbl.setConstraints(b1,c);
                app.p.add(b1);
                ThreeDBorder b2= new ThreeDBorder(app.Num_Particles);
                gbl.setConstraints(b2,c);
                app.p.add(b2);
                ThreeDBorder b3= new ThreeDBorder(app.Plot);
                gbl.setConstraints(b3,c);
                app.p.add(b3);
                ThreeDBorder b4= new ThreeDBorder(app.Plot2);
                gbl.setConstraints(b4,c);
                app.p.add(b4);
                ThreeDBorder b5= new ThreeDBorder(app.Fit);
                gbl.setConstraints(b5,c);
                app.p.add(b5);
                setLayout(new BorderLayout());
                add("North",new ThreeDBorder(app.p));
              //  app.c= new GraphCanvas(app);
              //  app.c.resize(300,200);
              //  add("South",new ThreeDBorder(app.c));
        }

        public Insets insets(){
                return new Insets(5,5,5,5);
        }

}

class GraphCanvas4 extends Canvas {

        OsmPressure4 app1;        

        static final int POSX = 35;
        static final int NEGX = 5;
        static final int POSY = 70;
        static final int NEGY = 90;
        
        int a,b,N,x1,x2,y1,y2;
        String Str0=" ";
        String Str1=" ";
        String Str2=" ";
        
        double xstep;
        double ystep;

        public GraphCanvas4 (Applet a){
                app1=(OsmPressure4)a;
        }
        
        public synchronized void paint (Graphics g) {
                Dimension dm=size();
                int x,y,i,j;
                xstep=(double)dm.width/(double)(POSX+NEGX);
                ystep=(double)dm.height/(double)(POSY+NEGY);
                int xax=(int)Math.round(ystep*POSY);
                int yax=(int)Math.round(xstep*NEGX);
                g.setColor(Color.blue);
                g.drawLine(yax,0,yax,dm.height);
                g.drawLine(0,xax,dm.width,xax);
                for(i=0;i<(NEGX+POSX+1);i++){
                        x=(int)Math.round(i*xstep);
                        g.drawLine(x,xax-2,x,xax+2);
                }
                for(i=0;i<(NEGY+POSY+1);i+=5){
                        y=(int)Math.round(xax+i*ystep);
                        g.drawLine(yax-2,y,yax+2,y);
                }         
                for(i=0;i<(POSY+1);i+=5){
                        y=(int)Math.round(xax-i*ystep);
                        g.drawLine(yax-2,y,yax+2,y);
                } 
                g.setColor(Color.black);
                Font f = new Font("TimesRoman",Font.PLAIN,11);
                g.setFont(f);
                g.drawString("0",yax+(int)Math.round(0*xstep)-2,xax+12);      
                g.drawString("2",yax+(int)Math.round(2*xstep)-2,xax+12);
                g.drawString("5",yax+(int)Math.round(5*xstep)-4,xax+12);
                g.drawString("10",yax+(int)Math.round(10*xstep)-9,xax+12);
                g.drawString("20",yax+(int)Math.round(20*xstep)-9,xax+12);
                g.drawString("30",yax+(int)Math.round(30*xstep)-9,xax+12);
                f = new Font("TimesRoman",Font.PLAIN+Font.BOLD,12);
                g.setFont(f);
                g.setColor(Color.magenta);
                g.drawString("C",dm.width-15,xax+15);
                g.drawString("P",yax-12,12);
                f = new Font("TimesRoman",Font.PLAIN,11);
                g.setFont(f);
                g.setColor(Color.black);
                for(i=50;i>(-NEGY);i-=20){
                        g.drawString(" "+i,yax-22,(int)Math.round((POSY-i)*ystep)+5);
                }
                if(app1.PLOT_AV == 1){                
                        a=app1.Particles[app1.Num_Particles.getSelectedIndex()];
                        b=(int)Math.round(app1.CDiff.total_sum/app1.CDiff.n_runs);
                        app1.xs[app1.is]=a;
                        app1.ys[app1.is]=b;
                        app1.zs[app1.is]=1;
                        app1.is++;
                        app1.PLOT_AV=0;
                }
                if(app1.PLOT_REAL == 1){
                        a=app1.Particles[app1.Num_Particles.getSelectedIndex()];
                        b=(int)Math.round(app1.CDiff.Rin-app1.CDiff.Rout);
                        app1.xs[app1.is]=a;
                        app1.ys[app1.is]=b;
                        app1.zs[app1.is]=2;
                        app1.is++;
                        app1.PLOT_REAL=0;
                }
                for(i=0;i<(app1.is);i++){
                        a=app1.xs[i];
                        b=app1.ys[i];
                        N=4;
                        int hit=0;
                        for(j=i;j>=0;j--){
                                if (a == app1.xs[j]){
                                        if(b == app1.ys[j]){
                                                if(j != i){
                                                        hit++;
                                                        if(hit == 5){
                                                                N++;
                                                                hit=0;
                                                        }        
                                                }       
                                        }
                                }
                        }
                        if(app1.zs[i] == 1){
                                g.setColor(Color.red);
                        }
                        if(app1.zs[i] == 2){ 
                                g.setColor(Color.green);
                        }
                        g.fillOval((int)Math.round((NEGX+a)*xstep-N/2),(int)Math.round((POSY-b)*ystep)-N/3,N,N);
                }
                if(app1.ERR == 1){
                        g.setColor(Color.black);
                        g.drawString("Not enough data",dm.width-75,15);
                        app1.ERR=0;        
                }
                if(app1.PLT == 1){
                        g.setColor(Color.black);
                        Str2=Str1;
                        Str1=Str0;
                        Str0="P="+(float)((float)Math.round(app1.a*100)/(float)100.0)+"+"+(float)((float)Math.round(app1.b*100)/(float)100.0)+"dC";
                        x1=(int)Math.round((NEGX-3)*xstep);
                        y1=(int)Math.round((POSY-(app1.a+app1.b*(-3)))*ystep);
                        x2=(int)Math.round((NEGX+30)*xstep);
                        y2=(int)Math.round((POSY-(app1.a+app1.b*30))*ystep);
                        g.drawLine(x1,y1,x2,y2);
                        app1.PLT=0;
                } 
                g.setColor(Color.black);
                g.drawString(Str0,dm.width-75,15);
                g.drawString(Str1,dm.width-75,30);
                g.drawString(Str2,dm.width-75,45);       
         }
}