version 1.2, 2001/09/04 02:45:33 |
version 1.6, 2001/10/09 01:36:25 |
|
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/parse/quote.c,v 1.5 2001/09/05 09:01:28 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; |
|
|
*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 *c) |
V v; |
{ |
QUOTE *c; |
QUOTE d,s,u; |
|
|
|
if ( UNIQ(COEF(dc)) ) { |
|
if ( DEG(dc) ) { |
|
if ( UNIQ(DEG(dc)) ) |
|
*c = v; |
|
else { |
|
objtoquote((Obj)DEG(dc),&d); |
|
pwrquote(CO,v,d,c); |
|
} |
|
} else |
|
objtoquote((Obj)ONE,c); |
|
} else { |
|
objtoquote((Obj)COEF(dc),&u); |
|
if ( DEG(dc) ) { |
|
if ( UNIQ(DEG(dc)) ) |
|
s = v; |
|
else { |
|
objtoquote((Obj)DEG(dc),&d); |
|
pwrquote(CO,v,d,&s); |
|
} |
|
mulquote(CO,u,s,c); |
|
} else |
|
*c = u; |
|
} |
|
} |
|
|
|
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; |