コード例 #1
0
ファイル: test.gp.c プロジェクト: dan13ram/simple_polynomials
void
init_test(void)	  /* void */
{
  GEN p1;
  GEN p2;	  /* vec */
  p = pol_x(fetch_user_var("p"));
  d = pol_x(fetch_user_var("d"));
  f = pol_x(fetch_user_var("f"));
  /*allocatemem(800000000); */
  p = stoi(65537);
  d = gadd(genrand(powis(gen_2, 7)), powis(gen_2, 7));
  p1 = gaddgs(d, 1);
  {
    long x;
    p2 = cgetg(gtos(p1)+1, t_VEC);
    for (x = 1; gcmpsg(x, p1) <= 0; ++x)
      gel(p2, x) = genrand(p);
  }
  f = p2;
  gel(f, 1) = gen_1;
  f = gtopoly(f, -1);
  return;
}
コード例 #2
0
ファイル: anal.c プロジェクト: jkeuffer/pari
GEN gadd1e(GEN *x)
{
  *x=typ(*x)==t_INT?addis(*x,1):gaddgs(*x,1);
  return *x;
}