#include <stdio.h>
#include <math.h>
#include <time.h>
#include <conio.h>
unsigned sqrt_fpu1(long l)
{
if (l == 0)
return 0;
double f_rslt = sqrt(l);
unsigned rslt = (int)f_rslt;
if (!(f_rslt - rslt < .5))
rslt ++;
return rslt;
}
unsigned sqrt_cpu1(long l)
{
long temp;
unsigned div, rslt = l;
if (l <= 0)
return 0;
else if (l & 0xFFFF0000L)
if (l & 0xFF000000L)
div = 0x3FFF;
else
div = 0x3FF;
else
if (l & 0x0FF00L)
div = 0x3F;
else
div = (l > 4) ? 0x7 : l;
while (1)
{
temp = l / div + div;
div = temp >> 1;
div += temp & 1;
if (rslt > div)
rslt = div;
else
{
if (l / rslt == rslt - 1 && l % rslt == 0)
rslt--;
break;
}
}
return (unsigned)rslt;
}
unsigned sqrt_cpu2(long l)
{
if (l <= 0)
return 0;
long rslt = l, div = l;
while (1)
{
div = (l / div + div) / 2;
if (rslt > div)
rslt = div;
else
break;
}
return (unsigned)rslt;
}
unsigned sqrt_cpu3(long l)
{
unsigned div = 1;
unsigned rslt = 0;
while (l > 0)
{
l-= div, div += 2;
rslt += l < 0 ? 0 : 1;
}
return rslt;
}
unsigned sqrt_cpu4(long l)
{
unsigned div = 1, rslt = 0;
while (l > 0)
{
l-= div, div += 2;
rslt += l < 0 ? 0 : 1;
}
if (l != 0)
rslt++;
return rslt;
}
#define steps 1000
#define range1 1000L
#define range2 10000L
#define range3 100000L
#define count 2000
double times[3][3];
void CalcTime()
{
long l;
int i;
time_t first, second;
printf("testing range1=[%lu..%lu]...\n", 0L, range1);
printf("fpu1...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = 0l; l < range1; l += range1 / steps)
sqrt_fpu1(l);
second = time(NULL);
times[0][0] = difftime(second, first);
printf("cpu1...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = 0l; l < range1; l += range1 / steps)
sqrt_cpu1(l);
second = time(NULL);
times[0][1] = difftime(second, first);
printf("cpu2...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = 0l; l < range1; l += range1 / steps)
sqrt_cpu2(l);
second = time(NULL);
times[0][2] = difftime(second, first);
printf("testing range2=[%lu..%lu]...\n", range1, range2);
printf("fpu1...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = range1; l < range2; l += (range2 - range1) / steps)
sqrt_fpu1(l);
second = time(NULL);
times[1][0] = difftime(second, first);
printf("cpu1...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = range1; l < range2; l += (range2 - range1) / steps)
sqrt_cpu1(l);
second = time(NULL);
times[1][1] = difftime(second, first);
printf("cpu2...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = range1; l < range2; l += (range2 - range1) / steps)
sqrt_cpu2(l);
second = time(NULL);
times[1][2] = difftime(second, first);
printf("testing range3=[%lu..%lu]...\n", range2, range3);
printf("fpu1...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = range2; l < range3; l += (range3 - range2) / steps)
sqrt_fpu1(l);
second = time(NULL);
times[2][0] = difftime(second, first);
printf("cpu1...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = range2; l < range3; l += (range3 - range2) / steps)
sqrt_cpu1(l);
second = time(NULL);
times[2][1] = difftime(second, first);
printf("cpu2...\n");
first = time(NULL);
for (i = 0; i < count; i++)
for (l = range2; l < range3; l += (range3 - range2) / steps)
sqrt_cpu2(l);
second = time(NULL);
times[2][2] = difftime(second, first);
printf("Done.\n");
printf("range\t\t fpu1\t cpu1\t cpu2\n");
printf(
"%lu\t\t%5.3f\t%5.3f\t%5.3f\n",
range1,
times[0][0],
times[0][1],
times[0][2]
);
printf(
"%lu\t\t%5.3f\t%5.3f\t%5.3f\n",
range2,
times[1][0],
times[1][1],
times[1][2]
);
printf(
"%lu\t\t%5.3f\t%5.3f\t%5.3f\n",
range3,
times[2][0],
times[2][1],
times[2][2]
);
}
typedef unsigned (*sqrt_func)(long L);
void ViewDifferents(long rang1, long rang2, unsigned step, sqrt_func fpsqrt)
{
unsigned long l;
unsigned rf, ri;
printf("testing with range[%lu..%lu]\n", rang1, rang2);
for (l = rang1; l < rang2; l += (rang2 - rang1) / step)
{
rf = sqrt_fpu1(l);
ri = (*fpsqrt)(l);
if (rf != ri)
printf("sqrt(%lu) %u %u %6.2f\n", l, rf, ri, sqrt(l));
}
printf("Done.\n");
}
void main()
{
ViewDifferents(0, 1000, 1000, sqrt_cpu1);
CalcTime();
while (!kbhit());
}
1 2
8 8 8
| 
*
Обратная связь