===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/plot/plotf.c,v
retrieving revision 1.11
retrieving revision 1.15
diff -u -p -r1.11 -r1.15
--- OpenXM_contrib2/asir2000/plot/plotf.c	2001/10/09 01:36:27	1.11
+++ OpenXM_contrib2/asir2000/plot/plotf.c	2002/08/02 08:59:48	1.15
@@ -45,7 +45,7 @@
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *
- * $OpenXM: OpenXM_contrib2/asir2000/plot/plotf.c,v 1.10 2001/08/22 09:19:21 noro Exp $ 
+ * $OpenXM: OpenXM_contrib2/asir2000/plot/plotf.c,v 1.14 2002/08/02 05:39:22 noro Exp $ 
 */
 #include "ca.h"
 #include "parse.h"
@@ -53,8 +53,9 @@
 #include "ifplot.h"
 
 void Pifplot(), Pconplot(), Pplotover(), Pplot(), Parrayplot(), Pdrawcircle();
+void Ppolarplot();
 void Pmemory_ifplot();
-void Popen_canvas(), Pclear_canvas(), Pdraw_obj();
+void Popen_canvas(), Pclear_canvas(), Pdraw_obj(), Pdraw_string();
 void Pox_rpc();
 void Pox_cmo_rpc();
 
@@ -63,19 +64,19 @@ struct ftab plot_tab[] = {
 	{"memory_ifplot",Pmemory_ifplot,-7},
 	{"conplot",Pconplot,-8},
 	{"plot",Pplot,-6},
+	{"polarplot",Ppolarplot,-6},
 	{"plotover",Pplotover,3},
 	{"drawcircle",Pdrawcircle,5},
 	{"open_canvas",Popen_canvas,-3},
 	{"clear_canvas",Pclear_canvas,2},
 	{"draw_obj",Pdraw_obj,-4},
+	{"draw_string",Pdraw_string,-5},
 /*
 	{"arrayplot",Parrayplot,2},
 */
 	{0,0,0},
 };
 
-int current_s;
-
 void Popen_canvas(NODE arg,Obj *rp)
 {
 	Q w300,s_id;
@@ -106,10 +107,7 @@ void Popen_canvas(NODE arg,Obj *rp)
 	   	wname=name (STRING)]
 	*/
 
-	if ( stream < 0 )
-		stream = current_s;
-	else
-		current_s = stream;
+	stream = validate_ox_plot_stream(stream);
 	STOQ(stream,s_id);
 	if ( !geom ) {
 		STOQ(300,w300);
@@ -120,7 +118,7 @@ void Popen_canvas(NODE arg,Obj *rp)
 	Pox_cmo_rpc(arg,rp);
 }
 
-void Pifplot(NODE arg,Obj *rp)
+void Pifplot(NODE arg,Q *rp)
 {
 	Q m2,p2,w300,s_id;
 	NODE defrange;
@@ -132,10 +130,12 @@ void Pifplot(NODE arg,Obj *rp)
 	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;
@@ -212,10 +212,7 @@ void Pifplot(NODE arg,Obj *rp)
 	   	wname=name (STRING)]
 	*/
 
-	if ( stream < 0 )
-		stream = current_s;
-	else
-		current_s = stream;
+	stream = validate_ox_plot_stream(stream);
 	STOQ(stream,s_id);
 	if ( !geom ) {
 		STOQ(300,w300);
@@ -223,10 +220,11 @@ void Pifplot(NODE arg,Obj *rp)
 	}
 	MKSTR(fname,"plot");
 	arg = mknode(8,s_id,fname,poly,xrange,yrange,0,geom,wname);
-	Pox_rpc(arg,rp);
+	Pox_rpc(arg,&t);
+	*rp = s_id;
 }
 
-void Pmemory_ifplot(NODE arg,Obj *rp)
+void Pmemory_ifplot(NODE arg,Q *rp)
 {
 	Q m2,p2,w300,s_id;
 	NODE defrange;
@@ -238,6 +236,7 @@ void Pmemory_ifplot(NODE arg,Obj *rp)
 	P var;
 	NODE n,n0;
 	STRING fname;
+	Obj t;
 
 	STOQ(-2,m2); STOQ(2,p2);
 	MKNODE(n,p2,0); MKNODE(defrange,m2,n); 
@@ -315,10 +314,7 @@ void Pmemory_ifplot(NODE arg,Obj *rp)
 	   	geom=[xsize,ysize] (LIST)]
 	*/
 
-	if ( stream < 0 )
-		stream = current_s;
-	else
-		current_s = stream;
+	stream = validate_ox_plot_stream(stream);
 	STOQ(stream,s_id);
 	if ( !geom ) {
 		STOQ(300,w300);
@@ -326,10 +322,11 @@ void Pmemory_ifplot(NODE arg,Obj *rp)
 	}
 	MKSTR(fname,"memory_plot");
 	arg = mknode(7,s_id,fname,poly,xrange,yrange,0,geom);
-	Pox_rpc(arg,rp);
+	Pox_rpc(arg,&t);
+	*rp = s_id;
 }
 
-void Pconplot(NODE arg,Obj *rp)
+void Pconplot(NODE arg,Q *rp)
 {
 	Q m2,p2,w300,s_id;
 	NODE defrange;
@@ -341,10 +338,12 @@ void Pconplot(NODE arg,Obj *rp)
 	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;
@@ -422,11 +421,6 @@ void Pconplot(NODE arg,Obj *rp)
 		default:
 			error("ifplot : cannot happen"); break;
 	}
-	if ( stream < 0 )
-		stream = current_s;
-	else
-		current_s = stream;
-
 	/* conplot in ox_plot requires 
 	   [s_id (Q),
 	   	formula (Obj),
@@ -437,10 +431,7 @@ void Pconplot(NODE arg,Obj *rp)
 	   	wname=name (STRING)]
 	*/
 
-	if ( stream < 0 )
-		stream = current_s;
-	else
-		current_s = stream;
+	stream = validate_ox_plot_stream(stream);
 	STOQ(stream,s_id);
 	if ( !geom ) {
 		STOQ(300,w300);
@@ -448,10 +439,11 @@ void Pconplot(NODE arg,Obj *rp)
 	}
 	MKSTR(fname,"plot");
 	arg = mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom,wname);
-	Pox_rpc(arg,rp);
+	Pox_rpc(arg,&t);
+	*rp = s_id;
 }
 
-void Pplot(NODE arg,Obj *rp)
+void Pplot(NODE arg,Q *rp)
 {
 	Q m2,p2,w300,s_id;
 	NODE defrange;
@@ -463,10 +455,12 @@ void Pplot(NODE arg,Obj *rp)
 	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] = 0;
 	for ( ; arg; arg = NEXT(arg) )
 		if ( !BDY(arg) )
 			stream = 0;
@@ -526,10 +520,7 @@ void Pplot(NODE arg,Obj *rp)
 	   	geom=[xsize,ysize] (LIST),
 	   	wname=name (STRING)]
 	*/
-	if ( stream < 0 )
-		stream = current_s;
-	else
-		current_s = stream;
+	stream = validate_ox_plot_stream(stream);
 	STOQ(stream,s_id);
 	if ( !geom ) {
 		STOQ(300,w300);
@@ -537,32 +528,119 @@ void Pplot(NODE arg,Obj *rp)
 	}
 	MKSTR(fname,"plot");
 	arg = mknode(8,s_id,fname,poly,xrange,0,0,geom,wname);
-	Pox_rpc(arg,rp);
+	Pox_rpc(arg,&t);
+	*rp = s_id;
 }
 
-void Pplotover(NODE arg,Obj *rp)
+#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");
+	arg = mknode(8,s_id,fname,poly,0,0,zrange,geom,wname);
+	Pox_rpc(arg,&t);
+	*rp = s_id;
+}
+
+void Pplotover(NODE arg,Q *rp)
+{
 	Q s_id,index;
 	P poly;
 	STRING fname;
+	Obj t;
 
 	poly = (P)ARG0(arg);
 	s_id = (Q)ARG1(arg);
 	index = (Q)ARG2(arg);
 	MKSTR(fname,"plotover");
 	arg = mknode(4,s_id,fname,index,poly);
-	Pox_rpc(arg,rp);
+	Pox_rpc(arg,&t);
+	*rp = s_id;
 }
 
 /* arg = [x,y,r,s_id,index] */
 
-void Pdrawcircle(NODE arg,Obj *rp)
+void Pdrawcircle(NODE arg,Q *rp)
 {
 	Q s_id,index;
 	Obj x,y,r;
 	STRING fname;
 	NODE n;
 	LIST pos;
+	Obj t;
 
 	x = (Obj)ARG0(arg);
 	y = (Obj)ARG1(arg);
@@ -572,15 +650,17 @@ void Pdrawcircle(NODE arg,Obj *rp)
 	MKSTR(fname,"drawcircle");
 	n = mknode(3,x,y,r); MKLIST(pos,n);
 	arg = mknode(4,s_id,fname,index,pos);
-	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,Obj *rp)
+void Pdraw_obj(NODE arg,Q *rp)
 {
 	static STRING fname;
 	Q s_id,index;
 	LIST obj;
+	Obj t;
 
 	if ( !fname ) {
 		MKSTR(fname,"draw_obj");
@@ -593,21 +673,49 @@ void Pdraw_obj(NODE arg,Obj *rp)
 		arg = mknode(5,s_id,fname,index,obj,ARG3(arg));
 	else
 		arg = mknode(4,s_id,fname,index,obj);
-	Pox_cmo_rpc(arg,rp);
+	Pox_cmo_rpc(arg,&t);
+	*rp = s_id;
 }
 
-void Pclear_canvas(NODE arg,Obj *rp)
+/* 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,rp);
+	Pox_cmo_rpc(arg,&t);
+	*rp = s_id;
 }
 
 #if 0