===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/plot/plotp.c,v
retrieving revision 1.6
retrieving revision 1.15
diff -u -p -r1.6 -r1.15
--- OpenXM_contrib2/asir2000/plot/plotp.c	2000/11/09 01:51:13	1.6
+++ OpenXM_contrib2/asir2000/plot/plotp.c	2011/08/11 06:25:06	1.15
@@ -45,7 +45,7 @@
  * 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.5 2000/11/07 06:06:40 noro Exp $ 
+ * $OpenXM: OpenXM_contrib2/asir2000/plot/plotp.c,v 1.14 2011/08/10 04:51:58 saito Exp $ 
 */
 #include "ca.h"
 #include "parse.h"
@@ -63,11 +63,73 @@
 static POINT oldpos;
 #endif
 
-if_print(display,tab,can)
-DISPLAY *display;
-double **tab;
-struct canvas *can;
+#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;
@@ -82,17 +144,41 @@ struct canvas *can;
 			iy<height-1 ;iy++, px++, px1++, px2++ )
 			if ( ((*px >= 0) && ((*px1 <= 0) || (*px2 <= 0))) ||
 				 ((*px <= 0) && ((*px1 >= 0) || (*px2 >= 0))) ) {
-				DRAWPOINT(display,pix,drawGC,ix,height-iy-1);
+				DRAWPOINT(display,pix,cdrawGC,ix,height-iy-1);
 				count_and_flush();
 			}
 	flush();
 }
 
-con_print(display,tab,can)
-DISPLAY *display;
-double **tab;
-struct canvas *can;
+#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;
@@ -132,17 +218,78 @@ struct canvas *can;
 	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,drawGC,XC(pa[ix]),YC(pa[ix]));
+			DRAWPOINT(display,pix,cdrawGC,XC(pa[ix]),YC(pa[ix]));
 			count_and_flush();
 		}
 	flush();
 }
 
-qif_print(display,tab,can)
-DISPLAY *display;
-char **tab;
-struct canvas *can;
+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;
@@ -152,17 +299,15 @@ struct canvas *can;
 	for( ix = 0; ix < width; ix++ )
 		for(iy = 0, px = tab[ix]; iy < height ;iy++, px++ )
 			if ( *px ) {
-				DRAWPOINT(display,pix,drawGC,ix,height-iy-1);
+				DRAWPOINT(display,pix,cdrawGC,ix,height-iy-1);
 				count_and_flush();
 			}
 	flush();
 }
 
-plot_print(display,can)
-DISPLAY *display;
-struct canvas *can;
+void plot_print(DISPLAY *display,struct canvas *can)
 {
-	int len,i;
+	int len;
 	POINT *pa;
 
 #if defined(VISUAL)	
@@ -180,11 +325,7 @@ struct canvas *can;
 #endif
 }
 
-draw_point(display,can,x,y,color)
-DISPLAY *display;
-struct canvas *can;
-int x,y;
-int color;
+void draw_point(DISPLAY *display,struct canvas *can,int x,int y,int color)
 {
 #if defined(VISUAL)
 	HDC dc;
@@ -194,17 +335,14 @@ int color;
 	SetPixel(dc,x,y,(COLORREF)color);
 	ReleaseDC(can->hwnd,dc);
 #else
-	DRAWPOINT(display,can->pix,drawGC,x,y);
-	DRAWPOINT(display,can->window,drawGC,x,y);
+	set_drawcolor(color);
+	DRAWPOINT(display,can->pix,cdrawGC,x,y);
+	DRAWPOINT(display,can->window,cdrawGC,x,y);
 	XFlush(display);
 #endif
 }
 
-draw_line(display,can,x,y,u,v,color)
-DISPLAY *display;
-struct canvas *can;
-int x,y,u,v;
-int color;
+void draw_line(DISPLAY *display,struct canvas *can,int x,int y,int u,int v,int color)
 {
 #if defined(VISUAL)
 	HDC dc;
@@ -230,20 +368,48 @@ int color;
 		ReleaseDC(can->hwnd,dc);
 	}
 #else
-	DRAWLINE(display,can->pix,drawGC,x,y,u,v);
-	DRAWLINE(display,can->window,drawGC,x,y,u,v);
+	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
 
-pline(display,can,d)
-DISPLAY *display;
-struct canvas *can;
-DRAWABLE d;
+void pline(DISPLAY *display,struct canvas *can,DRAWABLE d)
 {
-	double w,w1,k,e,n;
+	double w,w1,e,n;
 	int x0,y0,x,y,xadj,yadj;
 	char buf[BUFSIZ];
 	double adjust_scale();
@@ -275,7 +441,7 @@ DRAWABLE d;
 	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;
+		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));
@@ -286,7 +452,7 @@ DRAWABLE d;
 	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);
+		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 ) {
@@ -303,8 +469,7 @@ DRAWABLE d;
 	}
 }
 
-double adjust_scale(e,w)
-double e,w;
+double adjust_scale(double e,double w)
 {
 	switch ( (int)floor(w/e) ) {
 		case 1:
@@ -319,9 +484,7 @@ double e,w;
 	}
 }
 
-initmarker(can,message)
-struct canvas *can;
-char *message;
+void initmarker(struct canvas *can,char *message)
 {
 #if defined(VISUAL)
 	can->real_can->percentage = 0;
@@ -333,17 +496,15 @@ char *message;
 #endif
 }
 
-marker(can,dir,p)
-struct canvas *can;
-int dir,p;
+void marker(struct canvas *can,int dir,int p)
 {
 #if defined(VISUAL)
 	if ( dir == DIR_X )
-		can->real_can->percentage = ceil((float)p/(float)can->real_can->width*100);
+		can->real_can->percentage = (int)ceil((float)p/(float)can->real_can->width*100);
 	else if ( dir == DIR_Y )
-		can->real_can->percentage = ceil((float)p/(float)can->real_can->height*100);
+		can->real_can->percentage = (int)ceil((float)p/(float)can->real_can->height*100);
 	else
-		can->real_can->percentage = ceil((float)p/(float)can->real_can->nzstep*100);
+		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);
@@ -358,9 +519,7 @@ int dir,p;
 #endif
 }
 
-define_cursor(w,cur)
-WINDOW w;
-CURSOR cur;
+void define_cursor(WINDOW w,CURSOR cur)
 {
 #if !defined(VISUAL)
 	XDefineCursor(display,w,cur); flush();
@@ -374,14 +533,14 @@ static int flush_count;
 #define MAX_COUNT 32
 #endif
 
-count_and_flush() {
+void count_and_flush() {
 #if !defined(VISUAL)
 	if ( ++flush_count == MAX_COUNT )
 		flush();
 #endif
 }
 
-flush() {
+void flush() {
 #if !defined(VISUAL)
 	flush_count = 0;
 	XFlush(display);