GP0
w_l = z * (GP0 + t1) + w;
w = ((x - half) * ss + hln2pi) + z * (GP0 + zz * GP1 + (zz * zz) * GP2);