| version 1.4, 2000/12/05 01:24:52 |
version 1.8, 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/engine/pari.c,v 1.3 2000/08/22 05:04:06 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/pari.c,v 1.7 2005/10/10 15:16:38 saito Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| #if PARI |
#if defined(PARI) |
| #include "base.h" |
#include "base.h" |
| #include <math.h> |
#include <math.h> |
| #include "genpari.h" |
#include "genpari.h" |
| Line 62 extern long prec; |
|
| Line 62 extern long prec; |
|
| extern int paristack; |
extern int paristack; |
| |
|
| void risa_pari_init() { |
void risa_pari_init() { |
| char buf[BUFSIZ]; |
|
| int i; |
|
| |
|
| pari_init(paristack,2); |
pari_init(paristack,2); |
| prec = 4; |
prec = 4; |
| } |
} |
|
|
| if ( index > max_varn ) { |
if ( index > max_varn ) { |
| for ( i = max_varn+1; i <= index; i++ ) { |
for ( i = max_varn+1; i <= index; i++ ) { |
| sprintf(name,"x%d",i); |
sprintf(name,"x%d",i); |
| |
#if (PARI_VERSION_CODE < 131594) |
| fetch_named_var(name,0); |
fetch_named_var(name,0); |
| |
#else |
| |
fetch_named_var(name); |
| |
#endif |
| } |
} |
| max_varn = index; |
max_varn = index; |
| } |
} |
|
|
| } |
} |
| } |
} |
| |
|
| #if defined(LONG_IS_32BIT) |
#if SIZEOF_LONG == 4 |
| void ritopa_i(a,s,rp) |
void ritopa_i(a,s,rp) |
| N a; |
N a; |
| int s; |
int s; |
|
|
| for ( j = 0, b = (unsigned int *)BD(z); j < l; j++ ) |
for ( j = 0, b = (unsigned int *)BD(z); j < l; j++ ) |
| b[l-j-1] = ((unsigned int *)g)[j+2]; |
b[l-j-1] = ((unsigned int *)g)[j+2]; |
| } |
} |
| #endif |
|
| |
|
| #if defined(LONG_IS_64BIT) |
#elif SIZEOF_LONG == 8 |
| void ritopa_i(a,s,rp) |
void ritopa_i(a,s,rp) |
| N a; |
N a; |
| int s; |
int s; |