logint(2^99,2)
logint(2^100,2,&z)
z
logint(5^100,2,&z)
z
logint(5^100,3)
logint(5^100,4)
logint(5^100,5)
logint(2^64-1, 2^32)
logint(2^64-1, 2^16)
test(x,y)=
{
  my(e = logint(x, y), z);
  if (logint(x, y, &z) != e, error([x,y]));
  [e, z];
}
x=Pi^1000;
test(x, 2)
test(x, 3)
x=quadgen(5)^100;
test(x, 2)
x=(3/2)^10;
test(x,2)

\\ ERRORS
logint(-Pi, 2)
logint(-quadgen(5), 2)
logint(-1/2, 2)
