version 1.10, 2005/08/02 07:16:42 |
version 1.12, 2005/10/12 14:43:36 |
|
|
/* |
/* |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dalg.c,v 1.9 2005/07/11 00:24:02 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dalg.c,v 1.11 2005/08/24 06:28:39 noro Exp $ |
*/ |
*/ |
|
|
#include "ca.h" |
#include "ca.h" |
Line 651 NODE inv_or_split_dalg(DAlg a,DAlg *c) |
|
Line 651 NODE inv_or_split_dalg(DAlg a,DAlg *c) |
|
struct order_spec *current_spec; |
struct order_spec *current_spec; |
struct oEGT eg0,eg1; |
struct oEGT eg0,eg1; |
extern struct oEGT eg_le; |
extern struct oEGT eg_le; |
|
extern int DP_Print; |
|
|
if ( !(nf=current_numberfield) ) |
if ( !(nf=current_numberfield) ) |
error("invdalg : current_numberfield is not set"); |
error("invdalg : current_numberfield is not set"); |
Line 669 NODE inv_or_split_dalg(DAlg a,DAlg *c) |
|
Line 670 NODE inv_or_split_dalg(DAlg a,DAlg *c) |
|
simp = (DAlg *)ALLOCA(dim*sizeof(DAlg)); |
simp = (DAlg *)ALLOCA(dim*sizeof(DAlg)); |
current_spec = dp_current_spec; initd(nf->spec); |
current_spec = dp_current_spec; initd(nf->spec); |
for ( i = 0; i < dim; i++ ) { |
for ( i = 0; i < dim; i++ ) { |
|
if ( DP_Print ) { fprintf(asir_out,"."); fflush(asir_out); } |
m = mb[i]; |
m = mb[i]; |
for ( j = i-1; j >= 0; j-- ) |
for ( j = i-1; j >= 0; j-- ) |
if ( dp_redble(m,mb[j]) ) |
if ( dp_redble(m,mb[j]) ) |
break; |
break; |
if ( j >= 0 ) { |
if ( j >= 0 ) { |
dp_subd(m,mb[j],&d); |
dp_subd(m,mb[j],&d); |
muld(CO,d,simp[j]->nm,&u); |
if ( simp[j] ) { |
MKDAlg(u,simp[j]->dn,t); |
muld(CO,d,simp[j]->nm,&u); |
simpdalg(t,&simp[i]); |
MKDAlg(u,simp[j]->dn,t); |
|
simpdalg(t,&simp[i]); |
|
} else |
|
simp[i] = 0; |
} else { |
} else { |
MKDAlg(m,ONE,t); |
MKDAlg(m,ONE,t); |
muldalg(t,a0,&simp[i]); |
muldalg(t,a0,&simp[i]); |
} |
} |
gcdn(NM(simp[i]->dn),ln,&gn); divsn(ln,gn,&qn); |
if ( simp[i] ) { |
muln(NM(simp[i]->dn),qn,&ln); |
gcdn(NM(simp[i]->dn),ln,&gn); divsn(ln,gn,&qn); |
|
muln(NM(simp[i]->dn),qn,&ln); |
|
} |
} |
} |
initd(current_spec); |
initd(current_spec); |
NTOQ(ln,1,dn); |
NTOQ(ln,1,dn); |
Line 691 NODE inv_or_split_dalg(DAlg a,DAlg *c) |
|
Line 698 NODE inv_or_split_dalg(DAlg a,DAlg *c) |
|
mat = (Q **)BDY(mobj); |
mat = (Q **)BDY(mobj); |
mulq(dn,a->dn,&mat[0][dim]); |
mulq(dn,a->dn,&mat[0][dim]); |
for ( j = 0; j < dim; j++ ) { |
for ( j = 0; j < dim; j++ ) { |
divq(dn,simp[j]->dn,&mul); |
if ( simp[j] ) { |
for ( i = dim-1, mp = BDY(simp[j]->nm); mp && i >= 0; i-- ) |
divq(dn,simp[j]->dn,&mul); |
if ( dl_equal(n,BDY(mb[i])->dl,mp->dl) ) { |
for ( i = dim-1, mp = BDY(simp[j]->nm); mp && i >= 0; i-- ) |
mulq(mul,(Q)mp->c,&mat[i][j]); |
if ( dl_equal(n,BDY(mb[i])->dl,mp->dl) ) { |
mp = NEXT(mp); |
mulq(mul,(Q)mp->c,&mat[i][j]); |
} |
mp = NEXT(mp); |
|
} |
|
} |
} |
} |
get_eg(&eg0); |
get_eg(&eg0); |
rank = generic_gauss_elim_hensel(mobj,&sol,&dnsol,&rinfo,&cinfo); |
rank = generic_gauss_elim_hensel(mobj,&sol,&dnsol,&rinfo,&cinfo); |