version 1.3, 2000/03/31 01:57:14 |
version 1.23, 2012/10/04 04:56:39 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/plot/plotf.c,v 1.2 2000/03/28 06:32:23 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/plotf.c,v 1.22 2011/08/10 04:51:58 saito Exp $ |
|
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
#include "ox.h" |
#include "ox.h" |
#include "ifplot.h" |
#include "ifplot.h" |
|
|
void Pifplot(), Pconplot(), Pplotover(), Pplot(), Parrayplot(), Pdrawcircle(); |
void Pifplot(NODE, Obj *), Pconplot(NODE, Obj *), Pplotover(NODE, Q *); |
|
void Pplot(NODE, Obj *), Parrayplot(NODE, Q *), Pdrawcircle(NODE, Q *); |
|
#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 *); |
|
#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[] = { |
struct ftab plot_tab[] = { |
{"ifplot",Pifplot,-7}, |
{"ifplot",Pifplot,-7}, |
|
{"memory_ifplot",Pmemory_ifplot,-6}, |
|
#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}, |
|
#endif |
{"conplot",Pconplot,-8}, |
{"conplot",Pconplot,-8}, |
|
{"memory_conplot",Pmemory_conplot,-7}, |
{"plot",Pplot,-6}, |
{"plot",Pplot,-6}, |
{"plotover",Pplotover,3}, |
{"memory_plot",Pmemory_plot,-5}, |
{"drawcircle",Pdrawcircle,5}, |
{"polarplot",Ppolarplot,-6}, |
/* |
{"plotover",Pplotover,-4}, |
{"arrayplot",Parrayplot,2}, |
{"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}, |
{0,0,0}, |
}; |
}; |
|
|
int current_s; |
void Popen_canvas(NODE arg,Q *rp) |
|
{ |
|
Q w300,s_id; |
|
LIST geom; |
|
int stream; |
|
NODE n,n0; |
|
STRING fname,wname; |
|
|
void Pifplot(arg,rp) |
geom = 0; wname = 0; stream = -1; |
NODE arg; |
for ( ; arg; arg = NEXT(arg) ) |
Obj *rp; |
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); |
|
STOQ(stream,s_id); |
|
if ( !geom ) { |
|
STOQ(300,w300); |
|
MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); |
|
} |
|
MKSTR(fname,"open_canvas"); |
|
arg = mknode(4,s_id,fname,geom,wname); |
|
Pox_cmo_rpc(arg,rp); |
|
*rp = s_id; |
|
} |
|
|
|
void Pifplot(NODE arg,Obj *rp) |
{ |
{ |
Q m2,p2,w300,sid; |
ifplot_main(arg,0,rp); |
|
} |
|
|
|
void Pmemory_ifplot(NODE arg,Obj *rp) |
|
{ |
|
ifplot_main(arg,1,rp); |
|
} |
|
|
|
void ifplot_main(NODE arg,int is_memory, Obj *rp) |
|
{ |
|
Q m2,p2,w300,s_id; |
NODE defrange; |
NODE defrange; |
LIST xrange,yrange,range[2],list,geom; |
LIST xrange,yrange,range[2],list,geom; |
VL vl,vl0; |
VL vl,vl0; |
V v[2],av[2]; |
V v[2],av[2]; |
int stream,ri,id,i; |
int stream,ri,i,sign; |
P poly; |
P poly; |
P var; |
P var; |
NODE n,n0; |
NODE n,n0; |
STRING fname,wname; |
STRING fname,wname; |
|
Obj t; |
|
|
STOQ(-2,m2); STOQ(2,p2); |
STOQ(-2,m2); STOQ(2,p2); |
MKNODE(n,p2,0); MKNODE(defrange,m2,n); |
MKNODE(n,p2,0); MKNODE(defrange,m2,n); |
poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; |
poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; |
|
v[0] = v[1] = 0; |
for ( ; arg; arg = NEXT(arg) ) |
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("ifplot : 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("ifplot : 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("ifplot : invalid argument"); break; |
|
} |
|
if ( !poly ) 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); |
|
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]; |
|
} 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]; |
|
} else error("ifplot : invalid argument"); |
|
break; |
|
default: |
|
error("ifplot : cannot happen"); break; |
|
} |
|
/* ifplot 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); |
|
} |
|
if ( is_memory ) { |
|
MKSTR(fname,"memory_plot"); |
|
arg = mknode(8,s_id,fname,poly,xrange,yrange,0,geom); |
|
Pox_rpc(arg,&t); |
|
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); |
|
Pox_rpc(arg,&t); |
|
*rp = (Obj)s_id; |
|
} |
|
} |
|
|
|
#if defined(INTERVAL) |
|
void Pitvplot(NODE arg, Obj *rp) |
|
{ |
|
itvplot_main(arg, 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) |
|
{ |
|
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) ) |
if ( !BDY(arg) ) |
stream = 0; |
stream = 0; |
else |
else |
switch ( OID(BDY(arg)) ) { |
switch ( OID(BDY(arg)) ) { |
case O_P: |
case O_P: |
poly = (P)BDY(arg); |
poly = (P)BDY(arg); |
get_vars_recursive(poly,&vl); |
get_vars_recursive((Obj)poly,&vl); |
for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) |
for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) |
if ( vl0->v->attr == V_IND ) |
if ( vl0->v->attr == (pointer)V_IND ) |
if ( i >= 2 ) |
if ( i >= 2 ) |
error("ifplot : invalid argument"); |
error("itvplot : invalid argument"); |
else |
else |
v[i++] = vl0->v; |
v[i++] = vl0->v; |
break; |
break; |
|
|
list = (LIST)BDY(arg); |
list = (LIST)BDY(arg); |
if ( OID(BDY(BDY(list))) == O_P ) |
if ( OID(BDY(BDY(list))) == O_P ) |
if ( ri > 1 ) |
if ( ri > 1 ) |
error("ifplot : invalid argument"); |
error("itvplot : invalid argument"); |
else |
else |
range[ri++] = list; |
range[ri++] = list; |
else |
else |
|
|
case O_STR: |
case O_STR: |
wname = (STRING)BDY(arg); break; |
wname = (STRING)BDY(arg); break; |
default: |
default: |
error("ifplot : invalid argument"); break; |
error("itvplot : invalid argument"); break; |
} |
} |
if ( !poly ) |
if ( !poly ) |
error("ifplot : invalid argument"); |
error("itvplot : invalid argument"); |
switch ( ri ) { |
switch ( ri ) { |
case 0: |
case 0: |
if ( !v[1] ) |
if ( !v[1] ) |
error("ifplot : please specify all variables"); |
error("itvplot : please specify all variables"); |
MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); |
MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); |
MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); |
MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n); |
break; |
break; |
case 1: |
case 1: |
if ( !v[1] ) |
if ( !v[1] ) |
error("ifplot : please specify all variables"); |
error("itvplot : please specify all variables"); |
av[0] = VR((P)BDY(BDY(range[0]))); |
av[0] = VR((P)BDY(BDY(range[0]))); |
if ( v[0] == av[0] ) { |
if ( v[0] == av[0] ) { |
xrange = range[0]; |
xrange = range[0]; |
|
|
MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); |
MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n); |
yrange = range[0]; |
yrange = range[0]; |
} else |
} else |
error("ifplot : invalid argument"); |
error("itvplot : invalid argument"); |
break; |
break; |
case 2: |
case 2: |
av[0] = VR((P)BDY(BDY(range[0]))); |
av[0] = VR((P)BDY(BDY(range[0]))); |
|
|
((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { |
((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) { |
xrange = range[0]; yrange = range[1]; |
xrange = range[0]; yrange = range[1]; |
} else |
} else |
error("ifplot : invalid argument"); |
error("itvplot : invalid argument"); |
break; |
break; |
default: |
default: |
error("ifplot : cannot happen"); break; |
error("itvplot : cannot happen"); break; |
} |
} |
/* ifplot in ox_plot requires |
stream = validate_ox_plot_stream(stream); |
[sid (Q), |
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), |
formula (Obj), |
xrange=[x,xmin,xmax] (LIST), |
xrange=[x,xmin,xmax] (LIST), |
yrange=[y,ymin,ymax] (LIST), |
yrange=[y,ymin,ymax] (LIST), |
|
|
wname=name (STRING)] |
wname=name (STRING)] |
*/ |
*/ |
|
|
if ( stream < 0 ) |
stream = validate_ox_plot_stream(stream); |
stream = current_s; |
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 |
else |
current_s = stream; |
color = (Q)BDY(arg); |
STOQ(stream,sid); |
/* 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 ) { |
if ( !geom ) { |
STOQ(300,w300); |
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,"plot"); |
MKSTR(fname, "ineqn"); |
arg = mknode(8,sid,fname,poly,xrange,yrange,0,geom,wname); |
arg = mknode(8, s_id, fname, poly, color, xrange, yrange, geom, wname); |
Pox_rpc(arg,rp); |
Pox_cmo_rpc(arg,rp); |
|
*rp = (Obj)s_id; |
} |
} |
|
|
void Pconplot(arg,rp) |
void Pineqnor(NODE arg, Obj *rp) |
NODE arg; |
|
Obj *rp; |
|
{ |
{ |
Q m2,p2,w300,sid; |
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; |
NODE defrange; |
LIST xrange,yrange,zrange,range[3],list,geom; |
LIST xrange,yrange,zrange,range[3],list,geom; |
VL vl,vl0; |
VL vl,vl0; |
V v[2],av[2]; |
V v[2],av[2]; |
int stream,ri,id,i; |
int stream,ri,i; |
P poly; |
P poly; |
P var; |
P var; |
NODE n,n0; |
NODE n,n0; |
STRING fname,wname; |
STRING fname,wname; |
|
Obj t; |
|
|
STOQ(-2,m2); STOQ(2,p2); |
STOQ(-2,m2); STOQ(2,p2); |
MKNODE(n,p2,0); MKNODE(defrange,m2,n); |
MKNODE(n,p2,0); MKNODE(defrange,m2,n); |
poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; |
poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; |
|
v[0] = v[1] = 0; |
for ( ; arg; arg = NEXT(arg) ) |
for ( ; arg; arg = NEXT(arg) ) |
if ( !BDY(arg) ) |
if ( !BDY(arg) ) |
stream = 0; |
stream = 0; |
|
|
switch ( OID(BDY(arg)) ) { |
switch ( OID(BDY(arg)) ) { |
case O_P: |
case O_P: |
poly = (P)BDY(arg); |
poly = (P)BDY(arg); |
get_vars_recursive(poly,&vl); |
get_vars_recursive((Obj)poly,&vl); |
for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) |
for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) |
if ( vl0->v->attr == V_IND ) |
if ( vl0->v->attr == (pointer)V_IND ) |
if ( i >= 2 ) |
if ( i >= 2 ) |
error("ifplot : invalid argument"); |
error("ifplot : invalid argument"); |
else |
else |
|
|
default: |
default: |
error("ifplot : cannot happen"); break; |
error("ifplot : cannot happen"); break; |
} |
} |
if ( stream < 0 ) |
|
stream = current_s; |
|
else |
|
current_s = stream; |
|
|
|
/* conplot in ox_plot requires |
/* conplot in ox_plot requires |
[sid (Q), |
[s_id (Q), |
formula (Obj), |
formula (Obj), |
xrange=[x,xmin,xmax] (LIST), |
xrange=[x,xmin,xmax] (LIST), |
yrange=[y,ymin,ymax] (LIST), |
yrange=[y,ymin,ymax] (LIST), |
|
|
wname=name (STRING)] |
wname=name (STRING)] |
*/ |
*/ |
|
|
if ( stream < 0 ) |
stream = validate_ox_plot_stream(stream); |
stream = current_s; |
STOQ(stream,s_id); |
else |
|
current_s = stream; |
|
STOQ(stream,sid); |
|
if ( !geom ) { |
if ( !geom ) { |
STOQ(300,w300); |
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,"plot"); |
if ( is_memory ) { |
arg = mknode(8,sid,fname,poly,xrange,yrange,zrange,geom,wname); |
MKSTR(fname,"memory_plot"); |
Pox_rpc(arg,rp); |
arg = mknode(7,s_id,fname,poly,xrange,yrange,zrange,geom); |
|
Pox_rpc(arg,&t); |
|
arg = mknode(1,s_id); |
|
Pox_pop_cmo(arg,rp); |
|
} else { |
|
MKSTR(fname,"plot"); |
|
arg = mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom,wname); |
|
Pox_rpc(arg,&t); |
|
*rp = (Obj)s_id; |
|
} |
} |
} |
|
|
void Pplot(arg,rp) |
void Pplot(NODE arg,Obj *rp) |
NODE arg; |
|
Obj *rp; |
|
{ |
{ |
Q m2,p2,w300,sid; |
plot_main(arg,0,rp); |
|
} |
|
|
|
void Pmemory_plot(NODE arg,Obj *rp) |
|
{ |
|
plot_main(arg,1,rp); |
|
} |
|
|
|
void plot_main(NODE arg,int is_memory,Obj *rp) |
|
{ |
|
Q m2,p2,w300,s_id; |
NODE defrange; |
NODE defrange; |
LIST xrange,range[1],list,geom; |
LIST xrange,range[1],list,geom; |
VL vl,vl0; |
VL vl,vl0; |
V v[1],av[1]; |
V v[1],av[1]; |
int stream,ri,id,i; |
int stream,ri,i; |
P poly; |
P poly; |
P var; |
P var; |
NODE n,n0; |
NODE n,n0; |
STRING fname,wname; |
STRING fname,wname; |
|
Obj t; |
|
|
STOQ(-2,m2); STOQ(2,p2); |
STOQ(-2,m2); STOQ(2,p2); |
MKNODE(n,p2,0); MKNODE(defrange,m2,n); |
MKNODE(n,p2,0); MKNODE(defrange,m2,n); |
poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; |
poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0; |
|
v[0] = 0; |
for ( ; arg; arg = NEXT(arg) ) |
for ( ; arg; arg = NEXT(arg) ) |
if ( !BDY(arg) ) |
if ( !BDY(arg) ) |
stream = 0; |
stream = 0; |
|
|
switch ( OID(BDY(arg)) ) { |
switch ( OID(BDY(arg)) ) { |
case O_P: case O_R: |
case O_P: case O_R: |
poly = (P)BDY(arg); |
poly = (P)BDY(arg); |
get_vars_recursive(poly,&vl); |
get_vars_recursive((Obj)poly,&vl); |
for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) |
for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) ) |
if ( vl0->v->attr == V_IND ) |
if ( vl0->v->attr == (pointer)V_IND ) |
if ( i >= 1 ) |
if ( i >= 1 ) |
error("ifplot : invalid argument"); |
error("ifplot : invalid argument"); |
else |
else |
v[i++] = vl0->v; |
v[i++] = vl0->v; |
|
if ( i != 1 ) |
|
error("ifplot : invalid argument"); |
break; |
break; |
case O_LIST: |
case O_LIST: |
list = (LIST)BDY(arg); |
list = (LIST)BDY(arg); |
|
|
error("plot : cannot happen"); break; |
error("plot : cannot happen"); break; |
} |
} |
/* conplot in ox_plot requires |
/* conplot in ox_plot requires |
[sid (Q), |
[s_id (Q), |
formula (Obj), |
formula (Obj), |
xrange=[x,xmin,xmax] (LIST), |
xrange=[x,xmin,xmax] (LIST), |
yrange=0, |
yrange=0, |
|
|
geom=[xsize,ysize] (LIST), |
geom=[xsize,ysize] (LIST), |
wname=name (STRING)] |
wname=name (STRING)] |
*/ |
*/ |
if ( stream < 0 ) |
stream = validate_ox_plot_stream(stream); |
stream = current_s; |
STOQ(stream,s_id); |
else |
|
current_s = stream; |
|
STOQ(stream,sid); |
|
if ( !geom ) { |
if ( !geom ) { |
STOQ(300,w300); |
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); |
|
Pox_rpc(arg,&t); |
|
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); |
|
Pox_rpc(arg,&t); |
|
*rp = (Obj)s_id; |
|
} |
|
} |
|
|
|
#define Pi 3.14159265358979323846264 |
|
|
|
void Ppolarplot(NODE arg,Q *rp) |
|
{ |
|
Q m2,p2,w300,s_id; |
|
NODE defrange; |
|
LIST zrange,range[1],list,geom; |
|
VL vl,vl0; |
|
V v[1],av[1]; |
|
int stream,ri,i; |
|
P poly; |
|
P var; |
|
NODE n,n0; |
|
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); |
|
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 ) { |
|
STOQ(300,w300); |
|
MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n); |
|
} |
MKSTR(fname,"plot"); |
MKSTR(fname,"plot"); |
arg = mknode(8,sid,fname,poly,xrange,0,0,geom,wname); |
arg = mknode(8,s_id,fname,poly,NULLP,NULLP,zrange,geom,wname); |
Pox_rpc(arg,rp); |
Pox_rpc(arg,&t); |
|
*rp = s_id; |
} |
} |
|
|
void Pplotover(arg,rp) |
void Pplotover(NODE arg,Q *rp) |
NODE arg; |
|
Obj *rp; |
|
{ |
{ |
Q sid,index; |
Q s_id,index; |
P poly; |
P poly; |
STRING fname; |
STRING fname; |
|
Obj t; |
|
|
poly = (P)ARG0(arg); |
poly = (P)ARG0(arg); |
sid = (Q)ARG1(arg); |
s_id = (Q)ARG1(arg); |
index = (Q)ARG2(arg); |
index = (Q)ARG2(arg); |
MKSTR(fname,"plotover"); |
MKSTR(fname,"plotover"); |
arg = mknode(4,sid,fname,index,poly); |
if ( argc(arg) == 4 ) |
Pox_rpc(arg,rp); |
arg = mknode(5,s_id,fname,index,poly,(Q)ARG3(arg)); |
|
else |
|
arg = mknode(4,s_id,fname,index,poly); |
|
Pox_rpc(arg,&t); |
|
*rp = s_id; |
} |
} |
|
|
/* arg = [x,y,r,sid,index] */ |
/* arg = [x,y,r,s_id,index] */ |
|
|
void Pdrawcircle(arg,rp) |
void Pdrawcircle(NODE arg,Q *rp) |
NODE arg; |
|
Obj *rp; |
|
{ |
{ |
Q sid,index; |
Q s_id,index; |
Obj x,y,r; |
Obj x,y,r,c; |
STRING fname; |
STRING fname; |
NODE n; |
NODE n; |
LIST pos; |
LIST pos; |
|
Obj t; |
|
|
x = (Obj)ARG0(arg); |
x = (Obj)ARG0(arg); |
y = (Obj)ARG1(arg); |
y = (Obj)ARG1(arg); |
r = (Obj)ARG2(arg); |
r = (Obj)ARG2(arg); |
sid = (Q)ARG3(arg); |
c = (Obj)ARG3(arg); |
index = (Q)ARG4(arg); |
s_id = (Q)ARG4(arg); |
|
index = (Q)ARG5(arg); |
MKSTR(fname,"drawcircle"); |
MKSTR(fname,"drawcircle"); |
n = mknode(3,x,y,r); MKLIST(pos,n); |
n = mknode(3,x,y,r,c); MKLIST(pos,n); |
arg = mknode(4,sid,fname,index,pos); |
arg = mknode(5,s_id,fname,index,pos,c); |
Pox_rpc(arg,rp); |
Pox_rpc(arg,&t); |
|
*rp = s_id; |
} |
} |
|
|
|
/* draw_obj(s_id,cindex,point|line); point = [x,y], line = [xa,ya,xb,yb] */ |
|
void Pdraw_obj(NODE arg,Q *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); |
|
Pox_cmo_rpc(arg,&t); |
|
*rp = s_id; |
|
} |
|
|
|
/* draw_string(s_id,cindex,pos,string,[,color]); pos=[x,y] */ |
|
void Pdraw_string(NODE arg,Q *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); |
|
Pox_cmo_rpc(arg,&t); |
|
*rp = s_id; |
|
} |
|
|
|
void Pclear_canvas(NODE arg,Q *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); |
|
Pox_cmo_rpc(arg,&t); |
|
*rp = s_id; |
|
} |
|
|
#if 0 |
#if 0 |
void Parrayplot(arg,rp) |
void Parrayplot(NODE arg,Obj *rp) |
NODE arg; |
|
Obj *rp; |
|
{ |
{ |
int s; |
int s; |
int id; |
int id; |