===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2018/builtin/array.c,v
retrieving revision 1.4
retrieving revision 1.10
diff -u -p -r1.4 -r1.10
--- OpenXM_contrib2/asir2018/builtin/array.c	2018/10/19 23:27:38	1.4
+++ OpenXM_contrib2/asir2018/builtin/array.c	2022/01/13 08:15:02	1.10
@@ -45,7 +45,7 @@
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *
- * $OpenXM: OpenXM_contrib2/asir2018/builtin/array.c,v 1.3 2018/10/01 05:49:06 noro Exp $
+ * $OpenXM: OpenXM_contrib2/asir2018/builtin/array.c,v 1.9 2021/03/26 09:05:41 ohara Exp $
 */
 #include "ca.h"
 #include "base.h"
@@ -170,6 +170,13 @@ void solve_u(int *,ent **,int,int *,int);
 static int *ul,*ll;
 static ent **u,**l;
 static int modulus;
+#if defined(ANDROID)
+int getw(FILE *fp)
+{
+	int x;
+	return (fread((void *)&x, sizeof(x), 1, fp) == 1 ? x : EOF);
+}
+#endif
 
 void Plusolve_prep(NODE arg,Q *rp)
 {
@@ -615,7 +622,7 @@ void Pnewvect(NODE arg,VECT *rp)
       return;
     }
 #endif
-    for ( i = 0, tn = BDY(list), vb = BDY(vect); tn; i++, tn = NEXT(tn) )
+    for ( i = 0, tn = BDY(list), vb = BDY(vect); tn && i<len; i++, tn = NEXT(tn) )
       vb[i] = (pointer)BDY(tn);
   }
   *rp = vect;
@@ -1164,7 +1171,7 @@ void Pgeneric_gauss_elim(NODE arg,LIST *rp)
   if ( is_hensel )
     rank = generic_gauss_elim_hensel(m,&nm,&dn,&ri,&ci);
   else
-    rank = generic_gauss_elim64(m,&nm,&dn,&ri,&ci);
+    rank = generic_gauss_elim(m,&nm,&dn,&ri,&ci);
   t = col-rank;
   MKVECT(rind,rank);
   MKVECT(cind,t);
@@ -1243,9 +1250,12 @@ void Pgeneric_gauss_elim_mod64(NODE arg,LIST *rp)
   Z q;
   mp_limb_t md;
   int i,j,k,l,row,col,t,rank;
+  Obj val;
+  int asis = 0;
 
   asir_assert(ARG0(arg),O_MAT,"generic_gauss_elim_mod64");
   asir_assert(ARG1(arg),O_N,"generic_gauss_elim_mod64");
+  if ( get_opt("asis",&val) && val ) asis = 1;
   m = (MAT)ARG0(arg); md = ZTOS((Z)ARG1(arg));
   row = m->row; col = m->col; tmat = (Z **)m->body;
   wmat = (mp_limb_t **)almat64(row,col);
@@ -1258,6 +1268,16 @@ void Pgeneric_gauss_elim_mod64(NODE arg,LIST *rp)
     for ( j = 0; j < col; j++ )
       wmat[i][j] = remqi64((Q)tmat[i][j],md);
   rank = generic_gauss_elim_mod64(wmat,row,col,md,colstat);
+  if ( asis ) {
+    MKMAT(mat,row,col);
+    tmat = (Z **)mat->body;
+    for ( i = 0; i < rank; i++ )
+      for ( j = 0; j < col; j++ ) {
+        UTOZ(wmat[i][j],tmat[i][j]);
+      } 
+    *rp = (LIST)mat;
+    return;
+  }
 
   MKVECT(rnum,rank);
   rnb = (Z *)rnum->body;
@@ -1300,9 +1320,12 @@ void Pgeneric_gauss_elim_mod(NODE arg,LIST *rp)
   Z q;
   long mdl;
   int md,i,j,k,l,row,col,t,rank;
+  Obj val;
+  int asis = 0;
 
   asir_assert(ARG0(arg),O_MAT,"generic_gauss_elim_mod");
   asir_assert(ARG1(arg),O_N,"generic_gauss_elim_mod");
+  if ( get_opt("asis",&val) && val ) asis = 1;
 #if SIZEOF_LONG==8
   mdl = ZTOS((Z)ARG1(arg));
   if ( mdl >= ((mp_limb_t)1)<<32 ) {
@@ -1323,6 +1346,16 @@ void Pgeneric_gauss_elim_mod(NODE arg,LIST *rp)
     for ( j = 0; j < col; j++ )
       wmat[i][j] = remqi((Q)tmat[i][j],md);
   rank = generic_gauss_elim_mod(wmat,row,col,md,colstat);
+  if ( asis ) {
+    MKMAT(mat,row,col);
+    tmat = (Z **)mat->body;
+    for ( i = 0; i < rank; i++ )
+      for ( j = 0; j < col; j++ ) {
+        UTOZ(wmat[i][j],tmat[i][j]);
+      } 
+    *rp = (LIST)mat;
+    return;
+  }
 
   MKVECT(rnum,rank);
   rnb = (Z *)rnum->body;
@@ -1423,7 +1456,7 @@ int gauss_elim_mod(int **mat,int row,int col,int md)
     return -1;
 }
 
-struct oEGT eg_mod,eg_elim,eg_elim1,eg_elim2,eg_chrem,eg_gschk,eg_intrat,eg_symb;
+struct oEGT eg_mod,eg_elim,eg_elim1,eg_elim2,eg_chrem,eg_gschk,eg_intrat,eg_symb,eg_back,eg_fore;
 struct oEGT eg_conv;
 
 #if 0
@@ -1536,7 +1569,7 @@ void lu_dec_cr(MAT mat,MAT lu,Q *dn,int **perm)
 
 int f4_nocheck;
 
-#define ONE_STEP1  if ( zzz = *s ) { DMAR(zzz,hc,*tj,md,*tj) } tj++; s++;
+#define ONE_STEP1  if ( ( zzz = *s ) != 0 ) { DMAR(zzz,hc,*tj,md,*tj) } tj++; s++;
 
 void reduce_reducers_mod(int **mat,int row,int col,int md)
 {
@@ -1552,7 +1585,7 @@ void reduce_reducers_mod(int **mat,int row,int col,int
     ind[i] = j;
     for ( l = i-1; l >= 0; l-- ) {
       /* reduce mat[i] by mat[l] */
-      if ( hc = t[ind[l]] ) {
+      if ( ( hc = t[ind[l]] ) != 0 ) {
         /* mat[i] = mat[i]-hc*mat[l] */
         j = ind[l];
         s = mat[l]+j;
@@ -1578,7 +1611,7 @@ void reduce_reducers_mod(int **mat,int row,int col,int
           ONE_STEP1 ONE_STEP1 ONE_STEP1 ONE_STEP1
         }
         for ( ; k > 0; k-- ) {
-          if ( zzz = *s ) { DMAR(zzz,hc,*tj,md,*tj) } tj++; s++;
+          if ( ( zzz = *s ) != 0 ) { DMAR(zzz,hc,*tj,md,*tj) } tj++; s++;
         }
       }
     }
@@ -1613,7 +1646,7 @@ void pre_reduce_mod(int **mat,int row,int col,int nred
         DMAR(t[k],inv,0,md,t[k])
     for ( l = i-1; l >= 0; l-- ) {
       /* reduce mat[i] by mat[l] */
-      if ( hc = t[ind[l]] ) {
+      if ( ( hc = t[ind[l]] ) != 0 ) {
         /* mat[i] = mat[i]-hc*mat[l] */
         for ( k = ind[l], hc = md-hc, s = mat[l]+k, tk = t+k;
           k < col; k++, tk++, s++ )
@@ -1627,7 +1660,7 @@ void pre_reduce_mod(int **mat,int row,int col,int nred
     t = mat[i];
     for ( l = nred-1; l >= 0; l-- ) {
       /* reduce mat[i] by mat[l] */
-      if ( hc = t[ind[l]] ) {
+      if ( ( hc = t[ind[l]] ) != 0 ) {
         /* mat[i] = mat[i]-hc*mat[l] */
         for ( k = ind[l], hc = md-hc, s = mat[l]+k, tk = t+k;
           k < col; k++, tk++, s++ )
@@ -1651,14 +1684,14 @@ void reduce_sp_by_red_mod(int *sp,int **redmat,int *in
   /* reduce the spolys by redmat */
   for ( i = nred-1; i >= 0; i-- ) {
     /* reduce sp by redmat[i] */
-    if ( hc = sp[ind[i]] ) {
+    if ( ( hc = sp[ind[i]] ) != 0 ) {
       /* sp = sp-hc*redmat[i] */
       j = ind[i];
       hc = md-hc;
       s = redmat[i]+j;
       tj = sp+j;
       for ( k = col-j; k > 0; k-- ) {
-        if ( zzz = *s ) { DMAR(zzz,hc,*tj,md,*tj) } tj++; s++;
+        if ( ( zzz = *s ) != 0 ) { DMAR(zzz,hc,*tj,md,*tj) } tj++; s++;
       }
     }
   }
@@ -1756,7 +1789,7 @@ int generic_gauss_elim_mod64(mp_limb_t **mat,int row,i
         *pk = mulmod64(*pk,inv,md);
     for ( i = rank+1; i < row; i++ ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         red_by_vect64mod(md,t+j,pivot+j,md-a,col-j);
     }
     rank++;
@@ -1766,7 +1799,7 @@ int generic_gauss_elim_mod64(mp_limb_t **mat,int row,i
       pivot = mat[l];
       for ( i = 0; i < l; i++ ) {
         t = mat[i];
-        if ( a = t[j] )
+        if ( ( a = t[j] ) != 0 )
           red_by_vect64mod(md,t+j,pivot+j,md-a,col-j);
       }
       l--;
@@ -1937,7 +1970,7 @@ void reduce_sp_by_red_mod_compress (int *sp,CDP *redma
   for ( i = nred-1; i >= 0; i-- ) {
     /* reduce sp by redmat[i] */
     usp[ind[i]] %= md;
-    if ( hc = usp[ind[i]] ) {
+    if ( ( hc = usp[ind[i]] ) != 0 ) {
       /* sp = sp-hc*redmat[i] */
       hc = md-hc;
       ri = redmat[i];
@@ -1983,7 +2016,7 @@ int generic_gauss_elim_mod(int **mat0,int row,int col,
       }
     for ( i = rank+1; i < row; i++ ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         red_by_vect(md,t+j,pivot+j,md-a,col-j);
     }
     rank++;
@@ -1994,7 +2027,7 @@ int generic_gauss_elim_mod(int **mat0,int row,int col,
       for ( i = 0; i < l; i++ ) {
         t = mat[i];
         t[j] %= md;
-        if ( a = t[j] )
+        if ( ( a = t[j] ) != 0 )
           red_by_vect(md,t+j,pivot+j,md-a,col-j);
       }
       l--;
@@ -2043,7 +2076,7 @@ int generic_gauss_elim_mod2(int **mat0,int row,int col
       }
     for ( i = rank+1; i < row; i++ ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         red_by_vect(md,t+j,pivot+j,md-a,col-j);
     }
     rank++;
@@ -2054,7 +2087,7 @@ int generic_gauss_elim_mod2(int **mat0,int row,int col
       for ( i = 0; i < l; i++ ) {
         t = mat[i];
         t[j] %= md;
-        if ( a = t[j] )
+        if ( ( a = t[j] ) != 0 )
           red_by_vect(md,t+j,pivot+j,md-a,col-j);
       }
       l--;
@@ -2099,7 +2132,7 @@ int indep_rows_mod(int **mat0,int row,int col,int md,i
       }
     for ( i = rank+1; i < row; i++ ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         red_by_vect(md,t+j,pivot+j,md-a,col-j);
     }
     rank++;
@@ -2133,7 +2166,7 @@ int generic_gauss_elim_sf(int **mat0,int row,int col,i
         *pk = _mulsf(*pk,inv);
     for ( i = rank+1; i < row; i++ ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         red_by_vect_sf(md,t+j,pivot+j,_chsgnsf(a),col-j);
     }
     rank++;
@@ -2143,7 +2176,7 @@ int generic_gauss_elim_sf(int **mat0,int row,int col,i
       pivot = mat[l];
       for ( i = 0; i < l; i++ ) {
         t = mat[i];
-        if ( a = t[j] )
+        if ( ( a = t[j] ) != 0 )
           red_by_vect_sf(md,t+j,pivot+j,_chsgnsf(a),col-j);
       }
       l--;
@@ -2179,7 +2212,7 @@ int lu_gfmmat(GFMMAT mat,unsigned int md,int *perm)
     pivot[k] = inv = invm(pivot[k],md);
     for ( i = k+1; i < row; i++ ) {
       t = a[i];
-      if ( m = t[k] ) {
+      if ( ( m = t[k] ) != 0 ) {
         DMAR(inv,m,0,md,t[k])
         for ( j = k+1, m = md - t[k]; j < col; j++ )
           if ( pivot[j] ) {
@@ -2235,7 +2268,7 @@ int find_lhs_and_lu_mod(unsigned int **a,int row,int c
     pivot[k] = inv = invm(pivot[k],md);
     for ( i = d+1; i < row; i++ ) {
       t = a[i];
-      if ( m = t[k] ) {
+      if ( ( m = t[k] ) != 0 ) {
         DMAR(inv,m,0,md,t[k])
         for ( j = k+1, m = md - t[k]; j < col; j++ )
           if ( pivot[j] ) {
@@ -2270,7 +2303,7 @@ int lu_mod(unsigned int **a,int n,unsigned int md,int 
     inv = invm(pivot[k],md);
     for ( i = k+1; i < n; i++ ) {
       t = a[i];
-      if ( m = t[k] ) {
+      if ( ( m = t[k] ) != 0 ) {
         DMAR(inv,m,0,md,t[k])
         for ( j = k+1, m = md - t[k]; j < n; j++ )
           if ( pivot[j] ) {
@@ -2299,8 +2332,8 @@ void solve_by_lu_mod(int **a,int n,int md,int **b,int 
   int i,j,k;
   unsigned int t,m,m2;
 
-  y = (int *)MALLOC_ATOMIC(n*sizeof(int));
-  c = (int *)MALLOC_ATOMIC(n*sizeof(int));
+  y = (unsigned int *)MALLOC_ATOMIC(n*sizeof(int));
+  c = (unsigned int *)MALLOC_ATOMIC(n*sizeof(int));
   m2 = md>>1;
   for ( k = 0; k < l; k++ ) {
     /* copy b[.][k] to c */
@@ -2460,7 +2493,7 @@ int gauss_elim_geninv_mod(unsigned int **mat,int row,i
       pivot[k] = dmar(pivot[k],inv,0,md);
     for ( i = j+1; i < row; i++ ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0  )
         for ( k = j, a = md - a; k < m; k++ )
           t[k] = dmar(pivot[k],a,t[k],md);
     }
@@ -2469,7 +2502,7 @@ int gauss_elim_geninv_mod(unsigned int **mat,int row,i
     pivot = mat[j];
     for ( i = j-1; i >= 0; i-- ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         for ( k = j, a = md - a; k < m; k++ )
           t[k] = dmar(pivot[k],a,t[k],md);
     }
@@ -2687,7 +2720,7 @@ int gauss_elim_geninv_mod_swap(unsigned int **mat,int 
         pivot[k] = (unsigned int)dmar(pivot[k],inv,0,md);
     for ( i = j+1; i < row; i++ ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         for ( k = j, a = md - a; k < m; k++ )
           if ( pivot[k] )
             t[k] = dmar(pivot[k],a,t[k],md);
@@ -2697,7 +2730,7 @@ int gauss_elim_geninv_mod_swap(unsigned int **mat,int 
     pivot = mat[j];
     for ( i = j-1; i >= 0; i-- ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         for ( k = j, a = md - a; k < m; k++ )
           if ( pivot[k] )
             t[k] = dmar(pivot[k],a,t[k],md);
@@ -2733,7 +2766,7 @@ void Pgeninv_sf_swap(NODE arg,LIST *rp)
   for ( i = 0; i < row; i++ ) {
     bzero((char *)wmat[i],(col+row)*sizeof(int));
     for ( j = 0; j < col; j++ )
-      if ( q = (GFS)mat[i][j] )
+      if ( ( q = (GFS)mat[i][j] ) != 0 )
         wmat[i][j] = FTOIF(CONT(q));
     wmat[i][col+i] = _onesf();
   }
@@ -2744,7 +2777,7 @@ void Pgeninv_sf_swap(NODE arg,LIST *rp)
     MKMAT(mat1,col,col);
     for ( i = 0, tmat = (GFS **)mat1->body; i < col; i++ )
       for ( j = 0; j < col; j++ )
-        if ( t = invmat[i][j] ) {
+        if ( ( t = invmat[i][j] ) != 0 ) {
           MKGFS(IFTOF(t),tmat[i][j]);
         }
     MKVECT(vect1,row);
@@ -2781,7 +2814,7 @@ int gauss_elim_geninv_sf_swap(int **mat,int row,int co
         pivot[k] = _mulsf(pivot[k],inv);
     for ( i = j+1; i < row; i++ ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         for ( k = j, a = _chsgnsf(a); k < m; k++ )
           if ( pivot[k] ) {
             u = _mulsf(pivot[k],a);
@@ -2793,7 +2826,7 @@ int gauss_elim_geninv_sf_swap(int **mat,int row,int co
     pivot = mat[j];
     for ( i = j-1; i >= 0; i-- ) {
       t = mat[i];
-      if ( a = t[j] )
+      if ( ( a = t[j] ) != 0 )
         for ( k = j, a = _chsgnsf(a); k < m; k++ )
           if ( pivot[k] ) {
             u = _mulsf(pivot[k],a);
@@ -2966,7 +2999,7 @@ int generate_ONB_polynomial(UP2 *rp,int m,int type)
       for ( i = 0; i < w; i++ )
         f->b[i] = 0xffffffff;
       /* mask the top word if necessary */
-      if ( r = (m+1)&31 )
+      if ( ( r = (m+1)&31 ) != 0 )
         f->b[w-1] &= (1<<r)-1;
       return 0;
       break;