fixmean.cpp

这是我本人编写的一个32位定点小数运算的函数库。对于没有浮点运算器的场合

// FixMean.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"

#include <stdlib.h>
#include <math.h>
#include <stdio.h>
#include <conio.h>

long FixMean32(long* x,int n);
long FixMean32_C(long* x,int n);

#define N 1151

int main(int argc, char* argv[])
{
	double sum=0,maxd=0;
	for (int i=0;i<10000;i++) {
		long x[N];
		double y=0;
		for (int j=0;j<N;j++) {
			x[j]=(rand()<<17)+(rand()<<2)+(rand()&0x3);
			y+=(double)x[j]/(1<<20);
		}
		long x1=FixMean32(x,N),xx=FixMean32_C(x,N);
		double yy=y/N;

		double d=fabs(yy-(double)xx/(1<<20));
		if (d>maxd) maxd=d;
		sum+=d;
	}
	printf("Max Delta = %16.12lf\n",maxd);
	printf("Aver Delta = %16.12lf\n",sum/10000);
	getch();

	return 0;
}