version 1.3, 2004/12/02 13:48:43 |
version 1.4, 2004/12/03 07:16:34 |
|
|
/* |
/* |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dalg.c,v 1.2 2004/12/02 08:39:54 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dalg.c,v 1.3 2004/12/02 13:48:43 noro Exp $ |
*/ |
*/ |
|
|
#include "ca.h" |
#include "ca.h" |
Line 58 void setfield_dalg(NODE alist) |
|
Line 58 void setfield_dalg(NODE alist) |
|
mb[i] = (DP)BDY(t); |
mb[i] = (DP)BDY(t); |
} |
} |
|
|
|
void qtodalg(Q q,DAlg *r) |
|
{ |
|
NumberField nf; |
|
Q t; |
|
DP nm; |
|
|
|
if ( !(nf=current_numberfield) ) |
|
error("qtodalg : current_numberfield is not set"); |
|
if ( !q ) |
|
*r = 0; |
|
else if ( NID(q) == N_DA ) |
|
*r = (DAlg)q; |
|
else if ( NID(q) == N_Q ) { |
|
if ( INT(q) ) { |
|
muldc(CO,nf->one->nm,(P)q,&nm); |
|
MKDAlg(nm,ONE,*r); |
|
} else { |
|
NTOQ(NM(q),SGN(q),t); |
|
muldc(CO,nf->one->nm,(P)t,&nm); |
|
NTOQ(DN(q),1,t); |
|
MKDAlg(nm,t,*r); |
|
} |
|
} else |
|
error("qtodalg : invalid argument"); |
|
} |
|
|
void algtodalg(Alg a,DAlg *r) |
void algtodalg(Alg a,DAlg *r) |
{ |
{ |
P ap,p,p1; |
P ap,p,p1; |
Line 157 void adddalg(DAlg a,DAlg b,DAlg *c) |
|
Line 183 void adddalg(DAlg a,DAlg b,DAlg *c) |
|
NumberField nf; |
NumberField nf; |
Q dna,dnb,a1,b1,dn,g; |
Q dna,dnb,a1,b1,dn,g; |
N an,bn,gn; |
N an,bn,gn; |
|
DAlg t; |
DP ta,tb,nm; |
DP ta,tb,nm; |
struct order_spec *current_spec; |
struct order_spec *current_spec; |
|
|
Line 167 void adddalg(DAlg a,DAlg b,DAlg *c) |
|
Line 194 void adddalg(DAlg a,DAlg b,DAlg *c) |
|
else if ( !b ) |
else if ( !b ) |
*c = a; |
*c = a; |
else { |
else { |
|
qtodalg((Q)a,&t); a = t; qtodalg((Q)b,&t); b = t; |
dna = a->dn; |
dna = a->dn; |
dnb = b->dn; |
dnb = b->dn; |
gcdn(NM(dna),NM(dnb),&gn); |
gcdn(NM(dna),NM(dnb),&gn); |
Line 192 void subdalg(DAlg a,DAlg b,DAlg *c) |
|
Line 220 void subdalg(DAlg a,DAlg b,DAlg *c) |
|
Q dna,dnb,a1,b1,dn,g; |
Q dna,dnb,a1,b1,dn,g; |
N an,bn,gn; |
N an,bn,gn; |
DP ta,tb,nm; |
DP ta,tb,nm; |
|
DAlg t; |
struct order_spec *current_spec; |
struct order_spec *current_spec; |
|
|
if ( !(nf=current_numberfield) ) |
if ( !(nf=current_numberfield) ) |
Line 201 void subdalg(DAlg a,DAlg b,DAlg *c) |
|
Line 230 void subdalg(DAlg a,DAlg b,DAlg *c) |
|
else if ( !b ) |
else if ( !b ) |
*c = a; |
*c = a; |
else { |
else { |
|
qtodalg((Q)a,&t); a = t; qtodalg((Q)b,&t); b = t; |
dna = a->dn; |
dna = a->dn; |
dnb = b->dn; |
dnb = b->dn; |
gcdn(NM(dna),NM(dnb),&gn); |
gcdn(NM(dna),NM(dnb),&gn); |
Line 233 void muldalg(DAlg a,DAlg b,DAlg *c) |
|
Line 263 void muldalg(DAlg a,DAlg b,DAlg *c) |
|
if ( !a || !b ) |
if ( !a || !b ) |
*c = 0; |
*c = 0; |
else { |
else { |
|
qtodalg((Q)a,&t); a = t; qtodalg((Q)b,&t); b = t; |
current_spec = dp_current_spec; initd(nf->spec); |
current_spec = dp_current_spec; initd(nf->spec); |
muld(CO,a->nm,b->nm,&nm); |
muld(CO,a->nm,b->nm,&nm); |
initd(current_spec); |
initd(current_spec); |
Line 245 void muldalg(DAlg a,DAlg b,DAlg *c) |
|
Line 276 void muldalg(DAlg a,DAlg b,DAlg *c) |
|
|
|
void divdalg(DAlg a,DAlg b,DAlg *c) |
void divdalg(DAlg a,DAlg b,DAlg *c) |
{ |
{ |
DAlg inv; |
DAlg inv,t; |
|
|
if ( !current_numberfield ) |
if ( !current_numberfield ) |
error("divdalg : current_numberfield is not set"); |
error("divdalg : current_numberfield is not set"); |
Line 254 void divdalg(DAlg a,DAlg b,DAlg *c) |
|
Line 285 void divdalg(DAlg a,DAlg b,DAlg *c) |
|
if ( !a ) |
if ( !a ) |
c = 0; |
c = 0; |
else { |
else { |
|
qtodalg((Q)a,&t); a = t; qtodalg((Q)b,&t); b = t; |
invdalg(b,&inv); |
invdalg(b,&inv); |
muldalg(a,inv,c); |
muldalg(a,inv,c); |
} |
} |
Line 298 void invdalg(DAlg a,DAlg *c) |
|
Line 330 void invdalg(DAlg a,DAlg *c) |
|
error("invdalg : current_numberfield is not set"); |
error("invdalg : current_numberfield is not set"); |
if ( !a ) |
if ( !a ) |
error("invdalg : division by 0"); |
error("invdalg : division by 0"); |
|
else if ( NID(a) == N_Q ) { |
|
invq((Q)a,&dn); *c = (DAlg)dn; |
|
return; |
|
} |
dim = nf->dim; |
dim = nf->dim; |
mb = nf->mb; |
mb = nf->mb; |
n = nf->n; |
n = nf->n; |
Line 351 void invdalg(DAlg a,DAlg *c) |
|
Line 387 void invdalg(DAlg a,DAlg *c) |
|
void chsgndalg(DAlg a,DAlg *c) |
void chsgndalg(DAlg a,DAlg *c) |
{ |
{ |
DP nm; |
DP nm; |
|
Q t; |
|
|
if ( !a ) *c = 0; |
if ( !a ) *c = 0; |
else { |
else if ( NID(a) == N_Q ) { |
|
chsgnq((Q)a,&t); *c = (DAlg)t; |
|
} else { |
chsgnd(a->nm,&nm); |
chsgnd(a->nm,&nm); |
MKDAlg(nm,a->dn,*c); |
MKDAlg(nm,a->dn,*c); |
} |
} |
Line 363 void pwrdalg(DAlg a,Q e,DAlg *c) |
|
Line 402 void pwrdalg(DAlg a,Q e,DAlg *c) |
|
{ |
{ |
NumberField nf; |
NumberField nf; |
DAlg t,z,y; |
DAlg t,z,y; |
|
Q q; |
N en,qn; |
N en,qn; |
int r; |
int r; |
|
|
if ( !(nf=current_numberfield) ) |
if ( !(nf=current_numberfield) ) |
error("pwrdalg : current_numberfield is not set"); |
error("pwrdalg : current_numberfield is not set"); |
if ( !e ) |
if ( !a ) |
|
*c = !e ? (DAlg)ONE : 0; |
|
else if ( NID(a) == N_Q ) { |
|
pwrq((Q)a,e,&q); *c = (DAlg)q; |
|
} else if ( !e ) |
*c = nf->one; |
*c = nf->one; |
else if ( UNIQ(e) ) |
else if ( UNIQ(e) ) |
*c = a; |
*c = a; |