===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/builtin/array.c,v
retrieving revision 1.53
retrieving revision 1.58
diff -u -p -r1.53 -r1.58
--- OpenXM_contrib2/asir2000/builtin/array.c	2006/06/12 11:52:10	1.53
+++ OpenXM_contrib2/asir2000/builtin/array.c	2009/03/03 10:04:10	1.58
@@ -45,7 +45,7 @@
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *
- * $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.52 2006/05/30 07:35:30 noro Exp $
+ * $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.57 2009/02/03 00:39:23 noro Exp $
 */
 #include "ca.h"
 #include "base.h"
@@ -54,7 +54,9 @@
 
 #include <sys/types.h>
 #include <sys/stat.h>
+#if !defined(_MSC_VER)
 #include <unistd.h>
+#endif
 
 #define F4_INTRAT_PERIOD 8 
 
@@ -391,11 +393,13 @@ void Pnewvect(NODE arg,VECT *rp)
 	if ( argc(arg) == 2 ) {
 		list = (LIST)ARG1(arg);
 		asir_assert(list,O_LIST,"newvect");
+#if 0
 		for ( r = 0, tn = BDY(list); tn; r++, tn = NEXT(tn) );
 		if ( r > len ) {
 			*rp = vect;
 			return;
 		}
+#endif
 		for ( i = 0, tn = BDY(list), vb = BDY(vect); tn; i++, tn = NEXT(tn) )
 			vb[i] = (pointer)BDY(tn);
 	}
@@ -457,15 +461,29 @@ void Pnewbytearray(NODE arg,BYTEARRAY *rp)
 
 	ac = argc(arg);
 	if ( ac == 1 ) {
-		/* ARG0(arg) must be a filename */
-		asir_assert(ARG0(arg),O_STR,"newbytearray");
-		fname = BDY((STRING)ARG0(arg));
-		fp = fopen(fname,"rb");
-		if ( !fp ) error("newbytearray : fopen failed");
-		if ( stat(fname,&sbuf) < 0 ) error("newbytearray : stat failed");
-		len = sbuf.st_size;
-		MKBYTEARRAY(array,len);
-		fread(BDY(array),len,sizeof(char),fp);
+		if ( !OID((Obj)ARG0(arg)) ) error("newbytearray : invalid argument");
+		switch ( OID((Obj)ARG0(arg)) ) {
+			case O_STR:
+				fname = BDY((STRING)ARG0(arg));
+				fp = fopen(fname,"rb");
+				if ( !fp ) error("newbytearray : fopen failed");
+				if ( stat(fname,&sbuf) < 0 ) 
+					error("newbytearray : stat failed");
+				len = sbuf.st_size;
+				MKBYTEARRAY(array,len);
+				fread(BDY(array),len,sizeof(char),fp);
+				break;
+			case O_N:
+				if ( !RATN(ARG0(arg)) )
+					error("newbytearray : invalid argument");
+				len = QTOS((Q)ARG0(arg)); 
+				if ( len < 0 )
+					error("newbytearray : invalid size");
+				MKBYTEARRAY(array,len);
+				break;
+			default:
+				error("newbytearray : invalid argument");
+		}
 	} else if ( ac == 2 ) {
 		asir_assert(ARG0(arg),O_N,"newbytearray");
 		len = QTOS((Q)ARG0(arg)); 
@@ -1542,7 +1560,7 @@ int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn
 	}
 }
 
-int generic_gauss_elim_hensel_dalg(MAT mat,MAT *nmmat,Q *dn,int **rindp,int **cindp)
+int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT *nmmat,Q *dn,int **rindp,int **cindp)
 {
 	MAT bmat,xmat;
 	Q **a0,**a,**b,**x,**nm;
@@ -1566,12 +1584,9 @@ int generic_gauss_elim_hensel_dalg(MAT mat,MAT *nmmat,
 	N wn;
 	Q wq;
 	NumberField nf;
-	DP *mb;
 	DP m;
 	int col1;
 
-	nf = get_numberfield();
-	mb = nf->mb;
 	a0 = (Q **)mat->body;
 	row = mat->row; col = mat->col;
 	w = (int **)almat(row,col);
@@ -1818,13 +1833,10 @@ int intmtoratm(MAT mat,N md,MAT nm,Q *dn)
 		return 0;
 	row = mat->row; col = mat->col;
 	bshiftn(md,1,&t);
-	isqrt(t,&b);
-#if 0
 	isqrt(t,&s);
 	bshiftn(s,64,&b);
 	if ( !b )
 		b = ONEN;
-#endif
 	dn0 = ONE;
 	tmat = (N **)mat->body;
 	rmat = (Q **)nm->body;
@@ -2066,8 +2078,7 @@ void red_by_compress(int m,unsigned int *p,unsigned in
 
 void red_by_vect(int m,unsigned int *p,unsigned int *r,unsigned int hc,int len)
 {
-	register unsigned int up,lo;
-	unsigned int dmy;
+	unsigned int up,lo,dmy;
 
 	*p++ = 0; r++; len--;
 	for ( ; len; len--, r++, p++ )