version 1.106, 2005/12/10 14:14:15 |
version 1.108, 2005/12/11 07:21:43 |
|
|
* 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/strobj.c,v 1.105 2005/12/09 08:10:44 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/strobj.c,v 1.107 2005/12/11 05:27:30 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
Line 92 void Pquotetotex(); |
|
Line 92 void Pquotetotex(); |
|
void Pquotetotex_env(); |
void Pquotetotex_env(); |
void Pflatten_quote(); |
void Pflatten_quote(); |
|
|
void Pqt_is_integer(),Pqt_is_rational(),Pqt_is_number(); |
void Pqt_is_integer(),Pqt_is_rational(),Pqt_is_number(),Pqt_is_coef(); |
void Pqt_is_dependent(),Pqt_is_function(); |
void Pqt_is_dependent(),Pqt_is_function(),Pqt_is_var(); |
void Pqt_set_ord(); |
void Pqt_set_ord(),Pqt_set_coef(); |
void Pqt_normalize(); |
void Pqt_normalize(); |
void Pnqt_comp(); |
void Pnqt_comp(); |
void Pnqt_match(); |
void Pnqt_match(); |
Line 109 void Pnbm_hxky(), Pnbm_xky_rest(); |
|
Line 109 void Pnbm_hxky(), Pnbm_xky_rest(); |
|
void Pnbm_hv(), Pnbm_rest(); |
void Pnbm_hv(), Pnbm_rest(); |
|
|
void Pquote_to_funargs(),Pfunargs_to_quote(),Pget_function_name(); |
void Pquote_to_funargs(),Pfunargs_to_quote(),Pget_function_name(); |
void Pqt_match(),Pget_quote_id(),Pqt_match_rewrite(); |
void Pqt_match(),Pget_quote_id(); |
void Pqt_to_nary(),Pqt_to_bin(); |
void Pqt_to_nary(),Pqt_to_bin(); |
void fnode_do_assign(NODE arg); |
void fnode_do_assign(NODE arg); |
void do_assign(NODE arg); |
void do_assign(NODE arg); |
Line 127 FNODE partial_eval(FNODE), fnode_to_nary(FNODE), fnode |
|
Line 127 FNODE partial_eval(FNODE), fnode_to_nary(FNODE), fnode |
|
FNODE nfnode_add(FNODE a1,FNODE a2,int expand); |
FNODE nfnode_add(FNODE a1,FNODE a2,int expand); |
FNODE nfnode_mul(FNODE a1,FNODE a2,int expand); |
FNODE nfnode_mul(FNODE a1,FNODE a2,int expand); |
FNODE nfnode_pwr(FNODE a1,FNODE a2,int expand); |
FNODE nfnode_pwr(FNODE a1,FNODE a2,int expand); |
FNODE nfnode_mul_coef(Num c,FNODE f,int expand); |
FNODE nfnode_mul_coef(Obj c,FNODE f,int expand); |
FNODE fnode_expand_pwr(FNODE f,int n,int expand); |
FNODE fnode_expand_pwr(FNODE f,int n,int expand); |
FNODE to_narymul(FNODE f); |
FNODE to_narymul(FNODE f); |
FNODE to_naryadd(FNODE f); |
FNODE to_naryadd(FNODE f); |
FNODE fnode_node_to_nary(ARF op,NODE n); |
FNODE fnode_node_to_nary(ARF op,NODE n); |
void fnode_base_exp(FNODE f,FNODE *bp,FNODE *ep); |
void fnode_base_exp(FNODE f,FNODE *bp,FNODE *ep); |
void fnode_coef_body(FNODE f,Num *cp,FNODE *bp); |
void fnode_coef_body(FNODE f,Obj *cp,FNODE *bp); |
FNODE nfnode_match_rewrite(FNODE f,FNODE p,FNODE c,FNODE a,int mode); |
FNODE nfnode_match_rewrite(FNODE f,FNODE p,FNODE c,FNODE a,int mode); |
FNODE fnode_apply(FNODE f,FNODE (*func)(),int expand); |
FNODE fnode_apply(FNODE f,FNODE (*func)(),int expand); |
FNODE fnode_normalize(FNODE f,int expand); |
FNODE fnode_normalize(FNODE f,int expand); |
|
FNODE rewrite_fnode(FNODE f,NODE arg,int qarg); |
|
|
struct ftab str_tab[] = { |
struct ftab str_tab[] = { |
{"sprintf",Psprintf,-99999999}, |
{"sprintf",Psprintf,-99999999}, |
Line 154 struct ftab str_tab[] = { |
|
Line 155 struct ftab str_tab[] = { |
|
{"string_to_tb",Pstring_to_tb,1}, |
{"string_to_tb",Pstring_to_tb,1}, |
{"get_quote_id",Pget_quote_id,1}, |
{"get_quote_id",Pget_quote_id,1}, |
|
|
|
{"qt_is_var",Pqt_is_var,1}, |
|
{"qt_is_coef",Pqt_is_coef,1}, |
{"qt_is_number",Pqt_is_number,1}, |
{"qt_is_number",Pqt_is_number,1}, |
{"qt_is_rational",Pqt_is_rational,1}, |
{"qt_is_rational",Pqt_is_rational,1}, |
{"qt_is_integer",Pqt_is_integer,1}, |
{"qt_is_integer",Pqt_is_integer,1}, |
{"qt_is_function",Pqt_is_function,1}, |
{"qt_is_function",Pqt_is_function,1}, |
{"qt_is_dependent",Pqt_is_dependent,2}, |
{"qt_is_dependent",Pqt_is_dependent,2}, |
|
|
|
{"qt_set_coef",Pqt_set_coef,-1}, |
{"qt_set_ord",Pqt_set_ord,-1}, |
{"qt_set_ord",Pqt_set_ord,-1}, |
{"qt_normalize",Pqt_normalize,-2}, |
{"qt_normalize",Pqt_normalize,-2}, |
{"qt_match",Pqt_match,2}, |
{"qt_match",Pqt_match,2}, |
{"qt_match_rewrite",Pqt_match_rewrite,-4}, |
|
{"nqt_match_rewrite",Pnqt_match_rewrite,3}, |
{"nqt_match_rewrite",Pnqt_match_rewrite,3}, |
|
|
{"nqt_comp",Pnqt_comp,2}, |
{"nqt_comp",Pnqt_comp,2}, |
Line 607 void Pqt_to_bin(NODE arg,QUOTE *rp) |
|
Line 610 void Pqt_to_bin(NODE arg,QUOTE *rp) |
|
MKQUOTE(*rp,f); |
MKQUOTE(*rp,f); |
} |
} |
|
|
|
void Pqt_is_var(NODE arg,Q *rp) |
|
{ |
|
QUOTE q; |
|
int ret; |
|
|
|
q = (QUOTE)ARG0(arg); |
|
asir_assert(q,O_QUOTE,"qt_is_var"); |
|
ret = fnode_is_var(BDY(q)); |
|
STOQ(ret,*rp); |
|
} |
|
|
|
void Pqt_is_coef(NODE arg,Q *rp) |
|
{ |
|
QUOTE q; |
|
int ret; |
|
|
|
q = (QUOTE)ARG0(arg); |
|
asir_assert(q,O_QUOTE,"qt_is_coef"); |
|
ret = fnode_is_coef(BDY(q)); |
|
STOQ(ret,*rp); |
|
} |
|
|
void Pqt_is_number(NODE arg,Q *rp) |
void Pqt_is_number(NODE arg,Q *rp) |
{ |
{ |
QUOTE q; |
QUOTE q; |
Line 710 void Pnqt_match(NODE arg,Q *rp) |
|
Line 735 void Pnqt_match(NODE arg,Q *rp) |
|
*rp = 0; |
*rp = 0; |
} |
} |
|
|
|
|
FNODE rewrite_fnode(FNODE,NODE,int); |
|
|
|
extern Obj VOIDobj; |
|
|
|
void Pqt_match_rewrite(NODE arg,Obj *rp) |
|
{ |
|
FNODE f,g,h,c,value; |
|
Obj obj; |
|
QUOTE q; |
|
NODE r,t,s,s0,pair; |
|
int ret,ind,ac; |
|
|
|
obj = (Obj)ARG0(arg); |
|
ret = qt_match(obj,(Obj)ARG1(arg),&r); |
|
if ( ret ) { |
|
for ( t = r, s0 = 0; t; t = NEXT(t) ) { |
|
NEXTNODE(s0,s); |
|
pair = BDY((LIST)BDY(t)); |
|
ind = (int)FA0((FNODE)BDY((QUOTE)BDY(pair))); |
|
value = BDY((QUOTE)(BDY(NEXT(pair)))); |
|
BDY(s) = mknode(2,ind,value); |
|
} |
|
if ( s0 ) NEXT(s) = 0; |
|
switch ( ac = argc(arg) ) { |
|
case 3: |
|
h = rewrite_fnode(BDY((QUOTE)ARG2(arg)),s0,0); |
|
MKQUOTE(q,h); *rp = (Obj)q; |
|
break; |
|
case 4: |
|
c = rewrite_fnode(BDY((QUOTE)ARG2(arg)),s0,0); |
|
if ( eval(c) ) { |
|
h = rewrite_fnode(BDY((QUOTE)ARG3(arg)),s0,0); |
|
MKQUOTE(q,h); *rp = (Obj)q; |
|
} else |
|
*rp = VOIDobj; |
|
break; |
|
default: |
|
error("qt_match_rewrite : invalid argument"); |
|
} |
|
} else |
|
*rp = VOIDobj; |
|
} |
|
|
|
void Pnqt_match_rewrite(NODE arg,Obj *rp) |
void Pnqt_match_rewrite(NODE arg,Obj *rp) |
{ |
{ |
FNODE f,p,c,a,r; |
FNODE f,p,c,a,r; |
Line 800 FNODE nfnode_match_rewrite(FNODE f,FNODE p,FNODE c,FNO |
|
Line 781 FNODE nfnode_match_rewrite(FNODE f,FNODE p,FNODE c,FNO |
|
ARF op; |
ARF op; |
NODE arg,h0,t,h,valuen; |
NODE arg,h0,t,h,valuen; |
NODE r,s0,s,pair; |
NODE r,s0,s,pair; |
FNODE any,pany,head,tail,a1; |
FNODE any,pany,head,tail,a1,a2; |
QUOTE q; |
QUOTE q; |
int ret; |
int ret; |
FNODE value; |
FNODE value; |
Line 841 FNODE nfnode_match_rewrite(FNODE f,FNODE p,FNODE c,FNO |
|
Line 822 FNODE nfnode_match_rewrite(FNODE f,FNODE p,FNODE c,FNO |
|
ret = nfnode_match(f,p,&r) && eval(rewrite_fnode(c,r,1)); |
ret = nfnode_match(f,p,&r) && eval(rewrite_fnode(c,r,1)); |
|
|
if ( ret ) { |
if ( ret ) { |
a1 = rewrite_fnode(a,r,1); |
a1 = rewrite_fnode(a,r,0); |
return fnode_normalize(partial_eval(a1),mode); |
a2 = partial_eval(a1); |
|
return fnode_normalize(a2,mode); |
} else |
} else |
return 0; |
return 0; |
} |
} |
Line 2184 VL reordvars(VL vl0,NODE head) |
|
Line 2166 VL reordvars(VL vl0,NODE head) |
|
return vl; |
return vl; |
} |
} |
|
|
VL qt_current_ord; |
VL qt_current_ord, qt_current_coef; |
LIST qt_current_ord_obj; |
LIST qt_current_ord_obj,qt_current_coef_obj; |
|
|
void Pqt_set_ord(NODE arg,LIST *rp) |
void Pqt_set_ord(NODE arg,LIST *rp) |
{ |
{ |
Line 2206 void Pqt_set_ord(NODE arg,LIST *rp) |
|
Line 2188 void Pqt_set_ord(NODE arg,LIST *rp) |
|
} |
} |
} |
} |
|
|
|
void Pqt_set_coef(NODE arg,LIST *rp) |
|
{ |
|
NODE r0,r,n; |
|
VL vl0,vl; |
|
P v; |
|
|
|
if ( !argc(arg) ) |
|
*rp = qt_current_coef_obj; |
|
else { |
|
n = BDY((LIST)ARG0(arg)); |
|
for ( vl0 = 0, r0 = 0; n; n = NEXT(n) ) { |
|
NEXTNODE(r0,r); |
|
NEXTVL(vl0,vl); |
|
vl->v = VR((P)BDY(n)); |
|
MKV(vl->v,v); BDY(r) = v; |
|
} |
|
if ( r0 ) NEXT(r) = 0; |
|
if ( vl0 ) NEXT(vl) = 0; |
|
qt_current_coef = vl0; |
|
MKLIST(*rp,r0); |
|
qt_current_coef_obj = *rp; |
|
} |
|
} |
|
|
void Pqt_normalize(NODE arg,QUOTE *rp) |
void Pqt_normalize(NODE arg,QUOTE *rp) |
{ |
{ |
QUOTE q,r; |
QUOTE q,r; |
Line 2471 void Pnqt_comp(NODE arg,Q *rp) |
|
Line 2477 void Pnqt_comp(NODE arg,Q *rp) |
|
STOQ(r,*rp); |
STOQ(r,*rp); |
} |
} |
|
|
INLINE int fnode_is_number(FNODE f) |
int fnode_is_var(FNODE f) |
{ |
{ |
Obj obj; |
Obj obj; |
|
VL vl,t,s; |
|
DCP dc; |
|
|
|
if ( fnode_is_coef(f) ) return 0; |
switch ( f->id ) { |
switch ( f->id ) { |
|
case I_PAREN: |
|
return fnode_is_var(FA0(f)); |
|
|
|
case I_FORMULA: |
|
obj = FA0(f); |
|
if ( obj && OID(obj) == O_P ) { |
|
dc = DC((P)obj); |
|
if ( !cmpq(DEG(dc),ONE) && !NEXT(dc) |
|
&& !arf_comp(CO,(Obj)COEF(dc),(Obj)ONE) ) return 1; |
|
else return 0; |
|
} else return 0; |
|
|
|
default: |
|
return 0; |
|
} |
|
} |
|
|
|
int fnode_is_coef(FNODE f) |
|
{ |
|
Obj obj; |
|
VL vl,t,s; |
|
|
|
switch ( f->id ) { |
case I_MINUS: case I_PAREN: |
case I_MINUS: case I_PAREN: |
|
return fnode_is_coef(FA0(f)); |
|
|
|
case I_FORMULA: |
|
obj = FA0(f); |
|
if ( !obj ) return 1; |
|
else if ( OID(obj) == O_QUOTE ) |
|
return fnode_is_coef(BDY((QUOTE)obj)); |
|
else if ( NUM(obj) ) return 1; |
|
else if ( OID(obj) == O_P || OID(obj) == O_R) { |
|
get_vars_recursive(obj,&vl); |
|
for ( t = vl; t; t = NEXT(t) ) { |
|
for ( s = qt_current_coef; s; s = NEXT(s) ) |
|
if ( t->v == s->v ) break; |
|
if ( !s ) return 0; |
|
} |
|
return 1; |
|
} else return 0; |
|
|
|
case I_BOP: |
|
return fnode_is_coef(FA1(f)) && fnode_is_coef(FA2(f)); |
|
|
|
default: |
|
return 0; |
|
} |
|
} |
|
|
|
int fnode_is_number(FNODE f) |
|
{ |
|
Obj obj; |
|
|
|
switch ( f->id ) { |
|
case I_MINUS: case I_PAREN: |
return fnode_is_number(FA0(f)); |
return fnode_is_number(FA0(f)); |
|
|
case I_FORMULA: |
case I_FORMULA: |
Line 2631 FNODE fnode_normalize(FNODE f,int expand) |
|
Line 2695 FNODE fnode_normalize(FNODE f,int expand) |
|
break; |
break; |
|
|
case I_MINUS: |
case I_MINUS: |
r = nfnode_mul_coef((Num)q, |
r = nfnode_mul_coef((Obj)q, |
fnode_normalize(FA0(f),expand),expand); |
fnode_normalize(FA0(f),expand),expand); |
break; |
break; |
|
|
Line 2644 FNODE fnode_normalize(FNODE f,int expand) |
|
Line 2708 FNODE fnode_normalize(FNODE f,int expand) |
|
r = nfnode_add(a1,a2,expand); |
r = nfnode_add(a1,a2,expand); |
break; |
break; |
case '-': |
case '-': |
a2 = nfnode_mul_coef((Num)q,a2,expand); |
a2 = nfnode_mul_coef((Obj)q,a2,expand); |
r = nfnode_add(a1,a2,expand); |
r = nfnode_add(a1,a2,expand); |
break; |
break; |
case '*': |
case '*': |
Line 2731 FNODE nfnode_add(FNODE f1,FNODE f2,int expand) |
|
Line 2795 FNODE nfnode_add(FNODE f1,FNODE f2,int expand) |
|
NODE n1,n2,r0,r; |
NODE n1,n2,r0,r; |
FNODE b1,b2; |
FNODE b1,b2; |
int s; |
int s; |
Num c1,c2,c; |
Obj c1,c2,c; |
|
|
if ( IS_ZERO(f1) ) return f2; |
if ( IS_ZERO(f1) ) return f2; |
else if ( IS_ZERO(f2) ) return f1; |
else if ( IS_ZERO(f2) ) return f1; |
Line 2745 FNODE nfnode_add(FNODE f1,FNODE f2,int expand) |
|
Line 2809 FNODE nfnode_add(FNODE f1,FNODE f2,int expand) |
|
} else if ( s < 0 ) { |
} else if ( s < 0 ) { |
NEXTNODE(r0,r); BDY(r) = BDY(n2); n2 = NEXT(n2); |
NEXTNODE(r0,r); BDY(r) = BDY(n2); n2 = NEXT(n2); |
} else { |
} else { |
addnum(0,c1,c2,&c); |
arf_add(CO,c1,c2,&c); |
if ( c ) { |
if ( c ) { |
NEXTNODE(r0,r); BDY(r) = nfnode_mul_coef(c,b1,expand); |
NEXTNODE(r0,r); BDY(r) = nfnode_mul_coef(c,b1,expand); |
} |
} |
Line 2781 FNODE nfnode_mul(FNODE f1,FNODE f2,int expand) |
|
Line 2845 FNODE nfnode_mul(FNODE f1,FNODE f2,int expand) |
|
FNODE b1,b2,e1,e2,cc,t,t1; |
FNODE b1,b2,e1,e2,cc,t,t1; |
FNODE *m; |
FNODE *m; |
int s; |
int s; |
Num c1,c2,c,e; |
Obj c1,c2,c; |
|
Num e; |
int l1,l,i,j; |
int l1,l,i,j; |
|
|
if ( IS_ZERO(f1) || IS_ZERO(f2) ) return mkfnode(1,I_FORMULA,0); |
if ( IS_ZERO(f1) || IS_ZERO(f2) ) return mkfnode(1,I_FORMULA,0); |
else if ( fnode_is_number(f1) ) |
else if ( fnode_is_coef(f1) ) |
return nfnode_mul_coef((Num)eval(f1),f2,expand); |
return nfnode_mul_coef((Obj)eval(f1),f2,expand); |
else if ( fnode_is_number(f2) ) |
else if ( fnode_is_coef(f2) ) |
return nfnode_mul_coef((Num)eval(f2),f1,expand); |
return nfnode_mul_coef((Obj)eval(f2),f1,expand); |
|
|
if ( expand && IS_NARYADD(f1) ) { |
if ( expand && IS_NARYADD(f1) ) { |
t = mkfnode(1,I_FORMULA,0); |
t = mkfnode(1,I_FORMULA,0); |
Line 2808 FNODE nfnode_mul(FNODE f1,FNODE f2,int expand) |
|
Line 2873 FNODE nfnode_mul(FNODE f1,FNODE f2,int expand) |
|
} |
} |
|
|
fnode_coef_body(f1,&c1,&b1); fnode_coef_body(f2,&c2,&b2); |
fnode_coef_body(f1,&c1,&b1); fnode_coef_body(f2,&c2,&b2); |
mulnum(0,c1,c2,&c); |
arf_mul(CO,c1,c2,&c); |
if ( !c ) return mkfnode(1,I_FORMULA,0); |
if ( !c ) return mkfnode(1,I_FORMULA,0); |
|
|
|
|
Line 2842 FNODE nfnode_mul(FNODE f1,FNODE f2,int expand) |
|
Line 2907 FNODE nfnode_mul(FNODE f1,FNODE f2,int expand) |
|
FNODE nfnode_pwr(FNODE f1,FNODE f2,int expand) |
FNODE nfnode_pwr(FNODE f1,FNODE f2,int expand) |
{ |
{ |
FNODE b,b1,e1,e,cc,r,mf2,mone,inv; |
FNODE b,b1,e1,e,cc,r,mf2,mone,inv; |
Num c,c1,nf2; |
Obj c,c1; |
|
Num nf2; |
int ee; |
int ee; |
NODE arg,n,t0,t1; |
NODE arg,n,t0,t1; |
Q q; |
Q q; |
|
|
if ( IS_ZERO(f2) ) return mkfnode(1,I_FORMULA,ONE); |
if ( IS_ZERO(f2) ) return mkfnode(1,I_FORMULA,ONE); |
else if ( IS_ZERO(f1) ) return mkfnode(1,I_FORMULA,0); |
else if ( IS_ZERO(f1) ) return mkfnode(1,I_FORMULA,0); |
else if ( fnode_is_number(f1) ) { |
else if ( fnode_is_coef(f1) ) { |
if ( fnode_is_integer(f2) ) { |
if ( fnode_is_integer(f2) ) { |
if ( fnode_is_one(f2) ) return f1; |
if ( fnode_is_one(f2) ) return f1; |
else { |
else { |
pwrnum(0,(Num)eval(f1),(Num)eval(f2),&c); |
arf_pwr(CO,eval(f1),(Obj)eval(f2),&c); |
return mkfnode(1,I_FORMULA,c); |
return mkfnode(1,I_FORMULA,c); |
} |
} |
} else |
} else { |
|
f1 = mkfnode(1,I_FORMULA,eval(f1)); |
return mkfnode(3,I_BOP,pwrfs,f1,f2); |
return mkfnode(3,I_BOP,pwrfs,f1,f2); |
|
} |
} else if ( IS_BINARYPWR(f1) ) { |
} else if ( IS_BINARYPWR(f1) ) { |
b1 = FA1(f1); e1 = FA2(f1); |
b1 = FA1(f1); e1 = FA2(f1); |
e = nfnode_mul(e1,f2,expand); |
e = nfnode_mul(e1,f2,expand); |
Line 2869 FNODE nfnode_pwr(FNODE f1,FNODE f2,int expand) |
|
Line 2937 FNODE nfnode_pwr(FNODE f1,FNODE f2,int expand) |
|
&& fnode_is_integer(f2) ) { |
&& fnode_is_integer(f2) ) { |
fnode_coef_body(f1,&c1,&b1); |
fnode_coef_body(f1,&c1,&b1); |
nf2 = (Num)eval(f2); |
nf2 = (Num)eval(f2); |
pwrnum(0,(Num)c1,nf2,&c); |
arf_pwr(CO,c1,(Obj)nf2,&c); |
ee = QTOS((Q)nf2); |
ee = QTOS((Q)nf2); |
cc = mkfnode(1,I_FORMULA,c); |
cc = mkfnode(1,I_FORMULA,c); |
if ( fnode_is_nonnegative_integer(f2) ) |
if ( fnode_is_nonnegative_integer(f2) ) |
Line 2961 FNODE to_narymul(FNODE f) |
|
Line 3029 FNODE to_narymul(FNODE f) |
|
return r; |
return r; |
} |
} |
|
|
FNODE nfnode_mul_coef(Num c,FNODE f,int expand) |
FNODE nfnode_mul_coef(Obj c,FNODE f,int expand) |
{ |
{ |
FNODE b1,cc; |
FNODE b1,cc; |
Num c1,c2; |
Obj c1,c2; |
NODE n,r,r0; |
NODE n,r,r0; |
|
|
if ( !c ) |
if ( !c ) |
return mkfnode(I_FORMULA,0); |
return mkfnode(I_FORMULA,0); |
else { |
else { |
fnode_coef_body(f,&c1,&b1); |
fnode_coef_body(f,&c1,&b1); |
mulnum(0,c,c1,&c2); |
arf_mul(CO,c,c1,&c2); |
if ( UNIQ(c2) ) return b1; |
if ( UNIQ(c2) ) return b1; |
else { |
else { |
cc = mkfnode(1,I_FORMULA,c2); |
cc = mkfnode(1,I_FORMULA,c2); |
Line 2996 FNODE nfnode_mul_coef(Num c,FNODE f,int expand) |
|
Line 3064 FNODE nfnode_mul_coef(Num c,FNODE f,int expand) |
|
} |
} |
} |
} |
|
|
void fnode_coef_body(FNODE f,Num *cp,FNODE *bp) |
void fnode_coef_body(FNODE f,Obj *cp,FNODE *bp) |
{ |
{ |
FNODE c; |
FNODE c; |
|
|
if ( fnode_is_number(f) ) { |
if ( fnode_is_coef(f) ) { |
*cp = eval(f); *bp = mkfnode(1,I_FORMULA,ONE); |
*cp = (Obj)eval(f); *bp = mkfnode(1,I_FORMULA,ONE); |
} else if ( IS_NARYMUL(f) ) { |
} else if ( IS_NARYMUL(f) ) { |
c=(FNODE)BDY((NODE)FA1(f)); |
c=(FNODE)BDY((NODE)FA1(f)); |
if ( fnode_is_number(c) ) { |
if ( fnode_is_coef(c) ) { |
*cp = eval(c); |
*cp = (Obj)eval(c); |
*bp = fnode_node_to_nary(mulfs,NEXT((NODE)FA1(f))); |
*bp = fnode_node_to_nary(mulfs,NEXT((NODE)FA1(f))); |
} else { |
} else { |
*cp = (Num)ONE; *bp = f; |
*cp = (Obj)ONE; *bp = f; |
} |
} |
} else { |
} else { |
*cp = (Num)ONE; *bp = f; |
*cp = (Obj)ONE; *bp = f; |
} |
} |
} |
} |
|
|
Line 3023 int nfnode_comp(FNODE f1,FNODE f2) |
|
Line 3091 int nfnode_comp(FNODE f1,FNODE f2) |
|
int r,i1,i2,ret; |
int r,i1,i2,ret; |
char *nm1,*nm2; |
char *nm1,*nm2; |
FNODE b1,b2,e1,e2,g,a1,a2,fn1,fn2; |
FNODE b1,b2,e1,e2,g,a1,a2,fn1,fn2; |
Num ee,ee1,c1,c2; |
Num ee,ee1; |
|
Obj c1,c2; |
|
|
if ( IS_NARYADD(f1) || IS_NARYADD(f2) ) { |
if ( IS_NARYADD(f1) || IS_NARYADD(f2) ) { |
f1 = to_naryadd(f1); f2 = to_naryadd(f2); |
f1 = to_naryadd(f1); f2 = to_naryadd(f2); |
Line 3038 int nfnode_comp(FNODE f1,FNODE f2) |
|
Line 3107 int nfnode_comp(FNODE f1,FNODE f2) |
|
if ( IS_NARYMUL(f1) || IS_NARYMUL(f2) ) { |
if ( IS_NARYMUL(f1) || IS_NARYMUL(f2) ) { |
fnode_coef_body(f1,&c1,&b1); |
fnode_coef_body(f1,&c1,&b1); |
fnode_coef_body(f2,&c2,&b2); |
fnode_coef_body(f2,&c2,&b2); |
if ( !compfnode(b1,b2) ) return compnum(0,c1,c2); |
if ( !compfnode(b1,b2) ) return arf_comp(CO,c1,c2); |
b1 = to_narymul(b1); b2 = to_narymul(b2); |
b1 = to_narymul(b1); b2 = to_narymul(b2); |
n1 = (NODE)FA1(b1); n2 = (NODE)FA1(b2); |
n1 = (NODE)FA1(b1); n2 = (NODE)FA1(b2); |
while ( 1 ) { |
while ( 1 ) { |