csqrtl
long double complex csqrtl(long double complex);
w = clogl(z + csqrtl(z + 1) * csqrtl(z - 1));
z2 = csqrtl(zz);