=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/plot/plotf.c,v retrieving revision 1.23 retrieving revision 1.33 diff -u -p -r1.23 -r1.33 --- OpenXM_contrib2/asir2000/plot/plotf.c 2012/10/04 04:56:39 1.23 +++ OpenXM_contrib2/asir2000/plot/plotf.c 2017/09/04 01:57:53 1.33 @@ -45,185 +45,244 @@ * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. * - * $OpenXM: OpenXM_contrib2/asir2000/plot/plotf.c,v 1.22 2011/08/10 04:51:58 saito Exp $ + * $OpenXM: OpenXM_contrib2/asir2000/plot/plotf.c,v 1.32 2017/09/01 01:48:40 noro Exp $ */ #include "ca.h" #include "parse.h" #include "ox.h" #include "ifplot.h" -void Pifplot(NODE, Obj *), Pconplot(NODE, Obj *), Pplotover(NODE, Q *); -void Pplot(NODE, Obj *), Parrayplot(NODE, Q *), Pdrawcircle(NODE, Q *); +int validate_ox_plot_stream(int); +void ListCheck(char *,LIST); +void Pplot(NODE,Obj *); +void Ppolarplot(NODE,Obj *); +void Pobj_cp(NODE,Obj *),Parrayplot(NODE,Obj*),Pdrawcircle(NODE,Obj*); +void Pifplot(NODE,Obj *),Pconplot(NODE,Obj *),Pplotover(NODE,Obj *); +void Pmemory_ifplot(NODE,Obj *),Pmemory_conplot(NODE,Obj *); +void Pmemory_plot(NODE,Obj *); +void ifplot_main(NODE,int,char *,Obj *); +void plot_main(NODE,int,char *,Obj *); +void conplot_main(NODE,int,Obj *); + +void Popen_canvas(NODE,Obj *),Pclear_canvas(NODE,Obj *),Pdraw_obj(NODE,Obj *); +void Pdraw_string(NODE,Obj *); +void Pox_rpc(NODE,Obj *), Pox_cmo_rpc(NODE,Obj *); + +//NG #if defined(INTERVAL) -void Pitvplot1(NODE, Obj *), itvplot_main1(NODE, Obj *); //NORMAL -void Pitvplot2(NODE, Obj *), itvplot_main2(NODE, Obj *); //TRANSFER -void Pitvplot3(NODE, Obj *), itvplot_main3(NODE, Obj *); //RECURSION -void Pitvplot4(NODE, Obj *), itvplot_main4(NODE, Obj *); //RECURSION+TRANSFER -void Pitvplot(NODE, Obj *), itvplot_main(NODE, Obj *); //NEWTYPE -void Pobj_cp(NODE, Obj *); -void Pineqn(NODE, Obj *), Pineqnor(NODE, Obj *), Pineqnand(NODE, Obj *); -void Pineqnxor(NODE, Obj *), Pmemory_ineqn(NODE, Obj *); -void ineqn_main(NODE, int, Obj *); -void ineqn_main_op(NODE, int, int, Obj *); +void Pitvifplot(NODE, Obj *); #endif -void Ppolarplot(NODE, Q *); -void Pmemory_ifplot(NODE, Obj *), Pmemory_conplot(NODE, Obj *); -void Pmemory_plot(NODE, Obj *); -void Popen_canvas(NODE, Q *), Pclear_canvas(NODE, Q *), Pdraw_obj(NODE, Q *); -void Pdraw_string(NODE, Q *); -void Pox_rpc(), Pox_cmo_rpc(); -void ifplot_main(NODE, int, Obj *); -void conplot_main(NODE, int, Obj *); -void plot_main(NODE, int, Obj *); - -struct ftab plot_tab[] = { - {"ifplot",Pifplot,-7}, - {"memory_ifplot",Pmemory_ifplot,-6}, +void PifplotD(NODE,Obj *),PifplotQ(NODE,Obj *),PifplotB(NODE,Obj *); +void PineqnD(NODE,Obj *), PineqnQ(NODE,Obj *), PineqnB(NODE,Obj *); +void PineqnandD(NODE,Obj *), PineqnandQ(NODE,Obj *), PineqnandB(NODE,Obj *); +void PineqnorD(NODE,Obj *), PineqnorQ(NODE,Obj *), PineqnorB(NODE,Obj *); +void PineqnxorD(NODE,Obj *), PineqnxorQ(NODE,Obj *), PineqnxorB(NODE,Obj *); +void PconplotD(NODE,Obj *),PconplotQ(NODE,Obj *),PconplotB(NODE,Obj *); +void PpolarplotD(NODE,Obj *); +void PplotoverD(NODE,Obj *),PplotoverQ(NODE,Obj *),PplotoverB(NODE,Obj *); +void ifplot_mainNG(NODE,char *,Obj *); +void conplot_mainNG(NODE,char *,Obj *); +void plotover_mainNG(NODE,char *,Obj *); + +struct ftab plot_tab[]={ + {PLOT,Pplot,-7}, #if defined(INTERVAL) - {"itvplot1",Pitvplot1,-7}, - {"itvplot2",Pitvplot2,-7}, - {"itvplot3",Pitvplot3,-7}, - {"itvplot4",Pitvplot4,-7}, - {"itvplot",Pitvplot,-7}, - {"ineqn",Pineqn,-8}, - {"ineqnor",Pineqnor,-7}, - {"ineqnand",Pineqnand,-7}, - {"ineqnxor",Pineqnxor,-7}, - {"memory_ineqn",Pmemory_ineqn,-6}, - {"obj_cp",Pobj_cp,4}, + {ITVIFPLOT,Pitvifplot,-8}, #endif - {"conplot",Pconplot,-8}, - {"memory_conplot",Pmemory_conplot,-7}, - {"plot",Pplot,-6}, - {"memory_plot",Pmemory_plot,-5}, - {"polarplot",Ppolarplot,-6}, - {"plotover",Pplotover,-4}, - {"drawcircle",Pdrawcircle,6}, - {"open_canvas",Popen_canvas,-3}, - {"clear_canvas",Pclear_canvas,2}, - {"draw_obj",Pdraw_obj,-4}, - {"draw_string",Pdraw_string,-5}, + {OBJ_CP,Pobj_cp,4}, + {POLARPLOT,Ppolarplot,-6}, + {POLARPLOTD,PpolarplotD,-7}, + {IFPLOT,Pifplot,-7}, + {IFPLOTD,PifplotD,-8}, + {IFPLOTQ,PifplotQ,-8}, + {IFPLOTB,PifplotB,-8}, + {INEQN,PineqnD,-8}, + {INEQND,PineqnD,-8}, + {INEQNQ,PineqnQ,-8}, + {INEQNB,PineqnB,-8}, + {INEQNAND,PineqnandD,-4}, + {INEQNDAND,PineqnandD,-4}, + {INEQNQAND,PineqnandQ,-4}, + {INEQNBAND,PineqnandB,-4}, + {INEQNOR,PineqnorD,-4}, + {INEQNDOR,PineqnorD,-4}, + {INEQNQOR,PineqnorQ,-4}, + {INEQNBOR,PineqnorB,-4}, + {INEQNXOR,PineqnxorD,-4}, + {INEQNDXOR,PineqnxorD,-4}, + {INEQNQXOR,PineqnxorQ,-4}, + {INEQNBXOR,PineqnxorB,-4}, + {CONPLOT,Pconplot,-7}, + {CONPLOTD,PconplotD,-8}, + {CONPLOTB,PconplotB,-8}, + {CONPLOTQ,PconplotQ,-8}, + {PLOTOVER,Pplotover,-4}, + {PLOTOVERD,PplotoverD,-4}, + {PLOTOVERQ,PplotoverQ,-4}, + {PLOTOVERB,PplotoverB,-4}, + {MEMORY_IFPLOT,Pmemory_ifplot,-7}, + {MEMORY_CONPLOT,Pmemory_conplot,-7}, + {MEMORY_PLOT,Pmemory_plot,-7}, + {DRAWCIRCLE,Pdrawcircle,6}, + {OPEN_CANVAS,Popen_canvas,-3}, + {CLEAR_CANVAS,Pclear_canvas,2}, + {DRAW_OBJ,Pdraw_obj,-4}, + {DRAW_STRING,Pdraw_string,-5}, {0,0,0}, }; +char *pfn[]={ + IFPLOT,CONPLOT,PLOT,INTERACTIVE,POLARPLOT,PLOTOVER, + IFPLOTD,IFPLOTQ,IFPLOTB,INEQND,INEQNQ,INEQNB, + INEQNDAND,INEQNQAND,INEQNBAND, + INEQNDOR,INEQNQOR,INEQNBOR, + INEQNDXOR,INEQNQXOR,INEQNBXOR, + CONPLOTD,CONPLOTQ,CONPLOTB,ITVIFPLOT, + PLOTOVERD,PLOTOVERQ,PLOTOVERB, + MEMORY_IFPLOT,MEMORY_CONPLOT,MEMORY_PLOT,ARRAYPLOT,OPEN_CANVAS, + DRAWCIRCLE,DRAW_OBJ,DRAW_STRING,OBJ_CP,CLEAR_CANVAS,POLARPLOTD}; +/* + IFPLOT:0,CONPLOT:1,PLOT:2,INTERACTIVE:3,POLARPLOT:4,PLOTOVER:5, + IFPLOTD:6,IFPLOTQ:7,IFPLOTB:8,INEQND:9,INEQNQ:10,INEQNB:11, + INEQNDAND:12,INEQNQAND:13,INEQNBAND:14, + INEQNDOR:15,INEQNQOR:16,INEQNBOR:17, + INEQNDXOR:18,INEQNQXOR:19,INEQNBXOR:20, + CONPLOTD:21,CONPLOTQ:22,CONPLOTB:23,ITVIFPLOT:24, + PLOTOVERD:25,PLOTOVERQ:26,PLOTOVERB:27, + MEMORY_IFPLOT:28,MEMORY_CONPLOT:29,MEMORY_PLOT:30,ARRAYPLOT:31, + OPEN_CANVAS:32,DRAWCIRCLE:33,DRAW_OBJ:34,DRAW_STRING:35,OBJ_CP:36, + CLEAR_CANVAS:37,POLARPLOTD:38 +*/ +int modeNO(char *fn){ + int i; + char **z; + for(i=0,z=pfn;*z!=NULL;z++,i++)if(!strcmp(fn,*z))return i; + return -1; +} -void Popen_canvas(NODE arg,Q *rp) -{ +void Popen_canvas(NODE arg,Obj *rp){ Q w300,s_id; LIST geom; int stream; NODE n,n0; STRING fname,wname; - geom = 0; wname = 0; stream = -1; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) - stream = 0; - else - switch ( OID(BDY(arg)) ) { - case O_LIST: - geom = (LIST)BDY(arg); - break; - case O_N: - stream = QTOS((Q)BDY(arg)); break; - case O_STR: - wname = (STRING)BDY(arg); break; - default: - error("open_canvas : invalid argument"); break; + geom=0;wname=0;stream=-1; + for(;arg;arg=NEXT(arg)) + if(!BDY(arg)) stream=0; + else switch (OID(BDY(arg))){ + case O_LIST: + geom=(LIST)BDY(arg); + break; + case O_N: + stream=QTOS((Q)BDY(arg)); + break; + case O_STR: + wname=(STRING)BDY(arg); + break; + default: + error("open_canvas : invalid argument"); + break; } - stream = validate_ox_plot_stream(stream); + stream=validate_ox_plot_stream(stream); STOQ(stream,s_id); - if ( !geom ) { + if(!geom){ STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); + MKNODE(n0,w300,0); + MKNODE(n,w300,n0); + MKLIST(geom,n); } - MKSTR(fname,"open_canvas"); - arg = mknode(4,s_id,fname,geom,wname); + MKSTR(fname,OPEN_CANVAS); + arg=mknode(4,s_id,fname,geom,wname); Pox_cmo_rpc(arg,rp); - *rp = s_id; + *rp=(Obj)s_id; } -void Pifplot(NODE arg,Obj *rp) -{ - ifplot_main(arg,0,rp); -} +void Pifplot(NODE arg,Obj *rp){ifplot_main(arg,0,IFPLOT,rp);} +void Pmemory_ifplot(NODE arg,Obj *rp){ifplot_main(arg,1,IFPLOT,rp);} -void Pmemory_ifplot(NODE arg,Obj *rp) -{ - ifplot_main(arg,1,rp); -} - -void ifplot_main(NODE arg,int is_memory, Obj *rp) -{ +void ifplot_main(NODE arg,int is_memory,char *fn,Obj *rp){ Q m2,p2,w300,s_id; NODE defrange; - LIST xrange,yrange,range[2],list,geom; + LIST xrange,yrange,zrange,range[2],list,geom; VL vl,vl0; V v[2],av[2]; int stream,ri,i,sign; - P poly; + Obj poly; P var; NODE n,n0; STRING fname,wname; Obj t; + int found_f; - STOQ(-2,m2); STOQ(2,p2); - MKNODE(n,p2,0); MKNODE(defrange,m2,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = v[1] = 0; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) stream = 0; - else - switch ( OID(BDY(arg)) ) { + STOQ(-2,m2);STOQ(2,p2);MKNODE(n,p2,0);MKNODE(defrange,m2,n); + poly=0;vl=0;geom=0;wname=0;stream=-1;ri=0;xrange=0;yrange=0;zrange=0; + v[0]=v[1]=0; + found_f = 0; + for(;arg;arg=NEXT(arg)) + if(!BDY(arg)){ + if ( !found_f ) { + poly = 0; + found_f = 1; + } else stream=0; + } else + switch(OID(BDY(arg))){ case O_P: - poly = (P)BDY(arg); + poly=(Obj)BDY(arg); get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 2 ) error("ifplot : invalid argument"); - else v[i++] = vl0->v; + for(vl0=vl,i=0;vl0;vl0=NEXT(vl0)){ + if(vl0->v->attr==(pointer)V_IND){ + if(i>=2)error("ifplot : invalid argument"); + else v[i++]=vl0->v; + } + } + found_f = 1; break; case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 1 ) error("ifplot : invalid argument"); - else range[ri++] = list; - else geom = list; + list=(LIST)BDY(arg); + if(OID(BDY(BDY(list)))==O_P) + if(ri>1) error("ifplot : invalid argument"); + else range[ri++]=list; + else geom=list; break; case O_N: - stream = QTOS((Q)BDY(arg)); break; + if ( !found_f ) { + poly = (Obj)BDY(arg); + found_f = 1; + } else stream=QTOS((Q)BDY(arg)); + break; case O_STR: - wname = (STRING)BDY(arg); break; + wname=(STRING)BDY(arg);break; default: - error("ifplot : invalid argument"); break; + error("ifplot : invalid argument");break; } - if ( !poly ) error("ifplot : invalid argument"); - switch ( ri ) { + if(!found_f) error("ifplot : invalid argument"); + switch(ri){ case 0: - if ( !v[1] ) error("ifplot : please specify all variables"); - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); + if(!v[1]) error("ifplot : please specify all variables"); + MKV(v[0],var);MKNODE(n,var,defrange);MKLIST(xrange,n); + MKV(v[1],var);MKNODE(n,var,defrange);MKLIST(yrange,n); break; case 1: - if ( !v[1] ) error("ifplot : please specify all variables"); - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) { - xrange = range[0]; - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - } else if ( v[1] == av[0] ) { - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - yrange = range[0]; + if(!v[1]) error("ifplot : please specify all variables"); + av[0]=VR((P)BDY(BDY(range[0]))); + if(!poly || NUM(poly) || v[0]==av[0]){ + xrange=range[0]; + MKV(v[1],var);MKNODE(n,var,defrange);MKLIST(yrange,n); + } else if(v[1]==av[0]){ + MKV(v[0],var);MKNODE(n,var,defrange);MKLIST(xrange,n); + yrange=range[0]; } else error("ifplot : invalid argument"); break; case 2: - av[0] = VR((P)BDY(BDY(range[0]))); - av[1] = VR((P)BDY(BDY(range[1]))); - if ( ((v[0] == av[0]) && (!v[1] || v[1] == av[1])) || - ((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { - xrange = range[0]; yrange = range[1]; + av[0]=VR((P)BDY(BDY(range[0]))); + av[1]=VR((P)BDY(BDY(range[1]))); + if(!poly || NUM(poly) || (((v[0]==av[0])&&(!v[1]||v[1]==av[1]))|| + ((v[0]==av[1])&&(!v[1]||v[1]==av[0])))){ + xrange=range[0];yrange=range[1]; } else error("ifplot : invalid argument"); break; default: - error("ifplot : cannot happen"); break; + error("ifplot : cannot happen");break; } /* ifplot in ox_plot requires [s_id (Q), @@ -234,710 +293,31 @@ void ifplot_main(NODE arg,int is_memory, Obj *rp) geom=[xsize,ysize] (LIST), wname=name (STRING)] */ - - stream = validate_ox_plot_stream(stream); + stream=validate_ox_plot_stream(stream); STOQ(stream,s_id); - if ( !geom ) { - STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); + if(!geom){ + STOQ(300,w300);MKNODE(n0,w300,0);MKNODE(n,w300,n0);MKLIST(geom,n); } - if ( is_memory ) { - MKSTR(fname,"memory_plot"); - arg = mknode(8,s_id,fname,poly,xrange,yrange,0,geom); + if(is_memory){ + MKSTR(fname,MEMORY_PLOT); + arg=mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom); Pox_rpc(arg,&t); - arg = mknode(1,s_id); + arg=mknode(1,s_id); Pox_pop_cmo(arg,rp); } else { - MKSTR(fname,"plot"); - arg = mknode(8,s_id,fname,poly,xrange,yrange,0,geom,wname); + MKSTR(fname,fn); + arg=mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom,wname); Pox_rpc(arg,&t); - *rp = (Obj)s_id; + *rp=(Obj)s_id; } } -#if defined(INTERVAL) -void Pitvplot(NODE arg, Obj *rp) -{ - itvplot_main(arg, rp); -} +void Pconplot(NODE arg,Obj *rp){conplot_main(arg, 0, rp);} +void Pmemory_conplot(NODE arg,Obj *rp){conplot_main(arg,1,rp);} -void itvplot_main(NODE arg, Obj *rp) -{ - Q m2,p2,w300,s_id, itvsize; - NODE defrange; - LIST xrange,yrange,range[2],list,geom; - VL vl,vl0; - V v[2],av[2]; - int stream,ri,i,sign; - P poly; - P var; - NODE n,n0; - STRING fname,wname; - Obj t; - - STOQ(-2,m2); STOQ(2,p2); - MKNODE(n,p2,0); MKNODE(defrange,m2,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = v[1] = 0; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) stream = 0; - else - switch ( OID(BDY(arg)) ) { - case O_P: - poly = (P)BDY(arg); - get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 2 ) error("itvplot : invalid argument"); - else v[i++] = vl0->v; - break; - case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 1 ) error("itvplot : invalid argument"); - else range[ri++] = list; - else geom = list; - break; - case O_N: - stream = QTOS((Q)BDY(arg)); break; - case O_STR: - wname = (STRING)BDY(arg); break; - default: - error("itvplot : invalid argument"); break; - } - if ( !poly ) error("itvplot : invalid argument"); - switch ( ri ) { - case 0: - if ( !v[1] ) error("itvplot : please specify all variables"); - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - break; - case 1: - if ( !v[1] ) error("itvplot : please specify all variables"); - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) { - xrange = range[0]; - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - } else if ( v[1] == av[0] ) { - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - yrange = range[0]; - } else error("itvplot : invalid argument"); - break; - case 2: - av[0] = VR((P)BDY(BDY(range[0]))); - av[1] = VR((P)BDY(BDY(range[1]))); - if ( ((v[0] == av[0]) && (!v[1] || v[1] == av[1])) || - ((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { - xrange = range[0]; yrange = range[1]; - } else error("itvplot : invalid argument"); - break; - default: - error("itvplot : cannot happen"); break; - } - /* itvplot in ox_plot requires - [s_id (Q), - formula (Obj), - xrange=[x,xmin,xmax] (LIST), - yrange=[y,ymin,ymax] (LIST), - zrange=0, - geom=[xsize,ysize] (LIST), - wname=name (STRING)] - */ - stream = validate_ox_plot_stream(stream); - STOQ(stream,s_id); - if ( !geom ) { - STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); - } - MKSTR(fname,"itvifplot"); - STOQ(Itvplot, itvsize); - arg = mknode(9,s_id,fname,poly,xrange,yrange,0,geom,wname,itvsize); - Pox_rpc(arg,&t); - *rp = (Obj)s_id; -} - -void Pitvplot1(NODE arg, Obj *rp) -{ - itvplot_main1(arg, rp); -} - -void Pitvplot2(NODE arg, Obj *rp) -{ - itvplot_main2(arg, rp); -} - -void Pitvplot3(NODE arg, Obj *rp) -{ - itvplot_main3(arg, rp); -} - -void Pitvplot4(NODE arg, Obj *rp) -{ - itvplot_main4(arg, rp); -} - -// NORMAL type -void itvplot_main1(NODE arg, Obj *rp) -{ +void conplot_main(NODE arg,int is_memory,Obj *rp){ Q m2,p2,w300,s_id; NODE defrange; - LIST xrange,yrange,range[2],list,geom; - VL vl,vl0; - V v[2],av[2]; - int stream,ri,i,sign; - P poly; - P var; - NODE n,n0; - STRING fname,wname; - Obj t; - - STOQ(-2,m2); STOQ(2,p2); - MKNODE(n,p2,0); MKNODE(defrange,m2,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = v[1] = 0; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) stream = 0; - else - switch ( OID(BDY(arg)) ) { - case O_P: - poly = (P)BDY(arg); - get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 2 ) error("itvplot : invalid argument"); - else v[i++] = vl0->v; - break; - case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 1 ) error("itvplot : invalid argument"); - else range[ri++] = list; - else geom = list; - break; - case O_N: - stream = QTOS((Q)BDY(arg)); break; - case O_STR: - wname = (STRING)BDY(arg); break; - default: - error("itvplot : invalid argument"); break; - } - if ( !poly ) error("itvplot : invalid argument"); - switch ( ri ) { - case 0: - if ( !v[1] ) error("itvplot : please specify all variables"); - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - break; - case 1: - if ( !v[1] ) error("itvplot : please specify all variables"); - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) { - xrange = range[0]; - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - } else if ( v[1] == av[0] ) { - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - yrange = range[0]; - } else error("itvplot : invalid argument"); - break; - case 2: - av[0] = VR((P)BDY(BDY(range[0]))); - av[1] = VR((P)BDY(BDY(range[1]))); - if ( ((v[0] == av[0]) && (!v[1] || v[1] == av[1])) || - ((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { - xrange = range[0]; yrange = range[1]; - } else error("itvplot : invalid argument"); - break; - default: - error("itvplot : cannot happen"); break; - } - stream = validate_ox_plot_stream(stream); - STOQ(stream,s_id); - if ( !geom ) { - STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); - } - MKSTR(fname,"itvplot1"); - arg = mknode(8,s_id,fname,poly,xrange,yrange,0,geom,wname); - Pox_rpc(arg,&t); - *rp = (Obj)s_id; -} - -// TRANSFER TYPE -void itvplot_main2(NODE arg, Obj *rp) -{ - Q m2,p2,w300,s_id; - NODE defrange; - LIST xrange,yrange,range[2],list,geom; - VL vl,vl0; - V v[2],av[2]; - int stream,ri,i,sign; - P poly; - P var; - NODE n,n0; - STRING fname,wname; - Obj t; - - STOQ(-2,m2); STOQ(2,p2); - MKNODE(n,p2,0); MKNODE(defrange,m2,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = v[1] = 0; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) - stream = 0; - else - switch ( OID(BDY(arg)) ) { - case O_P: - poly = (P)BDY(arg); - get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 2 ) - error("itvplot : invalid argument"); - else - v[i++] = vl0->v; - break; - case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 1 ) - error("itvplot : invalid argument"); - else - range[ri++] = list; - else - geom = list; - break; - case O_N: - stream = QTOS((Q)BDY(arg)); break; - case O_STR: - wname = (STRING)BDY(arg); break; - default: - error("itvplot : invalid argument"); break; - } - if ( !poly ) - error("itvplot : invalid argument"); - switch ( ri ) { - case 0: - if ( !v[1] ) - error("itvplot : please specify all variables"); - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - break; - case 1: - if ( !v[1] ) - error("itvplot : please specify all variables"); - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) { - xrange = range[0]; - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - } else if ( v[1] == av[0] ) { - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - yrange = range[0]; - } else - error("itvplot : invalid argument"); - break; - case 2: - av[0] = VR((P)BDY(BDY(range[0]))); - av[1] = VR((P)BDY(BDY(range[1]))); - if ( ((v[0] == av[0]) && (!v[1] || v[1] == av[1])) || - ((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { - xrange = range[0]; yrange = range[1]; - } else - error("itvplot : invalid argument"); - break; - default: - error("itvplot : cannot happen"); break; - } - stream = validate_ox_plot_stream(stream); - STOQ(stream,s_id); - if ( !geom ) { - STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); - } - MKSTR(fname,"itvplot2"); - arg = mknode(8,s_id,fname,poly,xrange,yrange,0,geom,wname); - Pox_rpc(arg,&t); - *rp = (Obj)s_id; -} -// RECURSION TYPE -void itvplot_main3(NODE arg, Obj *rp) -{ - Q m2,p2,w300,s_id, itvsize; - NODE defrange; - LIST xrange,yrange,range[2],list,geom; - VL vl,vl0; - V v[2],av[2]; - int stream,ri,i,sign; - P poly; - P var; - NODE n,n0; - STRING fname,wname; - Obj t; - - STOQ(-2,m2); STOQ(2,p2); - MKNODE(n,p2,0); MKNODE(defrange,m2,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = v[1] = 0; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) - stream = 0; - else - switch ( OID(BDY(arg)) ) { - case O_P: - poly = (P)BDY(arg); - get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 2 ) - error("itvplot : invalid argument"); - else - v[i++] = vl0->v; - break; - case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 1 ) - error("itvplot : invalid argument"); - else - range[ri++] = list; - else - geom = list; - break; - case O_N: - stream = QTOS((Q)BDY(arg)); break; - case O_STR: - wname = (STRING)BDY(arg); break; - default: - error("itvplot : invalid argument"); break; - } - if ( !poly ) - error("itvplot : invalid argument"); - switch ( ri ) { - case 0: - if ( !v[1] ) - error("itvplot : please specify all variables"); - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - break; - case 1: - if ( !v[1] ) - error("itvplot : please specify all variables"); - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) { - xrange = range[0]; - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - } else if ( v[1] == av[0] ) { - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - yrange = range[0]; - } else - error("itvplot : invalid argument"); - break; - case 2: - av[0] = VR((P)BDY(BDY(range[0]))); - av[1] = VR((P)BDY(BDY(range[1]))); - if ( ((v[0] == av[0]) && (!v[1] || v[1] == av[1])) || - ((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { - xrange = range[0]; yrange = range[1]; - } else - error("itvplot : invalid argument"); - break; - default: - error("itvplot : cannot happen"); break; - } - stream = validate_ox_plot_stream(stream); - STOQ(stream,s_id); - if ( !geom ) { - STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); - } - MKSTR(fname,"itvplot3"); - STOQ(Itvplot, itvsize); - arg = mknode(9,s_id,fname,poly,xrange,yrange,0,geom,wname,itvsize); - Pox_rpc(arg,&t); - *rp = (Obj)s_id; -} -// RECURSION and TRANSFER TYPE -void itvplot_main4(NODE arg, Obj *rp) -{ - Q m2,p2,w300,s_id,itvsize; - NODE defrange; - LIST xrange,yrange,range[2],list,geom; - VL vl,vl0; - V v[2],av[2]; - int stream,ri,i,sign; - P poly; - P var; - NODE n,n0; - STRING fname,wname; - Obj t; - - STOQ(-2,m2); STOQ(2,p2); - MKNODE(n,p2,0); MKNODE(defrange,m2,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = v[1] = 0; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) - stream = 0; - else - switch ( OID(BDY(arg)) ) { - case O_P: - poly = (P)BDY(arg); - get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 2 ) - error("itvplot : invalid argument"); - else - v[i++] = vl0->v; - break; - case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 1 ) - error("itvplot : invalid argument"); - else - range[ri++] = list; - else - geom = list; - break; - case O_N: - stream = QTOS((Q)BDY(arg)); break; - case O_STR: - wname = (STRING)BDY(arg); break; - default: - error("itvplot : invalid argument"); break; - } - if ( !poly ) - error("itvplot : invalid argument"); - switch ( ri ) { - case 0: - if ( !v[1] ) - error("itvplot : please specify all variables"); - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - break; - case 1: - if ( !v[1] ) - error("itvplot : please specify all variables"); - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) { - xrange = range[0]; - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - } else if ( v[1] == av[0] ) { - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - yrange = range[0]; - } else - error("itvplot : invalid argument"); - break; - case 2: - av[0] = VR((P)BDY(BDY(range[0]))); - av[1] = VR((P)BDY(BDY(range[1]))); - if ( ((v[0] == av[0]) && (!v[1] || v[1] == av[1])) || - ((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { - xrange = range[0]; yrange = range[1]; - } else - error("itvplot : invalid argument"); - break; - default: - error("itvplot : cannot happen"); break; - } - /* itvplot in ox_plot requires - [s_id (Q), - formula (Obj), - xrange=[x,xmin,xmax] (LIST), - yrange=[y,ymin,ymax] (LIST), - zrange=0, - geom=[xsize,ysize] (LIST), - wname=name (STRING)] - */ - - stream = validate_ox_plot_stream(stream); - STOQ(stream,s_id); - if ( !geom ) { - STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); - } - MKSTR(fname,"itvplot4"); - STOQ(Itvplot, itvsize); - arg = mknode(9,s_id,fname,poly,xrange,yrange,0,geom,wname,itvsize); - Pox_rpc(arg,&t); - *rp = (Obj)s_id; -} - -void Pineqn(NODE arg, Obj *rp) -{ - ineqn_main(arg, 1, rp); -} - -void ineqn_main(NODE arg,int is_memory, Obj *rp) -{ - Q s_id, m2, p2, w300, color; - NODE defrange, n, n0; - P poly, var; - VL vl, vl0; - V v[2], av[2]; - LIST xrange, yrange, range[2], list, geom; - int stream, ri, i,sign; - STRING fname,wname; - - STOQ(-2,m2); STOQ(2,p2); - MKNODE(n,p2,0); MKNODE(defrange,m2,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = v[1] = 0; - /* get polynomial */ - if ( !(OID(BDY(arg)) == O_P || (OID(BDY(arg)) ==O_R)) ) - error("ineqn : first argument must be a polynormial"); - else { - poly = (P)BDY(arg); - /* get vars */ - get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 2 ) - error( - "ineqn : first argument must be a univariate or bivariate polynormial"); - else - v[i++] = vl0->v; - } - /* get color */ - arg = NEXT(arg); - if ( OID(BDY(arg)) != O_N ) - error("ineqn : second argument must be color code"); - else - color = (Q)BDY(arg); - /* other argument is optional */ - arg = NEXT(arg); - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) - stream = 0; - else - switch ( OID(BDY(arg)) ) { - case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 1 ) - error("ineqn : invalid list argument"); - else - range[ri++] = list; - else - geom = list; - break; - case O_N: - stream = QTOS((Q)BDY(arg)); break; - case O_STR: - wname = (STRING)BDY(arg); break; - default: - error("ineqn : invalid argument"); break; - } - switch ( ri ) { - case 0: - if ( !v[1] ) error("ineqn : please specify all variables"); - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - break; - case 1: - if ( !v[1] ) error("ineqn : please specify all variables"); - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) { - xrange = range[0]; - MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); - } else if ( v[1] == av[0] ) { - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); - yrange = range[0]; - } else error("ineqn : invalid argument"); - break; - case 2: - av[0] = VR((P)BDY(BDY(range[0]))); - av[1] = VR((P)BDY(BDY(range[1]))); - if ( ((v[0] == av[0]) && (!v[1] || v[1] == av[1])) || - ((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { - xrange = range[0]; yrange = range[1]; - } else error("ineqn : invalid argument"); - break; - default: - error("ineqn : cannot happen"); break; - } - - stream = validate_ox_plot_stream(stream); - STOQ(stream,s_id); - if ( !geom ) { - STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); - } - MKSTR(fname, "ineqn"); - arg = mknode(8, s_id, fname, poly, color, xrange, yrange, geom, wname); - Pox_cmo_rpc(arg,rp); - *rp = (Obj)s_id; -} - -void Pineqnor(NODE arg, Obj *rp) -{ - ineqn_main_op(arg, 7, 0, rp); -} - -void Pineqnand(NODE arg, Obj *rp) -{ - ineqn_main_op(arg, 1, 0, rp); -} - -void Pineqnxor(NODE arg, Obj *rp) -{ - ineqn_main_op(arg, 6, 0, rp); -} - -void Pmemory_ineqn(NODE arg, Obj *rp) -{ - ineqn_main(arg, 1, rp); -} - -void ineqn_main_op(NODE arg, int op, int is_memory, Obj *rp) -{ - Q s_id, index, color, op_code; - P poly; - STRING fname; - Obj t; - - poly = (P)ARG0(arg); - color = (Q)ARG1(arg); - s_id = (Q)ARG2(arg); - index = (Q)ARG3(arg); - STOQ( op, op_code); - MKSTR(fname, "ineqnover"); - arg = mknode(6, s_id, fname, index, poly, color, op_code); - Pox_rpc(arg, &t); - *rp = (Obj)s_id; -} - -void Pobj_cp(NODE arg, Obj *rp) -{ - Q sysid, index_A, index_B, op_code; - STRING fname; - Obj t; - - sysid = (Q)ARG0(arg); - index_A = (Q)ARG1(arg); - index_B = (Q)ARG2(arg); - op_code = (Q)ARG3(arg); - MKSTR(fname, "objcp"); - arg = mknode(5, sysid, fname, index_A, index_B, op_code); - Pox_rpc(arg, &t); - *rp = (Obj)sysid; -} -#endif - -void Pconplot(NODE arg,Obj *rp) -{ - conplot_main(arg, 0, rp); -} - -void Pmemory_conplot(NODE arg,Obj *rp) -{ - conplot_main(arg,1,rp); -} - -void conplot_main(NODE arg,int is_memory,Obj *rp) -{ - Q m2,p2,w300,s_id; - NODE defrange; LIST xrange,yrange,zrange,range[3],list,geom; VL vl,vl0; V v[2],av[2]; @@ -1052,90 +432,95 @@ void conplot_main(NODE arg,int is_memory,Obj *rp) arg = mknode(1,s_id); Pox_pop_cmo(arg,rp); } else { - MKSTR(fname,"plot"); + MKSTR(fname,CONPLOT); arg = mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom,wname); Pox_rpc(arg,&t); *rp = (Obj)s_id; } } -void Pplot(NODE arg,Obj *rp) -{ - plot_main(arg,0,rp); -} +void Pplot(NODE arg,Obj *rp){plot_main(arg,0,PLOT,rp);} +void Pmemory_plot(NODE arg,Obj *rp){plot_main(arg,1,PLOT,rp);} -void Pmemory_plot(NODE arg,Obj *rp) -{ - plot_main(arg,1,rp); -} +int plot_by_bigfloat; -void plot_main(NODE arg,int is_memory,Obj *rp) -{ +void plot_main(NODE arg,int is_memory,char *fn,Obj *rp){ Q m2,p2,w300,s_id; NODE defrange; LIST xrange,range[1],list,geom; VL vl,vl0; V v[1],av[1]; - int stream,ri,i; - P poly; + int stream,ri,i,found_f; + Obj func; P var; NODE n,n0; STRING fname,wname; Obj t; + Q prec; - STOQ(-2,m2); STOQ(2,p2); - MKNODE(n,p2,0); MKNODE(defrange,m2,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = 0; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) - stream = 0; + STOQ(-2,m2);STOQ(2,p2); + MKNODE(n,p2,0);MKNODE(defrange,m2,n); + func=0;vl=0;geom=0;wname=0;stream=-1;ri=0; + v[0]=0; + found_f = 0; + for(;arg;arg=NEXT(arg) ) { + if(!BDY(arg) ) + if ( !found_f ) { + makevar("x",&var); + v[0] = VR(var); + found_f = 1; + } else + stream=0; else - switch ( OID(BDY(arg)) ) { + switch ( OID(BDY(arg)) ){ case O_P: case O_R: - poly = (P)BDY(arg); - get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 1 ) - error("ifplot : invalid argument"); - else - v[i++] = vl0->v; - if ( i != 1 ) - error("ifplot : invalid argument"); + func = (Obj)BDY(arg); + get_vars_recursive(func,&vl); + for(vl0=vl, i=0;vl0;vl0=NEXT(vl0) ){ + if(vl0->v->attr==(pointer)V_IND ){ + if(i >= 1 ) error("plot : function must be univariate"); + else v[i++]=vl0->v; + } + } + found_f = 1; break; case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 0 ) - error("plot : invalid argument"); - else - range[ri++] = list; - else - geom = list; + list=(LIST)BDY(arg); + if(OID(BDY(BDY(list)))==O_P ){ + if(ri > 0 ) error("plot : too many intervals"); + else range[ri++]=list; + } else geom=list; break; case O_N: - stream = QTOS((Q)BDY(arg)); break; + if ( !found_f ) { + func = (Obj)BDY(arg); + makevar("x",&var); + v[0] = VR(var); + found_f = 1; + } else + stream=QTOS((Q)BDY(arg)); + break; case O_STR: - wname = (STRING)BDY(arg); break; + wname=(STRING)BDY(arg);break; default: - error("plot : invalid argument"); break; + error("plot : invalid argument");break; } - if ( !poly ) + } + if(!found_f ) error("plot : invalid argument"); - switch ( ri ) { + switch ( ri ){ case 0: - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); + MKV(v[0],var);MKNODE(n,var,defrange);MKLIST(xrange,n); break; case 1: - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) - xrange = range[0]; + av[0]=VR((P)BDY(BDY(range[0]))); + if(!func || NUM(func) || v[0]==av[0] ) + xrange=range[0]; else error("plot : invalid argument"); break; default: - error("plot : cannot happen"); break; + error("plot : cannot happen");break; } /* conplot in ox_plot requires [s_id (Q), @@ -1146,227 +531,581 @@ void plot_main(NODE arg,int is_memory,Obj *rp) geom=[xsize,ysize] (LIST), wname=name (STRING)] */ - stream = validate_ox_plot_stream(stream); + stream=validate_ox_plot_stream(stream); STOQ(stream,s_id); - if ( !geom ) { + if(!geom ){ STOQ(300,w300); - MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); + MKNODE(n0,w300,0);MKNODE(n,w300,n0);MKLIST(geom,n); } - if ( is_memory ) { - MKSTR(fname,"memory_plot"); - arg = mknode(7,s_id,fname,poly,xrange,NULLP,NULLP,geom); + if(plot_by_bigfloat) + STOQ(plot_by_bigfloat,prec); + else + prec = 0; + if(is_memory ){ + MKSTR(fname,MEMORY_PLOT); + arg=mknode(8,s_id,fname,func,xrange,NULLP,NULLP,geom,prec); Pox_rpc(arg,&t); - arg = mknode(1,s_id); + arg=mknode(1,s_id); Pox_pop_cmo(arg,rp); } else { - MKSTR(fname,"plot"); - arg = mknode(8,s_id,fname,poly,xrange,NULLP,NULLP,geom,wname); + MKSTR(fname,fn); + arg=mknode(9,s_id,fname,func,xrange,NULLP,NULLP,geom,wname,prec); Pox_rpc(arg,&t); - *rp = (Obj)s_id; + *rp=(Obj)s_id; } } #define Pi 3.14159265358979323846264 -void Ppolarplot(NODE arg,Q *rp) -{ +void Ppolarplot(NODE arg,Obj *rp){ Q m2,p2,w300,s_id; - NODE defrange; - LIST zrange,range[1],list,geom; + NODE defrange,n,n0; + LIST zrange,range[1],geom,list; VL vl,vl0; V v[1],av[1]; int stream,ri,i; - P poly; - P var; - NODE n,n0; + P poly,var; STRING fname,wname; Real pi2; Obj t; - MKReal(2*Pi,pi2); - MKNODE(n,pi2,0); MKNODE(defrange,0,n); - poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; - v[0] = 0; - for ( ; arg; arg = NEXT(arg) ) - if ( !BDY(arg) ) - stream = 0; - else - switch ( OID(BDY(arg)) ) { - case O_P: case O_R: - poly = (P)BDY(arg); - get_vars_recursive((Obj)poly,&vl); - for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) - if ( vl0->v->attr == (pointer)V_IND ) - if ( i >= 1 ) - error("polarplot : invalid argument"); - else - v[i++] = vl0->v; - if ( i != 1 ) - error("polarplot : invalid argument"); - break; - case O_LIST: - list = (LIST)BDY(arg); - if ( OID(BDY(BDY(list))) == O_P ) - if ( ri > 0 ) - error("polarplot : invalid argument"); - else - range[ri++] = list; - else - geom = list; - break; - case O_N: - stream = QTOS((Q)BDY(arg)); break; - case O_STR: - wname = (STRING)BDY(arg); break; - default: - error("polarplot : invalid argument"); break; - } - if ( !poly ) - error("polarplot : invalid argument"); - switch ( ri ) { - case 0: - MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(zrange,n); + MKReal(2*Pi,pi2);MKNODE(n,pi2,0); MKNODE(defrange,0,n); + poly=0;vl=0;geom=0;wname=0;stream=-1;ri=0;v[0]=0; + for(;arg;arg=NEXT(arg)){ + if(!BDY(arg)) stream=0; + else switch(OID(BDY(arg))){ + case O_P: case O_R://formula + poly=(P)BDY(arg); + get_vars_recursive((Obj)poly,&vl); + for(vl0=vl,i=0;vl0;vl0=NEXT(vl0)) + if(vl0->v->attr==(pointer)V_IND) + if(i>=1)error("polarplot : invalid argument"); + else v[i++]=vl0->v; + if(i!=1)error("polarplot : invalid argument"); break; - case 1: - av[0] = VR((P)BDY(BDY(range[0]))); - if ( v[0] == av[0] ) - zrange = range[0]; - else - error("polarplot : invalid argument"); + case O_LIST://range,geomerty + list=(LIST)BDY(arg); + if(OID(BDY(BDY(list)))==O_P) + if(ri>0)error("polarplot : invalid argument"); + else range[ri++]=list; + else geom=list; break; + case O_N: + stream=QTOS((Q)BDY(arg)); + break; + case O_STR://wname + wname=(STRING)BDY(arg); + break; default: - error("polarplot : cannot happen"); break; + error("polarplot : invalid argument"); + break; + } } - stream = validate_ox_plot_stream(stream); + //formular check + if(!poly)error("polarplot : invalid argument"); + switch (ri){ + case 0: + MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(zrange,n); + break; + case 1: + av[0]=VR((P)BDY(BDY(range[0]))); + if(v[0]==av[0]) zrange = range[0]; + else error("polarplot : invalid argument"); + break; + default: + error("polarplot : cannot happen"); + break; + } + stream=validate_ox_plot_stream(stream); STOQ(stream,s_id); - if ( !geom ) { + if(!geom){ STOQ(300,w300); MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); } - MKSTR(fname,"plot"); - arg = mknode(8,s_id,fname,poly,NULLP,NULLP,zrange,geom,wname); + MKSTR(fname,POLARPLOT); + arg=mknode(8,s_id,fname,poly,NULLP,NULLP,zrange,geom,wname); Pox_rpc(arg,&t); - *rp = s_id; + *rp=s_id; } -void Pplotover(NODE arg,Q *rp) -{ - Q s_id,index; +void Pplotover(NODE arg,Obj *rp){ + Q s_id,w_id,color; P poly; STRING fname; Obj t; - - poly = (P)ARG0(arg); - s_id = (Q)ARG1(arg); - index = (Q)ARG2(arg); - MKSTR(fname,"plotover"); - if ( argc(arg) == 4 ) - arg = mknode(5,s_id,fname,index,poly,(Q)ARG3(arg)); - else - arg = mknode(4,s_id,fname,index,poly); + poly=(P)ARG0(arg); + s_id=(Q)ARG1(arg); + w_id=(Q)ARG2(arg); + if(argc(arg)==4)color=(Q)ARG3(arg); + else color=0; + MKSTR(fname,PLOTOVER); + arg=mknode(5,s_id,fname,w_id,poly,color); Pox_rpc(arg,&t); - *rp = s_id; + *rp=(Obj)s_id; } -/* arg = [x,y,r,s_id,index] */ - -void Pdrawcircle(NODE arg,Q *rp) -{ +void Pdrawcircle(NODE arg,Obj *rp){ Q s_id,index; - Obj x,y,r,c; + Obj x,y,r,c,t; STRING fname; NODE n; LIST pos; - Obj t; - x = (Obj)ARG0(arg); - y = (Obj)ARG1(arg); - r = (Obj)ARG2(arg); - c = (Obj)ARG3(arg); - s_id = (Q)ARG4(arg); - index = (Q)ARG5(arg); - MKSTR(fname,"drawcircle"); - n = mknode(3,x,y,r,c); MKLIST(pos,n); - arg = mknode(5,s_id,fname,index,pos,c); + x=(Obj)ARG0(arg); + y=(Obj)ARG1(arg); + r=(Obj)ARG2(arg); + c=(Obj)ARG3(arg); + s_id=(Q)ARG4(arg); + index=(Q)ARG5(arg); + MKSTR(fname,DRAWCIRCLE); + n=mknode(3,x,y,r,c); + MKLIST(pos,n); + arg=mknode(5,s_id,fname,index,pos,c); Pox_rpc(arg,&t); - *rp = s_id; + *rp=(Obj)s_id; } -/* draw_obj(s_id,cindex,point|line); point = [x,y], line = [xa,ya,xb,yb] */ -void Pdraw_obj(NODE arg,Q *rp) -{ +void Pdraw_obj(NODE arg,Obj *rp){ static STRING fname; Q s_id,index; LIST obj; Obj t; - if ( !fname ) { - MKSTR(fname,"draw_obj"); - } - s_id = (Q)ARG0(arg); - index = (Q)ARG1(arg); - obj = (LIST)ARG2(arg); - /* ARG3(arg) = color */ - if ( argc(arg) == 4 ) - arg = mknode(5,s_id,fname,index,obj,ARG3(arg)); - else - arg = mknode(4,s_id,fname,index,obj); + if(!fname)MKSTR(fname,DRAW_OBJ); + s_id=(Q)ARG0(arg); + index=(Q)ARG1(arg); + obj=(LIST)ARG2(arg); + // ARG3(arg)=color + if(argc(arg)==4) arg=mknode(5,s_id,fname,index,obj,ARG3(arg)); + else arg=mknode(4,s_id,fname,index,obj); Pox_cmo_rpc(arg,&t); - *rp = s_id; + *rp=(Obj)s_id; } -/* draw_string(s_id,cindex,pos,string,[,color]); pos=[x,y] */ -void Pdraw_string(NODE arg,Q *rp) -{ +void Pdraw_string(NODE arg,Obj *rp){ static STRING fname; STRING str; Q s_id,index; LIST pos; Obj t; - if ( !fname ) { - MKSTR(fname,"draw_string"); - } - s_id = (Q)ARG0(arg); - index = (Q)ARG1(arg); - pos = (LIST)ARG2(arg); - str = (STRING)ARG3(arg); - /* ARG4(arg) = color */ - if ( argc(arg) == 5 ) - arg = mknode(6,s_id,fname,index,pos,str,ARG4(arg)); - else - arg = mknode(5,s_id,fname,index,pos,str); + if(!fname)MKSTR(fname,DRAW_STRING); + s_id=(Q)ARG0(arg); + index=(Q)ARG1(arg); + pos=(LIST)ARG2(arg); + str=(STRING)ARG3(arg); + // ARG4(arg)=color + if(argc(arg)==5) arg=mknode(6,s_id,fname,index,pos,str,ARG4(arg)); + else arg=mknode(5,s_id,fname,index,pos,str); Pox_cmo_rpc(arg,&t); - *rp = s_id; + *rp=(Obj)s_id; } -void Pclear_canvas(NODE arg,Q *rp) -{ +void Pclear_canvas(NODE arg,Obj *rp){ static STRING fname; Q s_id,index; Obj t; - if ( !fname ) { - MKSTR(fname,"clear_canvas"); - } - s_id = (Q)ARG0(arg); - index = (Q)ARG1(arg); - arg = mknode(3,s_id,fname,index); + if(!fname) MKSTR(fname,CLEAR_CANVAS); + s_id=(Q)ARG0(arg); + index=(Q)ARG1(arg); + arg=mknode(3,s_id,fname,index); Pox_cmo_rpc(arg,&t); - *rp = s_id; + *rp=(Obj)s_id; } +//****************************ifplotNG +/* + * name ifplot,ineqn,ineqnor,ineqnand,ineqnxor,conplot + * type + * D:sign character using double + * Q:sign character use rational + * B:Boundary character use sturm theorem + */ +void ListCheck(char * head,LIST list){ + int i; + NODE n; + if(!list){ + printf("%s zero \n",head); + return; + } + for(i=0,n=(NODE)BDY(list);n;i++,n=NEXT(n)); + printf("%s length %d\n",head,i); + for(i=0,n=(NODE)BDY(list);n;i++,n=NEXT(n)){ + if(!BDY(n))printf("%d 0\n",i); + else if(OID(BDY(n))==O_P) printf("%d poly\n",i); + else if(OID(BDY(n))==O_R) printf("%d real\n",i); + else if(OID(BDY(n))==O_N) printf("%d %d\n",i,QTOS((Q)BDY(n))); + } +} -#if 0 -void Parrayplot(NODE arg,Obj *rp) -{ - int s; - int id; +void PifplotD(NODE arg,Obj *rp){ifplot_mainNG(arg,IFPLOTD,rp);} +void PifplotQ(NODE arg,Obj *rp){ifplot_mainNG(arg,IFPLOTQ,rp);} +void PifplotB(NODE arg,Obj *rp){ifplot_mainNG(arg,IFPLOTB,rp);} - if ( ID((Obj)ARG0(arg)) == O_VECT && ID((Obj)ARG1(arg)) == O_LIST ) { - s = current_s; - gensend(s,C_APLOT,0); - gensend(s,C_OBJ,ARG0(arg)); gensend(s,C_OBJ,ARG1(arg)); genflush(s); +void PconplotD(NODE arg,Obj *rp){ifplot_mainNG(arg,CONPLOTD,rp);} +void PconplotQ(NODE arg,Obj *rp){ifplot_mainNG(arg,CONPLOTQ,rp);} +void PconplotB(NODE arg,Obj *rp){ifplot_mainNG(arg,CONPLOTB,rp);} + +void PineqnD(NODE arg,Obj *rp){ifplot_mainNG(arg,INEQND,rp);} +void PineqnQ(NODE arg,Obj *rp){ifplot_mainNG(arg,INEQNQ,rp);} +void PineqnB(NODE arg,Obj *rp){ifplot_mainNG(arg,INEQNB,rp);} + +#if defined(INTERVAL) +void Pitvifplot(NODE arg,Obj *rp){ifplot_mainNG(arg,ITVIFPLOT,rp);} +#endif + +void PineqnorD(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNDOR,rp);} +void PineqnorQ(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNQOR,rp);} +void PineqnorB(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNBOR,rp);} + +void PineqnandD(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNDAND,rp);} +void PineqnandQ(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNQAND,rp);} +void PineqnandB(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNBAND,rp);} + +void PineqnxorD(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNDXOR,rp);} +void PineqnxorQ(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNQXOR,rp);} +void PineqnxorB(NODE arg,Obj *rp){plotover_mainNG(arg,INEQNBXOR,rp);} + +void PplotoverD(NODE arg,Obj *rp){plotover_mainNG(arg,PLOTOVERD,rp);} +void PplotoverQ(NODE arg,Obj *rp){plotover_mainNG(arg,PLOTOVERQ,rp);} +void PplotoverB(NODE arg,Obj *rp){plotover_mainNG(arg,PLOTOVERB,rp);} + +void Pobj_cp(NODE arg,Obj *rp){ + //copy canvas + //1:and,3:copy,6:xor,7:or + Q sysid,index_A,index_B,op_code; + STRING fname; + Obj t; + sysid=(Q)ARG0(arg); + index_A=(Q)ARG1(arg); + index_B=(Q)ARG2(arg); + op_code=(Q)ARG3(arg); + MKSTR(fname,OBJ_CP); + arg=mknode(5,sysid,fname,index_A,index_B,op_code); + Pox_rpc(arg,&t); + *rp=(Obj)sysid; +} + +void PpolarplotD(NODE arg,Obj *rp){ + Q m2,p2,w300,defstep,s_id,color; + NODE defrange,n,n0,n1,n2; + LIST range,geom,list[2]; + VL vl,vl0; + V v[1],av; + int stream,ri,i,len,iNo,lNo,vNo,sNo,pfine,findG; + P poly,var; + STRING fname,wname; + Real pi2; + Obj t,frst,sec,thr; + char ebuf[BUFSIZ]; + + iNo=lNo=sNo=findG=0;pfine=TRUE; + poly=0;vl=0;geom=0;wname=0;color=0;stream=-1;ri=0;v[0]=0; + for(;arg;arg=NEXT(arg)){ + if(!BDY(arg)) iNo++; + else switch(OID(BDY(arg))){ + case O_P: case O_R://formular + poly=(P)BDY(arg); + get_vars_recursive((Obj)poly,&vl); + for(vl0=vl,vNo=0;vl0;vl0=NEXT(vl0)){ + if(vl0->v->attr==(pointer)V_IND){ + if(vNo>=1){ + sprintf(ebuf,"%s : invalaid argument",POLARPLOT); + error(ebuf); + } else v[vNo++]=vl0->v; + } + } + if(vNo!=1){ + sprintf(ebuf,"%s : only uni-variate formular",POLARPLOT); + error(ebuf); + } + break; + case O_N://color,id,idx,division + switch (iNo){ + case 0://color arg + color=(Q)BDY(arg); + iNo++; + break; + case 1://stream arg + stream=QTOS((Q)BDY(arg)); + iNo++; + break; + default://error + sprintf(ebuf,"%s : invalid number arguments",POLARPLOT); + error(ebuf); + break; + } + break; + case O_LIST://range,geomerty + if(lNo<2)list[lNo++]=(LIST)BDY(arg); + else { + sprintf(ebuf,"%s : invalid list argument",POLARPLOT); + error(ebuf); + } + break; + case O_STR://wname + wname=(STRING)BDY(arg); + sNo++; + break; + default: + break; + } } - genrecv(s,&id,rp); + //formular check + if(!poly){ + sprintf(ebuf,"%s : invalid plot argument",POLARPLOT); + error(ebuf); + } + //vars check + get_vars_recursive((Obj)poly,&vl); + for(vl0=vl,vNo=0;vl0;vNo++,vl0=NEXT(vl0)){ + if(vl0->v->attr==(pointer)V_IND){ + if(vNo>=2){ + sprintf(ebuf,"%s : invalid plot argument",POLARPLOT); + error(ebuf); + } else v[vNo]=vl0->v; + } + } + //list check + for(i=0;iv->attr==(pointer)V_IND){ + if(vNo>=2){ + sprintf(ebuf,"%s : invalid plot argument",fn); + error(ebuf); + } else v[vNo++]=vl0->v; + } + } + //list check + xrange=yrange=zrange=geom=0;frst=sec=thr=0; + for(i=0;i4){ + sprintf(ebuf,"%s : invalid list length",fn); + error(ebuf); + } + if(OID(BDY(BDY(list[i])))!=O_P){ + // First list value is number + for(len=0,n=(NODE)BDY(list[i]);n;len++,n=NEXT(n)){ + if(len==0)frst=BDY(n); + else if(len==1)sec=BDY(n); + else thr=BDY(n); + } + switch(len){ + case 2: + if(!strcmp(fn,CONPLOT)){ + if(thr==0)thr=(Obj)mxgc; + MKNODE(n,thr,0);MKNODE(n1,sec,n);MKNODE(n2,frst,n1);MKLIST(zrange,n2); + } else geom=list[i]; + break; + case 3: + zrange=list[i]; + break; + case 0: + case 1: + default: + sprintf(ebuf,"%s : invalid list length",fn); + error(ebuf); + break; + } + } else { + //xrange,yrange + av=VR((P)BDY(BDY(list[i]))); + if(v[0]==av)xrange=list[i]; + else if(v[1]==av)yrange=list[i]; + else { + MKLIST(zrange,NEXT(BDY(list[i]))); + } + } + } + //set default + if(!xrange){ + MKV(v[0],var);MKNODE(n,var,defrange);MKLIST(xrange,n); + } + if(!yrange){ + MKV(v[1],var);MKNODE(n,var,defrange);MKLIST(yrange,n); + } + if(!geom){ + STOQ(300,w300);MKNODE(n0,w300,0);MKNODE(n,w300,n0);MKLIST(geom,n); + } + if(!(strcmp(fn,CONPLOTD)&strcmp(fn,CONPLOTQ)&strcmp(fn,CONPLOTB))&!zrange){ + MKNODE(n,mxgc,0);MKNODE(n1,m2,n);MKNODE(n2,m2,n1);MKLIST(zrange,n2); + } + /*new ifplot in ox_plot requires + [s_id (Q), + formula (Obj), + color (Q), + geom=[xsize,ysize] (LIST optional), + xrange=[x,xmin,xmax] (LIST optional), + yrange=[y,ymin,ymax] (LIST optional), + zrange=[zmin,zmax,zstep] (LIST optional), + wname=name (STRING optional)], + itvstep (Q) if ITVIFPLOT */ + stream=validate_ox_plot_stream(stream); + STOQ(stream,s_id); + MKSTR(fname,fn); + arg=mknode(9,s_id,fname,poly,color,xrange,yrange,zrange,geom,wname); + Pox_rpc(arg,&t); + *rp=(Obj)s_id; +} + +void plotover_mainNG(NODE arg,char *fn,Obj *rp){ + Q s_id,w_id,color; + P poly; + STRING fname; + Obj t; + int iNo,pfine,sfine; + char ebuf[BUFSIZ]; + + pfine=sfine=TRUE; + iNo=0;poly=0;color=s_id=w_id=0; + for(;arg;arg=NEXT(arg)){ + if(!BDY(arg)) iNo++; + else switch(OID(BDY(arg))){ + case O_P://formular + if(pfine){ + poly=(P)BDY(arg); + pfine=FALSE; + } else { + sprintf(ebuf,"%s : to many plot arguments",fn); + error(ebuf); + } + break; + case O_N://color,s_id,w_id + switch(iNo){ + case 0://color arg + color=(Q)BDY(arg); + iNo++; + break; + case 1://stream arg + s_id=(Q)BDY(arg); + iNo++; + break; + case 2://window arg + w_id=(Q)BDY(arg); + iNo++; + break; + default://error + sprintf(ebuf,"%s : to many numbers",fn); + error(ebuf); + break; + } + break; + default: + sprintf(ebuf,"%s : arguments type miss match",fn); + error(ebuf); + } + } + MKSTR(fname,fn); + //[s_id (Q), w_id (Q), formula (Obj), color (Q)] + arg=mknode(5,s_id,fname,w_id,poly,color); + Pox_rpc(arg,&t); + *rp=(Obj)s_id; +}