如何从R源访问hidden_​​attribute函数？

Question

我想在R中使用以下功能。

double attribute_hidden chebyshev_eval(double x, const double *a, const int n)
{
    double b0, b1, b2, twox;
    int i;

    if (n < 1 || n > 1000) ML_ERR_return_NAN;

    if (x < -1.1 || x > 1.1) ML_ERR_return_NAN;

    twox = x * 2;
    b2 = b1 = 0;
    b0 = 0;
    for (i = 1; i <= n; i++) {
b2 = b1;
b1 = b0;
b0 = twox * b1 - b2 + a[n - i];
    }
    return (b0 - b2) * 0.5;
}

作为参考，这是切比雪夫多项式的R源代码的一部分。

我似乎无法找到一种方法来触及R中的hidden_attribute函数。我尝试手动为真实用例进行计算，但我遇到了错误：

gamcs=c(    +.8571195590989331421920062399942e-2,
+.4415381324841006757191315771652e-2,
+.5685043681599363378632664588789e-1,
-.4219835396418560501012500186624e-2,
+.1326808181212460220584006796352e-2,
-.1893024529798880432523947023886e-3,
+.3606925327441245256578082217225e-4,
-.6056761904460864218485548290365e-5,
+.1055829546302283344731823509093e-5,
-.1811967365542384048291855891166e-6,
+.3117724964715322277790254593169e-7,
-.5354219639019687140874081024347e-8,
+.9193275519859588946887786825940e-9,
-.1577941280288339761767423273953e-9,
+.2707980622934954543266540433089e-10,
-.4646818653825730144081661058933e-11,
+.7973350192007419656460767175359e-12,
-.1368078209830916025799499172309e-12,
+.2347319486563800657233471771688e-13,
-.4027432614949066932766570534699e-14,
+.6910051747372100912138336975257e-15,
-.1185584500221992907052387126192e-15,
+.2034148542496373955201026051932e-16,
-.3490054341717405849274012949108e-17,
+.5987993856485305567135051066026e-18,
-.1027378057872228074490069778431e-18,
+.1762702816060529824942759660748e-19,
-.3024320653735306260958772112042e-20,
+.5188914660218397839717833550506e-21,
-.8902770842456576692449251601066e-22,
+.1527474068493342602274596891306e-22,
-.2620731256187362900257328332799e-23,
+.4496464047830538670331046570666e-24,
-.7714712731336877911703901525333e-25,
+.1323635453126044036486572714666e-25,
-.2270999412942928816702313813333e-26,
+.3896418998003991449320816639999e-27,
-.6685198115125953327792127999999e-28,
+.1146998663140024384347613866666e-28,
-.1967938586345134677295103999999e-29,
+.3376448816585338090334890666666e-30,
-.5793070335782135784625493333333e-31)
ngam = 22
y = .8
x = y*2-1
b1 = b2 = 0
b0 = 0
twox = x * 2

n = ngam
a = gamcs

for (i in 1:n) {
b2 = b1
b1 = b0
b0 = twox * b1 - b2 + a[n - i]
}

(b0 - b2) * 0.5

（取自https://github.com/wch/r-source/blob/0b493ddbefbc3cc06f675841c4672b3197dbfbbb/src/nmath/gamma.c#L143和https://github.com/wch/r-source/blob/0b493ddbefbc3cc06f675841c4672b3197dbfbbb/src/nmath/chebyshev.c#L69。y的值是一个虚拟值，其属性为abs(y * 2 - 1)<=1.1}

手动计算不仅单调乏味，而且还会导致b0的零值，因此最终结果(b0 - b2) * 0.5不会返回值。

解决方法是直接访问该功能。这是一个问题，因为无法从R访问attribute_hidden函数。

我想知道是否有显示这些计算的功能，变通方法，命令等。例如。 ShowCalculations(chebyshev_eval(.6,gamcs,ngam))。

优选地，该函数还显示/复制每一行的计算，或for - 循环的每次迭代等。

注意：我清楚但不是非常了解Rprintf选项。我对C代码的了解是有限的，所以如果有人决定提出这个选项，我们将非常感谢一些指导。