数据结构上课发现
昨天开始看数据结构,姥姥的课上讲了一个测量循环的时间的方法。将传统算法和秦九韶算法对比,其中源代码为
#include<stdio.h>
#include<time.h>
#include<math.h>
clock_t start,stop;
double duration;
#define MAXN 10
double f1(int n,double a[],double x)
{
int i;
double p =a[0];
for(i=1;i<=n;i++)
p+=(a[i]*pow(x,i));
return p;
}
double f2(int n,double a[],double x)
{
int i;
double p=a[n];
for(i=n;i>0;i--)
p=a[i-1]+x*p;
return p;
}
int main()
{
int i,j;
double a[MAXN];
for(i=0;i<MAXN;i++)
a[i]=(double)i;
start = clock();
for(j=0;j<10000000;j++)
f1(MAXN-1,a,1.1);
stop = clock();
duration=((double)(stop-start))/CLK_TCK;
printf("ticks1=%f\n",(double)(stop-start));
printf("duration1=%3.2e\n",duration);
start = clock();
for(j=0;j<10000000;j++)
f2(MAXN-1,a,1.1);
stop = clock();
duration=((double)(stop-start))/CLK_TCK;
printf("ticks2=%f\n",(double)(stop-start));
printf("duration2=%3.2e\n",duration);
return 0;
}
然后姥姥说,mian()中有两段代码极度相似,可以改写,我想试试,不过遇到问题。
#include<stdio.h>
#include<time.h>
#include<math.h>
clock_t start,stop;
double duration;
#define MAXN 10
double f1(int n,double a[],double x)
{int i;
double p=0;
for(i=0;i<=n;i++)
p+=(a[i]*pow(x,i));
return p;
}
double f2(int n,double a[],double x)
{
int i;
double p=a[n];
for(i=n;i>0;i--)
p=(p*x+a[i-1]);
return p;
}
double f3(double (*a)(int,double[],double))
{start=clock();
(*a);
stop=clock();
duration=((double)(stop-start))/CLK_TCK;
printf("ticks=%f\n",(double)(stop-start));
return duration;
}
int main()
{
int i,j;
double (*p)(int,double[],double);
double a[MAXN];
for(i=0;i<MAXN;i++)
a[i]=(double)i;
p=f1;
for(j=0;j<10000000;j++)
(*p)(MAXN-1,a,1.1);
f3(p);
printf("duration1=%f\n",f3);
p=f2;
for(j=0;j<10000000;j++)
(*p)(MAXN-1,a,1.1);
f3(p);
printf("duration2=%f\n",f3);
}
可以看到这里用了一个函数,代替重复语句,然后函数中用函数指针调用了f1,f2,不过当使用这个函数测试后,发现两次测速时间几乎一致。
仔细思考原来是我把指针调用的时间测出来了,实际计算运行在外面已经运行好了。但是由于指针的原因,函数不可能在f3内部运行,注定不能在f3内部测速,如果把测速部分拿出来,把stop和strat拿出来作为变量,函数过于繁琐,而且节约的步骤很少,不值得这样做。
如有大佬有其他方法,请别吝啬批评。