version 1.15, 2017/08/31 04:21:48 |
version 1.17, 2018/03/29 01:32:50 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/builtin/bfaux.c,v 1.14 2017/03/29 01:15:14 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/builtin/bfaux.c,v 1.16 2018/03/28 05:27:22 noro Exp $ */ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
|
|
Line 14 void Prk_ratmat(); |
|
Line 14 void Prk_ratmat(); |
|
void mp_sin(),mp_cos(),mp_tan(),mp_asin(),mp_acos(),mp_atan(); |
void mp_sin(),mp_cos(),mp_tan(),mp_asin(),mp_acos(),mp_atan(); |
void mp_sinh(),mp_cosh(),mp_tanh(),mp_asinh(),mp_acosh(),mp_atanh(); |
void mp_sinh(),mp_cosh(),mp_tanh(),mp_asinh(),mp_acosh(),mp_atanh(); |
void mp_exp(),mp_log(),mp_pow(); |
void mp_exp(),mp_log(),mp_pow(); |
|
void mp_factorial(),mp_abs(); |
|
|
struct ftab bf_tab[] = { |
struct ftab bf_tab[] = { |
{"eval",Peval,-2}, |
{"eval",Peval,-2}, |
{"setprec",Psetprec,-1}, |
{"setprec",Psetprec,-1}, |
{"setbprec",Psetbprec,-1}, |
{"setbprec",Psetbprec,-1}, |
{"setround",Psetround,-1}, |
{"setround",Psetround,-1}, |
{"todouble",Ptodouble,1}, |
{"todouble",Ptodouble,1}, |
{"mpfr_sin",mp_sin,-2}, |
{"mpfr_sin",mp_sin,-2}, |
{"mpfr_cos",mp_cos,-2}, |
{"mpfr_cos",mp_cos,-2}, |
{"mpfr_tan",mp_tan,-2}, |
{"mpfr_tan",mp_tan,-2}, |
{"mpfr_asin",mp_asin,-2}, |
{"mpfr_asin",mp_asin,-2}, |
{"mpfr_acos",mp_acos,-2}, |
{"mpfr_acos",mp_acos,-2}, |
{"mpfr_atan",mp_atan,-2}, |
{"mpfr_atan",mp_atan,-2}, |
{"mpfr_sinh",mp_sinh,-2}, |
{"mpfr_sinh",mp_sinh,-2}, |
{"mpfr_cosh",mp_cosh,-2}, |
{"mpfr_cosh",mp_cosh,-2}, |
{"mpfr_tanh",mp_tanh,-2}, |
{"mpfr_tanh",mp_tanh,-2}, |
{"mpfr_asinh",mp_asinh,-2}, |
{"mpfr_asinh",mp_asinh,-2}, |
{"mpfr_acosh",mp_acosh,-2}, |
{"mpfr_acosh",mp_acosh,-2}, |
{"mpfr_atanh",mp_atanh,-2}, |
{"mpfr_atanh",mp_atanh,-2}, |
{"mpfr_exp",mp_exp,-2}, |
{"mpfr_exp",mp_exp,-2}, |
{"mpfr_log",mp_log,-2}, |
{"mpfr_log",mp_log,-2}, |
{"mpfr_pow",mp_pow,-3}, |
{"mpfr_pow",mp_pow,-3}, |
{"mpfr_ai",Pmpfr_ai,-2}, |
{"mpfr_ai",Pmpfr_ai,-2}, |
{"mpfr_zeta",Pmpfr_zeta,-2}, |
{"mpfr_zeta",Pmpfr_zeta,-2}, |
{"mpfr_j0",Pmpfr_j0,-2}, |
{"mpfr_j0",Pmpfr_j0,-2}, |
{"mpfr_j1",Pmpfr_j1,-2}, |
{"mpfr_j1",Pmpfr_j1,-2}, |
{"mpfr_y0",Pmpfr_y0,-2}, |
{"mpfr_y0",Pmpfr_y0,-2}, |
{"mpfr_y1",Pmpfr_y1,-2}, |
{"mpfr_y1",Pmpfr_y1,-2}, |
{"mpfr_eint",Pmpfr_eint,-2}, |
{"mpfr_eint",Pmpfr_eint,-2}, |
{"mpfr_erf",Pmpfr_erf,-2}, |
{"mpfr_erf",Pmpfr_erf,-2}, |
{"mpfr_erfc",Pmpfr_erfc,-2}, |
{"mpfr_erfc",Pmpfr_erfc,-2}, |
{"mpfr_li2",Pmpfr_li2,-2}, |
{"mpfr_li2",Pmpfr_li2,-2}, |
{"mpfr_gamma",Pmpfr_gamma,-2}, |
{"mpfr_gamma",Pmpfr_gamma,-2}, |
{"mpfr_lngamma",Pmpfr_gamma,-2}, |
{"mpfr_lngamma",Pmpfr_gamma,-2}, |
{"mpfr_digamma",Pmpfr_gamma,-2}, |
{"mpfr_digamma",Pmpfr_gamma,-2}, |
{"mpfr_floor",Pmpfr_floor,-2}, |
{"mpfr_floor",Pmpfr_floor,-2}, |
{"mpfr_ceil",Pmpfr_ceil,-2}, |
{"mpfr_ceil",Pmpfr_ceil,-2}, |
{"mpfr_round",Pmpfr_round,-2}, |
{"mpfr_round",Pmpfr_round,-2}, |
{"rk_ratmat",Prk_ratmat,7}, |
{"rk_ratmat",Prk_ratmat,7}, |
{0,0,0}, |
{0,0,0}, |
}; |
}; |
|
|
int mpfr_roundmode = MPFR_RNDN; |
int mpfr_roundmode = MPFR_RNDN; |
|
|
void todoublen(Num a,Num *rp) |
void todoublen(Num a,Num *rp) |
{ |
{ |
double r,i; |
double r,i; |
Real real,imag; |
Real real,imag; |
|
|
if ( !a ) { |
if ( !a ) { |
*rp = 0; |
*rp = 0; |
return; |
return; |
} |
} |
switch ( NID(a) ) { |
switch ( NID(a) ) { |
case N_R: case N_Q: case N_B: |
case N_R: case N_Q: case N_B: |
r = ToReal(a); |
r = ToReal(a); |
MKReal(r,real); |
MKReal(r,real); |
*rp = (Num)real; |
*rp = (Num)real; |
break; |
break; |
case N_C: |
case N_C: |
r = ToReal(((C)a)->r); |
r = ToReal(((C)a)->r); |
i = ToReal(((C)a)->i); |
i = ToReal(((C)a)->i); |
MKReal(r,real); |
MKReal(r,real); |
MKReal(i,imag); |
MKReal(i,imag); |
reimtocplx((Num)real,(Num)imag,rp); |
reimtocplx((Num)real,(Num)imag,rp); |
break; |
break; |
default: |
default: |
*rp = a; |
*rp = a; |
break; |
break; |
} |
} |
} |
} |
|
|
void todoublep(P a,P *rp) |
void todoublep(P a,P *rp) |
Line 119 void todoubler(R a,R *rp) |
|
Line 120 void todoubler(R a,R *rp) |
|
|
|
void todouble(Obj a,Obj *b) |
void todouble(Obj a,Obj *b) |
{ |
{ |
Obj t; |
Obj t; |
LIST l; |
LIST l; |
V v; |
V v; |
int row,col,len; |
int row,col,len; |
VECT vect; |
VECT vect; |
MAT mat; |
MAT mat; |
int i,j; |
int i,j; |
NODE n0,n,nd; |
NODE n0,n,nd; |
MP m,mp,mp0; |
MP m,mp,mp0; |
DP d; |
DP d; |
|
|
if ( !a ) { |
if ( !a ) { |
*b = 0; |
*b = 0; |
return; |
return; |
} |
} |
switch ( OID(a) ) { |
switch ( OID(a) ) { |
case O_N: |
case O_N: |
todoublen((Num)a,(Num *)b); |
todoublen((Num)a,(Num *)b); |
break; |
break; |
case O_P: |
case O_P: |
todoublep((P)a,(P *)b); |
todoublep((P)a,(P *)b); |
break; |
break; |
case O_R: |
case O_R: |
todoubler((R)a,(R *)b); |
todoubler((R)a,(R *)b); |
break; |
break; |
case O_LIST: |
case O_LIST: |
n0 = 0; |
n0 = 0; |
for ( nd = BDY((LIST)a); nd; nd = NEXT(nd) ) { |
for ( nd = BDY((LIST)a); nd; nd = NEXT(nd) ) { |
NEXTNODE(n0,n); |
NEXTNODE(n0,n); |
todouble((Obj)BDY(nd),(Obj *)&BDY(n)); |
todouble((Obj)BDY(nd),(Obj *)&BDY(n)); |
} |
} |
if ( n0 ) |
if ( n0 ) |
NEXT(n) = 0; |
NEXT(n) = 0; |
MKLIST(l,n0); |
MKLIST(l,n0); |
*b = (Obj)l; |
*b = (Obj)l; |
break; |
break; |
case O_VECT: |
case O_VECT: |
len = ((VECT)a)->len; |
len = ((VECT)a)->len; |
MKVECT(vect,len); |
MKVECT(vect,len); |
for ( i = 0; i < len; i++ ) { |
for ( i = 0; i < len; i++ ) { |
todouble((Obj)BDY((VECT)a)[i],(Obj *)&BDY(vect)[i]); |
todouble((Obj)BDY((VECT)a)[i],(Obj *)&BDY(vect)[i]); |
} |
} |
*b = (Obj)vect; |
*b = (Obj)vect; |
break; |
break; |
case O_MAT: |
case O_MAT: |
row = ((MAT)a)->row; |
row = ((MAT)a)->row; |
col = ((MAT)a)->col; |
col = ((MAT)a)->col; |
MKMAT(mat,row,col); |
MKMAT(mat,row,col); |
for ( i = 0; i < row; i++ ) |
for ( i = 0; i < row; i++ ) |
for ( j = 0; j < col; j++ ) { |
for ( j = 0; j < col; j++ ) { |
todouble((Obj)BDY((MAT)a)[i][j],(Obj *)&BDY(mat)[i][j]); |
todouble((Obj)BDY((MAT)a)[i][j],(Obj *)&BDY(mat)[i][j]); |
} |
} |
*b = (Obj)mat; |
*b = (Obj)mat; |
break; |
break; |
case O_DP: |
case O_DP: |
mp0 = 0; |
mp0 = 0; |
for ( m = BDY((DP)a); m; m = NEXT(m) ) { |
for ( m = BDY((DP)a); m; m = NEXT(m) ) { |
todouble(C(m),&t); |
todouble(C(m),&t); |
if ( t ) { |
if ( t ) { |
NEXTMP(mp0,mp); |
NEXTMP(mp0,mp); |
C(mp) = t; |
C(mp) = t; |
mp->dl = m->dl; |
mp->dl = m->dl; |
} |
} |
} |
} |
if ( mp0 ) { |
if ( mp0 ) { |
MKDP(NV((DP)a),mp0,d); |
MKDP(NV((DP)a),mp0,d); |
d->sugar = ((DP)a)->sugar; |
d->sugar = ((DP)a)->sugar; |
*b = (Obj)d; |
*b = (Obj)d; |
} else |
} else |
*b = 0; |
*b = 0; |
|
|
break; |
break; |
default: |
default: |
error("todouble : invalid argument"); |
error("todouble : invalid argument"); |
} |
} |
} |
} |
|
|
void Ptodouble(NODE arg,Obj *rp) |
void Ptodouble(NODE arg,Obj *rp) |
Line 205 void Peval(NODE arg,Obj *rp) |
|
Line 206 void Peval(NODE arg,Obj *rp) |
|
{ |
{ |
int prec; |
int prec; |
|
|
asir_assert(ARG0(arg),O_R,"eval"); |
asir_assert(ARG0(arg),O_R,"eval"); |
if ( argc(arg) == 2 ) { |
if ( argc(arg) == 2 ) { |
prec = QTOS((Q)ARG1(arg))*3.32193; |
prec = QTOS((Q)ARG1(arg))*3.32193; |
if ( prec < MPFR_PREC_MIN ) prec = MPFR_PREC_MIN; |
if ( prec < MPFR_PREC_MIN ) prec = MPFR_PREC_MIN; |
else if ( prec > MPFR_PREC_MAX ) prec = MPFR_PREC_MAX; |
else if ( prec > MPFR_PREC_MAX ) prec = MPFR_PREC_MAX; |
} else |
} else |
prec = 0; |
prec = 0; |
evalr(CO,(Obj)ARG0(arg),prec,rp); |
evalr(CO,(Obj)ARG0(arg),prec,rp); |
} |
} |
|
|
/* set/get decimal precision */ |
/* set/get decimal precision */ |
|
|
void Psetprec(NODE arg,Obj *rp) |
void Psetprec(NODE arg,Obj *rp) |
{ |
{ |
int p; |
int p; |
Q q; |
Q q; |
int prec,dprec; |
int prec,dprec; |
|
|
prec = mpfr_get_default_prec(); |
prec = mpfr_get_default_prec(); |
/* decimal precision */ |
/* decimal precision */ |
dprec = prec*0.30103; |
dprec = prec*0.30103; |
STOQ(dprec,q); *rp = (Obj)q; |
STOQ(dprec,q); *rp = (Obj)q; |
if ( arg ) { |
if ( arg ) { |
asir_assert(ARG0(arg),O_N,"setprec"); |
asir_assert(ARG0(arg),O_N,"setprec"); |
p = QTOS((Q)ARG0(arg))*3.32193; |
p = QTOS((Q)ARG0(arg))*3.32193; |
if ( p > 0 ) |
if ( p > 0 ) |
prec = p; |
prec = p; |
} |
} |
if ( prec < MPFR_PREC_MIN ) prec = MPFR_PREC_MIN; |
if ( prec < MPFR_PREC_MIN ) prec = MPFR_PREC_MIN; |
else if ( prec > MPFR_PREC_MAX ) prec = MPFR_PREC_MAX; |
else if ( prec > MPFR_PREC_MAX ) prec = MPFR_PREC_MAX; |
mpfr_set_default_prec(prec); |
mpfr_set_default_prec(prec); |
} |
} |
|
|
/* set/get bit precision */ |
/* set/get bit precision */ |
|
|
void Psetbprec(NODE arg,Obj *rp) |
void Psetbprec(NODE arg,Obj *rp) |
{ |
{ |
int p; |
int p; |
Q q; |
Q q; |
int prec; |
int prec; |
|
|
prec = mpfr_get_default_prec(); |
prec = mpfr_get_default_prec(); |
STOQ(prec,q); *rp = (Obj)q; |
STOQ(prec,q); *rp = (Obj)q; |
if ( arg ) { |
if ( arg ) { |
asir_assert(ARG0(arg),O_N,"setbprec"); |
asir_assert(ARG0(arg),O_N,"setbprec"); |
p = QTOS((Q)ARG0(arg)); |
p = QTOS((Q)ARG0(arg)); |
if ( p > 0 ) |
if ( p > 0 ) |
prec = p; |
prec = p; |
} |
} |
if ( prec < MPFR_PREC_MIN ) prec = MPFR_PREC_MIN; |
if ( prec < MPFR_PREC_MIN ) prec = MPFR_PREC_MIN; |
else if ( prec > MPFR_PREC_MAX ) prec = MPFR_PREC_MAX; |
else if ( prec > MPFR_PREC_MAX ) prec = MPFR_PREC_MAX; |
mpfr_set_default_prec(prec); |
mpfr_set_default_prec(prec); |
} |
} |
|
|
void Psetround(NODE arg,Q *rp) |
void Psetround(NODE arg,Q *rp) |
Line 265 void Psetround(NODE arg,Q *rp) |
|
Line 266 void Psetround(NODE arg,Q *rp) |
|
|
|
STOQ(mpfr_roundmode,*rp); |
STOQ(mpfr_roundmode,*rp); |
if ( arg ) { |
if ( arg ) { |
asir_assert(ARG0(arg),O_N,"setround"); |
asir_assert(ARG0(arg),O_N,"setround"); |
round = QTOS((Q)ARG0(arg)); |
round = QTOS((Q)ARG0(arg)); |
switch ( round ) { |
switch ( round ) { |
case 0: |
case 0: |
Line 301 Num tobf(Num a,int prec); |
|
Line 302 Num tobf(Num a,int prec); |
|
void mp_pi(NODE arg,BF *rp) |
void mp_pi(NODE arg,BF *rp) |
{ |
{ |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
prec = arg ? QTOS((Q)ARG0(arg)) : 0; |
prec = arg ? QTOS((Q)ARG0(arg)) : 0; |
NEWBF(r); |
NEWBF(r); |
prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body); |
prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body); |
mpfr_const_pi(r->body,mpfr_roundmode); |
mpfr_const_pi(r->body,mpfr_roundmode); |
if ( !cmpbf((Num)r,0) ) r = 0; |
if ( !cmpbf((Num)r,0) ) r = 0; |
*rp = r; |
*rp = r; |
} |
} |
Line 314 void mp_pi(NODE arg,BF *rp) |
|
Line 315 void mp_pi(NODE arg,BF *rp) |
|
void mp_e(NODE arg,BF *rp) |
void mp_e(NODE arg,BF *rp) |
{ |
{ |
int prec; |
int prec; |
mpfr_t one; |
mpfr_t one; |
BF r; |
BF r; |
|
|
prec = arg ? QTOS((Q)ARG0(arg)) : 0; |
prec = arg ? QTOS((Q)ARG0(arg)) : 0; |
NEWBF(r); |
NEWBF(r); |
prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body); |
prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body); |
mpfr_init(one); |
mpfr_init(one); |
mpfr_set_ui(one,1,mpfr_roundmode); |
mpfr_set_ui(one,1,mpfr_roundmode); |
mpfr_exp(r->body,one,mpfr_roundmode); |
mpfr_exp(r->body,one,mpfr_roundmode); |
if ( !cmpbf((Num)r,0) ) r = 0; |
if ( !cmpbf((Num)r,0) ) r = 0; |
*rp = r; |
*rp = r; |
} |
} |
|
|
void mpfr_or_mpc(NODE arg,int (*mpfr_f)(),int (*mpc_f)(),Num *rp) |
void mpfr_or_mpc(NODE arg,int (*mpfr_f)(),int (*mpc_f)(),Num *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r,re,im; |
BF r,re,im; |
C c; |
C c; |
mpc_t mpc,a1; |
mpc_t mpc,a1; |
|
|
prec = NEXT(arg) ? QTOS((Q)ARG1(arg)) : mpfr_get_default_prec(); |
prec = NEXT(arg) ? QTOS((Q)ARG1(arg)) : mpfr_get_default_prec(); |
a = tobf(ARG0(arg),prec); |
a = tobf(ARG0(arg),prec); |
if ( a && NID(a)==N_C ) { |
if ( a && NID(a)==N_C ) { |
mpc_init2(mpc,prec); mpc_init2(a1,prec); |
mpc_init2(mpc,prec); mpc_init2(a1,prec); |
re = (BF)((C)a)->r; im = (BF)((C)a)->i; |
re = (BF)((C)a)->r; im = (BF)((C)a)->i; |
Line 353 void mpfr_or_mpc(NODE arg,int (*mpfr_f)(),int (*mpc_f) |
|
Line 354 void mpfr_or_mpc(NODE arg,int (*mpfr_f)(),int (*mpc_f) |
|
*rp = (Num)c; |
*rp = (Num)c; |
} |
} |
} else { |
} else { |
NEWBF(r); |
NEWBF(r); |
mpfr_init2(r->body,prec); |
mpfr_init2(r->body,prec); |
(*mpfr_f)(r->body,((BF)a)->body,mpfr_roundmode); |
(*mpfr_f)(r->body,((BF)a)->body,mpfr_roundmode); |
if ( !cmpbf((Num)r,0) ) r = 0; |
if ( !cmpbf((Num)r,0) ) r = 0; |
*rp = (Num)r; |
*rp = (Num)r; |
} |
} |
Line 431 void mp_log(NODE arg,Num *rp) |
|
Line 432 void mp_log(NODE arg,Num *rp) |
|
mpfr_or_mpc(arg,mpfr_log,mpc_log,rp); |
mpfr_or_mpc(arg,mpfr_log,mpc_log,rp); |
} |
} |
|
|
|
void mp_abs(NODE arg,Num *rp) |
|
{ |
|
mpfr_or_mpc(arg,mpfr_abs,mpc_abs,rp); |
|
} |
|
|
|
void mp_factorial(NODE arg,Num *rp) |
|
{ |
|
struct oNODE arg0; |
|
Num a,a1; |
|
|
|
a = (Num)ARG0(arg); |
|
if ( !a ) *rp = (Num)ONE; |
|
else if ( INT(a) ) Pfac(arg,rp); |
|
else { |
|
addnum(0,a,(Num)ONE,&a1); |
|
arg0.body = (pointer)a1; |
|
arg0.next = arg->next; |
|
Pmpfr_gamma(&arg0,rp); |
|
} |
|
} |
|
|
void mp_pow(NODE arg,Num *rp) |
void mp_pow(NODE arg,Num *rp) |
{ |
{ |
Num a,e; |
Num a,e; |
int prec; |
int prec; |
BF r,re,im; |
BF r,re,im; |
C c; |
C c; |
mpc_t mpc,a1,e1; |
mpc_t mpc,a1,e1; |
|
|
prec = NEXT(NEXT(arg)) ? QTOS((Q)ARG2(arg)) : mpfr_get_default_prec(); |
prec = NEXT(NEXT(arg)) ? QTOS((Q)ARG2(arg)) : mpfr_get_default_prec(); |
a = tobf(ARG0(arg),prec); |
a = tobf(ARG0(arg),prec); |
e = tobf(ARG1(arg),prec); |
e = tobf(ARG1(arg),prec); |
if ( NID(a) == N_C || NID(e) == N_C || MPFR_SIGN(((BF)a)->body) < 0 ) { |
if ( NID(a) == N_C || NID(e) == N_C || MPFR_SIGN(((BF)a)->body) < 0 ) { |
mpc_init2(mpc,prec); mpc_init2(a1,prec); mpc_init2(e1,prec); |
mpc_init2(mpc,prec); mpc_init2(a1,prec); mpc_init2(e1,prec); |
if ( NID(a) == N_C ) { |
if ( NID(a) == N_C ) { |
Line 470 void mp_pow(NODE arg,Num *rp) |
|
Line 492 void mp_pow(NODE arg,Num *rp) |
|
*rp = (Num)c; |
*rp = (Num)c; |
} |
} |
} else { |
} else { |
NEWBF(r); |
NEWBF(r); |
mpfr_init2(r->body,prec); |
mpfr_init2(r->body,prec); |
mpfr_pow(r->body,((BF)a)->body,((BF)e)->body,mpfr_roundmode); |
mpfr_pow(r->body,((BF)a)->body,((BF)e)->body,mpfr_roundmode); |
*rp = (Num)r; |
*rp = (Num)r; |
} |
} |
} |
} |
Line 487 void mp_pow(NODE arg,Num *rp) |
|
Line 509 void mp_pow(NODE arg,Num *rp) |
|
|
|
void Pmpfr_gamma(NODE arg,BF *rp) |
void Pmpfr_gamma(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_gamma(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_gamma(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_lngamma(NODE arg,BF *rp) |
void Pmpfr_lngamma(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_lngamma(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_lngamma(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_digamma(NODE arg,BF *rp) |
void Pmpfr_digamma(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_digamma(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_digamma(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_zeta(NODE arg,BF *rp) |
void Pmpfr_zeta(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_zeta(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_zeta(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_eint(NODE arg,BF *rp) |
void Pmpfr_eint(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_eint(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_eint(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_erf(NODE arg,BF *rp) |
void Pmpfr_erf(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_erf(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_erf(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_erfc(NODE arg,BF *rp) |
void Pmpfr_erfc(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_erfc(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_erfc(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_j0(NODE arg,BF *rp) |
void Pmpfr_j0(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_j0(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_j0(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_j1(NODE arg,BF *rp) |
void Pmpfr_j1(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_j1(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_j1(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_y0(NODE arg,BF *rp) |
void Pmpfr_y0(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_y0(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_y0(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_y1(NODE arg,BF *rp) |
void Pmpfr_y1(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_y1(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_y1(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_li2(NODE arg,BF *rp) |
void Pmpfr_li2(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_li2(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_li2(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_ai(NODE arg,BF *rp) |
void Pmpfr_ai(NODE arg,BF *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
|
|
SETPREC |
SETPREC |
mpfr_ai(r->body,((BF)a)->body,mpfr_roundmode); |
mpfr_ai(r->body,((BF)a)->body,mpfr_roundmode); |
*rp = r; |
*rp = r; |
} |
} |
|
|
void Pmpfr_floor(NODE arg,Q *rp) |
void Pmpfr_floor(NODE arg,Q *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
mpz_t t; |
mpz_t t; |
GZ rz; |
GZ rz; |
|
|
SETPREC |
SETPREC |
mpfr_floor(r->body,((BF)a)->body); |
mpfr_floor(r->body,((BF)a)->body); |
mpz_init(t); |
mpz_init(t); |
mpfr_get_z(t,r->body,mpfr_roundmode); |
mpfr_get_z(t,r->body,mpfr_roundmode); |
MPZTOGZ(t,rz); |
MPZTOGZ(t,rz); |
*rp = gztoz(rz); |
*rp = gztoz(rz); |
} |
} |
|
|
void Pmpfr_ceil(NODE arg,Q *rp) |
void Pmpfr_ceil(NODE arg,Q *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
mpz_t t; |
mpz_t t; |
GZ rz; |
GZ rz; |
|
|
SETPREC |
SETPREC |
mpfr_ceil(r->body,((BF)a)->body); |
mpfr_ceil(r->body,((BF)a)->body); |
mpz_init(t); |
mpz_init(t); |
mpfr_get_z(t,r->body,mpfr_roundmode); |
mpfr_get_z(t,r->body,mpfr_roundmode); |
MPZTOGZ(t,rz); |
MPZTOGZ(t,rz); |
*rp = gztoz(rz); |
*rp = gztoz(rz); |
} |
} |
|
|
void Pmpfr_round(NODE arg,Q *rp) |
void Pmpfr_round(NODE arg,Q *rp) |
{ |
{ |
Num a; |
Num a; |
int prec; |
int prec; |
BF r; |
BF r; |
mpz_t t; |
mpz_t t; |
GZ rz; |
GZ rz; |
|
|
SETPREC |
SETPREC |
mpfr_round(r->body,((BF)a)->body); |
mpfr_round(r->body,((BF)a)->body); |
mpz_init(t); |
mpz_init(t); |
mpfr_get_z(t,r->body,mpfr_roundmode); |
mpfr_get_z(t,r->body,mpfr_roundmode); |
MPZTOGZ(t,rz); |
MPZTOGZ(t,rz); |
*rp = gztoz(rz); |
*rp = gztoz(rz); |
} |
} |
|
|
double **almat_double(int n) |
double **almat_double(int n) |