#include <windows.h>
#include <elecwnd.h>
#include <canvas.h>
#include <dlg.h>
#include <main.h>

#include <fstream.h>;
#include <string.h>;

 #include <math.h>
//****************************************

double ZVal(
	double X1, double Y1, double Z1,
	double X2, double Y2, double Z2,
	double X3, double Y3, double Z3,
	double X, double Y )

{
double denom = (X1*Y2 - X1*Y3 - X2*Y1 + X2*Y3 + X3*Y1 - X3*Y2);

if(denom == 0 ) { return 0; }

double Z =
(
  X*(Z1*Y2 - Z1*Y3 - Y1*Z2 + Y1*Z3 + Z2*Y3 - Y2*Z3)
+ Y*( - Z1*X2+ Z1*X3 + X1*Z2 - X1*Z3 + X2*Z3 - X3*Z2)
+ Z1*X2*Y3 - Z1*X3*Y2 - X1*Z2*Y3 + X1*Y2*Z3 - X2*Y1*Z3 + X3*Y1*Z2

)
/denom;
return Z;
}




double Round(double x)
       {
	double xx = floor(x);
	if( (x - xx)<0.5 )
	  { return xx;       }
	else
	  { return ceil(x);  }
       }

//****************************************

//ColorScale colscal(16, -100, 100);
//extern ColorScale colscal(16,-4000,4000);
extern  colscal;
TArray<Color>;
//***************** stream operators *************************

ostream&  operator << (ostream& stm,   const Color& c)
{
stm<<"COLOR"<<","<<c._r<<","<<c._g<<","<<c._b;
}


istream&  operator >> (istream& istm,  Color& c)
{
char tmpc;
int i;
//eating initial string
//      C    O     L      O     R     ,
istm>>tmpc>>tmpc>>tmpc>>tmpc>>tmpc>>tmpc;

istm >>c._r>>tmpc >>c._g>>tmpc >>c._b;

return istm;
}

//---------------
//friend
ostream&  operator << (ostream& stm,   const ColorScale& cs)
{
stm<<endl<<cs._cnum;
for(int i=0;i < cs._cnum; i++ )
{

stm<<endl<<cs.GetColor(i);
}

}



//friend
istream&  operator >> (istream& istm,  ColorScale& cs)
{
const long  s = cs.Size();
//cout<<endl<<"---- ColorScale operator >>  ----";
int size=-1;
//long size=-1;
istm>>size;
if( 0<size )
{
//cout<<endl<<size;
cs.Resize(size);
cs._cnum=size;
for( int i=0;i<size;i++)
//for( long i=0;i<size;i++)
   {

   if ( istm.good() )
      { istm>>cs.p()[i]; }
   else
      {}
   }

}
else
{
//cout<<endl<<"file reading error";
}



}




//*******************************************

ElectrodesWindow::Init(int maxHid, int maxVid)
{
_maxHid = maxHid;
_maxVid = maxVid;
_raster.ResizeRaster( _maxHid, _maxVid );
_status = EWND_OK;
}



void ElectrodesWindow::Register(HANDLE hInstance)
{
	  WNDCLASS    wndclass ;
	  wndclass.style         = CS_HREDRAW | CS_VREDRAW ;
	  wndclass.lpfnWndProc   = ElecWndProc ;
//wndclass.lpfnWndProc   = NULL;
	  wndclass.cbClsExtra    = 0 ;
	  wndclass.cbWndExtra    = 0 ;
	  wndclass.hInstance     = hInstance ;
	  wndclass.hIcon         = LoadIcon (NULL, IDI_APPLICATION) ;
	  wndclass.hCursor       = LoadCursor (NULL, IDC_ARROW) ;
	  wndclass.hbrBackground = GetStockObject (WHITE_BRUSH) ;
	  wndclass.lpszMenuName  = NULL ;
	  wndclass.lpszClassName = "EWindow" ;

	  RegisterClass (&wndclass) ;
}



ElectrodesWindow::ElectrodesWindow(HWND hwndParent, HANDLE hInstance, HANDLE hPrevInstance)
:_hwndParent(hwndParent),
_hInstance(hInstance),
_status(EWND_EMPTY), _maxHid(0), _maxVid(0), _cxClient(0), _cyClient(0),
_raster(0,0),
_pcol(0),
_balloon( CreateDC ("DISPLAY", NULL, NULL, NULL)  )
{


	  if (!hPrevInstance)
	  {
	      Register(hInstance);
	  }
	      //Register(hInstance);

     _hwnd =
     CreateWindow ("EWindow",               // window class name
		    "E Window",         // window caption
	  WS_CHILD |WS_OVERLAPPEDWINDOW|WS_VISIBLE,     // window style
		    0,           // initial x position
		    0,           // initial y position
		    100,           // initial x size
		    100,           // initial y size
		    hwndParent,                    // parent window handle
		    25,                    // window menu handle
		    hInstance,         // program instance handle
		    (LPSTR) this) ;	     // creation parameters

     //ShowWindow (hwnd, Main::nCmdShow) ;
     //UpdateWindow (hwnd) ;




//          _hwnd = CreateWindow ("STATIC", NULL,
  //		      WS_BORDER | WS_CHILD | WS_VISIBLE,
    //		      0, 0, 0, 0,
      //		      hwndParent, 25, hInstance, NULL) ;


}




ElectrodesWindow::ElectrodesWindow(HWND hwndParent,
	     HANDLE hInstance,
	     HANDLE hPrevInstance, ColorScale& colscal)
:
_hInstance(hInstance),
_hwndParent(hwndParent),
_status(EWND_EMPTY), _maxHid(0), _maxVid(0), _cxClient(0), _cyClient(0),
_raster(0,0),
_pcol(&colscal),
_balloon( CreateDC ("DISPLAY", NULL, NULL, NULL)  )
{

	  if (!hPrevInstance)
	  {
	      Register(hInstance);
	  }
	      //Register(hInstance);

     _hwnd =
     CreateWindow ("EWindow",               // window class name
		    "E Window",         // window caption
	  WS_CHILD |WS_OVERLAPPEDWINDOW|WS_VISIBLE,     // window style
		    0,           // initial x position
		    0,           // initial y position
		    100,           // initial x size
		    100,           // initial y size
		    hwndParent,                    // parent window handle
		    25,                    // window menu handle
		    hInstance,         // program instance handle
		    (LPSTR) this) ;	     // creation parameters

     //ShowWindow (hwnd, Main::nCmdShow) ;
     //UpdateWindow (hwnd) ;


}




ElectrodesWindow::Size(int cx, int cy,int width,int height)
{
      MoveWindow (_hwnd, cx, cy,width,height,TRUE) ;

_cxClient = width;
_cyClient = height;
}



ElectrodesWindow::PaintRaster(ColorScale &colscal)
{


    RECT rect;
    POINT p;
    ScreenCanvas s( hwnd() );

    GetClientRect(hwnd(),(LPRECT) &rect);
    //rect.right = rect.right-20;
    //rect.top = rect.top-20;
    SetViewportOrg( s.DC(), 0, rect.bottom) ;
    s.MapMode( MM_HIMETRIC ) ;

    s.Rect(10,10,100,100);
    s.Line(0,0,1000, 1000);


    DPtoLP(s.DC(), (LPPOINT) &rect,2);
    int nx = _raster.dimx();
    int ny = _raster.dimy();
    int rx = rect.right/nx;
    int ry = rect.top/ny;

//set random values in raster
//    for(int i1=1; i1<=nx; i1++)
//    {
//    for(int i2=1; i2<=ny; i2++)
//    {
//    double x =100.555*(rand() % 100) ;
//    Value tmpv = _raster.GetEntry(i1,i2) ;
//    tmpv.SetVal(x) ;
//    _raster.SetEntry(i1,i2,tmpv);
//    }
//    }

//ofstream fil("val.tst");

    for (int i=1; i<=nx; i++)
    {
    for (int ii=1; ii<=ny; ii++)
    {
       Value tmpv = _raster.GetEntry(i,ii) ;
       //Color col = colscal.GetColor(i);
       Color col = colscal.ValToCol(tmpv.u() );
       Brusch b(s.DC(), col.r(), col.g(), col.b() );
       //Brusch b(s.DC(), 255, 255, 255 );

       Value val = _raster.GetEntry(i,ii);
       //fil<<endl<<val;
       if( val.status() != 1 )
         {
	  s.Rect((i-1)*rx,(ii-1)*ry,i*rx,ii*ry);
         }
    }
    }




//ofstream ff("raster.tst");
//ff<<_raster;




}



ElectrodesWindow::PaintRaster()
{
if( _raster.dimx() != 0 && _raster.dimy() != 0)
  {

    RECT rect;
    POINT p;
    ScreenCanvas s( hwnd() );

    GetClientRect(hwnd(),(LPRECT) &rect);
    //rect.right = rect.right-20;
    //rect.top = rect.top-20;

    SetViewportOrg( s.DC(), 0, rect.bottom) ;
    s.MapMode( MM_HIMETRIC ) ;

    s.Rect(10,10,100,100);
    s.Line(0,0,1000, 1000);


    DPtoLP(s.DC(), (LPPOINT) &rect,2);
    int nx = _raster.dimx();
    int ny = _raster.dimy();
    int rx = rect.right/nx;
    int ry = rect.top/ny;

//set random values in raster
//    for(int i1=1; i1<=nx; i1++)
//    {
//    for(int i2=1; i2<=ny; i2++)
//    {
//    double x =100.555*(rand() % 100) ;
//    Value tmpv = _raster.GetEntry(i1,i2) ;
//    tmpv.SetVal(x) ;
//    _raster.SetEntry(i1,i2,tmpv);
//    }
//    }

ofstream fil("val.tst");

    Color col;
    for (int i=1; i<=nx; i++)
    {
    for (int ii=1; ii<=ny; ii++)
    {
       Value tmpv = _raster.GetEntry(i,ii) ;
       //Color col = colscal.GetColor(i);
       if( _pcol != 0 )
         {
	col = _pcol->ValToCol(tmpv.u() );
         }
	Brusch b(s.DC(), col.r(), col.g(), col.b() );
	//Brusch b(s.DC(), 255, 255, 255 );
	Value val = _raster.GetEntry(i,ii);

	fil<<endl<<val;

         if( val.status() != 1 )
           {
	  s.Rect((i-1)*rx,(ii-1)*ry,i*rx,ii*ry);
           }

    }
    }




ofstream ff("raster.tst");
ff<<_raster;

}


}


ElectrodesWindow::Paint()
{

    PAINTSTRUCT ps;
    RECT rect;
    PaintCanvas pc(_hwnd);



    GetClientRect( pc.hwnd(), (LPRECT) &rect );
    char buf[30];
    itoa(rect.right,buf,10);
    pc.Text(10,10,buf,4);


}



ElectrodesWindow::PaintRaster1(ColorScale &colscal)
{


    RECT rect;
    POINT p;
    PaintCanvas s( hwnd() );

    GetClientRect(hwnd(),(LPRECT) &rect);
    //rect.right = rect.right-20;
    //rect.top = rect.top-20;

    SetViewportOrg( s.DC(), 0, rect.bottom) ;
    s.MapMode( MM_HIMETRIC ) ;

    s.Rect(10,10,100,100);
    s.Line(0,0,1000, 1000);


    DPtoLP(s.DC(), (LPPOINT) &rect,2);
    int nx = _raster.dimx();
    int ny = _raster.dimy();
    int rx = rect.right/nx;
    int ry = rect.top/ny;

//set random values in raster
//    for(int i1=1; i1<=nx; i1++)
//    {
//    for(int i2=1; i2<=ny; i2++)
//    {
//    double x =100.555*(rand() % 100) ;
//    Value tmpv = _raster.GetEntry(i1,i2) ;
//    tmpv.SetVal(x) ;
//    _raster.SetEntry(i1,i2,tmpv);
//    }
//    }

//ofstream fil("val.tst");

    for (int i=1; i<=nx; i++)
    {
    for (int ii=1; ii<=ny; ii++)
    {
       Value tmpv = _raster.GetEntry(i,ii) ;
       //Color col = colscal.GetColor(i);
       Color col = colscal.ValToCol(tmpv.u() );
       Brusch b(s.DC(), col.r(), col.g(), col.b() );
       //Brusch b(s.DC(), 255, 255, 255 );

       Value val = _raster.GetEntry(i,ii);
       //fil<<endl<<val;
       if( val.status() != 1 )
         {
	  s.Rect((i-1)*rx,(ii-1)*ry,i*rx,ii*ry);
         }
    }
    }




//ofstream ff("raster.tst");
//ff<<_raster;




}


ElectrodesWindow::LButtonUp(WORD wParam, LONG lParam)
{
_balloon.HideBalloon();
}


ElectrodesWindow::LButtonDown(WORD wParam, LONG lParam)
{

   POINT p;
   p.x = LOWORD (lParam) ;
   p.y = HIWORD (lParam) ;

   RECT r;
   GetClientRect(hwnd(),(LPRECT)&r);

   int ixx, iyy;
   ixx =  1+( _raster.dimx()*p.x)/r.right ;
   iyy =  1+( _raster.dimy()*(r.bottom - p.y) )/r.bottom ;

if( ixx < 0 ||  _raster.dimx() < ixx) { MessageBeep(0); return 0; }
if( iyy < 0 ||  _raster.dimy() < iyy) { MessageBeep(0); return 0; }


Value v = _raster.GetEntry(ixx,iyy);
if( v.status() == 1 ) { MessageBeep(0); return 0; }

char bufx[30];
char bufy[30];
char bufv[30];
char buf[30]="electrode(";

gcvt(v.u(),5,bufv);
lstrcat(bufx,"");
itoa(ixx,bufx,10);
lstrcat(bufx,",");

itoa(iyy,bufy,10);
lstrcat(bufy,"");

lstrcat(buf,bufx);
lstrcat(buf,bufy);
lstrcat(buf,")");

lstrcat(buf,"\n");
lstrcat(buf,bufv);


//itoa(v.status(),buf,10);
//lstrcat(buf,"");


Font fff( _balloon.DC(), 15 );
Brusch b2( _balloon.DC(),255, 255, 0 );



     POINT pp = p ;
      ClientToScreen(hwnd(), (LPPOINT) &pp) ;

if( (p.x < r.right/2) && (p.y < r.bottom/2) )
_balloon.ShowBalloon(pp.x, pp.y, 30, buf, BLOWERRIGHT  );


if( (p.x > r.right/2) && (p.y < r.bottom/2)  )
_balloon.ShowBalloon(pp.x, pp.y, 30, buf, BLOWERLEFT  );


if( (p.x < r.right/2) && (p.y > r.bottom/2)  )
_balloon.ShowBalloon(pp.x, pp.y, 30, buf, BUPPERRIGHT  );

if( (p.x > r.right/2) && (p.y > r.bottom/2)  )
_balloon.ShowBalloon(pp.x, pp.y, 30, buf, BUPPERLEFT  );


//Neighbour investigation
double   dx =  ( _raster.dimx()*p.x)/r.right ;
double   dy =  ( _raster.dimy()*(r.bottom - p.y) )/r.bottom ;

double x1 = Round(dx);
double y1 = Round(dx);



}


ElectrodesWindow::RButtonDown(WORD wParam, LONG lParam, ColorScale &colscal )
{

   POINT p;
   p.x = LOWORD (lParam) ;
   p.y = HIWORD (lParam) ;

   RECT r;
   GetClientRect(hwnd(),(LPRECT)&r);

   int ixx, iyy;
   ixx =  1+( _raster.dimx()*p.x)/r.right ;
   iyy =  1+( _raster.dimy()*(r.bottom - p.y) )/r.bottom ;

if( ixx < 0 ||  _raster.dimx() < ixx) { MessageBeep(0); return 0; }
if( iyy < 0 ||  _raster.dimy() < iyy) { MessageBeep(0); return 0; }


Value v = _raster.GetEntry(ixx,iyy);
if( v.status() == 1 ) { MessageBeep(0); return 0; }

    int cnum = colscal.ValToColNum(v.u());

    Color col = colscal.GetColor(cnum);

    CtrlColorDialog ddd(1, Main::hInstance, col.r(),col.g(),col.b());
int dresult =  PerformDialog(hwnd(), &ddd);
if(dresult == CTRLCOLORDIALOG_OK)
{
MessageBeep(0);
colscal.SetColor( cnum, ddd.r(), ddd.g(), ddd.b() );
}



}




ElectrodesWindow::SaveRasterAsContourGraphics(char * path)
{
 ofstream file(path);

 file<<"(*^"<<
 endl<<"::[	frontEndVersion = "
     <<'"'<<"Microsoft Windows Mathematica Notebook Front End Version 2.1"<<'"'<<';'<<
 endl<<"	microsoftWindowsStandardFontEncoding;]";

 file<<endl<<':';
 file<<endl<<"Show[ContourGraphics[";


 if( _raster.dimx() != 0 && _raster.dimy() != 0)
  {
    file<<endl<<"{";
    for (int i=1; i<=_raster.dimy(); i++)
    {

    file<<"{";
    for (int ii=1; ii<=_raster.dimx(); ii++)
    {
       Value tmpv = _raster.GetEntry(ii,i) ;

	file<<tmpv.u();
	if( ii < _raster.dimx() ) { file<<","; }
    }
        file<<"}";
        if( i < _raster.dimy() ) { file<<","; }
        file<<endl;
    }
    file<<"}";

  }





 file<<"]]";
 file<<endl<<"^*)";


}