version 1.2, 2001/09/04 02:45:33 |
version 1.7, 2004/03/03 09:25:30 |
|
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/parse/quote.c,v 1.6 2001/10/09 01:36:25 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
|
|
void addquote(vl,a,b,c) |
void addquote(VL vl,QUOTE a,QUOTE b,QUOTE *c) |
VL vl; |
|
QUOTE a,b; |
|
QUOTE *c; |
|
{ |
{ |
FNODE fn; |
FNODE fn; |
|
|
|
|
MKQUOTE(*c,fn); |
MKQUOTE(*c,fn); |
} |
} |
|
|
void subquote(vl,a,b,c) |
void subquote(VL vl,QUOTE a,QUOTE b,QUOTE *c) |
VL vl; |
|
QUOTE a,b; |
|
QUOTE *c; |
|
{ |
{ |
FNODE fn; |
FNODE fn; |
|
|
|
|
MKQUOTE(*c,fn); |
MKQUOTE(*c,fn); |
} |
} |
|
|
void mulquote(vl,a,b,c) |
void mulquote(VL vl,QUOTE a,QUOTE b,QUOTE *c) |
VL vl; |
|
QUOTE a,b; |
|
QUOTE *c; |
|
{ |
{ |
FNODE fn; |
FNODE fn; |
|
|
|
|
MKQUOTE(*c,fn); |
MKQUOTE(*c,fn); |
} |
} |
|
|
void divquote(vl,a,b,c) |
void divquote(VL vl,QUOTE a,QUOTE b,QUOTE *c) |
VL vl; |
|
QUOTE a,b; |
|
QUOTE *c; |
|
{ |
{ |
FNODE fn; |
FNODE fn; |
|
|
|
|
MKQUOTE(*c,fn); |
MKQUOTE(*c,fn); |
} |
} |
|
|
void pwrquote(vl,a,b,c) |
void pwrquote(VL vl,QUOTE a,QUOTE b,QUOTE *c) |
VL vl; |
|
QUOTE a,b; |
|
QUOTE *c; |
|
{ |
{ |
FNODE fn; |
FNODE fn; |
|
|
|
|
MKQUOTE(*c,fn); |
MKQUOTE(*c,fn); |
} |
} |
|
|
void chsgnquote(a,c) |
void chsgnquote(QUOTE a,QUOTE *c) |
QUOTE a; |
|
QUOTE *c; |
|
{ |
{ |
FNODE fn; |
FNODE fn; |
|
|
|
|
MKQUOTE(*c,fn); |
MKQUOTE(*c,fn); |
} |
} |
|
|
void polytoquote(), dctoquote(), vartoquote(); |
void objtoquote(Obj a,QUOTE *c) |
|
|
void objtoquote(a,c) |
|
Obj a; |
|
QUOTE *c; |
|
{ |
{ |
QUOTE nm,dn; |
QUOTE nm,dn; |
|
NODE arg,t0,t,t1,t2,t3; |
|
FNODE fn; |
|
Obj obj; |
|
Obj *b; |
|
Obj **m; |
|
int i,j,len,row,col; |
|
Q q,qrow,qcol; |
|
FUNC f; |
|
|
if ( !a ) { |
if ( !a ) { |
MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a)); |
MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a)); |
|
|
} |
} |
switch ( OID(a) ) { |
switch ( OID(a) ) { |
case O_N: |
case O_N: |
|
if ( negative_number((Num)a) ) { |
|
arf_chsgn(a,&obj); |
|
MKQUOTE(*c,mkfnode(1,I_MINUS, |
|
mkfnode(1,I_FORMULA,(pointer)obj))); |
|
} else { |
|
MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a)); |
|
} |
|
break; |
|
case O_STR: |
MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a)); |
MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a)); |
break; |
break; |
case O_P: |
case O_P: |
|
|
polytoquote(DN((R)a),&dn); |
polytoquote(DN((R)a),&dn); |
divquote(CO,nm,dn,c); |
divquote(CO,nm,dn,c); |
break; |
break; |
|
case O_LIST: |
|
t0 = 0; |
|
for ( arg = BDY((LIST)a); arg; arg = NEXT(arg) ) { |
|
NEXTNODE(t0,t); |
|
objtoquote(BDY(arg),&nm); |
|
BDY(t) = BDY(nm); |
|
} |
|
if ( t0 ) |
|
NEXT(t) = 0; |
|
MKQUOTE(*c,mkfnode(1,I_LIST,t0)); |
|
break; |
|
case O_VECT: |
|
len = ((VECT)a)->len; |
|
b = (Obj *)BDY(((VECT)a)); |
|
t = 0; |
|
for ( i = len-1; i >= 0; i-- ) { |
|
objtoquote(b[i],&nm); |
|
MKNODE(t1,BDY(nm),t); |
|
t = t1; |
|
} |
|
STOQ(len,q); |
|
t = mknode(2,mkfnode(1,I_FORMULA,q),mkfnode(1,I_LIST,t)); |
|
gen_searchf("vector",&f); |
|
MKQUOTE(*c,mkfnode(2,I_FUNC,f,mkfnode(1,I_LIST,t))); |
|
break; |
|
case O_MAT: |
|
row = ((MAT)a)->row; |
|
col = ((MAT)a)->row; |
|
m = (Obj **)BDY(((MAT)a)); |
|
t2 = 0; |
|
for ( i = row-1; i >= 0; i-- ) { |
|
t = 0; |
|
for ( j = col-1; j >= 0; j-- ) { |
|
objtoquote(m[i][j],&nm); |
|
MKNODE(t1,BDY(nm),t); |
|
t = t1; |
|
} |
|
fn = mkfnode(1,I_LIST,t); |
|
MKNODE(t3,fn,t2); |
|
t2 = t3; |
|
} |
|
fn = mkfnode(1,I_LIST,t2); |
|
|
|
STOQ(row,qrow); |
|
STOQ(col,qcol); |
|
t = mknode(3, |
|
mkfnode(1,I_FORMULA,qrow),mkfnode(1,I_FORMULA,qcol),fn); |
|
gen_searchf("matrix",&f); |
|
MKQUOTE(*c,mkfnode(2,I_FUNC,f,mkfnode(1,I_LIST,t))); |
|
break; |
|
case O_DP: |
|
dptoquote((DP)a,c); |
|
break; |
case O_QUOTE: |
case O_QUOTE: |
*c = (QUOTE)a; |
*c = (QUOTE)a; |
break; |
break; |
|
|
} |
} |
} |
} |
|
|
void polytoquote(a,c) |
void polytoquote(P a,QUOTE *c) |
P a; |
|
QUOTE *c; |
|
{ |
{ |
DCP dc,t; |
DCP dc,t; |
DCP *dca; |
DCP *dca; |
int n,i; |
int n,i,sgn; |
QUOTE v,r,s,u; |
QUOTE v,r,s,u; |
|
|
if ( !a || (OID(a) == O_N) ) { |
if ( !a || (OID(a) == O_N) ) { |
MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a)); |
MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a)); |
return; |
return; |
} |
} |
dc = DC((P)a); |
|
vartoquote(VR((P)a),&v); |
vartoquote(VR((P)a),&v); |
for ( t = dc, n = 0; t; t = NEXT(t), n++ ); |
dc = DC((P)a); |
dca = (DCP *)ALLOCA(n*sizeof(DCP)); |
dctoquote(dc,v,&r,&sgn); |
for ( t = dc, i = 0; t; t = NEXT(t), i++ ) |
if ( sgn == -1 ) { |
dca[i] = t; |
MKQUOTE(u,mkfnode(1,I_MINUS,BDY(r))); |
dctoquote(dca[n-1],v,&r); |
|
for ( i = n-2; i >= 0; i-- ) { |
|
dctoquote(dca[i],v,&s); |
|
addquote(CO,s,r,&u); |
|
r = u; |
r = u; |
} |
} |
|
for (dc = NEXT(dc); dc; dc = NEXT(dc) ) { |
|
dctoquote(dc,v,&s,&sgn); |
|
if ( sgn == -1 ) |
|
subquote(CO,r,s,&u); |
|
else |
|
addquote(CO,r,s,&u); |
|
r = u; |
|
} |
*c = r; |
*c = r; |
} |
} |
|
|
void dctoquote(dc,v,c) |
void dptoquote(DP a,QUOTE *c) |
DCP dc; |
|
QUOTE v; |
|
QUOTE *c; |
|
{ |
{ |
QUOTE r,d,s,u; |
MP t; |
|
MP *m; |
|
int i,n,nv; |
|
QUOTE s,r,u; |
|
|
objtoquote(COEF(dc),&r); |
if ( !a ) { |
if ( DEG(dc) ) { |
MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a)); |
objtoquote(DEG(dc),&d); |
return; |
pwrquote(CO,v,d,&s); |
} |
mulquote(CO,r,s,&u); |
for ( t = BDY(a), n = 0; t; t = NEXT(t), n++ ); |
|
m = (MP *)ALLOCA(n*sizeof(MP)); |
|
for ( t = BDY(a), i = 0; t; t = NEXT(t), i++ ) |
|
m[i] = t; |
|
nv = NV(a); |
|
mptoquote(m[n-1],nv,&r); |
|
for ( i = n-2; i >= 0; i-- ) { |
|
mptoquote(m[i],nv,&s); |
|
addquote(CO,s,r,&u); |
r = u; |
r = u; |
} |
} |
*c = r; |
*c = r; |
} |
} |
|
|
void vartoquote(v,c) |
void dctoquote(DCP dc,QUOTE v,QUOTE *q,int *sgn) |
V v; |
{ |
QUOTE *c; |
QUOTE t,s,u,r; |
|
P c; |
|
Q d; |
|
|
|
if ( mmono(COEF(dc)) ) { |
|
/* -xyz... */ |
|
chsgnp(COEF(dc),&c); |
|
*sgn = -1; |
|
} else { |
|
c = COEF(dc); |
|
*sgn = 1; |
|
} |
|
d = DEG(dc); |
|
if ( UNIQ(c) ) { |
|
if ( d ) { |
|
if ( UNIQ(d) ) |
|
r = v; |
|
else { |
|
objtoquote((Obj)d,&t); |
|
pwrquote(CO,v,t,&r); |
|
} |
|
} else |
|
objtoquote((Obj)ONE,&r); |
|
} else { |
|
objtoquote((Obj)c,&u); |
|
if ( !NUM(c) && NEXT(DC(c)) ) { |
|
MKQUOTE(t,mkfnode(1,I_PAREN,BDY(u))); |
|
u = t; |
|
} |
|
if ( d ) { |
|
if ( UNIQ(d) ) |
|
s = v; |
|
else { |
|
objtoquote((Obj)d,&t); |
|
pwrquote(CO,v,t,&s); |
|
} |
|
mulquote(CO,u,s,&r); |
|
} else |
|
r = u; |
|
} |
|
*q = r; |
|
} |
|
|
|
void mptoquote(MP m,int n,QUOTE *c) |
|
{ |
|
QUOTE s,u; |
|
NODE t,t1; |
|
FNODE f; |
|
Q q; |
|
DL dl; |
|
int i; |
|
|
|
objtoquote((Obj)C(m),&s); |
|
dl = m->dl; |
|
for ( i = n-1, t = 0; i >= 0; i-- ) { |
|
STOQ(dl->d[i],q); |
|
f = mkfnode(1,I_FORMULA,q); |
|
MKNODE(t1,f,t); |
|
t = t1; |
|
} |
|
MKQUOTE(u,mkfnode(1,I_EV,t)); |
|
mulquote(CO,s,u,c); |
|
} |
|
|
|
void vartoquote(V v,QUOTE *c) |
{ |
{ |
P x; |
P x; |
PF pf; |
PF pf; |