===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/lib/weight,v
retrieving revision 1.12
retrieving revision 1.17
diff -u -p -r1.12 -r1.17
--- OpenXM_contrib2/asir2000/lib/weight	2003/12/08 07:18:14	1.12
+++ OpenXM_contrib2/asir2000/lib/weight	2004/01/06 17:23:23	1.17
@@ -163,7 +163,7 @@ def getgcd(A,B){
 	for(I=0;I<VarsNumA;I++){
 
 		for(J=0;J<VarsNumB;J++)
-			if(C[J]==A[I][0])
+			if(B[J]==A[I][0])
 				break$
 
 		if(J<VarsNumB)
@@ -204,24 +204,25 @@ def makeret(Res,Vars,FLAG){
 	VarsNum=length(Vars)$
 
 	ResVec=newvect(ResNum)$
-	for(M=0,I=0;I<ResNum;I++){
-		if(member(Res[I][0],Vars)){
-			ResVec[I]=Res[I][1]$
 
-			if(FLAG && type(ResVec[I])==1){
-				if(M==0)
-					M=ResVec[I]$
-				else
-					if(ResVec[I]<M)
-						M=ResVec[I]$
-			}
-		}
-	}		
-	
+        for(M=0,I=0;I<ResNum;I++){
+                if(member(Res[I][0],Vars)){
+                        ResVec[I]=Res[I][1]$
+
+                        if(FLAG && type(ResVec[I])==1){
+                                if(M==0)
+                                        M=ResVec[I]$
+                                else
+                                        if(ResVec[I]<M)
+                                                M=ResVec[I]$
+                        }
+                }
+        }               
+
 	if(M!=0)
 		ResVec=ResVec/M;
 
-	RET=newvect(VarsNum,Vars)$
+ 	RET=newvect(VarsNum,Vars)$
 
 	for(I=0;I<ResNum;I++){
 		for(J=0;J<VarsNum;J++)
@@ -249,7 +250,7 @@ def roundret(V){
 	VN=size(V)[0]$
 
 	RET0=V$
-	for(I=2;I<1000;I++){
+	for(I=1;I<1000;I++){
 		RET1=I*RET0$
 		for(J=0;J<VN;J++){
 			X=drint(RET1[J])$
@@ -448,9 +449,11 @@ def leastsq(NormMat,ExpMat,Vars,FLAG){
 	Res=getgcd(Res,Rea)$
 
 	if(nonposdegchk(Res)){
+
 		TMP1=makeret(Res,Vars,1)$
+
 		if(TMP1[0]==0){	
-			TMP=roundret(TMP1[1]*1.0)$
+			TMP=roundret(TMP1[1])$
 			if(TMP!=[])
 				RET=append(RET,wsort(TMP1[1],Vars,TMP,FLAG))$
 
@@ -559,14 +562,17 @@ def weightr(ExpMat,Vars,PolyListNum,OneMat,FLAG){
 			if(member(Rea[I],Vars))
 				TMP=cons(Rea[I],TMP)$
 	
-		TMP=cons(ExtVars[F],TMP)$
+		if(member(ExtVars[F],Vars))
+			TMP=cons(ExtVars[F],TMP)$
+
 		Res=getgcd(Res,TMP)$
 	
 		if(nonposdegchk(Res)){
 
 			TMP1=makeret(Res,Vars,1)$
+
 			if(TMP1[0]==0){	
-				TMP=roundret(TMP1[1]*1.0)$
+				TMP=roundret(TMP1[1])$
 				if(TMP!=[])
 					RET=append(RET,wsort(TMP1[1],Vars,
 						TMP,FLAG))$