__k_j1f
f = (float)__k_j1f(fx);
return (x * (u / v) + tpi * (__k_j1f(fx) * log(x) - one / x));
b = __k_j1f(fx);