__kernel_tandf
float __kernel_tandf(double,int);
t = __kernel_tandf(M_PI * x, 1);
t = -__kernel_tandf(M_PI * (0.5 - x), -1);