GP2
t2 = z4 * (GP2 + z4 * (GP4 + z4 * GP6));
w = ((x - half) * ss + hln2pi) + z * (GP0 + zz * GP1 + (zz * zz) * GP2);