===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/nd.c,v
retrieving revision 1.177
retrieving revision 1.180
diff -u -p -r1.177 -r1.180
--- OpenXM_contrib2/asir2000/engine/nd.c	2009/10/12 10:43:45	1.177
+++ OpenXM_contrib2/asir2000/engine/nd.c	2010/02/05 04:24:00	1.180
@@ -1,4 +1,4 @@
-/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.176 2009/09/24 07:13:00 noro Exp $ */
+/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.179 2010/02/03 07:23:06 noro Exp $ */
 
 #include "nd.h"
 
@@ -2845,7 +2845,7 @@ void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
     }
     nd_init_ord(ord);
     mrank = 0;
-    for ( t = BDY(f), max = 0; t; t = NEXT(t) )
+    for ( t = BDY(f), max = 1; t; t = NEXT(t) )
         for ( tv = vv; tv; tv = NEXT(tv) ) {
             if ( nd_module ) {
                 s = BDY((LIST)BDY(t));
@@ -2975,7 +2975,7 @@ void nd_gr_postproc(LIST f,LIST v,int m,struct order_s
         nvar += nalg;
     }
     nd_init_ord(ord);
-    for ( t = BDY(f), max = 0; t; t = NEXT(t) )
+    for ( t = BDY(f), max = 1; t; t = NEXT(t) )
         for ( tv = vv; tv; tv = NEXT(tv) ) {
             e = getdeg(tv->v,(P)BDY(t));
             max = MAX(e,max);
@@ -3087,7 +3087,7 @@ void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
     m = trace > 1 ? trace : get_lprime(mindex);
     nd_init_ord(ord);
     mrank = 0;
-    for ( t = BDY(f), max = 0; t; t = NEXT(t) )
+    for ( t = BDY(f), max = 1; t; t = NEXT(t) )
         for ( tv = vv; tv; tv = NEXT(tv) ) {
             if ( nd_module ) {
                 s = BDY((LIST)BDY(t));
@@ -3448,7 +3448,7 @@ void ndv_homogenize(NDV p,int obpe,int oadv,EPOS oepos
     NMV m,mr0,mr,t;
 
     len = p->len;
-    for ( m = BDY(p), i = 0, max = 0; i < len; NMV_OADV(m), i++ )
+    for ( m = BDY(p), i = 0, max = 1; i < len; NMV_OADV(m), i++ )
         max = MAX(max,TD(DL(m)));
     mr0 = nmv_adv>oadv?(NMV)REALLOC(BDY(p),len*nmv_adv):BDY(p);
     m = (NMV)((char *)mr0+(len-1)*oadv);
@@ -4052,7 +4052,9 @@ void weyl_mul_nm_nmv(int n,int mod,NM m0,NMV m1,NM *ta
     NEWNM(m); d = DL(m);
     if ( mod ) {
         c0 = CM(m0); c1 = CM(m1); DMAR(c0,c1,0,mod,c); CM(m) = c;
-    } else
+    } else if ( nd_vc )
+        mulp(nd_vc,CP(m0),CP(m1),&CP(m));
+	else
         mulq(CQ(m0),CQ(m1),&CQ(m));
     for ( i = 0; i < nd_wpd; i++ ) d[i] = 0;
     homo = n&1 ? 1 : 0;
@@ -4107,7 +4109,9 @@ void weyl_mul_nm_nmv(int n,int mod,NM m0,NMV m1,NM *ta
                         ndl_addto(DL(tab[u]),d);
                         if ( mod ) {
                             c0 = CM(tab[u]); DMAR(c0,c,0,mod,c1); CM(tab[u]) = c1;
-                        } else {
+                        } else if ( nd_vc )
+                            mulp(nd_vc,CP(tab[u]),(P)q,&CP(tab[u]));
+						else {
                             mulq(CQ(tab[u]),q,&q1); CQ(tab[u]) = q1;
                         }
                     }
@@ -4119,7 +4123,9 @@ void weyl_mul_nm_nmv(int n,int mod,NM m0,NMV m1,NM *ta
                         ndl_add(DL(tab[u]),d,DL(t));
                         if ( mod ) {
                             c0 = CM(tab[u]); DMAR(c0,c,0,mod,c1); CM(t) = c1;
-                        } else
+                        } if ( nd_vc )
+                            mulp(nd_vc,CP(tab[u]),(P)q,&CP(t));
+						else
                             mulq(CQ(tab[u]),q,&CQ(t));
                         *p = t;
                     }
@@ -4908,7 +4914,7 @@ void nd_nf_p(P f,LIST g,LIST v,int m,struct order_spec
     for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
 
     /* get the degree bound */
-    for ( t = BDY(g), max = 0; t; t = NEXT(t) )
+    for ( t = BDY(g), max = 1; t; t = NEXT(t) )
         for ( tv = vv; tv; tv = NEXT(tv) ) {
             e = getdeg(tv->v,(P)BDY(t));
             max = MAX(e,max);
@@ -6503,7 +6509,7 @@ void nd_det(int mod,MAT f,P *rp)
         m = (pointer **)w;
     }
 
-    for ( i = 0, max = 0; i < n; i++ )
+    for ( i = 0, max = 1; i < n; i++ )
         for ( j = 0; j < n; j++ )
             for ( tv = fv; tv; tv = NEXT(tv) ) {
                 e = getdeg(tv->v,(P)m[i][j]);
@@ -6511,7 +6517,7 @@ void nd_det(int mod,MAT f,P *rp)
             }
     nd_setup_parameters(nvar,max);
     dm = (NDV **)almat_pointer(n,n);
-    for ( i = 0, max = 0; i < n; i++ )
+    for ( i = 0, max = 1; i < n; i++ )
         for ( j = 0; j < n; j++ ) {
             dm[i][j] = ptondv(CO,fv,m[i][j]);
             if ( mod ) ndv_mod(mod,dm[i][j]);