p + q: Mean +- std dev: [main] 35.7 us +- 1.3 us -> [pr] 31.8 us +- 0.6 us: 1.12x faster
p * q: Mean +- std dev: [main] 37.3 us +- 1.1 us -> [pr] 33.9 us +- 1.1 us: 1.10x faster

Geometric mean: 1.11x faster

import pyperf
from numpy.polynomial import Polynomial

setup = """
from numpy.polynomial import Polynomial

p=Polynomial([.5,4,5])
q=Polynomial([1,0,5])
"""

runner = pyperf.Runner()
runner.timeit(name=f"p + q", stmt="_ = p + q", setup=setup)
runner.timeit(name=f"p * q", stmt="_ = p * q", setup=setup)