version 1.76, 2018/03/29 01:32:50 |
version 1.77, 2019/03/03 05:21:16 |
|
|
* 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; |