i would suggest to make `g` and `f`  - constant aliases  (`const poly&`)... 

Although, this may require to fix the declaration of `p_Diff` to take `const poly` input.