===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/nd.c,v
retrieving revision 1.147
retrieving revision 1.149
diff -u -p -r1.147 -r1.149
--- OpenXM_contrib2/asir2000/engine/nd.c	2006/11/27 07:54:41	1.147
+++ OpenXM_contrib2/asir2000/engine/nd.c	2006/12/04 01:40:51	1.149
@@ -1,4 +1,4 @@
-/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.146 2006/11/27 07:31:25 noro Exp $ */
+/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.148 2006/11/28 02:17:54 noro Exp $ */
 
 #include "nd.h"
 
@@ -52,6 +52,8 @@ void nd_det_reconstruct(NDV **dm,int n,int j,NDV d);
 int nd_monic(int m,ND *p);
 NDV plain_vect_to_ndv_q(Q *mat,int col,UINT *s0vect);
 
+extern int Denominator,DP_Multiple;
+
 void nd_free_private_storage()
 {
 	_nm_free_list = 0;
@@ -1685,6 +1687,7 @@ again:
 	else
 		for ( t = g; t; t = NEXT(t) )
 			BDY(t) = (pointer)nd_ps[(int)BDY(t)];
+	if ( !checkonly && DP_Print ) { printf("nd_gb done.\n"); fflush(stdout); }
 	return g;
 }
 
@@ -1864,6 +1867,7 @@ again:
 		print_eg("invdalg",&eg_invdalg);
 		print_eg("le",&eg_le);
 	}
+	if ( DP_Print ) { printf("nd_gb_trace done.\n"); fflush(stdout); }
 	return g;
 }
 
@@ -2411,6 +2415,8 @@ void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
 	if ( !m && Demand ) nd_demand = 1;
 	else nd_demand = 0;
 
+	if ( DP_Multiple )
+		nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1);
 #if 0
 	ndv_alloc = 0;
 #endif
@@ -2575,6 +2581,9 @@ void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
 	struct order_spec *ord1;
 	struct oEGT eg_check,eg0,eg1;
 
+	if ( DP_Multiple )
+		nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1);
+
 	get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
 	for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
 	switch ( ord->id ) {
@@ -3314,10 +3323,12 @@ ND nd_copy(ND p)
 
 int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
 {
-	NM m;
+	NM m1,m2;
 	NDV p1,p2;
 	ND t1,t2;
 	UINT *lcm;
+	P gp,tp;
+	Q g,t;
 	int td;
 
 	if ( !mod && nd_demand ) {
@@ -3330,23 +3341,29 @@ int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
 		}
 	}
 	lcm = LCM(p);
-	NEWNM(m);
-	CQ(m) = HCQ(p2);
-	ndl_sub(lcm,HDL(p1),DL(m));
-	if ( ndl_check_bound2(p->i1,DL(m)) )
-		return 0;
-	t1 = ndv_mul_nm(mod,m,p1);
-	if ( mod == -1 ) CM(m) = _chsgnsf(HCM(p1));
-	else if ( mod ) CM(m) = mod-HCM(p1); 
-	else chsgnp(HCP(p1),&CP(m));
-	ndl_sub(lcm,HDL(p2),DL(m));
-	if ( ndl_check_bound2(p->i2,DL(m)) ) {
-		nd_free(t1);
-		return 0;
+	NEWNM(m1); ndl_sub(lcm,HDL(p1),DL(m1));
+	if ( ndl_check_bound2(p->i1,DL(m1)) ) {
+		FREENM(m1); return 0;
 	}
-	t2 = ndv_mul_nm(mod,m,p2);
+	NEWNM(m2); ndl_sub(lcm,HDL(p2),DL(m2));
+	if ( ndl_check_bound2(p->i2,DL(m2)) ) {
+		FREENM(m1); FREENM(m2); return 0;
+	}
+
+	if ( mod == -1 ) {
+		CM(m1) = HCM(p2); CM(m2) = _chsgnsf(HCM(p1));
+	} else if ( mod ) {
+		CM(m1) = HCM(p2); CM(m2) = mod-HCM(p1); 
+	} else if ( nd_vc ) {
+		ezgcdpz(nd_vc,HCP(p1),HCP(p2),&gp);
+		divsp(nd_vc,HCP(p2),gp,&CP(m1));
+		divsp(nd_vc,HCP(p1),gp,&tp); chsgnp(tp,&CP(m2));
+	} else {
+		igcd_cofactor(HCQ(p1),HCQ(p2),&g,&t,&CQ(m1)); chsgnq(t,&CQ(m2));
+	}
+	t1 = ndv_mul_nm(mod,m1,p1); t2 = ndv_mul_nm(mod,m2,p2);
 	*rp = nd_add(mod,t1,t2);
-	FREENM(m);
+	FREENM(m1); FREENM(m2);
 	return 1;
 }