#include <complex.h> // for csqrtf, cacoshf, clogf, complex

float complex cacoshf(float complex z)
{
	return clogf(z + csqrtf(z + 1) * csqrtf(z - 1));
}