Fold phase factors into coefficients throughout. There seem to be precision issues in these routines