[BACK]Return to plotp.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / plot

Diff for /OpenXM_contrib2/asir2000/plot/plotp.c between version 1.1 and 1.9

version 1.1, 1999/12/03 07:39:13 version 1.9, 2001/10/09 01:36:28
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/asir99/plot/plotp.c,v 1.1.1.1 1999/11/10 08:12:34 noro Exp $ */  /*
    * 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.8 2001/08/22 09:19:21 noro Exp $
   */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
 #include "ox.h"  #include "ox.h"
 #include "ifplot.h"  #include "ifplot.h"
   #include <math.h>
   
 #if defined(sun) && !defined(__svr4__)  #if defined(sun) && !defined(__svr4__)
 #define EXP10(a) exp10(a)  #define EXP10(a) exp10(a)
Line 10 
Line 59 
 #define EXP10(a) pow(10.0,a)  #define EXP10(a) pow(10.0,a)
 #endif  #endif
   
 if_print(display,tab,can)  #if defined(VISUAL)
 DISPLAY *display;  static POINT oldpos;
 double **tab;  #endif
 struct canvas *can;  
   void if_print(DISPLAY *display,double **tab,struct canvas *can)
 {  {
         int ix,iy,width,height;          int ix,iy,width,height;
         double *px,*px1,*px2;          double *px,*px1,*px2;
Line 35  struct canvas *can;
Line 85  struct canvas *can;
         flush();          flush();
 }  }
   
 con_print(display,tab,can)  #define MEMORY_DRAWPOINT(a,len,x,y) (((a)[(len)*(y)+((x)>>3)]) |= (1<<((x)&7)))
 DISPLAY *display;  
 double **tab;  void memory_if_print(double **tab,struct canvas *can,BYTEARRAY *bytes)
 struct canvas *can;  
 {  {
           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;          int ix,iy,iz,width,height,pas,pai,len;
         double zstep,z;          double zstep,z;
         double *px,*px1,*px2;          double *px,*px1,*px2;
Line 85  struct canvas *can;
Line 159  struct canvas *can;
         flush();          flush();
 }  }
   
 qif_print(display,tab,can)  void memory_con_print(double **tab,struct canvas *can,BYTEARRAY *bytes)
 DISPLAY *display;  
 char **tab;  
 struct canvas *can;  
 {  {
           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 qif_print(DISPLAY *display,char **tab,struct canvas *can)
   {
         int ix,iy,width,height;          int ix,iy,width,height;
         char *px;          char *px;
         DRAWABLE pix;          DRAWABLE pix;
Line 105  struct canvas *can;
Line 222  struct canvas *can;
         flush();          flush();
 }  }
   
 plot_print(display,can)  void plot_print(DISPLAY *display,struct canvas *can)
 DISPLAY *display;  
 struct canvas *can;  
 {  {
         int len,i;          int len;
         POINT *pa;          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,          XDrawLines(display,can->pix,drawGC,
                 can->pa[0].pos,can->pa[0].length,CoordModeOrigin);                  can->pa[0].pos,can->pa[0].length,CoordModeOrigin);
         XFlush(display);          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
           if ( color ) {
                   set_drawcolor(color);
                   DRAWPOINT(display,can->pix,cdrawGC,x,y);
                   DRAWPOINT(display,can->window,cdrawGC,x,y);
           } else {
                   DRAWPOINT(display,can->pix,drawGC,x,y);
                   DRAWPOINT(display,can->window,drawGC,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
           if ( color ) {
                   set_drawcolor(color);
                   DRAWLINE(display,can->pix,cdrawGC,x,y,u,v);
                   DRAWLINE(display,can->window,cdrawGC,x,y,u,v);
           } else {
                   DRAWLINE(display,can->pix,drawGC,x,y,u,v);
                   DRAWLINE(display,can->window,drawGC,x,y,u,v);
           }
           XFlush(display);
   #endif
   }
   
 #define D 5  #define D 5
   
 pline(display,can,d)  void pline(DISPLAY *display,struct canvas *can,DRAWABLE d)
 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;          int x0,y0,x,y,xadj,yadj;
         char buf[BUFSIZ];          char buf[BUFSIZ];
         double adjust_scale();          double adjust_scale();
   
         if ( can->noaxis )          /* XXX : should be cleaned up */
           if ( can->noaxis || (can->mode == MODE_PLOT && !can->pa) )
                 return;                  return;
           if ( can->mode == MODE_INTERACTIVE )
                   return;
   
         xadj = yadj = 0;          xadj = yadj = 0;
         if ( (can->xmin < 0) && (can->xmax > 0) ) {          if ( (can->xmin < 0) && (can->xmax > 0) ) {
                 x0 = (int)((can->width-1)*(-can->xmin/(can->xmax-can->xmin)));                  x0 = (int)((can->width)*(-can->xmin/(can->xmax-can->xmin)));
                 DRAWLINE(display,d,dashGC,x0,0,x0,can->height);                  DRAWLINE(display,d,dashGC,x0,0,x0,can->height);
         } else if ( can->xmin >= 0 )          } else if ( can->xmin >= 0 )
                 x0 = 0;                  x0 = 0;
         else          else
                 x0 = can->width-1-D;                  x0 = can->width-D;
         if ( (can->ymin < 0) && (can->ymax > 0) ) {          if ( (can->ymin < 0) && (can->ymax > 0) ) {
                 y0 = (int)((can->height-1)*(can->ymax/(can->ymax-can->ymin)));                  y0 = (int)((can->height)*(can->ymax/(can->ymax-can->ymin)));
                 DRAWLINE(display,d,dashGC,0,y0,can->width,y0);                  DRAWLINE(display,d,dashGC,0,y0,can->width,y0);
         } else if ( can->ymin >= 0 )          } else if ( can->ymin >= 0 )
                 y0 = can->height-1;                  y0 = can->height;
         else          else
                 y0 = D;                  y0 = D;
   
           /* scale on x-axis */
         w = can->xmax-can->xmin;          w = can->xmax-can->xmin;
         w1 = w * DEFAULTWIDTH/can->width;          w1 = w * DEFAULTWIDTH/can->width;
         e = adjust_scale(EXP10(floor(log10(w1))),w1);          e = adjust_scale(EXP10(floor(log10(w1))),w1);
         for ( n = ceil(can->xmin/e); n*e<= can->xmax; n++ ) {          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);                  DRAWLINE(display,d,drawGC,x,y0,x,y0-D);
                 sprintf(buf,"%g",n*e);                  sprintf(buf,"%g",n*e);
                 DRAWSTRING(display,d,scaleGC,x+2,y0,buf,strlen(buf));                  DRAWSTRING(display,d,scaleGC,x+2,y0+2,buf,strlen(buf));
         }          }
   
           /* scale on y-axis */
         w = can->ymax-can->ymin;          w = can->ymax-can->ymin;
         w1 = w * DEFAULTHEIGHT/can->height;          w1 = w * DEFAULTHEIGHT/can->height;
         e = adjust_scale(EXP10(floor(log10(w1))),w1);          e = adjust_scale(EXP10(floor(log10(w1))),w1);
         for ( n = ceil(can->ymin/e); n*e<= can->ymax; n++ ) {          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);                  DRAWLINE(display,d,drawGC,x0,y,x0+D,y);
                 sprintf(buf,"%g",n*e);                  sprintf(buf,"%g",n*e);
                 if ( can->xmax <= 0 )                  if ( can->xmax <= 0 ) {
   #if !defined(VISUAL)
                         xadj = TEXTWIDTH(sffs,buf,strlen(buf));                          xadj = TEXTWIDTH(sffs,buf,strlen(buf));
                 DRAWSTRING(display,d,scaleGC,x0-xadj,y,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(e,w)  double adjust_scale(double e,double w)
 double e,w;  
 {  {
         switch ( (int)floor(w/e) ) {          switch ( (int)floor(w/e) ) {
                 case 1:                  case 1:
Line 185  double e,w;
Line 381  double e,w;
         }          }
 }  }
   
 initmarker(can,message)  void initmarker(struct canvas *can,char *message)
 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->xdone,0.0,0.0);
         XawScrollbarSetThumb(can->ydone,1.0,0.0);          XawScrollbarSetThumb(can->ydone,1.0,0.0);
         XFlush(display);          XFlush(display);
   #endif
 }  }
   
 marker(can,dir,p)  void marker(struct canvas *can,int dir,int p)
 struct canvas *can;  
 int dir,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 ) {          if ( dir == DIR_X ) {
                 XawScrollbarSetThumb(can->xdone,(float)p/(float)can->width,0.05);                  XawScrollbarSetThumb(can->xdone,(float)p/(float)can->width,0.05);
                 count_and_flush();                  count_and_flush();
Line 208  int dir,p;
Line 413  int dir,p;
                 XawScrollbarSetThumb(can->ydone,1.0-(float)p/(float)can->nzstep,0.05);                  XawScrollbarSetThumb(can->ydone,1.0-(float)p/(float)can->nzstep,0.05);
                 flush();                  flush();
         }          }
   #endif
 }  }
   
 define_cursor(w,cur)  void define_cursor(WINDOW w,CURSOR cur)
 WINDOW w;  
 CURSOR cur;  
 {  {
   #if !defined(VISUAL)
         XDefineCursor(display,w,cur); flush();          XDefineCursor(display,w,cur); flush();
   #endif
 }  }
   
 static int flush_count;  static int flush_count;
Line 224  static int flush_count;
Line 430  static int flush_count;
 #define MAX_COUNT 32  #define MAX_COUNT 32
 #endif  #endif
   
 count_and_flush() {  void count_and_flush() {
   #if !defined(VISUAL)
         if ( ++flush_count == MAX_COUNT )          if ( ++flush_count == MAX_COUNT )
                 flush();                  flush();
   #endif
 }  }
   
 flush() {  void flush() {
   #if !defined(VISUAL)
         flush_count = 0;          flush_count = 0;
         XFlush(display);          XFlush(display);
   #endif
 }  }
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>