=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/plot/if.c,v retrieving revision 1.23 retrieving revision 1.24 diff -u -p -r1.23 -r1.24 --- OpenXM_contrib2/asir2000/plot/if.c 2011/08/11 06:25:06 1.23 +++ OpenXM_contrib2/asir2000/plot/if.c 2013/12/19 05:48:24 1.24 @@ -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/if.c,v 1.22 2011/08/10 04:51:58 saito Exp $ + * $OpenXM: OpenXM_contrib2/asir2000/plot/if.c,v 1.23 2011/08/11 06:25:06 saito Exp $ */ #include "ca.h" #include "parse.h" @@ -53,7 +53,7 @@ #include "ifplot.h" #if defined(INTERVAL) -/* Time message and func*/ +// Time message and func #include #include #include @@ -65,7 +65,7 @@ double get_rtime(); void tstart() { get_eg(<ime); - r0 = get_rtime(); + r0=get_rtime(); } void tstop(struct canvas *can) @@ -77,9 +77,9 @@ void tstop(struct canvas *can) Widget warnshell,warndialog; get_eg(&egt1); - e = egt1.exectime - ltime.exectime; - g = egt1.gctime - ltime.gctime; - r = get_rtime() - r0; + 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); @@ -94,34 +94,31 @@ void tstop(struct canvas *can) extern JMP_BUF ox_env; -int open_canvas(NODE arg) -{ +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; + 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)))); + can->width=QTOS((Q)BDY(BDY(wsize))); + can->height=QTOS((Q)BDY(NEXT(BDY(wsize)))); } - if ( wname ) - can->wname = BDY(wname); - else - can->wname = ""; + if(wname)can->wname=BDY(wname); + else can->wname=""; create_canvas(can); return id; } -int plot(NODE arg) -{ + +int plot(NODE arg,int fn){ int id; NODE n; struct canvas *can; @@ -130,72 +127,79 @@ int plot(NODE arg) STRING wname; V v; - formula = (P)ARG0(arg); - xrange = (LIST)ARG1(arg); - yrange = (LIST)ARG2(arg); - zrange = (LIST)ARG3(arg); - wsize = (LIST)ARG4(arg); - wname = (STRING)ARG5(arg); - - 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); + formula=(P)ARG0(arg); + xrange=(LIST)ARG1(arg); + yrange=(LIST)ARG2(arg); + zrange=(LIST)ARG3(arg); + wsize=(LIST)ARG4(arg); + wname=(STRING)ARG5(arg); + 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(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; + can->mode=fn; + if(zrange){ + n=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=QTOS((Q)BDY(n)); else { - /* XXX */ - can->vx = VR((P)BDY(BDY(zrange))); - can->nzstep = n ? QTOS((Q)BDY(n)) : DEFAULTPOLARSTEP; + can->vx=VR((P)BDY(BDY(zrange))); + can->nzstep=QTOS((Q)BDY(n)); } } - - if ( !wsize ) { - can->width = DEFAULTWIDTH; can->height = DEFAULTHEIGHT; + if(!wsize){ + can->width=DEFAULTWIDTH; + can->height=DEFAULTHEIGHT; } else { - can->width = QTOS((Q)BDY(BDY(wsize))); - can->height = QTOS((Q)BDY(NEXT(BDY(wsize)))); + 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 ) { + if(wname) can->wname = BDY(wname); + else can->wname=""; + can->formula=formula; + if(can->mode==MODE_PLOT){ + //plot plotcalc(can); create_canvas(can); plot_print(display,can); - } else if ( can->mode == MODE_POLARPLOT ) { + } else if(can->mode==MODE_POLARPLOT){ + //polarplot polarplotcalc(can); create_canvas(can); plot_print(display,can); } else { + //ifplot,conplot create_canvas(can); - ifplotmain(can); + ifplotmainOld(can); } copy_to_canvas(can); return id; } -int memory_plot(NODE arg,LIST *bytes) -{ +void ifplotmainOld(struct canvas *can){ + int i,width,height; + double ** tabe; + + width=can->width;height=can->height; + tabe=(double **)ALLOCA((width+1)*sizeof(double *)); + for(i=0;iwindow,runningcur); + set_busy(can); set_selection(); + calc(tabe,can,0); + if_printOld(display,tabe,can); + reset_selection(); reset_busy(can); + define_cursor(can->window,normalcur); +} + +int memory_plot(NODE arg,LIST *bytes){ NODE n; struct canvas tmp_can; struct canvas *can; @@ -207,125 +211,115 @@ int memory_plot(NODE arg,LIST *bytes) 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); + 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); - 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 ( zrange ) { - n = NEXT(BDY(zrange)); - can->zmin = ToReal(BDY(n)); n = NEXT(n); can->zmax = ToReal(BDY(n)); - if ( n = NEXT(n) ) - can->nzstep = QTOS((Q)BDY(n)); - else - can->nzstep = MAXGC; - can->mode = MODE_CONPLOT; + can=&tmp_can; + 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( zrange ){ + n=NEXT(BDY(zrange)); + can->zmin=ToReal(BDY(n)); n=NEXT(n); can->zmax=ToReal(BDY(n)); + if(n=NEXT(n)) can->nzstep=QTOS((Q)BDY(n)); + else can->nzstep=MAXGC; + can->mode=MODE_CONPLOT; } else - can->mode = MODE_IFPLOT; + can->mode=MODE_IFPLOT; } else - can->mode = MODE_PLOT; - if ( !wsize ) { - can->width = DEFAULTWIDTH; can->height = DEFAULTHEIGHT; + can->mode=MODE_PLOT; + if( !wsize ){ + can->width=DEFAULTWIDTH; can->height=DEFAULTHEIGHT; } else { - can->width = QTOS((Q)BDY(BDY(wsize))); - can->height = QTOS((Q)BDY(NEXT(BDY(wsize)))); + can->width=QTOS((Q)BDY(BDY(wsize))); + can->height=QTOS((Q)BDY(NEXT(BDY(wsize)))); } - can->wname = ""; - can->formula = formula; - if ( can->mode == MODE_PLOT ) { + can->wname=""; + can->formula=formula; + if( can->mode == MODE_PLOT ){ plotcalc(can); memory_print(can,&barray); STOQ(can->width,qw); STOQ(can->height,qh); - n = mknode(3,qw,qh,barray); + n=mknode(3,qw,qh,barray); MKLIST(*bytes,n); } 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)); + width=can->width; height=can->height; + tabe=(double **)ALLOCA(width*sizeof(double *)); + for( i=0; icolor; - if ( !can->window ) - return -1; + id=QTOS((Q)ARG0(arg)); + formula=(P)ARG1(arg); + color=QTOS((Q)ARG2(arg)); + + can=canvas[id]; + orgcolor=can->color; + can->color=color; get_vars_recursive((Obj)formula,&vl); - for ( vl0 = vl; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( vl->v != can->vx && vl->v != can->vy ) - return -1; - if ( argc(arg) == 3 ) - can->color = QTOS((Q)ARG2(arg)); - else - can->color = 0; + for(vl0=vl;vl0;vl0=NEXT(vl0)) + if(vl0->v->attr==(pointer)V_IND) + if(vl->v!=can->vx && vl->v!=can->vy)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; - can->formula = formula; - if ( can->mode == MODE_PLOT ) { + current_can=can; + can->formula=formula; + if(can->mode==MODE_PLOT){ plotcalc(can); plot_print(display,can); - } else - ifplotmain(can); + } else ifplotmainOld(can); copy_to_canvas(can); - can->color = color; + can->color=color; #if !defined(VISUAL) set_drawcolor(can->color); #endif return index; } -int drawcircle(NODE arg) -{ +int drawcircle(NODE arg){ #if !defined(VISUAL) - int id; - int index; + int id,index,wx,wy,wr,c; pointer ptr; Q ret; LIST xyr; Obj x,y,r; - int wx,wy,wr,c; struct canvas *can; - index = QTOS((Q)ARG0(arg)); - xyr = (LIST)ARG1(arg); - x = (Obj)ARG0(BDY(xyr)); y = (Obj)ARG1(BDY(xyr)); r = (Obj)ARG2(BDY(xyr)); - c = QTOS((Q)ARG2(arg)); - can = canvas[index]; - if ( !can->window ) - return -1; + index=QTOS((Q)ARG0(arg)); + xyr=(LIST)ARG1(arg); + x=(Obj)ARG0(BDY(xyr)); y=(Obj)ARG1(BDY(xyr)); r=(Obj)ARG2(BDY(xyr)); + c=QTOS((Q)ARG2(arg)); + can=canvas[index]; + if(!can->window)return -1; else { - current_can = can; + current_can=can; set_drawcolor(c); - wx = (ToReal(x)-can->xmin)*can->width/(can->xmax-can->xmin); - wy = (can->ymax-ToReal(y))*can->height/(can->ymax-can->ymin); - wr = ToReal(r); + wx=(ToReal(x)-can->xmin)*can->width/(can->xmax-can->xmin); + wy=(can->ymax-ToReal(y))*can->height/(can->ymax-can->ymin); + wr=ToReal(r); XFillArc(display,can->pix,cdrawGC,wx-wr/2,wy-wr/2,wr,wr,0,360*64); copy_to_canvas(can); set_drawcolor(can->color); @@ -334,52 +328,47 @@ int drawcircle(NODE arg) #endif } -int draw_obj(NODE arg) -{ - int index; - int x,y,u,v,len,r; +int draw_obj(NODE arg){ + int index,color,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]; + 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) ) { + 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)); + 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; + 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)); + 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; + 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)); + 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; + can->history=n; break; default: set_lasterror("draw_obj : invalid request"); @@ -391,35 +380,31 @@ int draw_obj(NODE arg) return 0; } -int draw_string(NODE arg) -{ - int index,x,y; +int draw_string(NODE arg){ + int index,x,y,color; 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]; + 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) ) { + 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)); + 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); @@ -427,25 +412,22 @@ int draw_string(NODE arg) return 0; } -int clear_canvas(NODE arg) -{ +int clear_canvas(NODE arg){ int index; struct canvas *can; - index = QTOS((Q)ARG0(arg)); - can = canvas[index]; - if ( !can || !can->window ) - return -1; + 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; + // clear the history + can->history=0; } #define RealtoDbl(r) ((r)?BDY(r):0.0) -int arrayplot(NODE arg) -{ +int arrayplot(NODE arg){ int id,ix,w,h; VECT array; LIST xrange,wsize; @@ -456,61 +438,58 @@ int arrayplot(NODE arg) struct canvas *can; POINT *pa; - array = (VECT)ARG0(arg); - xrange = (LIST)ARG1(arg); - can = canvas[id = search_canvas()]; - 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 ( !wsize ) { - can->width = DEFAULTWIDTH; can->height = DEFAULTHEIGHT; + array=(VECT)ARG0(arg); + xrange=(LIST)ARG1(arg); + can=canvas[id=search_canvas()]; + 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(!wsize){ + can->width=DEFAULTWIDTH; + can->height=DEFAULTHEIGHT; } else { - can->width = QTOS((Q)BDY(BDY(wsize))); - can->height = QTOS((Q)BDY(NEXT(BDY(wsize)))); + can->width=QTOS((Q)BDY(BDY(wsize))); + can->height=QTOS((Q)BDY(NEXT(BDY(wsize)))); } - can->wname = wname; can->formula = 0; can->mode = MODE_PLOT; + can->wname=wname; can->formula=0; can->mode=MODE_PLOT; create_canvas(can); - w = array->len; - h = can->height; - tab = (Real *)BDY(array); - if ( can->ymax == can->ymin ) { - for ( ymax = ymin = RealtoDbl(tab[0]), ix = 1; ix < w; ix++ ) { - if ( RealtoDbl(tab[ix]) > ymax ) - ymax = RealtoDbl(tab[ix]); - if ( RealtoDbl(tab[ix]) < ymin ) - ymin = RealtoDbl(tab[ix]); + w=array->len; + h=can->height; + tab=(Real *)BDY(array); + if(can->ymax==can->ymin){ + for(ymax=ymin=RealtoDbl(tab[0]),ix=1; ixymax)ymax=RealtoDbl(tab[ix]); + if(RealtoDbl(tab[ix])ymax = ymax; can->ymin = ymin; + can->ymax=ymax; + can->ymin=ymin; } else { - ymax = can->ymax; ymin = can->ymin; + ymax=can->ymax; + ymin=can->ymin; } - dy = ymax-ymin; - can->pa = (struct pa *)MALLOC(sizeof(struct pa)); - can->pa[0].length = w; - can->pa[0].pos = pa = (POINT *)MALLOC(w*sizeof(POINT)); - xstep = (double)can->width/(double)(w-1); - for ( ix = 0; ix < w; ix++ ) { + dy=ymax-ymin; + can->pa=(struct pa *)MALLOC(sizeof(struct pa)); + can->pa[0].length=w; + can->pa[0].pos=pa=(POINT *)MALLOC(w*sizeof(POINT)); + xstep=(double)can->width/(double)(w-1); + for(ix=0;ix MAXSHORT ) - pa[ix].y = MAXSHORT; - else if ( t < -MAXSHORT ) - pa[ix].y = -MAXSHORT; - else - pa[ix].y = (long)t; + pa[ix].x=(int)(ix*xstep); + t=(h - 1)*(ymax - RealtoDbl(tab[ix]))/dy; + if(t>MAXSHORT)pa[ix].y=MAXSHORT; + else if(t<-MAXSHORT)pa[ix].y=-MAXSHORT; + else pa[ix].y=(long)t; } plot_print(display,can); copy_to_canvas(can); return id; } -void ifplot_resize(struct canvas *can,POINT spos,POINT epos) -{ +void ifplot_resize(struct canvas *can,POINT spos,POINT epos){ struct canvas *ncan; struct canvas fakecan; Q dx,dy,dx2,dy2,xmin,xmax,ymin,ymax,xmid,ymid; @@ -519,41 +498,57 @@ void ifplot_resize(struct canvas *can,POINT spos,POINT int new; int w,h,m; - if ( XC(spos) < XC(epos) && YC(spos) < YC(epos) ) { - if ( can->precise && !can->wide ) { - fakecan = *can; ncan = &fakecan; + if(XC(spos)precise && !can->wide){ + fakecan=*can; + ncan=&fakecan; } else { - new = search_canvas(); ncan = canvas[new]; + new=search_canvas(); + ncan=canvas[new]; } - ncan->mode = can->mode; - ncan->zmin = can->zmin; ncan->zmax = can->zmax; - ncan->nzstep = can->nzstep; - ncan->wname = can->wname; - ncan->vx = can->vx; ncan->vy = can->vy; - ncan->formula = can->formula; - w = XC(epos)-XC(spos); - h = YC(epos)-YC(spos); - m = MAX(can->width,can->height); - if ( can->precise ) { - ncan->width = w; ncan->height = h; - } else if ( w > h ) { - ncan->width = m; ncan->height = m * h/w; + ncan->mode=can->mode; + ncan->zmin=can->zmin; ncan->zmax=can->zmax; + ncan->nzstep=can->nzstep; + ncan->wname=can->wname; + ncan->vx=can->vx; ncan->vy=can->vy; + ncan->formula=can->formula; + w=XC(epos)-XC(spos); + h=YC(epos)-YC(spos); + m=MAX(can->width,can->height); + if(can->precise){ + ncan->width=w; + ncan->height=h; + } else if(w>h){ + ncan->width=m; + ncan->height=m*h/w; } else { - ncan->width = m * w/h; ncan->height = m; + ncan->width=m*w/h; + ncan->height=m; } - if ( can->wide ) { - STOQ(10,ten); STOQ(2,two); - subq(can->qxmax,can->qxmin,&t); mulq(t,ten,&dx); - subq(can->qymax,can->qymin,&t); mulq(t,ten,&dy); - addq(can->qxmax,can->qxmin,&t); divq(t,two,&xmid); - addq(can->qymax,can->qymin,&t); divq(t,two,&ymid); - divq(dx,two,&dx2); divq(dy,two,&dy2); - subq(xmid,dx2,&xmin); addq(xmid,dx2,&xmax); - subq(ymid,dy2,&ymin); addq(ymid,dy2,&ymax); + if(can->wide){ + STOQ(10,ten); + STOQ(2,two); + subq(can->qxmax,can->qxmin,&t); + mulq(t,ten,&dx); + subq(can->qymax,can->qymin,&t); + mulq(t,ten,&dy); + addq(can->qxmax,can->qxmin,&t); + divq(t,two,&xmid); + addq(can->qymax,can->qymin,&t); + divq(t,two,&ymid); + divq(dx,two,&dx2); + divq(dy,two,&dy2); + subq(xmid,dx2,&xmin); + addq(xmid,dx2,&xmax); + subq(ymid,dy2,&ymin); + addq(ymid,dy2,&ymax); } else { - subq(can->qxmax,can->qxmin,&dx); subq(can->qymax,can->qymin,&dy); - xmin = can->qxmin; xmax = can->qxmax; - ymin = can->qymin; ymax = can->qymax; + subq(can->qxmax,can->qxmin,&dx); + subq(can->qymax,can->qymin,&dy); + xmin=can->qxmin; + xmax=can->qxmax; + ymin=can->qymin; + ymax=can->qymax; } STOQ(XC(spos),sx); STOQ(YC(spos),sy); STOQ(XC(epos),ex); STOQ(YC(epos),ey); STOQ(can->width,cw); STOQ(can->height,ch); @@ -561,30 +556,27 @@ void ifplot_resize(struct canvas *can,POINT spos,POINT mulq(ex,dx,&t); divq(t,cw,&s); addq(xmin,s,&ncan->qxmax); mulq(ey,dy,&t); divq(t,ch,&s); subq(ymax,s,&ncan->qymin); mulq(sy,dy,&t); divq(t,ch,&s); subq(ymax,s,&ncan->qymax); - ncan->xmin = ToReal(ncan->qxmin); ncan->xmax = ToReal(ncan->qxmax); - ncan->ymin = ToReal(ncan->qymin); ncan->ymax = ToReal(ncan->qymax); - if ( can->precise && !can->wide ) { - current_can = can; + ncan->xmin=ToReal(ncan->qxmin); ncan->xmax=ToReal(ncan->qxmax); + ncan->ymin=ToReal(ncan->qymin); ncan->ymax=ToReal(ncan->qymax); + if(can->precise && !can->wide){ + current_can=can; alloc_pixmap(ncan); #if defined(VISUAL) - ncan->real_can = can; + ncan->real_can=can; #endif qifplotmain(ncan); copy_subimage(ncan,can,spos); copy_to_canvas(can); } else { create_canvas(ncan); - if ( can->precise ) - qifplotmain(ncan); - else - ifplotmain(ncan); + if( can->precise ) qifplotmain(ncan); + else ifplotmain(ncan); copy_to_canvas(ncan); } } } -void plot_resize(struct canvas *can,POINT spos,POINT epos) -{ +void plot_resize(struct canvas *can,POINT spos,POINT epos){ struct canvas *ncan; Q dx,dx2,xmin,xmax,xmid; double dy,ymin,ymax,ymid; @@ -593,45 +585,45 @@ void plot_resize(struct canvas *can,POINT spos,POINT e int new; int w,h,m; - if ( XC(spos) < XC(epos) && YC(spos) < YC(epos) ) { - new = search_canvas(); ncan = canvas[new]; - ncan->mode = can->mode; - ncan->zmin = can->zmin; ncan->zmax = can->zmax; - ncan->nzstep = can->nzstep; - ncan->wname = can->wname; - ncan->vx = can->vx; ncan->vy = can->vy; - ncan->formula = can->formula; - w = XC(epos)-XC(spos); - h = YC(epos)-YC(spos); - m = MAX(can->width,can->height); - if ( w > h ) { - ncan->width = m; ncan->height = m * h/w; + if( XC(spos)mode=can->mode; + ncan->zmin=can->zmin; ncan->zmax=can->zmax; + ncan->nzstep=can->nzstep; + ncan->wname=can->wname; + ncan->vx=can->vx; ncan->vy=can->vy; + ncan->formula=can->formula; + w=XC(epos)-XC(spos); + h=YC(epos)-YC(spos); + m=MAX(can->width,can->height); + if( w>h ){ + ncan->width=m; ncan->height=m * h/w; } else { - ncan->width = m * w/h; ncan->height = m; + ncan->width=m * w/h; ncan->height=m; } - if ( can->wide ) { + if( can->wide ){ STOQ(10,ten); STOQ(2,two); subq(can->qxmax,can->qxmin,&t); mulq(t,ten,&dx); addq(can->qxmax,can->qxmin,&t); divq(t,two,&xmid); divq(dx,two,&dx2); subq(xmid,dx2,&xmin); addq(xmid,dx2,&xmax); - dy = (can->ymax-can->ymin)*10; - ymid = (can->ymax+can->ymin)/2; - ymin = ymid-dy/2; ymax = ymid+dy/2; + dy=(can->ymax-can->ymin)*10; + ymid=(can->ymax+can->ymin)/2; + ymin=ymid-dy/2; ymax=ymid+dy/2; } else { subq(can->qxmax,can->qxmin,&dx); - xmin = can->qxmin; xmax = can->qxmax; + xmin=can->qxmin; xmax=can->qxmax; - dy = can->ymax-can->ymin; - ymin = can->ymin; ymax = can->ymax; + dy=can->ymax-can->ymin; + ymin=can->ymin; ymax=can->ymax; } STOQ(XC(spos),sx); STOQ(XC(epos),ex); STOQ(can->width,cw); mulq(sx,dx,&t); divq(t,cw,&s); addq(xmin,s,&ncan->qxmin); mulq(ex,dx,&t); divq(t,cw,&s); addq(xmin,s,&ncan->qxmax); - ncan->xmin = ToReal(ncan->qxmin); ncan->xmax = ToReal(ncan->qxmax); + ncan->xmin=ToReal(ncan->qxmin); ncan->xmax=ToReal(ncan->qxmax); - ncan->ymin = ymax-YC(epos)*dy/can->height; - ncan->ymax = ymax-YC(spos)*dy/can->height; + ncan->ymin=ymax-YC(epos)*dy/can->height; + ncan->ymax=ymax-YC(spos)*dy/can->height; create_canvas(ncan); plotcalc(ncan); @@ -640,134 +632,322 @@ void plot_resize(struct canvas *can,POINT spos,POINT e } } -void ifplotmain(struct canvas *can) +void qifplotmain(struct canvas *can) { int width,height; - double **tabe; + char **tabe,*tabeb; int i; -#if defined(INTERVAL) - tstart(); -#endif - width = can->width; height = can->height; - tabe = (double **)ALLOCA((width+1)*sizeof(double *)); - for ( i = 0; i < width; i++ ) - tabe[i] = (double *)ALLOCA((height+1)*sizeof(double)); + width=can->width; height=can->height; + tabe=(char **)ALLOCA(width*sizeof(char *)+width*height*sizeof(char)); + bzero((void *)tabe,width*sizeof(char *)+width*height*sizeof(char)); + for( i=0, tabeb=(char *)(tabe+width); iwindow,runningcur); set_busy(can); set_selection(); - calc(tabe,can,0); if_print(display,tabe,can); + qcalc(tabe,can); qif_print(display,tabe,can); reset_selection(); reset_busy(can); define_cursor(can->window,normalcur); -#if defined(INTERVAL) - tstop(can); -#endif } -void qifplotmain(struct canvas *can) -{ - int width,height; - char **tabe,*tabeb; - int i; +//*******************ifplotNG +int ifplotNG(NODE arg,char *func){ + int id,orgcolor,color,op_code; + NODE n; + struct canvas *can; + P formula; + LIST xrange,yrange,zrange,wsize; + STRING wname; - width = can->width; height = can->height; - tabe = (char **)ALLOCA(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++ ) - tabe[i] = tabeb + height*i; + formula=(P)ARG0(arg); + color=QTOS((Q)ARG1(arg)); + xrange=(LIST)ARG2(arg); + yrange=(LIST)ARG3(arg); + zrange=(LIST)ARG4(arg); + wsize=(LIST)ARG5(arg); + wname=(STRING)ARG6(arg); + + can=canvas[id=search_canvas()]; + orgcolor=can->color; + can->color=color; + can->division=0; + // set canvas data + 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(zrange){ + n=BDY(zrange); can->zmin=ToReal(BDY(n)); + n=NEXT(n); can->zmax=ToReal(BDY(n)); + n=NEXT(n); can->nzstep=QTOS((Q)BDY(n)); + } + 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; + set_drawcolor(color); + if(!strcmp(func,"ifplot"))can->mode=MODE_IFPLOT; + else if(!strcmp(func,"ifplotD"))can->mode=MODE_IFPLOTD; + else if(!strcmp(func,"ifplotQ"))can->mode=MODE_IFPLOTQ; + else if(!strcmp(func,"ifplotB"))can->mode=MODE_IFPLOTB; + else if(!strcmp(func,"ineqnD"))can->mode=MODE_INEQND; + else if(!strcmp(func,"ineqnQ"))can->mode=MODE_INEQNQ; + else if(!strcmp(func,"ineqnB"))can->mode=MODE_INEQNB; + else if(!strcmp(func,"conplotD"))can->mode=MODE_CONPLOTD; + else if(!strcmp(func,"conplotQ"))can->mode=MODE_CONPLOTQ; + else if(!strcmp(func,"conplotB"))can->mode=MODE_CONPLOTB; + else if(!strcmp(func,"itvifplot")){ + can->mode=MODE_ITVIFPLOT; + can->division=QTOS((Q)ARG7(arg)); + } + else can->mode=MODE_IFPLOTD; + create_canvas(can); + ifplotmain(can); + set_drawcolor(orgcolor); + copy_to_canvas(can); + can->color=orgcolor; + return id; +} + +int ifplotOP(NODE arg,char *func){ + //ineqnor[D,Q,B],ineqnand[D,Q,B],ineqnxor[D,Q,b] + int index,orgcolor,color,op_code; + NODE n; + struct canvas *can; + P formula; + //s_id,fname,w_id,poly,color + index=QTOS((Q)ARG0(arg)); + formula=(P)ARG1(arg); + color=QTOS((Q)ARG2(arg)); + can=canvas[index]; + orgcolor=can->color; + can->color=color; + // set canvas data + can->formula=formula; + if(!strcmp(func,"ineqnandD"))can->mode=MODE_INEQNANDD; + else if(!strcmp(func,"ineqnandQ"))can->mode=MODE_INEQNANDQ; + else if(!strcmp(func,"ineqnandB"))can->mode=MODE_INEQNANDB; + else if(!strcmp(func,"ineqnorD"))can->mode=MODE_INEQNORD; + else if(!strcmp(func,"ineqnorQ"))can->mode=MODE_INEQNORQ; + else if(!strcmp(func,"ineqnorB"))can->mode=MODE_INEQNORB; + else if(!strcmp(func,"ineqnxorD"))can->mode=MODE_INEQNXORD; + else if(!strcmp(func,"ineqnxorQ"))can->mode=MODE_INEQNXORQ; + else if(!strcmp(func,"ineqnxorB"))can->mode=MODE_INEQNXORB; + else if(!strcmp(func,"plotoverD"))can->mode=MODE_PLOTOVERD; + else if(!strcmp(func,"plotoverQ"))can->mode=MODE_PLOTOVERQ; + else if(!strcmp(func,"plotoverB"))can->mode=MODE_PLOTOVERB; + else can->mode=MODE_IFPLOTD; + set_drawcolor(color); + ifplotmain(can); + set_drawcolor(orgcolor); + copy_to_canvas(can); + can->color=orgcolor; + return index; +} + +void ifplotmain(struct canvas *can){ + int width,height,i,j,ix,iy,**mask; + double **tabe; + + width=can->width; height=can->height; + tabe=(double **)ALLOCA((width+1)*sizeof(double *)); + for(i=0;iwindow,runningcur); set_busy(can); set_selection(); - qcalc(tabe,can); qif_print(display,tabe,can); + set_drawcolor(can->color); + switch(can->mode){ + case MODE_IFPLOTD: + calc(tabe,can,0); + if_print(display,tabe,can,1); + break; + case MODE_IFPLOTQ: + calcq(tabe,can,0); + if_print(display,tabe,can,1); + break; + case MODE_IFPLOTB: + calcb(tabe,can,0); + if_print(display,tabe,can,0); + break; + case MODE_INEQND: + calc(tabe,can,0); + area_print(display,tabe,can,0); + break; + case MODE_INEQNQ: + calcq(tabe,can,0); + area_print(display,tabe,can,0); + break; + case MODE_INEQNB: + calcb(tabe,can,0); + area_print(display,tabe,can,0); + break; + case MODE_INEQNANDD: + calc(tabe,can,0); + area_print(display,tabe,can,0); + break; + case MODE_INEQNANDQ: + calcq(tabe,can,0); + area_print(display,tabe,can,2); + break; + case MODE_INEQNANDB: + calcb(tabe,can,0); + area_print(display,tabe,can,2); + break; + case MODE_INEQNORD: + calc(tabe,can,0); + area_print(display,tabe,can,3); + break; + case MODE_INEQNORQ: + calcq(tabe,can,0); + area_print(display,tabe,can,3); + break; + case MODE_INEQNORB: + calcb(tabe,can,0); + area_print(display,tabe,can,3); + break; + case MODE_INEQNXORD: + calc(tabe,can,0); + area_print(display,tabe,can,4); + break; + case MODE_INEQNXORQ: + calcq(tabe,can,0); + area_print(display,tabe,can,4); + break; + case MODE_INEQNXORB: + calcb(tabe,can,0); + area_print(display,tabe,can,4); + break; + case MODE_CONPLOTD: + calc(tabe,can,0); + con_print(display,tabe,can); + break; + case MODE_CONPLOTB: + calcb(tabe,can,0); + con_print(display,tabe,can); + break; + case MODE_CONPLOTQ: + calcq(tabe,can,0); + con_print(display,tabe,can); + break; + case MODE_PLOTOVERD: + calc(tabe,can,0); + over_print(display,tabe,can,0); + break; + case MODE_PLOTOVERQ: + calcq(tabe,can,0); + over_print(display,tabe,can,0); + break; + case MODE_PLOTOVERB: + calcb(tabe,can,0); + over_print(display,tabe,can,0); + break; +#if defined(INTERVAL) + case MODE_ITVIFPLOT: + itvcalc(tabe,can,1); + if_print(display,tabe,can,1); + break; +#endif + } + set_drawcolor(can->color); reset_selection(); reset_busy(can); define_cursor(can->window,normalcur); } -#if defined(INTERVAL) -int objcp(NODE arg) -{ +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]; + 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) -{ +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); + 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) { + 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;idata = imgsrc->data; + imgtrg->data=imgsrc->data; break; case 6:/* xor case */ for(i=0;icolor; - 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; -} - +//ifplotNG +/* int ineqnover(NODE arg) { - int id; - struct canvas *can; - int orgcolor, op_code; + 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; + 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; + 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)); + 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); + define_cursor(can->window,runningcur); + set_busy(can); set_selection(); +// ineqncalc(tbl, can, 1); + calc(tab,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; + set_drawcolor(orgcolor); + can->color = orgcolor; #endif - copy_to_canvas(can); + copy_to_canvas(can); } - -int itvifplot(NODE arg) -{ - int id, op_code, orgcolor; - struct canvas *can; - LIST xrange, yrange, zrange, geom; +*/ +int polarplotNG(NODE arg){ + int id,color,orgcolor; 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; + P formula; + LIST xrange,wsize; STRING wname; - double **tbl; + V v; - 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; -} + formula=(P)ARG0(arg); + xrange=(LIST)ARG1(arg); + color=QTOS((Q)ARG2(arg)); + wsize=(LIST)ARG3(arg); + wname=(STRING)ARG4(arg); -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; + 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); } - - 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 = ""; + can->mode=MODE_POLARPLOT; + 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; + orgcolor=can->color; + can->color=color; + polarplotcalc(can); create_canvas(can); - itvplotmain2(can); -#if !defined(VISUAL) - set_drawcolor(orgcolor); - can->color = orgcolor; -#endif + plot_print(display,can); + can->color=orgcolor; 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; +int plotoverD(NODE arg){ + int index, color; + P formula; struct canvas *can; - LIST xrange, yrange, zrange, geom; - NODE n; - STRING wname; - double **tbl; - int itvsize; + VL vl,vl0; - 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); + index=QTOS((Q)ARG0(arg)); + formula=(P)ARG1(arg); + can=canvas[index]; + color=can->color; + if(!can->window)return -1; + get_vars_recursive((Obj)formula,&vl); + for(vl0=vl;vl0;vl0=NEXT(vl0)) + if(vl0->v->attr==(pointer)V_IND) + if(vl->v!=can->vx && vl->v!=can->vy) return -1; + if(argc(arg)== 3)can->color=QTOS((Q)ARG2(arg)); + else can->color=0; #if !defined(VISUAL) - set_drawcolor(orgcolor); - can->color = orgcolor; + set_drawcolor(can->color); #endif + current_can=can; + can->formula=formula; + if(can->mode==MODE_PLOT){ + calcq(can); + plot_print(display,can); + } else ifplotmain(can); 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); + can->color=color; #if !defined(VISUAL) - set_drawcolor(orgcolor); - can->color = orgcolor; + set_drawcolor(can->color); #endif - copy_to_canvas(can); - return id; + return index; } - -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