=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/plot/if.c,v retrieving revision 1.1 retrieving revision 1.22 diff -u -p -r1.1 -r1.22 --- OpenXM_contrib2/asir2000/plot/if.c 1999/12/03 07:39:13 1.1 +++ OpenXM_contrib2/asir2000/plot/if.c 2011/08/10 04:51:58 1.22 @@ -1,11 +1,125 @@ -/* $OpenXM: OpenXM/src/asir99/plot/if.c,v 1.1.1.1 1999/11/10 08:12:34 noro Exp $ */ +/* + * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED + * All rights reserved. + * + * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited, + * non-exclusive and royalty-free license to use, copy, modify and + * redistribute, solely for non-commercial and non-profit purposes, the + * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and + * conditions of this Agreement. For the avoidance of doubt, you acquire + * only a limited right to use the SOFTWARE hereunder, and FLL or any + * third party developer retains all rights, including but not limited to + * copyrights, in and to the SOFTWARE. + * + * (1) FLL does not grant you a license in any way for commercial + * purposes. You may use the SOFTWARE only for non-commercial and + * non-profit purposes only, such as academic, research and internal + * business use. + * (2) The SOFTWARE is protected by the Copyright Law of Japan and + * international copyright treaties. If you make copies of the SOFTWARE, + * with or without modification, as permitted hereunder, you shall affix + * to all such copies of the SOFTWARE the above copyright notice. + * (3) An explicit reference to this SOFTWARE and its copyright owner + * shall be made on your publication or presentation in any form of the + * results obtained by use of the SOFTWARE. + * (4) In the event that you modify the SOFTWARE, you shall notify FLL by + * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification + * for such modification or the source code of the modified part of the + * SOFTWARE. + * + * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL + * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND + * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS + * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES' + * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY + * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY. + * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT, + * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY + * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL + * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES + * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES + * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY + * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF + * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART + * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY + * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, + * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. + * + * $OpenXM: OpenXM_contrib2/asir2000/plot/if.c,v 1.21 2006/11/09 15:54:35 saito Exp $ +*/ #include "ca.h" #include "parse.h" #include "ox.h" #include "ifplot.h" -extern jmp_buf ox_env; +#if defined(INTERVAL) +/* Time message and func*/ +#include +#include +#include +static struct oEGT ltime; +static double r0; +double get_rtime(); +#if defined(ITV_TIME_CHECK) +void tstart() +{ + get_eg(<ime); + 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; + +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 id; @@ -14,6 +128,7 @@ int plot(NODE arg) P formula; LIST xrange,yrange,zrange,wsize; STRING wname; + V v; formula = (P)ARG0(arg); xrange = (LIST)ARG1(arg); @@ -23,6 +138,83 @@ int plot(NODE 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 ( 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); can->qxmin = (Q)BDY(n); n = NEXT(n); can->qxmax = (Q)BDY(n); can->xmin = ToReal(can->qxmin); can->xmax = ToReal(can->qxmax); @@ -48,79 +240,208 @@ int plot(NODE arg) 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->wname = ""; can->formula = formula; - create_canvas(can); if ( can->mode == MODE_PLOT ) { plotcalc(can); - plot_print(display,can); - } else - ifplotmain(can); - copy_to_canvas(can); - return id; + 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; + 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 index; + int index, color; P formula; struct canvas *can; - struct canvas fakecan; VL vl,vl0; index = QTOS((Q)ARG0(arg)); formula = (P)ARG1(arg); can = canvas[index]; + color = can->color; if ( !can->window ) return -1; - get_vars_recursive(formula,&vl); + get_vars_recursive((Obj)formula,&vl); 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 ) 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; - fakecan = *can; fakecan.formula = formula; + can->formula = formula; if ( can->mode == MODE_PLOT ) { - plotcalc(&fakecan); - plot_print(display,&fakecan); + plotcalc(can); + plot_print(display,can); } else - ifplotmain(&fakecan); - copy_to_canvas(&fakecan); + ifplotmain(can); + copy_to_canvas(can); + can->color = color; +#if !defined(VISUAL) + set_drawcolor(can->color); +#endif return index; } int drawcircle(NODE arg) { +#if !defined(VISUAL) int id; int index; pointer ptr; Q ret; LIST xyr; Obj x,y,r; - int wx,wy,wr; + int wx,wy,wr,c; struct canvas *can; - struct canvas fakecan; 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; + 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); - 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); + set_drawcolor(can->color); 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) int arrayplot(NODE arg) @@ -130,7 +451,6 @@ int arrayplot(NODE arg) LIST xrange,wsize; char *wname; NODE n; - Q ret; double ymax,ymin,dy,xstep; Real *tab; struct canvas *can; @@ -182,16 +502,14 @@ int arrayplot(NODE arg) else if ( t < -MAXSHORT ) pa[ix].y = -MAXSHORT; else - pa[ix].y = t; + pa[ix].y = (long)t; } plot_print(display,can); copy_to_canvas(can); return id; } -ifplot_resize(can,spos,epos) -struct canvas *can; -POINT spos,epos; +void ifplot_resize(struct canvas *can,POINT spos,POINT epos) { struct canvas *ncan; struct canvas fakecan; @@ -200,7 +518,6 @@ POINT spos,epos; Q s,t; int new; int w,h,m; - Q ret; if ( XC(spos) < XC(epos) && YC(spos) < YC(epos) ) { if ( can->precise && !can->wide ) { @@ -249,6 +566,9 @@ POINT spos,epos; if ( can->precise && !can->wide ) { current_can = can; alloc_pixmap(ncan); +#if defined(VISUAL) + ncan->real_can = can; +#endif qifplotmain(ncan); copy_subimage(ncan,can,spos); copy_to_canvas(can); @@ -263,15 +583,12 @@ POINT spos,epos; } } -plot_resize(can,spos,epos) -struct canvas *can; -POINT spos,epos; +void plot_resize(struct canvas *can,POINT spos,POINT epos) { struct canvas *ncan; Q dx,dx2,xmin,xmax,xmid; - double dy,dy2,ymin,ymax,ymid; + double dy,ymin,ymax,ymid; Q sx,ex,cw,ten,two; - double sy,ey; Q s,t; int new; int w,h,m; @@ -323,26 +640,30 @@ POINT spos,epos; } } -ifplotmain(can) -struct canvas *can; +void ifplotmain(struct canvas *can) { int width,height; - double **tabe,*tabeb; + double **tabe; int i; +#if defined(INTERVAL) + tstart(); +#endif 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++ ) - tabe[i] = (double *)ALLOCA(height*sizeof(double)); + tabe[i] = (double *)ALLOCA((height+1)*sizeof(double)); define_cursor(can->window,runningcur); 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); define_cursor(can->window,normalcur); +#if defined(INTERVAL) + tstop(can); +#endif } -qifplotmain(can) -struct canvas *can; +void qifplotmain(struct canvas *can) { int width,height; char **tabe,*tabeb; @@ -350,7 +671,7 @@ struct canvas *can; width = can->width; height = can->height; 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++ ) tabe[i] = tabeb + height*i; define_cursor(can->window,runningcur); @@ -359,3 +680,551 @@ struct canvas *can; reset_selection(); reset_busy(can); 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;idata = imgsrc->data; + break; + case 6:/* xor case */ + for(i=0;ipix, 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