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

Diff for /OpenXM_contrib2/asir2000/plot/if.c between version 1.2 and 1.21

version 1.2, 2000/08/21 08:31:50 version 1.21, 2006/11/09 15:54:35
Line 23 
Line 23 
  * shall be made on your publication or presentation in any form of the   * shall be made on your publication or presentation in any form of the
  * results obtained by use of the SOFTWARE.   * results obtained by use of the SOFTWARE.
  * (4) In the event that you modify the SOFTWARE, you shall notify FLL by   * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
  * e-mail at risa-admin@flab.fujitsu.co.jp of the detailed specification   * 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   * for such modification or the source code of the modified part of the
  * SOFTWARE.   * SOFTWARE.
  *   *
Line 45 
Line 45 
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *   *
  * $OpenXM: OpenXM_contrib2/asir2000/plot/if.c,v 1.1.1.1 1999/12/03 07:39:13 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/plot/if.c,v 1.20 2006/11/09 03:41:47 saito 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"
   
 extern jmp_buf ox_env;  extern JMP_BUF ox_env;
   
   int open_canvas(NODE arg)
   {
           int id;
           struct canvas *can;
           LIST wsize;
           STRING wname;
   
           wsize = (LIST)ARG0(arg);
           wname = (STRING)ARG1(arg);
   
           can = canvas[id = search_canvas()];
           can->mode = MODE_INTERACTIVE;
           if ( !wsize ) {
                   can->width = DEFAULTWIDTH; can->height = DEFAULTHEIGHT;
           } else {
                   can->width = QTOS((Q)BDY(BDY(wsize)));
                   can->height = QTOS((Q)BDY(NEXT(BDY(wsize))));
           }
           if ( wname )
                   can->wname = BDY(wname);
           else
                   can->wname = "";
           create_canvas(can);
           return id;
   }
   
 int plot(NODE arg)  int plot(NODE arg)
 {  {
         int id;          int id;
Line 62  int plot(NODE arg)
Line 88  int plot(NODE arg)
         P formula;          P formula;
         LIST xrange,yrange,zrange,wsize;          LIST xrange,yrange,zrange,wsize;
         STRING wname;          STRING wname;
           V v;
   
         formula = (P)ARG0(arg);          formula = (P)ARG0(arg);
         xrange = (LIST)ARG1(arg);          xrange = (LIST)ARG1(arg);
Line 71  int plot(NODE arg)
Line 98  int plot(NODE arg)
         wname = (STRING)ARG5(arg);          wname = (STRING)ARG5(arg);
   
         can = canvas[id = search_canvas()];          can = canvas[id = search_canvas()];
           if ( xrange ) {
                   n = BDY(xrange); can->vx = VR((P)BDY(n)); n = NEXT(n);
                   can->qxmin = (Q)BDY(n); n = NEXT(n); can->qxmax = (Q)BDY(n);
                   can->xmin = ToReal(can->qxmin); can->xmax = ToReal(can->qxmax);
           }
           if ( yrange ) {
                   n = BDY(yrange); can->vy = VR((P)BDY(n)); n = NEXT(n);
                   can->qymin = (Q)BDY(n); n = NEXT(n); can->qymax = (Q)BDY(n);
                   can->ymin = ToReal(can->qymin); can->ymax = ToReal(can->qymax);
           }
           if ( xrange && yrange )
                   can->mode = zrange ? MODE_CONPLOT : MODE_IFPLOT;
           else
                   can->mode = xrange ? MODE_PLOT : MODE_POLARPLOT;
   
           if ( zrange ) {
                   n = NEXT(BDY(zrange));
                   can->zmin = ToReal(BDY(n));
                   n = NEXT(n); can->zmax = ToReal(BDY(n));
                   n = NEXT(n);
                   if ( can->mode == MODE_CONPLOT )
                           can->nzstep = n ? QTOS((Q)BDY(n)) : MAXGC;
                   else {
                           /* XXX */
                           can->vx = VR((P)BDY(BDY(zrange)));
                           can->nzstep = n ? QTOS((Q)BDY(n)) : DEFAULTPOLARSTEP;
                   }
           }
   
           if ( !wsize ) {
                   can->width = DEFAULTWIDTH; can->height = DEFAULTHEIGHT;
           } else {
                   can->width = QTOS((Q)BDY(BDY(wsize)));
                   can->height = QTOS((Q)BDY(NEXT(BDY(wsize))));
           }
           if ( wname )
                   can->wname = BDY(wname);
           else
                   can->wname = "";
           can->formula = formula;
           if ( can->mode == MODE_PLOT ) {
                   plotcalc(can);
                   create_canvas(can);
                   plot_print(display,can);
           } else if ( can->mode == MODE_POLARPLOT ) {
                   polarplotcalc(can);
                   create_canvas(can);
                   plot_print(display,can);
           } else {
                   create_canvas(can);
                   ifplotmain(can);
           }
           copy_to_canvas(can);
           return id;
   }
   
   int memory_plot(NODE arg,LIST *bytes)
   {
           NODE n;
           struct canvas tmp_can;
           struct canvas *can;
           P formula;
           LIST xrange,yrange,zrange,wsize;
           int width,height;
           double **tabe;
           int i;
           BYTEARRAY barray;
           Q qw,qh;
   
           formula = (P)ARG0(arg);
           xrange = (LIST)ARG1(arg);
           yrange = (LIST)ARG2(arg);
           zrange = (LIST)ARG3(arg);
           wsize = (LIST)ARG4(arg);
   
           bzero((char *)&tmp_can,sizeof(tmp_can));
           can = &tmp_can;
         n = BDY(xrange); can->vx = VR((P)BDY(n)); n = NEXT(n);          n = BDY(xrange); can->vx = VR((P)BDY(n)); n = NEXT(n);
         can->qxmin = (Q)BDY(n); n = NEXT(n); can->qxmax = (Q)BDY(n);          can->qxmin = (Q)BDY(n); n = NEXT(n); can->qxmax = (Q)BDY(n);
         can->xmin = ToReal(can->qxmin); can->xmax = ToReal(can->qxmax);          can->xmin = ToReal(can->qxmin); can->xmax = ToReal(can->qxmax);
Line 96  int plot(NODE arg)
Line 200  int plot(NODE arg)
                 can->width = QTOS((Q)BDY(BDY(wsize)));                  can->width = QTOS((Q)BDY(BDY(wsize)));
                 can->height = QTOS((Q)BDY(NEXT(BDY(wsize))));                  can->height = QTOS((Q)BDY(NEXT(BDY(wsize))));
         }          }
         if ( wname )          can->wname = "";
                 can->wname = BDY(wname);  
         else  
                 can->wname = "";  
         can->formula = formula;          can->formula = formula;
         create_canvas(can);  
         if ( can->mode == MODE_PLOT ) {          if ( can->mode == MODE_PLOT ) {
                 plotcalc(can);                  plotcalc(can);
                 plot_print(display,can);                  memory_print(can,&barray);
         } else                  STOQ(can->width,qw); STOQ(can->height,qh);
                 ifplotmain(can);                  n = mknode(3,qw,qh,barray);
         copy_to_canvas(can);                  MKLIST(*bytes,n);
         return id;          } else {
                   width = can->width; height = can->height;
                   tabe = (double **)ALLOCA(width*sizeof(double *));
                   for ( i = 0; i < width; i++ )
                           tabe[i] = (double *)ALLOCA(height*sizeof(double));
                   calc(tabe,can,1);
                   memory_if_print(tabe,can,&barray);
                   STOQ(width,qw); STOQ(height,qh);
                   n = mknode(3,qw,qh,barray);
                   MKLIST(*bytes,n);
           }
 }  }
   
 int plotover(NODE arg)  int plotover(NODE arg)
 {  {
         int index;          int index, color;
         P formula;          P formula;
         struct canvas *can;          struct canvas *can;
         struct canvas fakecan;  
         VL vl,vl0;          VL vl,vl0;
   
         index = QTOS((Q)ARG0(arg));          index = QTOS((Q)ARG0(arg));
         formula = (P)ARG1(arg);          formula = (P)ARG1(arg);
         can = canvas[index];          can = canvas[index];
           color = can->color;
         if ( !can->window )          if ( !can->window )
                 return -1;                  return -1;
         get_vars_recursive(formula,&vl);          get_vars_recursive((Obj)formula,&vl);
         for ( vl0 = vl; vl0; vl0 = NEXT(vl0) )          for ( vl0 = vl; vl0; vl0 = NEXT(vl0) )
                 if ( vl0->v->attr == V_IND )                  if ( vl0->v->attr == (pointer)V_IND )
                         if ( vl->v != can->vx && vl->v != can->vy )                          if ( vl->v != can->vx && vl->v != can->vy )
                                 return -1;                                  return -1;
           if ( argc(arg) == 3 )
                   can->color = QTOS((Q)ARG2(arg));
           else
                   can->color = 0;
   #if !defined(VISUAL)
           set_drawcolor(can->color);
   #endif
         current_can = can;          current_can = can;
         fakecan = *can; fakecan.formula = formula;          can->formula = formula;
         if ( can->mode == MODE_PLOT ) {          if ( can->mode == MODE_PLOT ) {
                 plotcalc(&fakecan);                  plotcalc(can);
                 plot_print(display,&fakecan);                  plot_print(display,can);
         } else          } else
                 ifplotmain(&fakecan);                  ifplotmain(can);
         copy_to_canvas(&fakecan);          copy_to_canvas(can);
           can->color = color;
   #if !defined(VISUAL)
           set_drawcolor(can->color);
   #endif
         return index;          return index;
 }  }
   
 int drawcircle(NODE arg)  int drawcircle(NODE arg)
 {  {
   #if !defined(VISUAL)
         int id;          int id;
         int index;          int index;
         pointer ptr;          pointer ptr;
         Q ret;          Q ret;
         LIST xyr;          LIST xyr;
         Obj x,y,r;          Obj x,y,r;
         int wx,wy,wr;          int wx,wy,wr,c;
         struct canvas *can;          struct canvas *can;
         struct canvas fakecan;  
   
         index = QTOS((Q)ARG0(arg));          index = QTOS((Q)ARG0(arg));
         xyr = (LIST)ARG1(arg);          xyr = (LIST)ARG1(arg);
         x = (Obj)ARG0(BDY(xyr)); y = (Obj)ARG1(BDY(xyr)); r = (Obj)ARG2(BDY(xyr));          x = (Obj)ARG0(BDY(xyr)); y = (Obj)ARG1(BDY(xyr)); r = (Obj)ARG2(BDY(xyr));
           c = QTOS((Q)ARG2(arg));
         can = canvas[index];          can = canvas[index];
         if ( !can->window )          if ( !can->window )
                 return -1;                  return -1;
         else {          else {
                 current_can = can;                  current_can = can;
                   set_drawcolor(c);
                 wx = (ToReal(x)-can->xmin)*can->width/(can->xmax-can->xmin);                  wx = (ToReal(x)-can->xmin)*can->width/(can->xmax-can->xmin);
                 wy = (can->ymax-ToReal(y))*can->height/(can->ymax-can->ymin);                  wy = (can->ymax-ToReal(y))*can->height/(can->ymax-can->ymin);
                 wr = ToReal(r);                  wr = ToReal(r);
                 XFillArc(display,can->pix,colorGC,wx-wr/2,wy-wr/2,wr,wr,0,360*64);                  XFillArc(display,can->pix,cdrawGC,wx-wr/2,wy-wr/2,wr,wr,0,360*64);
                 copy_to_canvas(can);                  copy_to_canvas(can);
                   set_drawcolor(can->color);
                 return index;                  return index;
         }          }
   #endif
 }  }
   
   int draw_obj(NODE arg)
   {
           int index;
           int x,y,u,v,len,r;
           NODE obj,n;
           RealVect *vect;
           struct canvas *can;
           int color;
   
           index = QTOS((Q)ARG0(arg));
           can = canvas[index];
           if ( !can && closed_canvas[index] ) {
                   canvas[index] = closed_canvas[index];
                   closed_canvas[index] = 0;
                   can = canvas[index];
                   popup_canvas(index);
                   current_can = can;
           } else if ( !can || (can && !can->window) ) {
                   set_lasterror("draw_obj : canvas does not exist");
                   return -1;
           }
   
           obj = BDY((LIST)ARG1(arg));
           if ( argc(arg) == 3 )
                   color = QTOS((Q)ARG2(arg));
           else
                   color = 0; /* black */
           switch ( len = length(obj) ) {
                   case 2: /* point */
                           x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));
                           draw_point(display,can,x,y,color);
                           MKRVECT3(vect,x,y,color); MKNODE(n,vect,can->history);
                           can->history = n;
                           break;
                   case 3: /* circle */
                           x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));
                           r = (int)ToReal((Q)ARG2(obj));
                           MKRVECT4(vect,x,y,r,color); MKNODE(n,vect,can->history);
                           can->history = n;
                           break;
                   case 4: /* line */
                           x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));
                           u = (int)ToReal((Q)ARG2(obj)); v = (int)ToReal((Q)ARG3(obj));
                           draw_line(display,can,x,y,u,v,color);
                           MKRVECT5(vect,x,y,u,v,color); MKNODE(n,vect,can->history);
                           can->history = n;
                           break;
                   default:
                           set_lasterror("draw_obj : invalid request");
                           return -1;
           }
   #if !defined(VISUAL)
           set_drawcolor(can->color);
   #endif
           return 0;
   }
   
   int draw_string(NODE arg)
   {
           int index,x,y;
           char *str;
           NODE pos;
           struct canvas *can;
           int color;
   
           index = QTOS((Q)ARG0(arg));
           can = canvas[index];
           if ( !can && closed_canvas[index] ) {
                   canvas[index] = closed_canvas[index];
                   closed_canvas[index] = 0;
                   can = canvas[index];
                   popup_canvas(index);
                   current_can = can;
           } else if ( !can || (can && !can->window) ) {
                   set_lasterror("draw_obj : canvas does not exist");
                   return -1;
           }
   
           pos = BDY((LIST)ARG1(arg));
           str = BDY((STRING)ARG2(arg));
           if ( argc(arg) == 4 )
                   color = QTOS((Q)ARG3(arg));
           else
                   color = 0; /* black */
           x = (int)ToReal((Q)ARG0(pos));
           y = (int)ToReal((Q)ARG1(pos));
           draw_character_string(display,can,x,y,str,color);
   #if !defined(VISUAL)
           set_drawcolor(can->color);
   #endif
           return 0;
   }
   
   int clear_canvas(NODE arg)
   {
           int index;
           struct canvas *can;
   
           index = QTOS((Q)ARG0(arg));
           can = canvas[index];
           if ( !can || !can->window )
                   return -1;
           clear_pixmap(can);
           copy_to_canvas(can);
           /* clear the history */
           can->history = 0;
   }
   
 #define RealtoDbl(r) ((r)?BDY(r):0.0)  #define RealtoDbl(r) ((r)?BDY(r):0.0)
   
 int arrayplot(NODE arg)  int arrayplot(NODE arg)
Line 178  int arrayplot(NODE arg)
Line 411  int arrayplot(NODE arg)
         LIST xrange,wsize;          LIST xrange,wsize;
         char *wname;          char *wname;
         NODE n;          NODE n;
         Q ret;  
         double ymax,ymin,dy,xstep;          double ymax,ymin,dy,xstep;
         Real *tab;          Real *tab;
         struct canvas *can;          struct canvas *can;
Line 230  int arrayplot(NODE arg)
Line 462  int arrayplot(NODE arg)
                 else if ( t < -MAXSHORT )                  else if ( t < -MAXSHORT )
                         pa[ix].y = -MAXSHORT;                          pa[ix].y = -MAXSHORT;
                 else                  else
                         pa[ix].y = t;                          pa[ix].y = (long)t;
         }          }
         plot_print(display,can);          plot_print(display,can);
         copy_to_canvas(can);          copy_to_canvas(can);
         return id;          return id;
 }  }
   
 ifplot_resize(can,spos,epos)  void ifplot_resize(struct canvas *can,POINT spos,POINT epos)
 struct canvas *can;  
 POINT spos,epos;  
 {  {
         struct canvas *ncan;          struct canvas *ncan;
         struct canvas fakecan;          struct canvas fakecan;
Line 248  POINT spos,epos;
Line 478  POINT spos,epos;
         Q s,t;          Q s,t;
         int new;          int new;
         int w,h,m;          int w,h,m;
         Q ret;  
   
         if ( XC(spos) < XC(epos) && YC(spos) < YC(epos) ) {          if ( XC(spos) < XC(epos) && YC(spos) < YC(epos) ) {
                 if ( can->precise && !can->wide ) {                  if ( can->precise && !can->wide ) {
Line 297  POINT spos,epos;
Line 526  POINT spos,epos;
                 if ( can->precise && !can->wide ) {                  if ( can->precise && !can->wide ) {
                         current_can = can;                          current_can = can;
                         alloc_pixmap(ncan);                          alloc_pixmap(ncan);
   #if defined(VISUAL)
                           ncan->real_can = can;
   #endif
                         qifplotmain(ncan);                          qifplotmain(ncan);
                         copy_subimage(ncan,can,spos);                          copy_subimage(ncan,can,spos);
                         copy_to_canvas(can);                          copy_to_canvas(can);
Line 311  POINT spos,epos;
Line 543  POINT spos,epos;
         }          }
 }  }
   
 plot_resize(can,spos,epos)  void plot_resize(struct canvas *can,POINT spos,POINT epos)
 struct canvas *can;  
 POINT spos,epos;  
 {  {
         struct canvas *ncan;          struct canvas *ncan;
         Q dx,dx2,xmin,xmax,xmid;          Q dx,dx2,xmin,xmax,xmid;
         double dy,dy2,ymin,ymax,ymid;          double dy,ymin,ymax,ymid;
         Q sx,ex,cw,ten,two;          Q sx,ex,cw,ten,two;
         double sy,ey;  
         Q s,t;          Q s,t;
         int new;          int new;
         int w,h,m;          int w,h,m;
Line 371  POINT spos,epos;
Line 600  POINT spos,epos;
         }          }
 }  }
   
 ifplotmain(can)  void ifplotmain(struct canvas *can)
 struct canvas *can;  
 {  {
         int width,height;          int width,height;
         double **tabe,*tabeb;          double **tabe;
         int i;          int i;
   
         width = can->width; height = can->height;          width = can->width; height = can->height;
Line 384  struct canvas *can;
Line 612  struct canvas *can;
                 tabe[i] = (double *)ALLOCA(height*sizeof(double));                  tabe[i] = (double *)ALLOCA(height*sizeof(double));
         define_cursor(can->window,runningcur);          define_cursor(can->window,runningcur);
         set_busy(can); set_selection();          set_busy(can); set_selection();
         calc(tabe,can); if_print(display,tabe,can);          calc(tabe,can,0); if_print(display,tabe,can);
         reset_selection(); reset_busy(can);          reset_selection(); reset_busy(can);
         define_cursor(can->window,normalcur);          define_cursor(can->window,normalcur);
 }  }
   
 qifplotmain(can)  void qifplotmain(struct canvas *can)
 struct canvas *can;  
 {  {
         int width,height;          int width,height;
         char **tabe,*tabeb;          char **tabe,*tabeb;
Line 398  struct canvas *can;
Line 625  struct canvas *can;
   
         width = can->width; height = can->height;          width = can->width; height = can->height;
         tabe = (char **)ALLOCA(width*sizeof(char *)+width*height*sizeof(char));          tabe = (char **)ALLOCA(width*sizeof(char *)+width*height*sizeof(char));
         bzero(tabe,width*sizeof(char *)+width*height*sizeof(char));          bzero((void *)tabe,width*sizeof(char *)+width*height*sizeof(char));
         for ( i = 0, tabeb = (char *)(tabe+width); i < width; i++ )          for ( i = 0, tabeb = (char *)(tabe+width); i < width; i++ )
                 tabe[i] = tabeb + height*i;                  tabe[i] = tabeb + height*i;
         define_cursor(can->window,runningcur);          define_cursor(can->window,runningcur);
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.21
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>