csqrtf
float complex csqrtf(float complex);
w = clogf(z + csqrtf(z + 1) * csqrtf(z - 1));
z2 = csqrtf(zz);