/*
 * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED 
 * All rights reserved.
 * 
 * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited,
 * non-exclusive and royalty-free license to use, copy, modify and
 * redistribute, solely for non-commercial and non-profit purposes, the
 * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and
 * conditions of this Agreement. For the avoidance of doubt, you acquire
 * only a limited right to use the SOFTWARE hereunder, and FLL or any
 * third party developer retains all rights, including but not limited to
 * copyrights, in and to the SOFTWARE.
 * 
 * (1) FLL does not grant you a license in any way for commercial
 * purposes. You may use the SOFTWARE only for non-commercial and
 * non-profit purposes only, such as academic, research and internal
 * business use.
 * (2) The SOFTWARE is protected by the Copyright Law of Japan and
 * international copyright treaties. If you make copies of the SOFTWARE,
 * with or without modification, as permitted hereunder, you shall affix
 * to all such copies of the SOFTWARE the above copyright notice.
 * (3) An explicit reference to this SOFTWARE and its copyright owner
 * shall be made on your publication or presentation in any form of the
 * results obtained by use of the SOFTWARE.
 * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
 * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
 * for such modification or the source code of the modified part of the
 * SOFTWARE.
 * 
 * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL
 * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND
 * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS
 * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES'
 * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY
 * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY.
 * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT,
 * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY
 * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL
 * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES
 * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES
 * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY
 * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF
 * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART
 * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY
 * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
 * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
 *
 * $OpenXM: OpenXM_contrib2/asir2000/plot/plotp.c,v 1.15 2011/08/11 06:25:06 saito Exp $ 
*/
#include "ca.h"
#include "parse.h"
#include "ox.h"
#include "ifplot.h"
#include <math.h>

#if defined(sun) && !defined(__svr4__)
#define EXP10(a) exp10(a)
#else
#define EXP10(a) pow(10.0,a)
#endif

#if defined(VISUAL)
static POINT oldpos;
#endif

#if defined(INTERVAL)
unsigned long GetColor(Display *, char *);

unsigned long GetColor(Display *dis, char *color_name)
{
	Colormap cmap;
	XColor near_color, true_color;

	cmap = DefaultColormap( dis, 0 );
	XAllocNamedColor( dis, cmap, color_name, &near_color, &true_color );
	return( near_color.pixel );
}

void area_print(DISPLAY *display, int **mask, struct canvas *can, int GXcode)
{
	int ix, iy;
	XImage *image;
	unsigned long color, black, white;
	int wc;

	flush();
	black = GetColor(display, "black");
	white = GetColor(display, "white");
	image = XGetImage(display, can->pix, 0, 0, can->width, can->height,
		-1, ZPixmap);
	for(iy = 0; iy < can->height; iy++){
		for(ix = 0; ix < can->width; ix++){
			color = XGetPixel(image, ix, iy);
			if (color == white) wc = -1;
			else if (color == black) wc = 0;
			else wc = 1;
			if ( wc != 0 ) {//XPutPixel(image,ix,iy,black);
				if ( mask[iy][ix] == 0 ) XPutPixel(image, ix, iy, black);
				else {
					switch (GXcode) {
					case 3: /* copy case */
						if (mask[iy][ix] == 1 ) XPutPixel(image, ix, iy, can->color);
						else XPutPixel(image, ix, iy, BackPixel);
						break;
					case 1: /* and case */
						if ((mask[iy][ix] == 1) && (wc == 1))
							XPutPixel(image, ix, iy, can->color);
						else XPutPixel(image, ix, iy, BackPixel);
						break;
					case 7: /* or case */
						if ((mask[iy][ix] == 1) || (wc == 1))
							XPutPixel(image, ix, iy, can->color);
						else XPutPixel(image, ix, iy, BackPixel);
						break;
					case 6: /* xor case */
						if ((mask[iy][ix] == 1) ^ (wc == 1))
							XPutPixel(image, ix, iy, can->color);
						else XPutPixel(image, ix, iy, BackPixel);
						break;
					}
				}
			}
		}
	}
	XPutImage(display,can->pix,drawGC,image,0,0,0,0,can->width, can->height);
	count_and_flush();
	flush();
}
#endif

void if_print(DISPLAY *display,double **tab,struct canvas *can)
{
	int ix,iy,width,height;
	double *px,*px1,*px2;
	DRAWABLE pix;

	if ( can->mode == MODE_CONPLOT ) {
		con_print(display,tab,can); return;
	}
	flush();
	width = can->width; height = can->height; pix = can->pix;
	for( ix=0; ix<width-1; ix++ )
		for(iy=0, px=tab[ix], px1 = tab[ix+1], px2 = px+1;
			iy<height-1 ;iy++, px++, px1++, px2++ )
			if ( ((*px >= 0) && ((*px1 <= 0) || (*px2 <= 0))) ||
				 ((*px <= 0) && ((*px1 >= 0) || (*px2 >= 0))) ) {
				DRAWPOINT(display,pix,cdrawGC,ix,height-iy-1);
				count_and_flush();
			}
	flush();
}

#define MEMORY_DRAWPOINT(a,len,x,y) (((a)[(len)*(y)+((x)>>3)]) |= (1<<((x)&7)))

void memory_if_print(double **tab,struct canvas *can,BYTEARRAY *bytes)
{
	int ix,iy,width,height;
	double *px,*px1,*px2;
	unsigned char *array;
	int scan_len;

	if ( can->mode == MODE_CONPLOT ) {
		memory_con_print(tab,can,bytes); return;
	}
	width = can->width; height = can->height;

	/* scan_len = byte length of the scan line */
	scan_len = (width+7)/8;
	MKBYTEARRAY(*bytes,scan_len*height);
	array = BDY(*bytes);
	for( ix=0; ix<width-1; ix++ )
		for(iy=0, px=tab[ix], px1 = tab[ix+1], px2 = px+1;
			iy<height-1 ;iy++, px++, px1++, px2++ )
			if ( ((*px >= 0) && ((*px1 <= 0) || (*px2 <= 0))) ||
				 ((*px <= 0) && ((*px1 >= 0) || (*px2 >= 0))) ) {
				MEMORY_DRAWPOINT(array,scan_len,ix,height-iy-1);
			}
}

void con_print(DISPLAY *display,double **tab,struct canvas *can)
{
	int ix,iy,iz,width,height,pas,pai,len;
	double zstep,z;
	double *px,*px1,*px2;
	DRAWABLE pix;
	POINT *pa,*pa1;
	struct pa *parray;

	
	width = can->width; height = can->height; pix = can->pix;
	zstep = (can->zmax-can->zmin)/can->nzstep;
	can->pa = (struct pa *)MALLOC((can->nzstep+1)*sizeof(struct pa));
	pas = width;
	pa = (POINT *)ALLOCA(pas*sizeof(POINT));
	initmarker(can,"Drawing...");
	for ( z = can->zmin, iz = 0; z <= can->zmax; z += zstep, iz++ ) {
		marker(can,DIR_Z,iz);
		pai = 0;	
		for( ix=0; ix<width-1; ix++ )
			for(iy=0, px=tab[ix], px1 = tab[ix+1], px2 = px+1;
				iy<height-1 ;iy++, px++, px1++, px2++ )
				if ( ((*px >= z) && ((*px1 <= z) || (*px2 <= z))) ||
				 	 ((*px <= z) && ((*px1 >= z) || (*px2 >= z))) ) {
					if ( pai == pas ) {
						pa1 = (POINT *)ALLOCA(2*pas*sizeof(POINT));
						bcopy((char *)pa,(char *)pa1,pas*sizeof(POINT)); pa = pa1;
						pas += pas;
					}
					XC(pa[pai]) = ix; YC(pa[pai]) = height-iy-1; pai++;
				}
		can->pa[iz].length = pai;
		if ( pai ) {
			pa1 = (POINT *)MALLOC(pai*sizeof(POINT));
			bcopy((char *)pa,(char *)pa1,pai*sizeof(POINT));
			can->pa[iz].pos = pa1;
		}
	}
	flush();
	for ( parray = can->pa, iz = 0; iz <= can->nzstep; iz++, parray++ )
		for ( pa = parray->pos, len = parray->length, ix = 0; ix < len; ix++ ) {
			DRAWPOINT(display,pix,cdrawGC,XC(pa[ix]),YC(pa[ix]));
			count_and_flush();
		}
	flush();
}

void memory_con_print(double **tab,struct canvas *can,BYTEARRAY *bytes)
{
	int ix,iy,iz,width,height,pas,pai,len;
	double zstep,z;
	double *px,*px1,*px2;
	POINT *pa,*pa1;
	struct pa *parray;
	unsigned char *array;
	int scan_len;

	width = can->width; height = can->height;
	zstep = (can->zmax-can->zmin)/can->nzstep;
	can->pa = (struct pa *)MALLOC((can->nzstep+1)*sizeof(struct pa));
	pas = width;
	pa = (POINT *)ALLOCA(pas*sizeof(POINT));
	for ( z = can->zmin, iz = 0; z <= can->zmax; z += zstep, iz++ ) {
		pai = 0;	
		for( ix=0; ix<width-1; ix++ )
			for(iy=0, px=tab[ix], px1 = tab[ix+1], px2 = px+1;
				iy<height-1 ;iy++, px++, px1++, px2++ )
				if ( ((*px >= z) && ((*px1 <= z) || (*px2 <= z))) ||
				 	 ((*px <= z) && ((*px1 >= z) || (*px2 >= z))) ) {
					if ( pai == pas ) {
						pa1 = (POINT *)ALLOCA(2*pas*sizeof(POINT));
						bcopy((char *)pa,(char *)pa1,pas*sizeof(POINT)); pa = pa1;
						pas += pas;
					}
					XC(pa[pai]) = ix; YC(pa[pai]) = height-iy-1; pai++;
				}
		can->pa[iz].length = pai;
		if ( pai ) {
			pa1 = (POINT *)MALLOC(pai*sizeof(POINT));
			bcopy((char *)pa,(char *)pa1,pai*sizeof(POINT));
			can->pa[iz].pos = pa1;
		}
	}
	/* scan_len = byte length of the scan line */
	scan_len = (width+7)/8;
	MKBYTEARRAY(*bytes,scan_len*height);
	array = BDY(*bytes);
	for ( parray = can->pa, iz = 0; iz <= can->nzstep; iz++, parray++ )
		for ( pa = parray->pos, len = parray->length, ix = 0; ix < len; ix++ ) {
			MEMORY_DRAWPOINT(array,scan_len,XC(pa[ix]),YC(pa[ix]));
		}
}

void memory_print(struct canvas *can,BYTEARRAY *bytes)
{
	int len,scan_len,i;
	POINT *pa;
	char *array;

	/* scan_len = byte length of the scan line */
	scan_len = (can->width+7)/8;
	MKBYTEARRAY(*bytes,scan_len*can->height);
	array = BDY(*bytes);

	len = can->pa[0].length;
	pa = can->pa[0].pos;
	for ( i = 0; i < len; i++ ) {
		MEMORY_DRAWPOINT(array,scan_len,pa[i].x,pa[i].y);
	}
}

void qif_print(DISPLAY *display,char **tab,struct canvas *can)
{
	int ix,iy,width,height;
	char *px;
	DRAWABLE pix;

	flush();
	width = can->width; height = can->height; pix = can->pix;
	for( ix = 0; ix < width; ix++ )
		for(iy = 0, px = tab[ix]; iy < height ;iy++, px++ )
			if ( *px ) {
				DRAWPOINT(display,pix,cdrawGC,ix,height-iy-1);
				count_and_flush();
			}
	flush();
}

void plot_print(DISPLAY *display,struct canvas *can)
{
	int len;
	POINT *pa;

#if defined(VISUAL)	
	len = can->pa[0].length;
	pa = can->pa[0].pos;

	Polyline(can->pix,pa,len);
//	for ( i = 0; i < len-1; i++ ) {
//		DRAWLINE(display,can->pix,drawGC,pa[i].x,pa[i].y,pa[i+1].x,pa[i+1].y);
//	}
#else
	XDrawLines(display,can->pix,drawGC,
		can->pa[0].pos,can->pa[0].length,CoordModeOrigin);
	XFlush(display);
#endif
}

void draw_point(DISPLAY *display,struct canvas *can,int x,int y,int color)
{
#if defined(VISUAL)
	HDC dc;

	SetPixel(can->pix,x,y,(COLORREF)color);
	dc = GetDC(can->hwnd);
	SetPixel(dc,x,y,(COLORREF)color);
	ReleaseDC(can->hwnd,dc);
#else
	set_drawcolor(color);
	DRAWPOINT(display,can->pix,cdrawGC,x,y);
	DRAWPOINT(display,can->window,cdrawGC,x,y);
	XFlush(display);
#endif
}

void draw_line(DISPLAY *display,struct canvas *can,int x,int y,int u,int v,int color)
{
#if defined(VISUAL)
	HDC dc;
	HPEN pen,oldpen;

	if ( color ) {
		pen = CreatePen(PS_SOLID,1,color);
		oldpen = SelectObject(can->pix,pen);
		DRAWLINE(display,can->pix,drawGC,x,y,u,v);
		SelectObject(can->pix,oldpen);

		dc = GetDC(can->hwnd);
		oldpen = SelectObject(dc,pen);
		DRAWLINE(display,dc,drawGC,x,y,u,v);
		SelectObject(dc,oldpen);
		ReleaseDC(can->hwnd,dc);

		DeleteObject(pen);
	} else {
		DRAWLINE(display,can->pix,drawGC,x,y,u,v);
		dc = GetDC(can->hwnd);
		DRAWLINE(display,dc,drawGC,x,y,u,v);
		ReleaseDC(can->hwnd,dc);
	}
#else
	set_drawcolor(color);
	DRAWLINE(display,can->pix,cdrawGC,x,y,u,v);
	DRAWLINE(display,can->window,cdrawGC,x,y,u,v);
	XFlush(display);
#endif
}

void draw_character_string(DISPLAY *display,struct canvas *can,int x,int y,char *str,int color)
{
#if defined(VISUAL)
	HDC dc;
	COLORREF oldcolor;

	if ( color ) {
		oldcolor = SetTextColor(can->pix,color);
		DRAWSTRING(display,can->pix,drawGC,x,y,str,strlen(str));
		SetTextColor(can->pix,oldcolor);

		dc = GetDC(can->hwnd);
		oldcolor = SetTextColor(dc,color);
		DRAWSTRING(display,dc,drawGC,x,y,str,strlen(str));
		SetTextColor(dc,oldcolor);
		ReleaseDC(can->hwnd,dc);
	} else {
		DRAWSTRING(display,can->pix,drawGC,x,y,str,strlen(str));
		dc = GetDC(can->hwnd);
		DRAWSTRING(display,dc,drawGC,x,y,str,strlen(str));
		ReleaseDC(can->hwnd,dc);
	}
#else
	set_drawcolor(color);
	DRAWSTRING(display,can->pix,cdrawGC,x,y,str,strlen(str));
	DRAWSTRING(display,can->window,cdrawGC,x,y,str,strlen(str));
	XFlush(display);
#endif
}

#define D 5

void pline(DISPLAY *display,struct canvas *can,DRAWABLE d)
{
	double w,w1,e,n;
	int x0,y0,x,y,xadj,yadj;
	char buf[BUFSIZ];
	double adjust_scale();

	/* XXX : should be cleaned up */
	if ( can->noaxis || (can->mode == MODE_PLOT && !can->pa) )
		return;
	if ( can->mode == MODE_INTERACTIVE )
		return;

	xadj = yadj = 0;
	if ( (can->xmin < 0) && (can->xmax > 0) ) {
		x0 = (int)((can->width)*(-can->xmin/(can->xmax-can->xmin)));
		DRAWLINE(display,d,dashGC,x0,0,x0,can->height);
	} else if ( can->xmin >= 0 )
		x0 = 0;
	else
		x0 = can->width-D;
	if ( (can->ymin < 0) && (can->ymax > 0) ) {
		y0 = (int)((can->height)*(can->ymax/(can->ymax-can->ymin)));
		DRAWLINE(display,d,dashGC,0,y0,can->width,y0);
	} else if ( can->ymin >= 0 )
		y0 = can->height;
	else
		y0 = D;

	/* scale on x-axis */
	w = can->xmax-can->xmin; 
	w1 = w * DEFAULTWIDTH/can->width;
	e = adjust_scale(EXP10(floor(log10(w1))),w1);
	for ( n = ceil(can->xmin/e); n*e<= can->xmax; n++ ) {
		x = (int)(can->width*(n*e-can->xmin)/w);
		DRAWLINE(display,d,drawGC,x,y0,x,y0-D);
		sprintf(buf,"%g",n*e);
		DRAWSTRING(display,d,scaleGC,x+2,y0+2,buf,strlen(buf));
	}

	/* scale on y-axis */
	w = can->ymax-can->ymin;
	w1 = w * DEFAULTHEIGHT/can->height;
	e = adjust_scale(EXP10(floor(log10(w1))),w1);
	for ( n = ceil(can->ymin/e); n*e<= can->ymax; n++ ) {
		y = (int)(can->height*(1-(n*e-can->ymin)/w));
		DRAWLINE(display,d,drawGC,x0,y,x0+D,y);
		sprintf(buf,"%g",n*e);
		if ( can->xmax <= 0 ) {
#if !defined(VISUAL)
			xadj = TEXTWIDTH(sffs,buf,strlen(buf));
#else
			SIZE size;

			GetTextExtentPoint32(d,buf,strlen(buf),&size);
			xadj = size.cx;
#endif
		}
		DRAWSTRING(display,d,scaleGC,x0+2-xadj,y+2,buf,strlen(buf));
	}
}

double adjust_scale(double e,double w)
{
	switch ( (int)floor(w/e) ) {
		case 1:
			return e/4; break;
		case 2: case 3:
			return e/2; break;
		case 4: case 5: case 6: case 7:
			return e; break;
			break;
		case 8: case 9: case 10: default:
			return 2*e; break;
	}
}

void initmarker(struct canvas *can,char *message)
{
#if defined(VISUAL)
	can->real_can->percentage = 0;
	can->real_can->prefix = message;
#else
	XawScrollbarSetThumb(can->xdone,0.0,0.0);
	XawScrollbarSetThumb(can->ydone,1.0,0.0);
	XFlush(display);
#endif
}

void marker(struct canvas *can,int dir,int p)
{
#if defined(VISUAL)
	if ( dir == DIR_X )
		can->real_can->percentage = (int)ceil((float)p/(float)can->real_can->width*100);
	else if ( dir == DIR_Y )
		can->real_can->percentage = (int)ceil((float)p/(float)can->real_can->height*100);
	else
		can->real_can->percentage = (int)ceil((float)p/(float)can->real_can->nzstep*100);
#else
	if ( dir == DIR_X ) {
		XawScrollbarSetThumb(can->xdone,(float)p/(float)can->width,0.05);
		count_and_flush();
	} else if ( dir == DIR_Y ) {
		XawScrollbarSetThumb(can->ydone,1.0-(float)p/(float)can->height,0.05);
		count_and_flush();
	} else {
		XawScrollbarSetThumb(can->ydone,1.0-(float)p/(float)can->nzstep,0.05);
		flush();
	}
#endif
}

void define_cursor(WINDOW w,CURSOR cur)
{
#if !defined(VISUAL)
	XDefineCursor(display,w,cur); flush();
#endif
}

static int flush_count;
#if defined(linux) || defined(__FreeBSD__) || defined(__NetBSD__)
#define MAX_COUNT 64
#else
#define MAX_COUNT 32
#endif

void count_and_flush() {
#if !defined(VISUAL)
	if ( ++flush_count == MAX_COUNT )
		flush();
#endif
}

void flush() {
#if !defined(VISUAL)
	flush_count = 0;
	XFlush(display);
#endif
}