[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.76 and 1.77

version 1.76, 2018/03/29 01:32:50 version 1.77, 2019/03/03 05:21:16
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.75 2017/09/17 02:34:02 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.76 2018/03/29 01:32:50 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "base.h"  #include "base.h"
Line 1872  int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT 
Line 1872  int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT 
       fprintf(asir_out,"LU decomposition.."); fflush(asir_out);        fprintf(asir_out,"LU decomposition.."); fflush(asir_out);
     }      }
     rank = find_lhs_and_lu_mod((unsigned int **)w,row,col,md,&rinfo,&cinfo);      rank = find_lhs_and_lu_mod((unsigned int **)w,row,col,md,&rinfo,&cinfo);
       printf("\n");
       for ( i = 0; i < row; i++ ) {
         for ( j = 0; j < col; j++ )
           printf("%d ",w[i][j]);
         printf("\n");
       }
     if ( DP_Print ) {      if ( DP_Print ) {
       fprintf(asir_out,"done.\n"); fflush(asir_out);        fprintf(asir_out,"done.\n"); fflush(asir_out);
     }      }
Line 1956  int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT 
Line 1962  int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT 
             } else              } else
               b[i][j] = 0;                b[i][j] = 0;
           }            }
           printf("\n");
           for ( i = 0; i < rank; i++ ) {
             for ( j = 0; j < ri; j++ ) {
               printexpr(CO,b[i][j]); printf(" ");
             }
             printf("\n");
           }
         get_eg(&tmp1);          get_eg(&tmp1);
         add_eg(&eg_mul,&tmp0,&tmp1);          add_eg(&eg_mul,&tmp0,&tmp1);
         /* q = q*md */          /* q = q*md */
Line 1996  int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT 
Line 2009  int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT 
                 fflush(asir_out);                  fflush(asir_out);
               }                }
               return rank;                return rank;
             }              } else
                 goto reset;
           } else {            } else {
   reset:
             period = period*3/2;              period = period*3/2;
             count = 0;              count = 0;
             nsize += period;              nsize += period;
Legend:
Removed from v.1.76  
changed lines
  Added in v.1.77
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>