1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
#define IN_INIT_PARI
#define NEED_MP_H
#ifndef STANDALONE
#include "include.h"
#endif
#ifdef GMP
/* static void */
/* init_pari(void) */
/* { */
/* ; */
/* } */
#else
GEN gnil,gzero,gun,gdeux,ghalf,gi;
plong lontyp[30]={0,0x10000,0x10000,1,1,1,1,2,1,0,2,2,1,1,1,0,1,1,1,1};
unsigned plong hiremainder,overflow;
#ifdef STANDALONE
#define FEerror printf
#define make_si_sfun(a,b,c)
#endif
#define INITIAL_PARI_STACK 400
char initial_pari_stack[400];
our_ulong bot= (our_ulong) initial_pari_stack;
our_ulong top = (our_ulong)(initial_pari_stack+INITIAL_PARI_STACK);
/* not initted */
our_ulong avma= 0;
void
err(s)
int s;
{ switch (s) {
case errpile:
FEerror("Out of bignum stack space, (si::MULTIPLY-BIGNUM-STACK n) to grow",0);
case dvmer1:
case diver4:
case diver2:
case diver1:
FEerror("Divide by zero",0);
case muler1:
FEerror("Multiply overflow",0);
case moder1:
FEerror("Mod by 0",0);
default:
FEerror("Integer Arithmetic error",0);
}}
multiply_bignum_stack(n)
int n;
{ int parisize = n* (top - bot);
in_saved_avma = 0;
if (n> 1)
{ if (bot != (our_ulong)initial_pari_stack) free(bot);
set_pari_stack(parisize);
}
return parisize;
}
set_pari_stack(parisize)
int parisize;
{
bot=(plong)malloc(parisize);
top = avma = bot + parisize;
}
static
init_pari()
{
if (avma==0)
{
make_si_sfun("MULTIPLY-BIGNUM-STACK",multiply_bignum_stack,
ARGTYPE1(f_fixnum) | RESTYPE(f_fixnum));
avma = top;
}
/* room for the permanent things */
gnil = cgeti(2);gnil[1]=2; setpere(gnil,255);
gzero = cgeti(2);gzero[1]=2; setpere(gzero, 255);
gun = stoi(1); setpere(gun, 255);
gdeux = stoi(2); setpere(gdeux, 255);
ghalf = cgetg(3,4);ghalf[1]=un;ghalf[2]=deux; setpere(ghalf, 255);
gi = cgetg(3,6); gi[1] = zero; gi[2] = un; setpere(gi, 255);
/* set_pari_stack(BIGNUM_STACK_SIZE);*/
}
#endif
|
