===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/builtin/subst.c,v
retrieving revision 1.5
retrieving revision 1.6
diff -u -p -r1.5 -r1.6
--- OpenXM_contrib2/asir2000/builtin/subst.c	2001/10/09 01:36:07	1.5
+++ OpenXM_contrib2/asir2000/builtin/subst.c	2003/11/27 08:28:40	1.6
@@ -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/builtin/subst.c,v 1.4 2001/09/05 04:43:58 noro Exp $ 
+ * $OpenXM: OpenXM_contrib2/asir2000/builtin/subst.c,v 1.5 2001/10/09 01:36:07 noro Exp $ 
 */
 #include "ca.h"
 #include "parse.h"
@@ -67,27 +67,102 @@ Obj *rp;
 	Obj a,b,t;
 	LIST l;
 	V v;
+	int row,col,len;
+	VECT vect;
+	MAT mat;
+	int i,j;
+	NODE n0,n,nd;
+	struct oNODE arg0;
+	MP m,mp,mp0;
+	DP d;
 
 	if ( !arg ) {
 		*rp = 0; return;
 	}
-	asir_assert(ARG0(arg),O_R,"subst");
-	reductr(CO,(Obj)ARG0(arg),&a);
-/*	a = (Obj)ARG0(arg); */
-	arg = NEXT(arg);
-	if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )
-		arg = BDY(l);
-	while ( arg ) {
-		asir_assert(BDY(arg),O_P,"subst");
-		v = VR((P)BDY(arg)); arg = NEXT(arg);
-		if ( !arg )
-			error("subst : invalid argument");
-		asir_assert(ARG0(arg),O_R,"subst");
-		reductr(CO,(Obj)BDY(arg),&b); arg = NEXT(arg);
-/*		b = (Obj)BDY(arg); arg = NEXT(arg); */
-		substr(CO,0,a,v,b,&t); a = t;
+	a = (Obj)ARG0(arg);
+	if ( !a ) {
+		*rp = 0;
+		return;
 	}
-	*rp = a;
+	switch ( OID(a) ) {
+		case O_N: case O_P: case O_R:
+			reductr(CO,(Obj)ARG0(arg),&a);
+			arg = NEXT(arg);
+			if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )
+				arg = BDY(l);
+			while ( arg ) {
+				asir_assert(BDY(arg),O_P,"subst");
+				v = VR((P)BDY(arg)); arg = NEXT(arg);
+				if ( !arg )
+					error("subst : invalid argument");
+				asir_assert(ARG0(arg),O_R,"subst");
+				reductr(CO,(Obj)BDY(arg),&b); arg = NEXT(arg);
+		/*		b = (Obj)BDY(arg); arg = NEXT(arg); */
+				substr(CO,0,a,v,b,&t); a = t;
+			}
+			*rp = a;
+			break;
+		case O_LIST:
+			n0 = 0;
+			for ( nd = BDY((LIST)a); nd; nd = NEXT(nd) ) {
+				NEXTNODE(n0,n);
+				arg0.body = (pointer)BDY(nd);
+				arg0.next = NEXT(arg);
+				Psubst(&arg0,&b);
+				BDY(n) = (pointer)b;
+			}
+			if ( n0 )
+				NEXT(n) = 0;
+			MKLIST(l,n0);
+			*rp = (Obj)l;
+			break;
+		case O_VECT:
+			len = ((VECT)a)->len;
+			MKVECT(vect,len);
+			for ( i = 0; i < len; i++ ) {
+				arg0.body = (pointer)BDY((VECT)a)[i];
+				arg0.next = NEXT(arg);
+				Psubst(&arg0,&b);
+				BDY(vect)[i] = (pointer)b;
+			}
+			*rp = (Obj)vect;
+			break;
+		case O_MAT:
+			row = ((MAT)a)->row;
+			col = ((MAT)a)->col;
+			MKMAT(mat,row,col);
+			for ( i = 0; i < row; i++ )
+				for ( j = 0; j < col; j++ ) {
+					arg0.body = (pointer)BDY((MAT)a)[i][j];
+					arg0.next = NEXT(arg);
+					Psubst(&arg0,&b);
+					BDY(mat)[i][j] = (pointer)b;
+				}
+			*rp = (Obj)mat;
+			break;
+		case O_DP:
+			mp0 = 0;
+			for ( m = BDY((DP)a); m; m = NEXT(m) ) {
+				arg0.body = (pointer)C(m);
+				arg0.next = NEXT(arg);
+				Psubst(&arg0,&b);
+				if ( b ) {
+					NEXTMP(mp0,mp);
+					C(mp) = (P)b;
+					mp->dl = m->dl;
+				}
+			}
+			if ( mp0 ) {
+				MKDP(NV((DP)a),mp0,d);
+				d->sugar = ((DP)a)->sugar;
+				*rp = (Obj)d;
+			} else
+				*rp = 0;
+
+			break;
+		default:
+			error("subst invalid argument");
+	}
 }
 
 FNODE subst_in_fnode();