version 1.41, 2004/12/04 09:39:27 |
version 1.42, 2004/12/13 23:04: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.40 2004/12/02 13:53:31 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.41 2004/12/04 09:39:27 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
Line 1224 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
Line 1224 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
init_eg(&eg_check); init_eg(&eg_intrat); |
init_eg(&eg_check); init_eg(&eg_intrat); |
period = F4_INTRAT_PERIOD; |
period = F4_INTRAT_PERIOD; |
for ( q = ONE, count = 0; ; count++ ) { |
for ( q = ONE, count = 0; ; count++ ) { |
if ( DP_Print ) |
if ( DP_Print > 3 ) |
fprintf(stderr,"o"); |
fprintf(stderr,"o"); |
/* wc = -b mod md */ |
/* wc = -b mod md */ |
for ( i = 0; i < rank; i++ ) |
for ( i = 0; i < rank; i++ ) |
Line 1278 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
Line 1278 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
ret = gensolve_check(mat,*nmmat,*dn,rind,cind); |
ret = gensolve_check(mat,*nmmat,*dn,rind,cind); |
get_eg(&tmp1); add_eg(&eg_check,&tmp0,&tmp1); |
get_eg(&tmp1); add_eg(&eg_check,&tmp0,&tmp1); |
if ( ret ) { |
if ( ret ) { |
if ( DP_Print ) { |
if ( DP_Print > 3 ) { |
fprintf(stderr,"\n"); |
fprintf(stderr,"\n"); |
print_eg("INV",&eg_inv); |
print_eg("INV",&eg_inv); |
print_eg("MUL",&eg_mul); |
print_eg("MUL",&eg_mul); |