wt1.c

来自「< C语言数值算法程序大全>>配套程序」· C语言 代码 · 共 13 行

void wt1(float a[], unsigned long n, int isign,
	void (*wtstep)(float [], unsigned long, int))
{
	unsigned long nn;

	if (n < 4) return;
	if (isign >= 0) {
		for (nn=n;nn>=4;nn>>=1) (*wtstep)(a,nn,isign);
	} else {
		for (nn=4;nn<=n;nn<<=1) (*wtstep)(a,nn,isign);
	}
}