/*
    Copyright (C) 2019 Alessandro Languasco
*/

/*-*- compile-command: "/usr/bin/gcc -c -o precT.gp.o -O3 -Wall -fno-strict-aliasing -fomit-frame-pointer -fPIC -I\"/usr/local/include\" precT.gp.c && /usr/bin/gcc -o precT.gp.so -bundle -undefined dynamic_lookup -O3 -Wall -fno-strict-aliasing -fomit-frame-pointer -fPIC precT.gp.o "; -*-*/
#include <pari/pari.h>
/*
 GP;install("my_init_precT","v","init_precT","./precT.gp.so");
 GP;install("my_DilcherTtab","D0,G,D0,G,p","DilcherTtab","./precT.gp.so");
 GP;install("my_precT","vD0,G,D0,G,p","precT","./precT.gp.so");
 */
void my_init_precT(void);
GEN my_DilcherTtab(GEN my_x, GEN my_tabstandalone, long prec);
void my_precT(GEN my_x, GEN my_y, long prec);
/*End of prototype*/


int main() {
  long a, b;
  GEN x, y;
  /* pari_init(40000000,0);   memory size, primes precomputed*/
  pari_init(4000000000,0);  /* for cluster computation; memory size, primes precomputed*/
  
  my_init_precT();
    
  printf("x = ");
  scanf("%ld",&a);
  printf("y = ");
  scanf("%ld",&b);
  x = stoi(a);
  y = stoi(b);

  my_precT(x,y,ndec2prec(38));
    /* ndec2prec(38): converts a number of decimal digits to a number of codewords */
 
  pari_close();
  return 0;
}

/*
MAC
 /usr/bin/gcc m_precT.c precT.gp.o -g -O3 -Wall -fno-strict-aliasing -fomit-frame-pointer -fPIC -I\"/usr/local/include\" -L/usr/local/lib -lpari -o m_precT.exe

pc ufficio
/usr/bin/gcc m_precT.c precT.gp.o -g -O3 -Wall -Wno-unused-result -fno-strict-aliasing -fomit-frame-pointer -fPIC -I\"/usr/local/include\"-lgmp -L/usr/local/lib -lpari -o m_precT.exe 
 */