===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/nd.c,v
retrieving revision 1.87
retrieving revision 1.90
diff -u -p -r1.87 -r1.90
--- OpenXM_contrib2/asir2000/engine/nd.c	2003/11/05 08:02:45	1.87
+++ OpenXM_contrib2/asir2000/engine/nd.c	2004/03/13 06:49:15	1.90
@@ -1,4 +1,4 @@
-/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.86 2003/10/28 08:47:01 noro Exp $ */
+/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.89 2004/02/03 23:31:57 noro Exp $ */
 
 #include "ca.h"
 #include "parse.h"
@@ -146,12 +146,11 @@ static int nd_psn,nd_pslen;
 static RHist *nd_red;
 
 static int nd_found,nd_create,nd_notfirst;
-static int nm_adv;
 static int nmv_adv;
 static int nd_dcomp;
 static int nd_demand;
 
-extern struct order_spec dp_current_spec;
+extern struct order_spec *dp_current_spec;
 extern char *Demand;
 extern VL CO;
 extern int Top,Reverse,DP_Print,dp_nelim,do_weyl,NoSugar;
@@ -179,6 +178,8 @@ extern int *current_weyl_weight_vector;
 #define GET_EXP_OLD(d,a) (((d)[oepos[a].i]>>oepos[a].s)&omask0)
 #define PUT_EXP_OLD(r,a,e) ((r)[oepos[a].i] |= ((e)<<oepos[a].s))
 
+#define ROUND_FOR_ALIGN(s) ((((s)+sizeof(void *)-1)/sizeof(void *))*sizeof(void *))
+
 /* macros for term comparison */
 #define TD_DL_COMPARE(d1,d2)\
 (TD(d1)>TD(d2)?1:(TD(d1)<TD(d2)?-1:ndl_lex_compare(d1,d2)))
@@ -2406,6 +2407,14 @@ void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
 	ndv_alloc = 0;
 	get_vars((Obj)f,&fv); pltovl(v,&vv);
 	for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
+	switch ( ord->id ) {
+		case 1:
+			if ( ord->nv != nvar )
+				error("nd_{gr,f4} : invalid order specification");
+			break;
+		default:
+			break;
+	}
 	nd_init_ord(ord);
 	for ( t = BDY(f), max = 0; t; t = NEXT(t) )
 		for ( tv = vv; tv; tv = NEXT(tv) ) {
@@ -2438,7 +2447,7 @@ void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
 
 void nd_gr_trace(LIST f,LIST v,int trace,int homo,struct order_spec *ord,LIST *rp)
 {
-	struct order_spec ord1;
+	struct order_spec *ord1;
 	VL tv,fv,vv,vc;
 	NODE fd,fd0,in0,in,r,r0,t,s,cand;
 	int m,nocheck,nvar,mindex,e,max;
@@ -2450,6 +2459,14 @@ void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru
 
 	get_vars((Obj)f,&fv); pltovl(v,&vv);
 	for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
+	switch ( ord->id ) {
+		case 1:
+			if ( ord->nv != nvar )
+				error("nd_gr_trace : invalid order specification");
+			break;
+		default:
+			break;
+	}
 	nocheck = 0;
 	mindex = 0;
 
@@ -2489,7 +2506,7 @@ void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru
 				wmax = MAX(TD(DL(a)),wmax);
 		}
 		homogenize_order(ord,nvar,&ord1);
-		nd_init_ord(&ord1);
+		nd_init_ord(ord1);
 		nd_setup_parameters(nvar+1,wmax);
 		for ( t = fd0; t; t = NEXT(t) )
 			ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos);
@@ -2528,7 +2545,7 @@ void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru
 			m = get_lprime(++mindex);
 			/* reset the parameters */
 			if ( !ishomo && homo ) {
-				nd_init_ord(&ord1);
+				nd_init_ord(ord1);
 				nd_setup_parameters(nvar+1,wmax);
 			} else {
 				nd_init_ord(ord);
@@ -2759,7 +2776,7 @@ void ndv_dehomogenize(NDV p,struct order_spec *ord)
 	for ( m = BDY(p), i = 0; i < len; NMV_ADV(m), i++ )
 		ndl_dehomogenize(DL(m));
 	if ( newwpd != nd_wpd ) {
-		newadv = sizeof(struct oNMV)+(newwpd-1)*sizeof(UINT);
+		newadv = ROUND_FOR_ALIGN(sizeof(struct oNMV)+(newwpd-1)*sizeof(UINT));
 		for ( m = r = BDY(p), i = 0; i < len; NMV_ADV(m), NDV_NADV(r), i++ ) {
 			CQ(r) = CQ(m);
 			for ( j = 0; j < newexporigin; j++ ) DL(r)[j] = DL(m)[j];
@@ -2955,8 +2972,7 @@ void nd_setup_parameters(int nvar,int max) {
 		nd_mask[nd_epw-i-1] = (nd_mask0<<(i*nd_bpe));
 		nd_mask1 |= (1<<(nd_bpe-1))<<(i*nd_bpe);
 	}
-	nm_adv = sizeof(struct oNM)+(nd_wpd-1)*sizeof(UINT);
-	nmv_adv = sizeof(struct oNMV)+(nd_wpd-1)*sizeof(UINT);
+	nmv_adv = ROUND_FOR_ALIGN(sizeof(struct oNMV)+(nd_wpd-1)*sizeof(UINT));
 	nd_epos = nd_create_epos(nd_ord);
 	nd_blockmask = nd_create_blockmask(nd_ord);
 }
@@ -4417,7 +4433,7 @@ NODE nd_f4_red_dist(int m,ND_pairs sp0,UINT *s0vect,in
 	nd_send_int(nd_wpd);
 	nd_send_int(nmv_adv);
 
-	saveobj(nd_write,dp_current_spec.obj); fflush(nd_write);
+	saveobj(nd_write,dp_current_spec->obj); fflush(nd_write);
 
 	nd_send_int(nd_psn);
 	for ( i = 0; i < nd_psn; i++ ) nd_send_ndv(nd_ps[i]);
@@ -4463,7 +4479,7 @@ void nd_exec_f4_red_dist()
 	ND_pairs *sp0;
 	int *colstat;
 	int a,sprow,rank;
-	struct order_spec ord;
+	struct order_spec *ord;
 	Obj ordspec;
 	ND spol;
 	int maxrs;
@@ -4478,8 +4494,8 @@ void nd_exec_f4_red_dist()
 	nmv_adv = nd_recv_int();
 
 	loadobj(nd_read,&ordspec);
-	create_order_spec(ordspec,&ord);
-	nd_init_ord(&ord);
+	create_order_spec(0,ordspec,&ord);
+	nd_init_ord(ord);
 	nd_setup_parameters(nd_nvar,0);
 
 	nd_psn = nd_recv_int();