===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/lib/weight,v
retrieving revision 1.14
retrieving revision 1.17
diff -u -p -r1.14 -r1.17
--- OpenXM_contrib2/asir2000/lib/weight	2003/12/11 15:12:35	1.14
+++ 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++)
@@ -448,7 +449,9 @@ 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])$
 			if(TMP!=[])
@@ -559,12 +562,15 @@ 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])$
 				if(TMP!=[])