- moved polynomial interpolation in heur_gcd() to its own function