GP1
t1 = z2 * (GP1 + z4 * (GP3 + z4 * (GP5 + z4 * GP7)));
w = ((x - half) * ss + hln2pi) + z * (GP0 + zz * GP1 + (zz * zz) * GP2);