| version 1.8, 2009/03/16 16:43:02 |
version 1.9, 2009/03/27 14:42:29 |
|
|
| * 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/engine/pari.c,v 1.7 2005/10/10 15:16:38 saito Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/pari.c,v 1.8 2009/03/16 16:43:02 ohara Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| |
|
| #if defined(PARI) |
#if defined(PARI) |
| #include "base.h" |
#include "base.h" |
| #include <math.h> |
#include <math.h> |
| Line 58 void patori_i(GEN,N *); |
|
| Line 59 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 *); |
| |
|
| |
// PARI_VERSION(2,2,12) == 131596 |
| |
#if PARI_VERSION_CODE >= 131596 |
| |
#define prec precreal |
| |
#endif |
| extern long prec; |
extern long prec; |
| extern int paristack; |
extern int paristack; |
| |
|
| |
long get_pariprec() { |
| |
return prec; |
| |
} |
| |
void set_pariprec(long p) { |
| |
prec = p; |
| |
} |
| |
|
| void risa_pari_init() { |
void risa_pari_init() { |
| pari_init(paristack,2); |
pari_init(paristack,2); |
| prec = 4; |
set_pariprec(4); |
| } |
} |
| |
|
| void create_pari_variable(index) |
void create_pari_variable(index) |
|
|
| GEN a; |
GEN a; |
| { |
{ |
| return lg(a); |
return lg(a); |
| |
} |
| |
|
| |
void gpui_ri(Obj a, Obj e, Obj *c) |
| |
{ |
| |
GEN pa,pe,z; |
| |
int ltop,lbot; |
| |
|
| |
ltop = avma; ritopa(a,&pa); ritopa(e,&pe); lbot = avma; |
| |
z = gerepile(ltop,lbot,gpui(pa,pe,get_pariprec())); |
| |
patori(z,c); cgiv(z); |
| } |
} |
| |
|
| void ritopa(a,rp) |
void ritopa(a,rp) |