let pgcd u v =
    let r =
      if is_zero v
      then u
      else gcd_big_int u v in
      r