version 1.39, 2004/12/01 12:55:19 |
version 1.40, 2004/12/02 13:53:31 |
|
|
* 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.38 2004/09/21 05:23:13 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.39 2004/12/01 12:55:19 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
Line 1222 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
Line 1222 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
init_eg(&eg_mul); init_eg(&eg_inv); |
init_eg(&eg_mul); init_eg(&eg_inv); |
period = F4_INTRAT_PERIOD; |
period = F4_INTRAT_PERIOD; |
for ( q = ONE, count = 0; ; count++ ) { |
for ( q = ONE, count = 0; ; count++ ) { |
fprintf(stderr,"."); |
if ( DP_Print ) |
|
fprintf(stderr,"."); |
/* wc = -b mod md */ |
/* wc = -b mod md */ |
for ( i = 0; i < rank; i++ ) |
for ( i = 0; i < rank; i++ ) |
for ( j = 0, bi = b[i], wi = wc[i]; j < ri; j++ ) |
for ( j = 0, bi = b[i], wi = wc[i]; j < ri; j++ ) |
Line 1269 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
Line 1270 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
else |
else |
cind[l++] = j; |
cind[l++] = j; |
if ( gensolve_check(mat,*nmmat,*dn,rind,cind) ) { |
if ( gensolve_check(mat,*nmmat,*dn,rind,cind) ) { |
fprintf(stderr,"\n"); |
if ( DP_Print ) { |
print_eg("INV",&eg_inv); |
fprintf(stderr,"\n"); |
print_eg("MUL",&eg_mul); |
print_eg("INV",&eg_inv); |
fflush(asir_out); |
print_eg("MUL",&eg_mul); |
|
fflush(asir_out); |
|
} |
return rank; |
return rank; |
} |
} |
} else |
} else |