[BACK]Return to array.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / builtin

Diff for /OpenXM_contrib2/asir2000/builtin/array.c between version 1.64 and 1.67

version 1.64, 2013/11/05 02:55:02 version 1.67, 2015/08/08 14:19:41
Line 45 
Line 45 
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *   *
  * $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.63 2013/09/09 07:29:25 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.66 2015/03/04 08:44:26 ohara Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "base.h"  #include "base.h"
Line 231  int lu_gauss(int *ul,ent **u,int *ll,ent **l,int n,int
Line 231  int lu_gauss(int *ul,ent **u,int *ll,ent **l,int n,int
         for ( i = 0; i < n; i++ ) ll2[i] = 0;          for ( i = 0; i < n; i++ ) ll2[i] = 0;
         for ( i = 0; i < n; i++ ) {          for ( i = 0; i < n; i++ ) {
                 fprintf(stderr,"i=%d\n",i);                  fprintf(stderr,"i=%d\n",i);
   #if defined(__MINGW32__) || defined(__MINGW64__)
                   fflush(stderr);
   #endif
                 inv = invm(u[i][0].e,mod);                  inv = invm(u[i][0].e,mod);
                 for ( k = i+1; k < n; k++ )                  for ( k = i+1; k < n; k++ )
                         if ( u[k][0].j == n-i ) {                          if ( u[k][0].j == n-i ) {
Line 895  void Pvtol(NODE arg,LIST *rp)
Line 898  void Pvtol(NODE arg,LIST *rp)
         pointer *a;          pointer *a;
         int len,i;          int len,i;
   
           if ( OID(ARG0(arg)) == O_LIST ) {
                   *rp = ARG0(arg);
                   return;
           }
         asir_assert(ARG0(arg),O_VECT,"vtol");          asir_assert(ARG0(arg),O_VECT,"vtol");
         v = (VECT)ARG0(arg); len = v->len; a = BDY(v);          v = (VECT)ARG0(arg); len = v->len; a = BDY(v);
         for ( i = len - 1, n = 0; i >= 0; i-- ) {          for ( i = len - 1, n = 0; i >= 0; i-- ) {
Line 906  void Pvtol(NODE arg,LIST *rp)
Line 913  void Pvtol(NODE arg,LIST *rp)
 void Pltov(NODE arg,VECT *rp)  void Pltov(NODE arg,VECT *rp)
 {  {
         NODE n;          NODE n;
         VECT v;          VECT v,v0;
         int len,i;          int len,i;
   
           if ( OID(ARG0(arg)) == O_VECT ) {
                   v0 = (VECT)ARG0(arg); len = v0->len;
                   MKVECT(v,len);
                   for ( i = 0; i < len; i++ ) {
                           BDY(v)[i] = BDY(v0)[i];
                   }
                   *rp = v;
                   return;
           }
         asir_assert(ARG0(arg),O_LIST,"ltov");          asir_assert(ARG0(arg),O_LIST,"ltov");
         n = (NODE)BDY((LIST)ARG0(arg));          n = (NODE)BDY((LIST)ARG0(arg));
         len = length(n);          len = length(n);
Line 1762  int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn
Line 1778  int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn
                                 }                                  }
                         }                          }
         }          }
   #if defined(__MINGW32__) || defined(__MINGW64__)
           fflush(stderr);
   #endif
 }  }
   
 int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,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)
Line 1947  int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT 
Line 1966  int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT 
                                 }                                  }
                         }                          }
         }          }
   #if defined(__MINGW32__) || defined(__MINGW64__)
           fflush(stderr);
   #endif
 }  }
   
 int f4_nocheck;  int f4_nocheck;
Line 2353  int generic_gauss_elim_mod(int **mat0,int row,int col,
Line 2375  int generic_gauss_elim_mod(int **mat0,int row,int col,
                         colstat[j] = 1;                          colstat[j] = 1;
                 if ( i != rank ) {                  if ( i != rank ) {
                         t = mat[i]; mat[i] = mat[rank]; mat[rank] = t;                          t = mat[i]; mat[i] = mat[rank]; mat[rank] = t;
                   }
                   pivot = mat[rank];
                   inv = invm(pivot[j],md);
                   for ( k = j, pk = pivot+k; k < col; k++, pk++ )
                           if ( *pk ) {
                                   if ( *pk >= (unsigned int)md )
                                           *pk %= md;
                                   DMAR(*pk,inv,0,md,*pk)
                           }
                   for ( i = rank+1; i < row; i++ ) {
                           t = mat[i];
                           if ( a = t[j] )
                                   red_by_vect(md,t+j,pivot+j,md-a,col-j);
                   }
                   rank++;
           }
           for ( j = col-1, l = rank-1; j >= 0; j-- )
                   if ( colstat[j] ) {
                           pivot = mat[l];
                           for ( i = 0; i < l; i++ ) {
                                   t = mat[i];
                                   t[j] %= md;
                                   if ( a = t[j] )
                                           red_by_vect(md,t+j,pivot+j,md-a,col-j);
                           }
                           l--;
                   }
           for ( j = 0, l = 0; l < rank; j++ )
                   if ( colstat[j] ) {
                           t = mat[l];
                           for ( k = j; k < col; k++ )
                                   if ( t[k] >= (unsigned int)md )
                                           t[k] %= md;
                           l++;
                   }
           return rank;
   }
   
   int generic_gauss_elim_mod2(int **mat0,int row,int col,int md,int *colstat,int *rowstat)
   {
           int i,j,k,l,inv,a,rank;
           unsigned int *t,*pivot,*pk;
           unsigned int **mat;
   
           for ( i = 0; i < row; i++ ) rowstat[i] = i;
           mat = (unsigned int **)mat0;
           for ( rank = 0, j = 0; j < col; j++ ) {
                   for ( i = rank; i < row; i++ )
                           mat[i][j] %= md;
                   for ( i = rank; i < row; i++ )
                           if ( mat[i][j] )
                                   break;
                   if ( i == row ) {
                           colstat[j] = 0;
                           continue;
                   } else
                           colstat[j] = 1;
                   if ( i != rank ) {
                           t = mat[i]; mat[i] = mat[rank]; mat[rank] = t;
                           k = rowstat[i]; rowstat[i] = rowstat[rank]; rowstat[rank] = k;
                 }                  }
                 pivot = mat[rank];                  pivot = mat[rank];
                 inv = invm(pivot[j],md);                  inv = invm(pivot[j],md);
Legend:
Removed from v.1.64  
changed lines
  Added in v.1.67
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>