version 1.8, 2005/10/05 08:57:25 |
version 1.12, 2015/08/14 13:51:56 |
|
|
* 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/parse/puref.c,v 1.7 2004/08/18 00:50:37 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/parse/puref.c,v 1.11 2015/08/08 14:19:42 fujimoto Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
Line 284 void derivr(VL vl,Obj a,V v,Obj *b) |
|
Line 284 void derivr(VL vl,Obj a,V v,Obj *b) |
|
} |
} |
} |
} |
|
|
|
void simple_derivr(VL vl,Obj a,V v,Obj *b) |
|
{ |
|
Obj r,s,t,u,nm,dn; |
|
|
|
if ( !a || NUM(a) ) |
|
*b = 0; |
|
else |
|
switch ( OID(a) ) { |
|
case O_P: |
|
pderivr(vl,a,v,b); break; |
|
case O_R: |
|
nm = (Obj)NM((R)a); dn = (Obj)DN((R)a); |
|
/* (nm/dn)' = nm'/dn - dn'/dn*nm/dn */ |
|
pderivr(vl,nm,v,&s); divr(vl,s,dn,&u); reductr(vl,u,&t); |
|
pderivr(vl,dn,v,&s); divr(vl,s,dn,&u); reductr(vl,u,&s); mulr(vl,s,a,&u); |
|
subr(vl,t,u,&s); reductr(vl,s,b); |
|
break; |
|
default: |
|
error("simple_derivr : invalid argument"); |
|
} |
|
} |
|
|
int obj_is_dependent(Obj a,V v) |
int obj_is_dependent(Obj a,V v) |
{ |
{ |
if ( !a || OID(a) <= O_N ) return 0; |
if ( !a || OID(a) <= O_N ) return 0; |
Line 628 void devalins(PFINS ins,Obj *rp) |
|
Line 650 void devalins(PFINS ins,Obj *rp) |
|
} |
} |
MKReal(d,r); *rp = (Obj)r; |
MKReal(d,r); *rp = (Obj)r; |
} |
} |
|
} |
|
|
|
extern int evalef; |
|
|
|
void simplify_elemfunc_ins(PFINS ins,Obj *rp) |
|
{ |
|
V v; |
|
P t; |
|
|
|
if ( evalef ) |
|
evalins(ins,0,rp); |
|
else { |
|
instov(ins,&v); MKV(v,t); *rp = (Obj)t; |
|
} |
} |
} |
|
|
void simplify_ins(PFINS ins,Obj *rp) |
void simplify_ins(PFINS ins,Obj *rp) |