sub bitsize ($)
{
  return pari2num(floor(Math::Pari::log(shift)/Matth::Pari::log(2)) + 1);
}