| version 1.10, 2001/10/09 01:36:06 |
version 1.15, 2009/03/16 16:43:02 |
|
|
| * 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/parif.c,v 1.9 2001/10/03 01:47:30 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/parif.c,v 1.14 2003/02/14 22:29:07 ohara Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| #include "parse.h" |
#include "parse.h" |
| |
|
| #if PARI |
#if defined(PARI) |
| #include "genpari.h" |
#include "genpari.h" |
| |
|
| extern long prec; |
extern long prec; |
| Line 60 void patori_i(GEN,N *); |
|
| Line 60 void patori_i(GEN,N *); |
|
| void ritopa(Obj,GEN *); |
void ritopa(Obj,GEN *); |
| void ritopa_i(N,int,GEN *); |
void ritopa_i(N,int,GEN *); |
| |
|
| |
void Ptodouble(); |
| void Peval(),Psetprec(),p_pi(),p_e(),p_mul(),p_gcd(); |
void Peval(),Psetprec(),p_pi(),p_e(),p_mul(),p_gcd(); |
| void asir_cgiv(GEN); |
void asir_cgiv(GEN); |
| |
|
| Line 70 void Psetprecword(); |
|
| Line 71 void Psetprecword(); |
|
| struct ftab pari_tab[] = { |
struct ftab pari_tab[] = { |
| {"eval",Peval,-2}, |
{"eval",Peval,-2}, |
| {"setprec",Psetprec,-1}, |
{"setprec",Psetprec,-1}, |
| |
{"todouble",Ptodouble,1}, |
| #if defined(INTERVAL) || 1 |
#if defined(INTERVAL) || 1 |
| {"setprecword",Psetprecword,-1}, |
{"setprecword",Psetprecword,-1}, |
| #endif |
#endif |
| Line 94 void f(ar,rp) NODE ar; Obj *rp;\ |
|
| Line 96 void f(ar,rp) NODE ar; Obj *rp;\ |
|
| |
|
| #if defined(INTERVAL) |
#if defined(INTERVAL) |
| #define PREC_CONV pariK1 |
#define PREC_CONV pariK1 |
| #else |
#elif SIZEOF_LONG == 4 |
| #if defined(LONG_IS_32BIT) |
|
| #define PREC_CONV 0.103810253 |
#define PREC_CONV 0.103810253 |
| #endif |
#elif SIZEOF_LONG == 8 |
| #if defined(LONG_IS_64BIT) |
|
| #define PREC_CONV 0.051905126 |
#define PREC_CONV 0.051905126 |
| #endif |
#endif |
| #endif |
|
| |
|
| /* XXX : we should be more careful when we free PARI pointers. */ |
/* XXX : we should be more careful when we free PARI pointers. */ |
| |
|
|
|
| return (int)(p*PREC_CONV+3); |
return (int)(p*PREC_CONV+3); |
| } |
} |
| |
|
| |
void Ptodouble(arg,rp) |
| |
NODE arg; |
| |
Num *rp; |
| |
{ |
| |
double r,i; |
| |
Real real,imag; |
| |
Num num; |
| |
|
| |
asir_assert(ARG0(arg),O_N,"todouble"); |
| |
num = (Num)ARG0(arg); |
| |
if ( !num ) { |
| |
*rp = 0; |
| |
return; |
| |
} |
| |
switch ( NID(num) ) { |
| |
case N_R: case N_Q: case N_B: |
| |
r = ToReal(num); |
| |
MKReal(r,real); |
| |
*rp = (Num)real; |
| |
break; |
| |
case N_C: |
| |
r = ToReal(((C)num)->r); |
| |
i = ToReal(((C)num)->i); |
| |
MKReal(r,real); |
| |
MKReal(i,imag); |
| |
reimtocplx((Num)real,(Num)imag,rp); |
| |
break; |
| |
default: |
| |
*rp = num; |
| |
break; |
| |
} |
| |
} |
| |
|
| void Peval(arg,rp) |
void Peval(arg,rp) |
| NODE arg; |
NODE arg; |
| Obj *rp; |
Obj *rp; |
| Line 351 struct pariftab pariftab[] = { |
|
| Line 383 struct pariftab pariftab[] = { |
|
| {"indsort",indexsort,1}, |
{"indsort",indexsort,1}, |
| {"initalg",initalg,1}, |
{"initalg",initalg,1}, |
| {"isfund",gisfundamental,1}, |
{"isfund",gisfundamental,1}, |
| {"isprime",gisprime,1}, |
{"isprime",gisprime,2}, |
| {"ispsp",gispsp,1}, |
{"ispsp",gispsp,1}, |
| {"isqrt",racine,1}, |
{"isqrt",racine,1}, |
| {"issqfree",gissquarefree,1}, |
{"issqfree",gissquarefree,1}, |
| Line 362 struct pariftab pariftab[] = { |
|
| Line 394 struct pariftab pariftab[] = { |
|
| {"keri",keri,1}, |
{"keri",keri,1}, |
| {"kerint",kerint,1}, |
{"kerint",kerint,1}, |
| {"kerintg1",kerint1,1}, |
{"kerintg1",kerint1,1}, |
| {"kerint2",kerint2,1}, |
|
| {"length",(GEN(*)())glength,1}, |
{"length",(GEN(*)())glength,1}, |
| {"lexsort",lexsort,1}, |
{"lexsort",lexsort,1}, |
| {"lift",lift,1}, |
{"lift",lift,1}, |
| {"lindep",lindep,1}, |
{"lindep",lindep,1}, |
| {"lll",lll,1}, |
{"lll",lll,1}, |
| {"lllg1",lll1,1}, |
|
| {"lllgen",lllgen,1}, |
{"lllgen",lllgen,1}, |
| {"lllgram",lllgram,1}, |
{"lllgram",lllgram,1}, |
| {"lllgramg1",lllgram1,1}, |
|
| {"lllgramgen",lllgramgen,1}, |
{"lllgramgen",lllgramgen,1}, |
| {"lllgramint",lllgramint,1}, |
{"lllgramint",lllgramint,1}, |
| {"lllgramkerim",lllgramkerim,1}, |
{"lllgramkerim",lllgramkerim,1}, |
| Line 379 struct pariftab pariftab[] = { |
|
| Line 408 struct pariftab pariftab[] = { |
|
| {"lllint",lllint,1}, |
{"lllint",lllint,1}, |
| {"lllkerim",lllkerim,1}, |
{"lllkerim",lllkerim,1}, |
| {"lllkerimgen",lllkerimgen,1}, |
{"lllkerimgen",lllkerimgen,1}, |
| {"lllrat",lllrat,1}, |
|
| {"lngamma",glngamma,1}, |
{"lngamma",glngamma,1}, |
| {"logagm",glogagm,1}, |
{"logagm",glogagm,1}, |
| {"mat",gtomat,1}, |
{"mat",gtomat,1}, |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to parif.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / builtin