在 CUDA 项目中使用 sqrtl function ？

Question

我需要使用称为“intersectbump”的设备 function 内部的一些“long double”类型变量执行计算。

这是计算：

long double q = 0.002344;
long double s = 0.0266;
long double v = 0.0744;
long double a = -q * pow(PDCx, 2);
long double b = s * PDCx - 2 * q*PCLx*PDCx - PDCz;
long double c = (-1.0)*(PCLz + q * pow(PCLx, 2) - s * PCLx + v);
long double d = b * b - 4.0 * a*c;

t = (-b + sqrtl(d))/(2.0*a);

我收到错误：

“不允许从设备函数（“intersectbump”）调用主机函数（“sqrtl”）。”

和

“标识符“sqrtl”在设备代码中未定义”

我知道 long double 在 CUDA 代码中被视为 double，但是我的计算取决于 long double 的使用。 如何在 CUDA 项目中使用 sqrtl 和 powl 等函数？