===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/nd.c,v
retrieving revision 1.200
retrieving revision 1.203
diff -u -p -r1.200 -r1.203
--- OpenXM_contrib2/asir2000/engine/nd.c	2012/12/17 07:20:44	1.200
+++ OpenXM_contrib2/asir2000/engine/nd.c	2013/01/31 01:13:47	1.203
@@ -1,10 +1,11 @@
-/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.199 2012/08/27 05:38:00 noro Exp $ */
+/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.202 2013/01/30 08:03:18 noro Exp $ */
 
 #include "nd.h"
 
 struct oEGT eg_search;
 
 int diag_period = 6;
+int weight_check = 1;
 int (*ndl_compare_function)(UINT *a1,UINT *a2);
 int nd_dcomp;
 NM _nm_free_list;
@@ -3952,6 +3953,31 @@ void nd_setup_parameters(int nvar,int max) {
         else if ( max < 65536 ) nd_bpe = 16;
         else nd_bpe = 32;
     }
+    if ( !do_weyl && weight_check && (current_dl_weight_vector || nd_matrix) ) {
+        UINT t;
+		int st;
+        int *v;
+	/* t = max(weights) */
+        t = 0;
+        if ( current_dl_weight_vector )
+            for ( i = 0, t = 0; i < nd_nvar; i++ ) {
+                if ( (st=current_dl_weight_vector[i]) < 0 ) st = -st;
+                if ( t < st ) t = st;
+            }
+        if ( nd_matrix )
+            for ( i = 0; i < nd_matrix_len; i++ )
+                for ( j = 0, v = nd_matrix[i]; j < nd_nvar; j++ ) {
+                    if ( (st=v[j]) < 0 ) st = -st;
+                    if ( t < st ) t = st;
+				}
+        /* i = bitsize of t */
+        for ( i = 0; t; t >>=1, i++ );
+        /* i += bitsize of nd_nvar */
+        for ( t = nd_nvar; t; t >>=1, i++);
+        /* nd_bpe+i = bitsize of max(weights)*max(exp)*nd_nvar */
+        if ( (nd_bpe+i) >= 31 )
+            error("nd_setup_parameters : too large weight");
+    }
     nd_epw = (sizeof(UINT)*8)/nd_bpe;
     elen = nd_nvar/nd_epw+(nd_nvar%nd_epw?1:0);
     nd_exporigin = nd_get_exporigin(nd_ord);
@@ -4931,6 +4957,8 @@ NODE ndv_reducebase(NODE x,int *perm)
 void nd_init_ord(struct order_spec *ord)
 {
     nd_module = (ord->id >= 256);
+	nd_matrix = 0;
+	nd_matrix_len = 0;
     switch ( ord->id ) {
         case 0:
             switch ( ord->ord.simple ) {