version 1.38, 2004/09/21 05:23:13 |
version 1.39, 2004/12/01 12:55:19 |
|
|
* 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.37 2004/09/15 01:43:32 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.38 2004/09/21 05:23:13 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
Line 1172 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
Line 1172 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
int *rind,*cind; |
int *rind,*cind; |
int count; |
int count; |
struct oEGT eg_mul,eg_inv,tmp0,tmp1; |
struct oEGT eg_mul,eg_inv,tmp0,tmp1; |
|
int period; |
|
|
a0 = (Q **)mat->body; |
a0 = (Q **)mat->body; |
row = mat->row; col = mat->col; |
row = mat->row; col = mat->col; |
Line 1219 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
Line 1220 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
*cindp = cind = (int *)MALLOC_ATOMIC((ri)*sizeof(int)); |
*cindp = cind = (int *)MALLOC_ATOMIC((ri)*sizeof(int)); |
|
|
init_eg(&eg_mul); init_eg(&eg_inv); |
init_eg(&eg_mul); init_eg(&eg_inv); |
|
period = F4_INTRAT_PERIOD; |
for ( q = ONE, count = 0; ; count++ ) { |
for ( q = ONE, count = 0; ; count++ ) { |
fprintf(stderr,"."); |
fprintf(stderr,"."); |
/* wc = -b mod md */ |
/* wc = -b mod md */ |
Line 1259 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
Line 1261 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
add_eg(&eg_mul,&tmp0,&tmp1); |
add_eg(&eg_mul,&tmp0,&tmp1); |
/* q = q*md */ |
/* q = q*md */ |
mulq(q,mdq,&u); q = u; |
mulq(q,mdq,&u); q = u; |
if ( !(count % F4_INTRAT_PERIOD) && intmtoratm_q(xmat,NM(q),*nmmat,dn) ) { |
if ( !(count % period) ) |
for ( j = k = l = 0; j < col; j++ ) |
if ( intmtoratm_q(xmat,NM(q),*nmmat,dn) ) { |
if ( cinfo[j] ) |
for ( j = k = l = 0; j < col; j++ ) |
rind[k++] = j; |
if ( cinfo[j] ) |
else |
rind[k++] = j; |
cind[l++] = j; |
else |
if ( gensolve_check(mat,*nmmat,*dn,rind,cind) ) { |
cind[l++] = j; |
fprintf(stderr,"\n"); |
if ( gensolve_check(mat,*nmmat,*dn,rind,cind) ) { |
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); |
return rank; |
fflush(asir_out); |
} |
return rank; |
} |
} |
|
} else |
|
period *=2; |
} |
} |
} |
} |
} |
} |