[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.13 and 1.22

version 1.13, 2002/08/02 09:29:17 version 1.22, 2011/08/10 04:51:58
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.12 2002/08/02 08:59:47 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/plot/if.c,v 1.21 2006/11/09 15:54:35 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"
   
   #if defined(INTERVAL)
   /* Time message and func*/
   #include <sys/types.h>
   #include <sys/resource.h>
   #include <sys/time.h>
   
   static struct oEGT ltime;
   static double r0;
   double get_rtime();
   #if defined(ITV_TIME_CHECK)
   void tstart()
   {
           get_eg(&ltime);
           r0 = get_rtime();
   }
   
   void tstop(struct canvas *can)
   {
           struct oEGT egt1;
           double e, g, r;
           char ts[100];
           void popdown_warning();
           Widget warnshell,warndialog;
   
           get_eg(&egt1);
           e = egt1.exectime - ltime.exectime;
           g = egt1.gctime - ltime.gctime;
           r = get_rtime() - r0;
           sprintf(ts,"(%8.6f + gc %8.6f) total %8.6f \n",e,g,r);
           create_popup(can->shell,"Message",&ts,&warnshell,&warndialog);
           XawDialogAddButton(warndialog,"OK",popdown_warning,warnshell);
           XtPopup(warnshell,XtGrabNone);
           SetWM_Proto(warnshell);
   }
   #else
   #define tstart()
   #define tstop(a)
   #endif
   #endif
   
 extern JMP_BUF ox_env;  extern JMP_BUF ox_env;
   
 int open_canvas(NODE arg)  int open_canvas(NODE arg)
Line 114  int plot(NODE arg)
Line 154  int plot(NODE arg)
                 can->mode = xrange ? MODE_PLOT : MODE_POLARPLOT;                  can->mode = xrange ? MODE_PLOT : MODE_POLARPLOT;
   
         if ( zrange ) {          if ( zrange ) {
                 n = BDY(zrange); v = VR((P)BDY(n)); n = NEXT(n);                  n = NEXT(BDY(zrange));
                 can->zmin = ToReal(BDY(n));                  can->zmin = ToReal(BDY(n));
                 n = NEXT(n); can->zmax = ToReal(BDY(n));                  n = NEXT(n); can->zmax = ToReal(BDY(n));
                 n = NEXT(n);                  n = NEXT(n);
Line 122  int plot(NODE arg)
Line 162  int plot(NODE arg)
                         can->nzstep = n ? QTOS((Q)BDY(n)) : MAXGC;                          can->nzstep = n ? QTOS((Q)BDY(n)) : MAXGC;
                 else {                  else {
                         /* XXX */                          /* XXX */
                         can->vx = v;                          can->vx = VR((P)BDY(BDY(zrange)));
                         can->nzstep = n ? QTOS((Q)BDY(n)) : DEFAULTPOLARSTEP;                          can->nzstep = n ? QTOS((Q)BDY(n)) : DEFAULTPOLARSTEP;
                 }                  }
         }          }
Line 138  int plot(NODE arg)
Line 178  int plot(NODE arg)
         else          else
                 can->wname = "";                  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);
                   create_canvas(can);
                 plot_print(display,can);                  plot_print(display,can);
         } else if ( can->mode == MODE_POLARPLOT ) {          } else if ( can->mode == MODE_POLARPLOT ) {
                 polarplotcalc(can);                  polarplotcalc(can);
                   create_canvas(can);
                 plot_print(display,can);                  plot_print(display,can);
         } else  
                 ifplotmain(can);  
         copy_to_canvas(can);  
         return id;  
 }  
   
 int polarplot(NODE arg)  
 {  
         int id;  
         NODE n;  
         struct canvas *can;  
         P formula;  
         LIST xrange,yrange,zrange,wsize;  
         STRING wname;  
   
         formula = (P)ARG0(arg);  
         xrange = (LIST)ARG1(arg);  
         wsize = (LIST)ARG4(arg);  
         wname = (STRING)ARG5(arg);  
   
         can = canvas[id = search_canvas()];  
         n = BDY(xrange); can->vx = VR((P)BDY(n)); n = NEXT(n);  
         /* XXX  use zmin, zmax to hold parameter range */  
         can->zmin = ToReal((Num)BDY(n)); n = NEXT(n);  
         can->zmax = ToReal((Num)BDY(n));  
         can->mode = MODE_POLARPLOT;  
         if ( !wsize ) {  
                 can->width = DEFAULTWIDTH; can->height = DEFAULTHEIGHT;  
         } else {          } else {
                 can->width = QTOS((Q)BDY(BDY(wsize)));                  create_canvas(can);
                 can->height = QTOS((Q)BDY(NEXT(BDY(wsize))));                  ifplotmain(can);
         }          }
         if ( wname )  
                 can->wname = BDY(wname);  
         else  
                 can->wname = "";  
         can->formula = formula;  
         polarplotcalc(can);  
         /* xmax, xmin etc are determined after calc */  
         create_canvas(can);  
         plot_print(display,can);  
         copy_to_canvas(can);          copy_to_canvas(can);
         return id;          return id;
 }  }
Line 209  int memory_plot(NODE arg,LIST *bytes)
Line 213  int memory_plot(NODE arg,LIST *bytes)
         zrange = (LIST)ARG3(arg);          zrange = (LIST)ARG3(arg);
         wsize = (LIST)ARG4(arg);          wsize = (LIST)ARG4(arg);
   
           bzero((char *)&tmp_can,sizeof(tmp_can));
         can = &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);
Line 237  int memory_plot(NODE arg,LIST *bytes)
Line 242  int memory_plot(NODE arg,LIST *bytes)
         }          }
         can->wname = "";          can->wname = "";
         can->formula = formula;          can->formula = formula;
         if ( can->mode == MODE_PLOT )          if ( can->mode == MODE_PLOT ) {
                 plotcalc(can);                  plotcalc(can);
         else {                  memory_print(can,&barray);
                   STOQ(can->width,qw); STOQ(can->height,qh);
                   n = mknode(3,qw,qh,barray);
                   MKLIST(*bytes,n);
           } else {
                 width = can->width; height = can->height;                  width = can->width; height = can->height;
                 tabe = (double **)ALLOCA(width*sizeof(double *));                  tabe = (double **)ALLOCA(width*sizeof(double *));
                 for ( i = 0; i < width; i++ )                  for ( i = 0; i < width; i++ )
Line 254  int memory_plot(NODE arg,LIST *bytes)
Line 263  int memory_plot(NODE arg,LIST *bytes)
   
 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((Obj)formula,&vl);          get_vars_recursive((Obj)formula,&vl);
Line 270  int plotover(NODE arg)
Line 279  int plotover(NODE arg)
                 if ( vl0->v->attr == (pointer)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;
 }  }
   
Line 290  int drawcircle(NODE arg)
Line 310  int drawcircle(NODE arg)
         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  #endif
Line 343  int draw_obj(NODE arg)
Line 365  int draw_obj(NODE arg)
                 case 2: /* point */                  case 2: /* point */
                         x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));                          x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));
                         draw_point(display,can,x,y,color);                          draw_point(display,can,x,y,color);
                         MKRVECT2(vect,x,y); MKNODE(n,vect,can->history); can->history = n;                          MKRVECT3(vect,x,y,color); MKNODE(n,vect,can->history);
                           can->history = n;
                         break;                          break;
                 case 3: /* circle */                  case 3: /* circle */
                         x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));                          x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));
                         r = (int)ToReal((Q)ARG2(obj));                          r = (int)ToReal((Q)ARG2(obj));
                         MKRVECT3(vect,x,y,r); MKNODE(n,vect,can->history); can->history = n;                          MKRVECT4(vect,x,y,r,color); MKNODE(n,vect,can->history);
                           can->history = n;
                         break;                          break;
                 case 4: /* line */                  case 4: /* line */
                         x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));                          x = (int)ToReal((Q)ARG0(obj)); y = (int)ToReal((Q)ARG1(obj));
                         u = (int)ToReal((Q)ARG2(obj)); v = (int)ToReal((Q)ARG3(obj));                          u = (int)ToReal((Q)ARG2(obj)); v = (int)ToReal((Q)ARG3(obj));
                         draw_line(display,can,x,y,u,v,color);                          draw_line(display,can,x,y,u,v,color);
                         MKRVECT4(vect,x,y,u,v); MKNODE(n,vect,can->history); can->history = n;                          MKRVECT5(vect,x,y,u,v,color); MKNODE(n,vect,can->history);
                           can->history = n;
                         break;                          break;
                 default:                  default:
                         set_lasterror("draw_obj : invalid request");                          set_lasterror("draw_obj : invalid request");
                         return -1;                          return -1;
         }          }
   #if !defined(VISUAL)
           set_drawcolor(can->color);
   #endif
         return 0;          return 0;
 }  }
   
Line 393  int draw_string(NODE arg)
Line 421  int draw_string(NODE arg)
         x = (int)ToReal((Q)ARG0(pos));          x = (int)ToReal((Q)ARG0(pos));
         y = (int)ToReal((Q)ARG1(pos));          y = (int)ToReal((Q)ARG1(pos));
         draw_character_string(display,can,x,y,str,color);          draw_character_string(display,can,x,y,str,color);
   #if !defined(VISUAL)
           set_drawcolor(can->color);
   #endif
         return 0;          return 0;
 }  }
   
Line 615  void ifplotmain(struct canvas *can)
Line 646  void ifplotmain(struct canvas *can)
         double **tabe;          double **tabe;
         int i;          int i;
   
   #if defined(INTERVAL)
           tstart();
   #endif
         width = can->width; height = can->height;          width = can->width; height = can->height;
         tabe = (double **)ALLOCA(width*sizeof(double *));          tabe = (double **)ALLOCA((width+1)*sizeof(double *));
         for ( i = 0; i < width; i++ )          for ( i = 0; i < width; i++ )
                 tabe[i] = (double *)ALLOCA(height*sizeof(double));                  tabe[i] = (double *)ALLOCA((height+1)*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,0); 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);
   #if defined(INTERVAL)
           tstop(can);
   #endif
 }  }
   
 void qifplotmain(struct canvas *can)  void qifplotmain(struct canvas *can)
Line 643  void qifplotmain(struct canvas *can)
Line 680  void qifplotmain(struct canvas *can)
         reset_selection(); reset_busy(can);          reset_selection(); reset_busy(can);
         define_cursor(can->window,normalcur);          define_cursor(can->window,normalcur);
 }  }
   
   #if defined(INTERVAL)
   int objcp(NODE arg)
   {
           int idsrc, idtrg, op_code;
           struct canvas *cansrc, *cantrg;
   
           idsrc = QTOS((Q)ARG0(arg));
           idtrg = QTOS((Q)ARG1(arg));
           op_code = QTOS((Q)ARG2(arg));
           cansrc = canvas[idsrc];
           cantrg = canvas[idtrg];
           obj_op(cansrc, cantrg, op_code);
           return idsrc;
   }
   
   void obj_op(struct canvas *cansrc, struct canvas *cantrg, int op)
   {
           XImage *imgsrc, *imgtrg;
           int width, height, i, j;
           unsigned long src, trg, black, white;
   
           width = cansrc->width; height = cansrc->height;
           imgsrc = XGetImage(display, cansrc->pix, 0, 0, width, height, -1, ZPixmap);
           imgtrg = XGetImage(display, cantrg->pix, 0, 0, width, height, -1, ZPixmap);
           black=GetColor(display, "black");
           white=GetColor(display, "white");
           flush();
           define_cursor(cantrg->window,runningcur);
           set_busy(cantrg); set_selection();
           cantrg->precise = cansrc->precise;
           cantrg->noaxis = cansrc->noaxis;
           cantrg->noaxisb = cansrc->noaxisb;
           cantrg->vx = cansrc->vx;
           cantrg->vy = cansrc->vy;
           cantrg->formula = cansrc->formula;
           cantrg->width = cansrc->width;
           cantrg->height = cansrc->height;
           cantrg->xmin = cansrc->xmin;
           cantrg->xmax = cansrc->xmax;
           cantrg->ymin = cansrc->ymin;
           cantrg->ymax = cansrc->ymax;
           cantrg->zmin = cansrc->zmin;
           cantrg->zmax = cansrc->zmax;
           cantrg->nzstep = cansrc->nzstep;
           cantrg->qxmin = cansrc->qxmin;
           cantrg->qxmax = cansrc->qxmax;
           cantrg->qymin = cansrc->qymin;
           cantrg->qymax = cansrc->qymax;
           cantrg->pa = cansrc->pa;
           switch (op) {
                   case 1:/* and case */
                           for(i=0;i<width;i++)for(j=0;j<height;j++){
                                   src = XGetPixel(imgsrc,i,j);
                                   trg = XGetPixel(imgtrg,i,j);
                                   if ( (src == black) || (trg == black) )
                                           XPutPixel(imgtrg,i,j,black);
                                   else if ( (src == white) || (trg == white) )
                                           XPutPixel(imgtrg,i,j,white);
                                   else XPutPixel(imgtrg,i,j,(src & trg));
                           }
                           break;
                   case 3:/* copy case */
                           imgtrg->data = imgsrc->data;
                           break;
                   case 6:/* xor case */
                           for(i=0;i<width;i++)for(j=0;j<height;j++){
                                   src = XGetPixel(imgsrc,i,j);
                                   trg = XGetPixel(imgtrg,i,j);
                                   if ( (src == black) || (trg == black) )
                                           XPutPixel(imgtrg,i,j,black);
                                   else if ( (src == white) && (trg == white) )
                                           XPutPixel(imgtrg,i,j,trg|src);
                                   else if ( (src != white) && (trg != white) )
                                           XPutPixel(imgtrg,i,j,white);
                                   else if ( src == white )
                                           XPutPixel(imgtrg,i,j,src);
                           }
                           break;
                   case 7:/* or case */
                           for(i=0;i<width;i++)for(j=0;j<height;j++){
                                   src = XGetPixel(imgsrc,i,j);
                                   trg = XGetPixel(imgtrg,i,j);
                                   if ( (src == black) || (trg == black) )
                                           XPutPixel(imgtrg,i,j,black);
                                   else if (src == white)
                                           XPutPixel(imgtrg,i,j,trg);
                                   else if (trg == white)
                                           XPutPixel(imgtrg,i,j,src);
                           }
                           break;
                   default:
                           break;
           }
           XPutImage(display, cantrg->pix, drawGC, imgtrg, 0, 0, 0, 0, width, height);
           reset_selection(); reset_busy(cantrg);
           define_cursor(cantrg->window,normalcur);
           copy_to_canvas(cantrg);
           count_and_flush();
           flush();
   }
   
   int ineqn(NODE arg)
   {
           int id, op_code, orgcolor;
           struct canvas *can;
           LIST xrange, yrange, geom;
           NODE n;
           STRING wname;
           double **tbl;
   
           can = canvas[id = search_canvas()];
           orgcolor = can->color;
           can->formula = (P)ARG0(arg);
           can->color = QTOS((Q)ARG1(arg));
           xrange = (LIST)ARG2(arg);
           yrange = (LIST)ARG3(arg);
           geom   = (LIST)ARG4(arg);
           wname  = (STRING)ARG5(arg);
           op_code = 3;
           /* set canvas data */
           can->mode = MODE_INEQNP;
           can->width = QTOS((Q)BDY(BDY(geom)));
           can->height = QTOS((Q)BDY(NEXT(BDY(geom))));
           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);
           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);
           can->mode = MODE_INEQNP;
           if ( wname )
                   can->wname = BDY(wname);
           else
                   can->wname = "";
           create_canvas(can);
           ineqnmain(can, orgcolor, op_code);
           return id;
   }
   
   int ineqnover(NODE arg)
   {
           int id;
           struct canvas *can;
           int orgcolor, op_code;
   
           id = QTOS((Q)ARG0(arg));
           can = canvas[id];
           orgcolor = can->color;
           can->formula = (P)ARG1(arg);
           can->color   = QTOS((Q)ARG2(arg));
           op_code      = QTOS((Q)ARG3(arg));
           can->mode    = MODE_INEQNP;
           ineqnmain(can, orgcolor, op_code);
           return id;
   }
   
   void ineqnmain(struct canvas *can, int orgcolor, int op_code)
   {
           int **mask;
           double **tbl;
           int i,j;
   
           current_can = can;
           tbl = (double **)ALLOCA((can->height+1)*sizeof(double *));
           for ( i = 0; i <= can->height; i++ )
                   tbl[i] = (double *)ALLOCA((can->width+1)*sizeof(double));
           mask = (int **)ALLOCA(can->height*sizeof(int *));
           for ( i = 0; i < can->height; i++)
                   mask[i] = (int *)ALLOCA(can->width*sizeof(int));
   
           define_cursor(can->window,runningcur);
           set_busy(can); set_selection();
           ineqncalc(tbl, can, 1);
           for (j = 0; j < can->height; j++){
                   for (i = 0; i < can->width; i++){
                           if ( tbl[j][i] >= 0 ){
                                   if ( (tbl[j+1][i] <= 0 ) ||
                                           (tbl[j][i+1] <= 0) ||
                                           (tbl[j+1][i+1] <= 0) ) mask[j][i] = 0;
                                   else mask[j][i] = 1;
                           } else {
                                   if( (tbl[j+1][i] >= 0) ||
                                           (tbl[j][i+1] >= 0) ||
                                           (tbl[j+1][i+1] >= 0) ) mask[j][i] = 0;
                                   else mask[j][i] = -1;
                           }
                   }
           }
           area_print(display, mask, can, op_code);
           reset_selection();
           reset_busy(can);
           define_cursor(can->window,normalcur);
   #if !defined(VISUAL)
           set_drawcolor(orgcolor);
           can->color = orgcolor;
   #endif
           copy_to_canvas(can);
   }
   
   #if defined(INTERVAL)
   int itvifplot(NODE arg)
   {
           int id, op_code, orgcolor;
           struct canvas *can;
           LIST xrange, yrange, zrange, geom;
           NODE n;
           STRING wname;
           double **tbl;
           int itvsize;
   
           can = canvas[id = search_canvas()];
           orgcolor = can->color;
           can->formula = (P)ARG0(arg);
           xrange = (LIST)ARG1(arg);
           yrange = (LIST)ARG2(arg);
           zrange = (LIST)ARG3(arg);
           geom   = (LIST)ARG4(arg);
           wname  = (STRING)ARG5(arg);
           itvsize = QTOS((Q)ARG6(arg));
           /* set canvas data */
           can->mode = MODE_INEQNP;
           can->width = QTOS((Q)BDY(BDY(geom)));
           can->height = QTOS((Q)BDY(NEXT(BDY(geom))));
           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);
           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);
           can->mode = MODE_INEQNP;
           if ( wname )
                   can->wname = BDY(wname);
           else
                   can->wname = "";
           create_canvas(can);
           itvplotmain(can, itvsize);
   #if !defined(VISUAL)
           set_drawcolor(orgcolor);
           can->color = orgcolor;
   #endif
           copy_to_canvas(can);
           return id;
   }
   
   void itvplotmain(struct canvas *can, int itvsize)
   {
           int **mask;
           double **tbl;
           int i,j;
           int op_code;
           pointer *prp;
   
           tstart(); /* time calc */
           op_code=3;
           current_can = can;
           can->color=0xff00;
           mask = (int **)ALLOCA(can->height*sizeof(int *));
           for ( i = 0; i < can->height; i++){
                   mask[i] = (int *)ALLOCA(can->width*sizeof(int));
                   for (j = 0; j< can->width; j++) mask[i][j] = -1;
           }
           define_cursor(can->window,runningcur);
           set_busy(can); set_selection();
           itvcalc(mask, can, 1, itvsize);
           area_print(display, mask, can, op_code);
           reset_selection();
           reset_busy(can);
           define_cursor(can->window,normalcur);
           tstop(can); /* time calc */
   }
   
   // NORMAL type
   int itvplot1(NODE arg)
   {
           int id, op_code, orgcolor;
           struct canvas *can;
           LIST xrange, yrange, zrange, geom;
           NODE n;
           STRING wname;
           double **tbl;
   
           can = canvas[id = search_canvas()];
           orgcolor = can->color;
           can->formula = (P)ARG0(arg);
           xrange = (LIST)ARG1(arg);
           yrange = (LIST)ARG2(arg);
           zrange = (LIST)ARG3(arg);
           geom   = (LIST)ARG4(arg);
           wname  = (STRING)ARG5(arg);
           /* set canvas data */
           can->mode = MODE_INEQNP;
           can->width = QTOS((Q)BDY(BDY(geom)));
           can->height = QTOS((Q)BDY(NEXT(BDY(geom))));
           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);
           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);
           can->mode = MODE_INEQNP;
           if ( wname )
                   can->wname = BDY(wname);
           else
                   can->wname = "";
           create_canvas(can);
           itvplotmain1(can);
   #if !defined(VISUAL)
           set_drawcolor(orgcolor);
           can->color = orgcolor;
   #endif
           copy_to_canvas(can);
           return id;
   }
   
   void itvplotmain1(struct canvas *can)
   {
           int **mask;
           double **tbl;
           int i,j;
           int op_code;
   
           op_code=3;
           current_can = can;
           mask = (int **)ALLOCA(can->height*sizeof(int *));
           for ( i = 0; i < can->height; i++){
                   mask[i] = (int *)ALLOCA(can->width*sizeof(int));
                   for (j = 0; j< can->width; j++) mask[i][j] = -1;
           }
   
           define_cursor(can->window,runningcur);
           set_busy(can); set_selection();
           itvcalc1(mask, can, 1);
           area_print(display, mask, can, op_code);
           reset_selection();
           reset_busy(can);
           define_cursor(can->window,normalcur);
   }
   
   // TRANSFER type
   int itvplot2(NODE arg)
   {
           int id, op_code, orgcolor;
           struct canvas *can;
           LIST xrange, yrange, zrange, geom;
           NODE n;
           STRING wname;
           double **tbl;
   
           can = canvas[id = search_canvas()];
           orgcolor = can->color;
           can->formula = (P)ARG0(arg);
           xrange = (LIST)ARG1(arg);
           yrange = (LIST)ARG2(arg);
           zrange = (LIST)ARG3(arg);
           geom   = (LIST)ARG4(arg);
           wname  = (STRING)ARG5(arg);
           /* set canvas data */
           can->mode = MODE_INEQNP;
           can->width = QTOS((Q)BDY(BDY(geom)));
           can->height = QTOS((Q)BDY(NEXT(BDY(geom))));
           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);
           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);
           can->mode = MODE_INEQNP;
           if ( wname )
                   can->wname = BDY(wname);
           else
                   can->wname = "";
           create_canvas(can);
           itvplotmain2(can);
   #if !defined(VISUAL)
           set_drawcolor(orgcolor);
           can->color = orgcolor;
   #endif
           copy_to_canvas(can);
           return id;
   }
   
   void itvplotmain2(struct canvas *can)
   {
           int **mask;
           double **tbl;
           int i,j;
           int op_code;
   
           op_code=3;
           current_can = can;
   
           mask = (int **)ALLOCA(can->height*sizeof(int *));
           for ( i = 0; i < can->height; i++){
                   mask[i] = (int *)ALLOCA(can->width*sizeof(int));
                   for (j = 0; j< can->width; j++) mask[i][j] = -1;
           }
   
           define_cursor(can->window,runningcur);
           set_busy(can); set_selection();
           itvcalc2(mask, can, 1);
           area_print(display, mask, can, op_code);
           reset_selection();
           reset_busy(can);
           define_cursor(can->window,normalcur);
   }
   
   // RECURSION type
   int itvplot3(NODE arg)
   {
           int id, op_code, orgcolor;
           struct canvas *can;
           LIST xrange, yrange, zrange, geom;
           NODE n;
           STRING wname;
           double **tbl;
           int itvsize;
   
           can = canvas[id = search_canvas()];
           orgcolor = can->color;
           can->formula = (P)ARG0(arg);
           xrange = (LIST)ARG1(arg);
           yrange = (LIST)ARG2(arg);
           zrange = (LIST)ARG3(arg);
           geom   = (LIST)ARG4(arg);
           wname  = (STRING)ARG5(arg);
           itvsize = QTOS((Q)ARG6(arg));
           /* set canvas data */
           can->mode = MODE_INEQNP;
           can->width = QTOS((Q)BDY(BDY(geom)));
           can->height = QTOS((Q)BDY(NEXT(BDY(geom))));
           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);
           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);
           can->mode = MODE_INEQNP;
           if ( wname )
                   can->wname = BDY(wname);
           else
                   can->wname = "";
           create_canvas(can);
           itvplotmain3(can, itvsize);
   #if !defined(VISUAL)
           set_drawcolor(orgcolor);
           can->color = orgcolor;
   #endif
           copy_to_canvas(can);
           return id;
   }
   
   void itvplotmain3(struct canvas *can, int itvsize)
   {
           int **mask;
           double **tbl;
           int i,j;
           int op_code;
   
           op_code=3;
           current_can = can;
   
           mask = (int **)ALLOCA(can->height*sizeof(int *));
           for ( i = 0; i < can->height; i++)
                   mask[i] = (int *)ALLOCA(can->width*sizeof(int));
   
           define_cursor(can->window,runningcur);
           set_busy(can); set_selection();
           itvcalc3(mask, can, 1, itvsize);
           area_print(display, mask, can, op_code);
           reset_selection();
           reset_busy(can);
           define_cursor(can->window,normalcur);
   }
   
   // RECURSION and TRANSFER type
   int itvplot4(NODE arg)
   {
           int id, op_code, orgcolor;
           struct canvas *can;
           LIST xrange, yrange, zrange, geom;
           NODE n;
           STRING wname;
           double **tbl;
           int itvsize;
   
           can = canvas[id = search_canvas()];
           orgcolor = can->color;
           can->formula = (P)ARG0(arg);
           xrange = (LIST)ARG1(arg);
           yrange = (LIST)ARG2(arg);
           zrange = (LIST)ARG3(arg);
           geom   = (LIST)ARG4(arg);
           wname  = (STRING)ARG5(arg);
           itvsize = QTOS((Q)ARG6(arg));
           /* set canvas data */
           can->mode = MODE_INEQNP;
           can->width = QTOS((Q)BDY(BDY(geom)));
           can->height = QTOS((Q)BDY(NEXT(BDY(geom))));
           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);
           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);
           can->mode = MODE_INEQNP;
           if ( wname )
                   can->wname = BDY(wname);
           else
                   can->wname = "";
           create_canvas(can);
           itvplotmain4(can, itvsize);
   #if !defined(VISUAL)
           set_drawcolor(orgcolor);
           can->color = orgcolor;
   #endif
           copy_to_canvas(can);
           return id;
   }
   
   void itvplotmain4(struct canvas *can, int itvsize)
   {
           int **mask;
           double **tbl;
           int i,j;
           int op_code;
   
           tstart();/* time calc */
           op_code=3;
           current_can = can;
   
           mask = (int **)ALLOCA(can->height*sizeof(int *));
           for ( i = 0; i < can->height; i++){
                   mask[i] = (int *)ALLOCA(can->width*sizeof(int));
                   for (j = 0; j< can->width; j++) mask[i][j] = -1;
           }
   
           define_cursor(can->window,runningcur);
           set_busy(can); set_selection();
           itvcalc4(mask, can, 1, itvsize);
           area_print(display, mask, can, op_code);
           reset_selection();
           reset_busy(can);
           define_cursor(can->window,normalcur);
           tstop(can); /* time calc */
   }
   #endif
   #endif

Legend:
Removed from v.1.13  
changed lines
  Added in v.1.22

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>