mp2dfr.c

来自「适合大型数值计算代码 现在网络上已经找不到了 购买需要20$」· C语言 代码 · 共 18 行

#define IAZ 48

void mp2dfr(unsigned char a[], unsigned char s[], int n, int *m)
{
	void mplsh(unsigned char u[], int n);
	void mpsmu(unsigned char w[], unsigned char u[], int n, int iv);
	int j;

	*m=(int) (2.408*n);
	for (j=1;j<=(*m);j++) {
		mpsmu(a,a,n,10);
		s[j]=a[1]+IAZ;
		mplsh(a,n);
	}
}
#undef IAZ