===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/plot/smoothing.c,v
retrieving revision 1.4
retrieving revision 1.9
diff -u -p -r1.4 -r1.9
--- OpenXM_contrib2/asir2000/plot/smoothing.c	2000/11/08 05:49:35	1.4
+++ OpenXM_contrib2/asir2000/plot/smoothing.c	2007/01/25 16:19:41	1.9
@@ -1,4 +1,7 @@
-/* $OpenXM: OpenXM_contrib2/asir2000/plot/smoothing.c,v 1.3 2000/10/24 01:53:53 takayama Exp $ */
+/* $OpenXM: OpenXM_contrib2/asir2000/plot/smoothing.c,v 1.8 2002/07/20 02:28:08 noro Exp $ */
+#include <stdio.h>
+#include <time.h>
+#include <math.h>
 #include "ca.h"
 #include "parse.h"
 #include "ox.h"
@@ -8,397 +11,99 @@
 #define MAG 1
 #define PRINT_XOFFSET   100
 #define PRINT_YOFFSET   100
+#define DPI 300
 
-struct polyLine {
-  int numberOfSegments;
-  int limit;
-  int *x;
-  int *y;
-};
-
 static FILE *Fp = NULL;
-
-static void *gcmalloc(int a);
-static struct polyLine *polyLine_new(void);
-static struct polyLine *polyLine_addNewSegment(struct polyLine *pl,
-											   int x,int y);
-static int polyLine_lastY(struct polyLine *pl);
-static void polyLine_outputPS(struct polyLine *pl);
-static void polyLine_output_bezier_PS(struct polyLine *pl);
-static void polyLine_outputProlog(int xmin, int ymin, int xmax, int ymax);
-static void polyLine_outputEpilog(void);
-static void polyLine_error(char *s);
-#define translateX(x)  (x*MAG+PRINT_XOFFSET)
-#define translateY(y)  (y*MAG+PRINT_YOFFSET)
-/* #define translateY(y) ((Ysize-y)*MAG+PRINT_YOFFSET) */
-#define IS_DOT    != 0
-#define IS_POLYLINE  > -1
-#define YES 1
-#define NO  0
-
-
-static int Xsize = 0;
-static int Ysize = 0;
-static int updatePolyLine(struct polyLine *pl[], int plSize,
-				   int prev[], int Curr[], int x);
-
+static void PSProlog(int, int, int, int);
+static void PSDashedLine(int, int, int, int);
 static void polyLine_outputPS_dashed_line(int x0,int y0,int x1,int y1);
-static int polyLine_pline(struct canvas *can);
+
 static int Strategy_generate_PS = 0;
 
-static void *gcmalloc(a)  {
-  void *m;
-  /* BUG:  interruption must be locked. */
-  m = (void *) GC_malloc(a);
-  if (m == NULL) polyLine_error("no memory");
-  return(m);
-}
+PSFromImage(FILE *fp,XImage *image,struct canvas *can) {
+  int x, y, ir, ig, ib;
+  float r, g, b;
+  int srcc, trg,grux, gruy, ix, iy;
+  float xwd, ywd, t;
 
-#define MINIMAL_POLYLINE_SIZE  10
-static struct polyLine *polyLine_new_with_size(int size) {
-  struct polyLine *pl;
-  pl = (struct polyLine *)gcmalloc(sizeof(struct polyLine));
-  pl-> limit = (size > 0? size: 1);
-  pl-> numberOfSegments = 0;
-  pl->x = (int *)gcmalloc(sizeof(int)*(pl->limit));
-  pl->y = (int *)gcmalloc(sizeof(int)*(pl->limit));
-  return(pl);
-}
-
-static struct polyLine *polyLine_new(void) {
-  return( polyLine_new_with_size(MINIMAL_POLYLINE_SIZE));
-}
-
-static struct polyLine *polyLine_addNewSegment(struct polyLine *pl,
-												int x,int y) {
-  int n,limit,i;
-  struct polyLine *new_pl;
-  limit = pl->limit;
-  n = pl->numberOfSegments;
-  if (n >= limit) {
-	new_pl = polyLine_new_with_size( 2*limit );
-	new_pl -> numberOfSegments = pl->numberOfSegments;
-	for (i=0; i<pl->numberOfSegments; i++) {
-	  new_pl->x[i] = pl->x[i];
-	  new_pl->y[i] = pl->y[i];
-	}
-	pl = new_pl;
+  Fp = fp;
+  srcc = 0xffffff;
+  PSProlog(0,0,can->width,can->height);
+  for(y=0;y< can->height;y++)
+    for(x=0;x < can->width; x++){
+      trg = XGetPixel(image,x,y) & 0xffffff;
+      if ( trg != 0xffffff ){
+        if ( trg != srcc ){
+          srcc = trg;
+          r = ir = (srcc >>16) & 0xff;
+          r = r/256;
+          g = ig = (srcc >> 8) & 0xff;
+          g = g/256;
+          b = ib = (srcc) & 0xff;
+          b = b/256;
+          if(ir == 0 && ig ==0 && ib == 0)
+            fprintf(Fp,"0 0 0 C\n");
+          else
+            fprintf(Fp,"%f %f %f C\n",r,g,b);
+        }
+        fprintf(Fp,"%d %d R\n", x,can->height - y);
+      }
+    }
+  fprintf(Fp,"0 0 0 C\n");
+  if ( can->noaxis == 0 ){
+    grux = ceil(can->xmax);
+    xwd = can->width/(can->xmax - can->xmin);
+    for(ix=ceil(can->xmin); ix < grux; ix++){
+      t = ix;
+      x = (t - can->xmin) * xwd;
+      PSDashedLine(x, 0, x, can->height);
+    }
+    gruy = ceil(can->ymax);
+    ywd = can->height/(can->ymax - can->ymin);
+    for(iy=can->ymin; iy < gruy; iy++){
+      t = iy;
+      y = (t - can->ymin ) * ywd;
+      PSDashedLine(0, y, can->width, y);
+    }
   }
-  pl->numberOfSegments++;
-  pl->x[n] = x;
-  pl->y[n] = y;
-  return(pl);
+  fprintf(Fp,"[] 0 setdash 0 0 moveto ");
+  fprintf(Fp,"%d 0 lineto %d %d lineto 0 %d lineto closepath stroke\n",
+    can->width, can->width, can->height, can->height);
+  fprintf(Fp,"showpage grestore\n");
+  fflush(Fp);
 }
 
-static int polyLine_lastY(struct polyLine *pl) {
-  int n;
-  n = pl->numberOfSegments;
-  if (n == 0) {
-	polyLine_error("polyLine_lastY:  empty polyLine.");
-  }
-  return(pl->y[n-1]);
-}
+static void PSProlog(int xmin, int ymin,int xmax, int ymax) {
+  int width, height;
+  struct tm *systime;
+  time_t t;
 
-
-static void polyLine_outputPS(struct polyLine *pl) {
-  int n,i;
-  n = pl->numberOfSegments;
-  if (n == 1) {
-	fprintf(Fp," %d %d ifplot_putpixel\n",translateX(pl->x[0]),translateY(pl->y[0]));
-  }else if (n > 1) {
-	fprintf(Fp," newpath ");
-	for (i=0; i<n; i++) {
-	  fprintf(Fp," %d %d ",translateX(pl->x[i]),translateY(pl->y[i]));
-	  if (i==0) fprintf(Fp," moveto ");
-	  else fprintf(Fp," lineto ");
-	}
-	fprintf(Fp," stroke\n");
-  }
-  fflush(Fp);
-}
-static void polyLine_output_bezier_PS(struct polyLine *pl) {
-  int n,i,m;
-  n = pl->numberOfSegments;
-  if (n == 1) {
-	fprintf(Fp," %d %d ifplot_putpixel\n",translateX(pl->x[0]),translateY(pl->y[0]));
-  }else if (n > 1) {
-	fprintf(Fp," newpath ");
-    i=0;
-	m = n-1-((n-1)/3)*3;
-    fprintf(Fp," %d %d ",translateX(pl->x[i]),translateY(pl->y[i]));
-	fprintf(Fp," moveto \n ");
-	for (i=1; i<= ((n-1)/3)*3; i++) {
-	  fprintf(Fp," %d %d ",translateX(pl->x[i]),translateY(pl->y[i]));
-	  if ( i%3 == 0) fprintf(Fp," curveto ");
-	}
-	for (i=n-m; i<n; i++) {
-	  fprintf(Fp," %d %d lineto ",translateX(pl->x[i]),translateY(pl->y[i]));
-	}
-	fprintf(Fp," stroke\n");
-  }
-  fflush(Fp);
-}
-static void polyLine_outputProlog(int xmin, int ymin,int xmax, int ymax) {
-  fprintf(Fp,"%%!PS-Adobe-1.0\n");
-  fprintf(Fp,"%%%%BoundingBox: %d %d %d %d \n",
-		  PRINT_XOFFSET+xmin*MAG,PRINT_YOFFSET+ymin*MAG,
-		  PRINT_XOFFSET+xmax*MAG, PRINT_YOFFSET+ymax*MAG);
+  t = time(NULL);
+  systime = localtime(&t);
+  width = 72*(xmax - xmin-1)/300;
+  height = 72*(ymax - ymin-1)/300;
+  fprintf(Fp,"%%!PS-Adobe-2.0 EPSF-2.0\n");
+  fprintf(Fp,"%%%%BoundingBox: -3 -2 %d %d \n", width*MAG+3, height*MAG+3);
   fprintf(Fp,"%%%%Creator: This is generated by ifplot\n");
   fprintf(Fp,"%%%%Title: ifplot\n");
+  fprintf(Fp,"%%%%CreationDate: %.4d-%.2d-%.2d\n",
+    systime->tm_year+1900,systime->tm_mon+1,systime->tm_mday);
+  fprintf(Fp,"%%%%Pages: 0\n");
   fprintf(Fp,"%%%%EndComments: \n");
-  fprintf(Fp,"0.1 setlinewidth \n");
-  fprintf(Fp,"2 setlinecap \n");
-  fprintf(Fp,"2 setlinejoin \n");
-  fprintf(Fp,"/ifplot_putpixel {  \n");
-  fprintf(Fp,"    /yyy 2 1 roll def /xxx 2 1 roll def \n");
-  fprintf(Fp,"    gsave newpath xxx yyy .5 0 360 arc \n");
-  fprintf(Fp,"    fill grestore \n");
-  fprintf(Fp,"} def \n");
+  fprintf(Fp,"%%%%EndProlog\n%%%%Page: one 1\n%%!\n");
+  fprintf(Fp,"/dpi %d def\n", DPI);
+  fprintf(Fp,"/SCALE 72 dpi div def\n");
+  fprintf(Fp,"/R {moveto 1 0 rlineto 0 1 rlineto -1 0 rlineto closepath");
+  fprintf(Fp," fill} bind def\n");
+  fprintf(Fp,"/C {setrgbcolor} bind def\n");
+  fprintf(Fp,"gsave newpath\n");
+  fprintf(Fp,"0 0 moveto SCALE SCALE scale 1 setlinewidth\n");
   fflush(Fp);
 }
-static void polyLine_outputEpilog(void) {
-  fprintf(Fp,"showpage \n"); fflush(Fp);
-}
-static void polyLine_error(char *s) {
-  fprintf(stderr,"Error in smoothing: %s\n",s);
-  exit(-1);
-}
-static int updatePolyLine(struct polyLine *pl[], int plSize,
-						  int prev[], int curr[],int x) {
-  /*  1 <= y <= Ysize ,  curr[0], curr[Ysize+1]  = 0 ,  i=y .
-	  prev = {-1, -1, -1, 0, 1, -1, 2, -1}  -1: no polyLine. 0:pl[0], 1:pl[1],
-	  curr = {0,  1,  0,  1, 1, 1,  0, 0}    0: no dot, 1: dot.
-  */
-  /* Return Value:  plSize, pl, prev */
-  int i;
-  static int *prevHasNext = NULL;
-  static struct polyLine **tmpPl = NULL;
-  static struct polyLine **newPl = NULL;
-  int tmpPlSize = 0;
-  static int *tmpPrev = NULL;
-  int newPlSize = 0;
-  extern int Ysize;
-  if (prevHasNext == NULL) prevHasNext = (int *)gcmalloc(sizeof(int)*(Ysize+2));
-  if (tmpPl == NULL)
-	tmpPl = (struct polyLine **)gcmalloc(sizeof(struct polyLine *)*(Ysize+2));
-  if (newPl == NULL)
-	newPl = (struct polyLine **)gcmalloc(sizeof(struct polyLine *)*(Ysize+2));
-  if (tmpPrev == NULL) tmpPrev = (int *)gcmalloc(sizeof(int)*(Ysize+2));
-  for (i=0; i<Ysize+2; i++) {
-	prevHasNext[i] = NO; newPl[i] = NULL; tmpPl[i] = NULL;
-	tmpPrev[i] = -1;
-  }
-  prev[0] = prev[Ysize+1] = -1;
-  for (i=1; i <= Ysize; i++) {
-	if (curr[i] IS_DOT) {
-	  if (prev[i-1] IS_POLYLINE && !prevHasNext[i-1]) {
-		pl[prev[i-1]] = polyLine_addNewSegment(pl[prev[i-1]],x,i);
-		prevHasNext[i-1] = YES;
-	  }else if (prev[i] IS_POLYLINE && !prevHasNext[i]) {
-		pl[prev[i]] = polyLine_addNewSegment(pl[prev[i]],x,i);
-		prevHasNext[i] = YES;
-	  }else if (prev[i+1] IS_POLYLINE && !prevHasNext[i+1]) {
-		pl[prev[i+1]] = polyLine_addNewSegment(pl[prev[i+1]],x,i);
-		prevHasNext[i+1] = YES;
-	  }else{
-		newPl[newPlSize] = polyLine_new();
-		newPl[newPlSize] = polyLine_addNewSegment(newPl[newPlSize],x,i);
-		newPlSize++;
-	  }
-	}
-  }
-  for (i=1; i<=Ysize; i++) {
-	if (prevHasNext[i] == NO && prev[i] IS_POLYLINE) {
-	  /* polyLine_outputPS(pl[prev[i]]); */
-	  polyLine_output_bezier_PS(pl[prev[i]]);
-	}
-  }
 
-  /* generate new pl and new prev */
-  tmpPlSize = 0; tmpPrev[0] = tmpPrev[Ysize+1] = -1;
-  for (i = 1; i<=Ysize; i++) {
-	if (prevHasNext[i] == YES) {
-	  tmpPl[tmpPlSize] = pl[prev[i]];
-	  tmpPrev[polyLine_lastY(pl[prev[i]])] = tmpPlSize;
-	  tmpPlSize++;
-	}
-  }
-  for (i=0; i<newPlSize; i++) {
-	tmpPl[tmpPlSize] = newPl[i];
-	tmpPrev[polyLine_lastY(newPl[i])] = tmpPlSize;
-	tmpPlSize++;
-  }
-
-  /* copy and update for the next step.  */
-  for (i=0; i< Ysize+2; i++) {
-	prev[i] = tmpPrev[i];
-  }
-  for (i=0; i<tmpPlSize; i++) {
-	pl[i] = tmpPl[i];
-  }
-  return(tmpPlSize);
-}
-
-/*
-#define TEST_YSIZE 12
-#define TEST_XSIZE 5
-generatePS_from_image(FILE *fp) {
-  int image_test[TEST_XSIZE][TEST_YSIZE] =
-                          {{0,0,1,0,0,0,0,0,0,0,0,0},
-                           {0,1,0,1,0,0,0,0,0,1,0,0},
-                           {0,1,0,1,1,0,0,0,1,0,1,0},
-                           {0,0,1,0,0,0,0,0,0,1,0,1},
-                           {0,0,0,0,0,0,1,0,0,1,0,0}};
-  struct polyLine **pl;
-  int plSize = 0;
-  int *prev;
-  int *curr;
-  int i,x,y;
-
-  Xsize = TEST_XSIZE;
-  Ysize = TEST_YSIZE;
-  pl = (struct polyLine **)gcmalloc(sizeof(struct polyLine *)*(Ysize+2));
-  prev = (int *)gcmalloc(sizeof(int)*(Ysize+2));
-  curr = (int *)gcmalloc(sizeof(int)*(Ysize+2));
-  Fp = fp;
-  polyLine_outputProlog(0,0,Xsize,Ysize);
-  for (i=0; i<= Ysize+1; i++) {
-	prev[i] = -1;
-  }
-  for (x=0; x<Xsize; x++) {
-	curr[0] = curr[Ysize+1] = 0;
-	for (y=0; y<Ysize; y++) {
-	  if (image_test[x][y]) curr[y+1]=1;
-	  else curr[y+1] = 0;
-	}
-	plSize = updatePolyLine(pl,plSize,prev,curr,x);
-  }
-  for (y=0; y<Ysize+2; y++) {
-	curr[y] = 0;
-  }
-  plSize = updatePolyLine(pl,plSize,prev,curr,Xsize);
-  polyLine_outputEpilog();
-}
-main() {
-  generatePS_from_image(stdout);
-}
-
-*/
-static void polyLine_outputPS_dashed_line(int x0,int y0,int x1,int y1) {
+static void PSDashedLine(int x0,int y0,int x1,int y1) {
   extern FILE *Fp;
-  fprintf(Fp," gsave [3] 0 setdash newpath \n");
+  fprintf(Fp,"0.1 setlinewidth [3 3] 0 setdash\n");
   fprintf(Fp," %d %d moveto %d %d lineto stroke \n",x0,y0,x1,y1);
-  fprintf(Fp," stroke grestore \n");
 }
-
-#define D 5
-static int polyLine_pline(can)
-struct canvas *can;
-{
-	double w,w1,k,e,n;
-	int x0,y0,x,y,xadj,yadj;
-	char buf[BUFSIZ];
-	double adjust_scale();
-
-	if ( can->noaxis )
-		return;
-
-	xadj = yadj = 0;
-	if ( (can->xmin < 0) && (can->xmax > 0) ) {
-		x0 = (int)((can->width-1)*(-can->xmin/(can->xmax-can->xmin)));
-		polyLine_outputPS_dashed_line(translateX(x0),translateY(0),
-									  translateX(x0),translateY(can->height));
-	} else if ( can->xmin >= 0 )
-		x0 = 0;
-	else
-		x0 = can->width-1-D;
-	if ( (can->ymin < 0) && (can->ymax > 0) ) {
-		y0 = (int)((can->height-1)*(can->ymax/(can->ymax-can->ymin)));
-		polyLine_outputPS_dashed_line(translateX(0),translateY(y0),
-									  translateX(can->width),translateY(y0));
-	} else if ( can->ymin >= 0 )
-		y0 = can->height-1;
-	else
-		y0 = D;
-	/* BUG:  not written yet a code for PS.
-	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,buf,strlen(buf));
-	}
-	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 )
-			xadj = TEXTWIDTH(sffs,buf,strlen(buf));
-		DRAWSTRING(display,d,scaleGC,x0-xadj,y,buf,strlen(buf));
-	}
-	*/
-}
-
-generatePS_from_image(FILE *fp,XImage *image,int xsize, int ysize,
-					  int color[],int colorSize,
-					  struct canvas *can) {
-  struct polyLine **pl;
-  int plSize = 0;
-  int *prev;
-  int *curr;
-  int i,x,y,c;
-  extern int Strategy_generate_PS;
-
-  Xsize = xsize;
-  Ysize = ysize;
-  pl = (struct polyLine **)gcmalloc(sizeof(struct polyLine *)*(Ysize+2));
-  prev = (int *)gcmalloc(sizeof(int)*(Ysize+2));
-  curr = (int *)gcmalloc(sizeof(int)*(Ysize+2));
-  Fp = fp;
-  polyLine_outputProlog(0,0,Xsize,Ysize);
-  switch(Strategy_generate_PS) {
-  case 1:
-	for (c=0; c<colorSize; c++) {
-	  /* Set color if necessary */
-	  for (i=0; i<= Ysize+1; i++) {
-		prev[i] = -1;
-	  }
-	  for (x=0; x<Xsize; x++) {
-		curr[0] = curr[Ysize+1] = 0;
-		for (y=0; y<Ysize; y++) {
-		  if ((int) XGetPixel(image,x,y) == color[c]) curr[y+1]=1;
-		  else curr[y+1] = 0;
-		}
-		plSize = updatePolyLine(pl,plSize,prev,curr,x);
-	  }
-	  for (y=0; y<Ysize+2; y++) {
-		curr[y] = 0;
-	  }
-	  plSize = updatePolyLine(pl,plSize,prev,curr,Xsize);
-	}
-	break;
-  default:
-	for (c=0; c<colorSize; c++) {
-	  /* Set color if necessary */
-	  for (x=0; x<Xsize; x++) {
-		for (y=0; y<Ysize; y++) {
-		  if ((int) XGetPixel(image,x,y) == color[c]){
-			fprintf(Fp," %d %d ", translateX(x),translateY(y) );
-			fprintf(Fp," ifplot_putpixel\n");
-		  }
-		}
-	  }
-	}
-	break;
-  }
-  polyLine_pline(can);
-  polyLine_outputEpilog();
-}
-
-