/*
    Copyright (C) 2019 Alessandro Languasco
*/

/*-*- compile-command: "/usr/bin/gcc -c -o precpsi.gp.o -O3 -Wall -fno-strict-aliasing -fomit-frame-pointer -fPIC -I\"/usr/local/include\" precpsi.gp.c && /usr/bin/gcc -o precpsi.gp.so -bundle -undefined dynamic_lookup -O3 -Wall -fno-strict-aliasing -fomit-frame-pointer -fPIC precpsi.gp.o "; -*-*/
#include <pari/pari.h>
/*
 GP;install("my_init_precpsi","v","init_precpsi","./precpsi.gp.so");
 GP;install("my_precpsi","vD0,G,D0,G,p","precpsi","./precpsi.gp.so");
 */
void my_init_precpsi(void);
void my_precpsi(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_precpsi();
    
  printf("x = ");
  scanf("%ld",&a);
  printf("y = ");
  scanf("%ld",&b);
  x = stoi(a);
  y = stoi(b);

  my_precpsi(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_precpsi.c precpsi.gp.o -g -O3 -Wall -fno-strict-aliasing -fomit-frame-pointer -fPIC -I\"/usr/local/include\" -L/usr/local/lib -lpari -o m_precpsi.exe

pc ufficio
/usr/bin/gcc m_precpsi.c precpsi.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_precpsi.exe 

 */