要调试我的R代码,我试图跟踪从fft()
库到R的stats
函数(快速傅立叶变换)的调用。我可以立即从R呼叫C_fft
,fftmx()
和fft_work()
吗?
来自文件src/library/stats/R/fft.R
:
fft <- function(z, inverse=FALSE) .Call(C_fft, z, inverse)
C_fft
在文件src/library/stats/src/fourier.c
中定义:
SEXP fft(SEXP z, SEXP inverse)
该逻辑在src/library/stats/src/fft.c
中定义的fftmx()
和fft_work()
中实现:
static void fftmx(double *a, double *b, int ntot, int n, int nspan, int isn,
int m, int kt, double *at, double *ck, double *bt, double *sk,
int *np, int *nfac)
...
到目前为止,我已经尝试.Call(C_fft, z, inverse)
和dyn.load("stats.dll")
无济于事。