void
help(const char *s)
{
entree *ep = is_entry(s);
if (ep && ep->help)
pari_printf("%s\n",ep->help);
else
pari_printf("Function %s not found\n",s);
}
int
main()
{
GEN x, y, d, u, v;
pari_init(1000000,2);
printf("x = "); x = gp_read_stream(stdin);
printf("y = "); y = gp_read_stream(stdin);
d = extgcd(x, y, &u, &v);
pari_printf("gcd = %Ps\nu = %Ps\nv = %Ps\n", d, u, v);
pari_close();
return 0;
}
int
main(void)
{
GEN M,N1,N2, F1,F2,D;
struct pari_thread pth[MAXTHREADS];
int numth = omp_get_max_threads(), i;
/* Initialise the main PARI stack and global objects (gen_0, etc.) */
pari_init(4000000,500000);
if (numth > MAXTHREADS)
{
numth = MAXTHREADS;
omp_set_num_threads(numth);
}
/* Compute in the main PARI stack */
N1 = addis(int2n(256), 1); /* 2^256 + 1 */
N2 = subis(int2n(193), 1); /* 2^193 - 1 */
M = mathilbert(80);
/*Allocate pari thread structures */
for (i = 1; i < numth; i++) pari_thread_alloc(&pth[i],4000000,NULL);
#pragma omp parallel
{
int this_th = omp_get_thread_num();
if (this_th) (void)pari_thread_start(&pth[this_th]);
#pragma omp sections
{
#pragma omp section
{
F1 = factor(N1);
}
#pragma omp section
{
F2 = factor(N2);
}
#pragma omp section
{
D = det(M);
}
} /* omp sections */
if (this_th) pari_thread_close();
} /* omp parallel */
pari_printf("F1=%Ps\nF2=%Ps\nlog(D)=%Ps\n", F1, F2, glog(D,3));
for (i = 1; i < numth; i++) pari_thread_free(&pth[i]);
return 0;
}
int do_factor(GEN n, long prec)
{
pari_sp ltop;
GEN sq = gfloor(gsqrt(n, prec));
GEN q = stoi(2);
ltop = avma;
for (;;) {
if (cmpii(q, sq) > 0)
return -1;
if (equalsi(0, gmod(n, q))) {
pari_printf("%Ps = %Ps * %Ps\n", n, q, gdiv(n, q));
return 0;
}
gaddz(gen_1, q, q);
avma = ltop;
}
}
void
print_all_user_fun(int member)
{
pari_sp av = avma;
long iL = 0, lL = 1024;
GEN L = cgetg(lL+1, t_VECSMALL);
entree *ep;
int i;
for (i = 0; i < functions_tblsz; i++)
for (ep = functions_hash[i]; ep; ep = ep->next)
{
const char *f;
int is_member;
if (EpVALENCE(ep) != EpVAR || typ((GEN)ep->value)!=t_CLOSURE) continue;
f = ep->name;
is_member = (f[0] == '_' && f[1] == '.');
if (member != is_member) continue;
if (iL >= lL)
{
GEN oL = L;
long j;
lL *= 2; L = cgetg(lL+1, t_VECSMALL);
for (j = 1; j <= iL; j++) gel(L,j) = gel(oL,j);
}
L[++iL] = (long)ep;
}
if (iL)
{
setlg(L, iL+1);
L = gen_sort(L, NULL, &cmp_epname);
for (i = 1; i <= iL; i++)
{
ep = (entree*)L[i];
pari_printf("%s =\n %Ps\n\n", ep->name, ep->value);
}
}
avma = av;
}
void
pariplot(GEN a, GEN b, GEN code, GEN ysmlu,GEN ybigu, long prec)
{
const char BLANK = ' ', YY = '|', XX_UPPER = '\'', XX_LOWER = '.';
long jz, j, i, sig;
pari_sp av = avma;
int jnew, jpre = 0; /* for lint */
GEN x, dx;
double diff, dyj, ysml, ybig, y[ISCR+1];
screen scr;
char buf[80], z;
sig=gcmp(b,a); if (!sig) return;
if (sig<0) { x=a; a=b; b=x; }
x = gtofp(a, prec); push_lex(x, code);
dx = divru(gtofp(gsub(b,a),prec), ISCR-1);
for (j=1; j<=JSCR; j++) scr[1][j]=scr[ISCR][j]=YY;
for (i=2; i<ISCR; i++)
{
scr[i][1] = XX_LOWER;
scr[i][JSCR]= XX_UPPER;
for (j=2; j<JSCR; j++) scr[i][j] = BLANK;
}
ysml = ybig = 0.; /* -Wall */
for (i=1; i<=ISCR; i++)
{
pari_sp av2 = avma;
y[i] = gtodouble( READ_EXPR(code,x) );
avma = av2;
if (i == 1)
ysml = ybig = y[1];
else
{
if (y[i] < ysml) ysml = y[i];
if (y[i] > ybig) ybig = y[i];
}
x = addrr(x,dx);
}
avma = av;
if (ysmlu) ysml = gtodouble(ysmlu);
if (ybigu) ybig = gtodouble(ybigu);
diff = ybig - ysml;
if (!diff) { ybig += 1; diff= 1.; }
dyj = ((JSCR-1)*3+2) / diff;
/* work around bug in gcc-4.8 (32bit): plot(x=-5,5,sin(x)))) */
jz = 3 - (long)(ysml*dyj + 0.5); /* 3 - DTOL(ysml*dyj) */
z = PICTZERO(jz); jz /= 3;
for (i=1; i<=ISCR; i++)
{
if (0<=jz && jz<=JSCR) scr[i][jz]=z;
j = 3 + DTOL((y[i]-ysml)*dyj);
jnew = j/3;
if (i > 1) fill_gap(scr, i, jnew, jpre);
if (0<=jnew && jnew<=JSCR) scr[i][jnew] = PICT(j);
jpre = jnew;
}
pari_putc('\n');
pari_printf("%s ", dsprintf9(ybig, buf));
for (i=1; i<=ISCR; i++) pari_putc(scr[i][JSCR]);
pari_putc('\n');
for (j=(JSCR-1); j>=2; j--)
{
pari_puts(" ");
for (i=1; i<=ISCR; i++) pari_putc(scr[i][j]);
pari_putc('\n');
}
pari_printf("%s ", dsprintf9(ysml, buf));
for (i=1; i<=ISCR; i++) pari_putc(scr[i][1]);
pari_putc('\n');
{
char line[10 + 32 + 32 + ISCR - 9];
sprintf(line, "%10s%-9.7g%*.7g\n"," ",todbl(a),ISCR-9,todbl(b));
pari_printf(line);
}
pop_lex(1);
}
static void
help(const char *s0, int flag)
{
const long long_help = flag & h_LONG;
long n;
entree *ep;
char *s = get_sep(s0);
if (isdigit((int)*s)) { digit_help(s,flag); return; }
if (flag & h_APROPOS) { external_help(s,-1); return; }
/* Get meaningful answer on '\ps 5' (e.g. from <F1>) */
if (*s == '\\') { char *t = s+1; pari_skip_alpha(&t); *t = '\0'; }
if (isalpha((int)*s))
{
if (!strncmp(s, "default", 7))
{ /* special-case ?default(dft_name), e.g. default(log) */
char *t = s+7;
pari_skip_space(&t);
if (*t == '(')
{
t++; pari_skip_space(&t);
cut_trailing_garbage(t);
if (pari_is_default(t)) { default_help(t,flag); return; }
}
}
cut_trailing_garbage(s);
}
if (long_help && (n = ok_external_help(&s))) { external_help(s,n); return; }
switch (*s)
{
case '*' : commands(-1); return;
case '\0': menu_commands(); return;
case '\\': slash_commands(); return;
case '.' : member_commands(); return;
}
ep = is_entry(s);
if (!ep)
{
if (pari_is_default(s))
default_help(s,flag);
else if (long_help)
external_help(s,3);
else if (!cb_pari_whatnow || !cb_pari_whatnow(pariOut, s,1))
simple_help(s,"unknown identifier");
return;
}
if (EpVALENCE(ep) == EpALIAS)
{
pari_printf("%s is aliased to:\n\n",s);
ep = do_alias(ep);
}
switch(EpVALENCE(ep))
{
case EpVAR:
if (!ep->help)
{
if (typ((GEN)ep->value)!=t_CLOSURE)
simple_help(s, "user defined variable");
else
{
GEN str = closure_get_text((GEN)ep->value);
if (typ(str) == t_VEC)
pari_printf("%s =\n %Ps\n", ep->name, ep->value);
}
return;
}
break;
case EpINSTALL:
if (!ep->help) { simple_help(s, "installed function"); return; }
break;
case EpNEW:
if (!ep->help) { simple_help(s, "new identifier"); return; };
break;
default: /* built-in function */
if (!ep->help) pari_err_BUG("gp_help (no help found)"); /*paranoia*/
if (long_help) { external_help(ep->name,3); return; }
}
print_text(ep->help);
}
static void
simple_help(const char *s1, const char *s2) { pari_printf("%s: %s\n", s1, s2); }
