OpenXM_contrib2/asir2000/builtin/pf.c - diff

Return to pf.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / builtin

Diff for /OpenXM_contrib2/asir2000/builtin/pf.c between version 1.18 and 1.23

version 1.18, 2013/06/14 04:47:17

version 1.23, 2018/03/27 06:29:19

Line 45

* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,

* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.

* $OpenXM: OpenXM_contrib2/asir2000/builtin/pf.c,v 1.17 2011/09/14 06:41:20 noro Exp $

* $OpenXM: OpenXM_contrib2/asir2000/builtin/pf.c,v 1.22 2015/08/14 13:51:54 fujimoto Exp $

#include "ca.h"

#include "math.h"

Line 90 struct ftab puref_tab[] = {

{0,0,0},

};

#if defined(PARI)

int mp_pi(),mp_e();

int p_pi(),p_e();

int mp_exp(), mp_log(), mp_pow();

int p_log(),p_exp(),p_pow();

int mp_sin(),mp_cos(),mp_tan(),mp_asin(),mp_acos(),mp_atan();

int p_sin(),p_cos(),p_tan(),p_asin(),p_acos(),p_atan();

int mp_sinh(),mp_cosh(),mp_tanh(),mp_asinh(),mp_acosh(),mp_atanh();

int p_sinh(),p_cosh(),p_tanh(),p_asinh(),p_acosh(),p_atanh();

#else

int p_pi,p_e;

int p_log,p_exp,p_pow;

int p_sin,p_cos,p_tan,p_asin,p_acos,p_atan;

int p_sinh,p_cosh,p_tanh,p_asinh,p_acosh,p_atanh;

#endif

static V *uarg,*darg;

static P x,y;

static PF pidef,edef;

Line 116 static PF asinhdef,acoshdef,atanhdef;

Line 110 static PF asinhdef,acoshdef,atanhdef;

double const_pi() { return 3.14159265358979323846264338327950288; }

double const_e() { return 2.718281828459045235360287471352662497; }

int simplify_elemfunc_ins();

void pf_init() {

uarg = (V *)CALLOC(1,sizeof(V));

uarg[0] = &oVAR[26]; MKV(uarg[0],x);

Line 124 void pf_init() {

Line 120 void pf_init() {

darg[0] = &oVAR[26];

darg[1] = &oVAR[27]; MKV(darg[1],y);

mkpf("@pi",0,0,0,(int (*)())p_pi,const_pi,0,&pidef);

mkpf("@pi",0,0,0,(int (*)())mp_pi,const_pi,simplify_elemfunc_ins,&pidef);

mkpf("@e",0,0,0,(int (*)())p_e,const_e,0,&edef);

mkpf("@e",0,0,0,(int (*)())mp_e,const_e,simplify_elemfunc_ins,&edef);

mkpf("log",0,1,uarg,(int (*)())p_log,log,0,&logdef);

mkpf("log",0,1,uarg,(int (*)())mp_log,log,simplify_elemfunc_ins,&logdef);

mkpf("exp",0,1,uarg,(int (*)())p_exp,exp,0,&expdef);

mkpf("exp",0,1,uarg,(int (*)())mp_exp,exp,simplify_elemfunc_ins,&expdef);

mkpf("pow",0,2,darg,(int (*)())p_pow,pow,(int (*)())simplify_pow,&powdef);

mkpf("pow",0,2,darg,(int (*)())mp_pow,pow,(int (*)())simplify_pow,&powdef);

mkpf("sin",0,1,uarg,(int (*)())p_sin,sin,0,&sindef);

mkpf("sin",0,1,uarg,(int (*)())mp_sin,sin,simplify_elemfunc_ins,&sindef);

mkpf("cos",0,1,uarg,(int (*)())p_cos,cos,0,&cosdef);

mkpf("cos",0,1,uarg,(int (*)())mp_cos,cos,simplify_elemfunc_ins,&cosdef);

mkpf("tan",0,1,uarg,(int (*)())p_tan,tan,0,&tandef);

mkpf("tan",0,1,uarg,(int (*)())mp_tan,tan,simplify_elemfunc_ins,&tandef);

mkpf("asin",0,1,uarg,(int (*)())p_asin,asin,0,&asindef);

mkpf("asin",0,1,uarg,(int (*)())mp_asin,asin,simplify_elemfunc_ins,&asindef);

mkpf("acos",0,1,uarg,(int (*)())p_acos,acos,0,&acosdef);

mkpf("acos",0,1,uarg,(int (*)())mp_acos,acos,simplify_elemfunc_ins,&acosdef);

mkpf("atan",0,1,uarg,(int (*)())p_atan,atan,0,&atandef);

mkpf("atan",0,1,uarg,(int (*)())mp_atan,atan,simplify_elemfunc_ins,&atandef);

mkpf("sinh",0,1,uarg,(int (*)())p_sinh,sinh,0,&sinhdef);

mkpf("sinh",0,1,uarg,(int (*)())mp_sinh,sinh,simplify_elemfunc_ins,&sinhdef);

mkpf("cosh",0,1,uarg,(int (*)())p_cosh,cosh,0,&coshdef);

mkpf("cosh",0,1,uarg,(int (*)())mp_cosh,cosh,simplify_elemfunc_ins,&coshdef);

mkpf("tanh",0,1,uarg,(int (*)())p_tanh,tanh,0,&tanhdef);

mkpf("tanh",0,1,uarg,(int (*)())mp_tanh,tanh,simplify_elemfunc_ins,&tanhdef);

#if !defined(VISUAL)

#if !defined(VISUAL) && !defined(__MINGW32__)

mkpf("asinh",0,1,uarg,(int (*)())p_asinh,asinh,0,&asinhdef);

mkpf("asinh",0,1,uarg,(int (*)())mp_asinh,asinh,simplify_elemfunc_ins,&asinhdef);

mkpf("acosh",0,1,uarg,(int (*)())p_acosh,acosh,0,&acoshdef);

mkpf("acosh",0,1,uarg,(int (*)())mp_acosh,acosh,simplify_elemfunc_ins,&acoshdef);

mkpf("atanh",0,1,uarg,(int (*)())p_atanh,atanh,0,&atanhdef);

mkpf("atanh",0,1,uarg,(int (*)())mp_atanh,atanh,simplify_elemfunc_ins,&atanhdef);

#endif

make_exp();

make_tri();

make_itri();

make_hyp();

#if !defined(VISUAL)

#if !defined(VISUAL) && !defined(__MINGW32__)

make_ihyp();

#endif

}

Line 280 Obj *r;

Line 276 Obj *r;

simplify_ins(ins,r);

}

extern int evalef;

void simplify_pow(ins,rp)

PFINS ins;

Obj *rp;

Line 290 Obj *rp;

Line 288 Obj *rp;

V v;

P t;

if ( evalef ) {

simplify_elemfunc_ins(ins,rp);

return;

}

pf = ins->pf; ad = ins->ad; a0 = ad[0].arg; a1 = ad[1].arg;

if ( !a1 )

*rp = (Obj)ONE;

Line 299 Obj *rp;

Line 301 Obj *rp;

else if ( RATN(a1) && SGN((Q)a1) < 0 )

error("simplify_pow : division by 0");

else {

instov(ins,&v); MKV(v,t); *rp = (Obj)t;

instoobj(ins,rp);

}

} else if ( NUM(a1) && INT(a1) )

arf_pwr(CO,a0,a1,rp);

else {

instov(ins,&v); MKV(v,t); *rp = (Obj)t;

instoobj(ins,rp);

}

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>