dotnet下时间测量(续):进行纳秒级测量

还能不能得到更精确的时间呢?那就用汇编吧,通过rdtsc直接取时钟周期数。在Feng Yuan的《Windows图形编程》上找到获取时钟周期的函数,在网上搜索到获取本机CPU主频函数,凑在一起,得到如下代码:

 1 #include <windows.h>
 2 
 3 extern "C"
 4 {
 5     __declspec(dllexport) unsigned __int64 GetCycleCount(void)
 6     {
 7         _asm    _emit 0x0F;
 8         _asm    _emit 0x31;
 9     }
10 
11     __declspec(dllexport) int GetFrequency(void//MHz
12     {
13         LARGE_INTEGER CurrTicks, TicksCount;
14         __int64 iStartCounter, iStopCounter;
15 
16         DWORD dwOldProcessP = GetPriorityClass(GetCurrentProcess());
17         DWORD dwOldThreadP = GetThreadPriority(GetCurrentThread());
18 
19         SetPriorityClass(GetCurrentProcess(), REALTIME_PRIORITY_CLASS);
20         SetThreadPriority(GetCurrentThread(), THREAD_PRIORITY_TIME_CRITICAL);
21 
22         QueryPerformanceFrequency(&TicksCount);
23         QueryPerformanceCounter(&CurrTicks);
24 
25         TicksCount.QuadPart /= 16;
26         TicksCount.QuadPart += CurrTicks.QuadPart;
27 
28         _asm rdtsc
29         _asm mov DWORD PTR iStartCounter, EAX
30         _asm mov DWORD PTR (iStartCounter+4), EDX
31 
32         while(CurrTicks.QuadPart<TicksCount.QuadPart)
33             QueryPerformanceCounter(&CurrTicks);
34 
35         _asm rdtsc
36         _asm mov DWORD PTR iStopCounter, EAX
37         _asm mov DWORD PTR (iStopCounter + 4), EDX
38 
39         SetThreadPriority(GetCurrentThread(), dwOldThreadP);
40         SetPriorityClass(GetCurrentProcess(), dwOldProcessP);
41 
42         return (int)((iStopCounter-iStartCounter)/62500);
43     }
44 }

编译为TimeStamp.dll,通过DllImport调用,代码如下(为测试而作,实用得再加几个方法):

1dotnet下时间测量(续):进行纳秒级测量    public sealed class TimeStamp
2dotnet下时间测量(续):进行纳秒级测量    {
3dotnet下时间测量(续):进行纳秒级测量        [DllImport("TimeStamp.dll")]
4dotnet下时间测量(续):进行纳秒级测量        public static extern long GetCycleCount();
5dotnet下时间测量(续):进行纳秒级测量
6dotnet下时间测量(续):进行纳秒级测量        [DllImport("TimeStamp.dll")]
7dotnet下时间测量(续):进行纳秒级测量        public static extern int GetFrequency();
8dotnet下时间测量(续):进行纳秒级测量    }

先测试TimeStamp.GetCycleCount()的执行时间,也既两次GetCycleCount()运行的时间差。测试代码如下:

 1dotnet下时间测量(续):进行纳秒级测量        static void TestAbsoluteError()
 2dotnet下时间测量(续):进行纳秒级测量        {
 3dotnet下时间测量(续):进行纳秒级测量            long start,end;
 4dotnet下时间测量(续):进行纳秒级测量            int frequency = TimeStamp.GetFrequency();
 5dotnet下时间测量(续):进行纳秒级测量
 6dotnet下时间测量(续):进行纳秒级测量            using (StreamWriter sw = new StreamWriter("Results.txt")) 
 7dotnet下时间测量(续):进行纳秒级测量            {
 8dotnet下时间测量(续):进行纳秒级测量
 9dotnet下时间测量(续):进行纳秒级测量                for(int i=0;i<1000;i++)
10dotnet下时间测量(续):进行纳秒级测量                {
11dotnet下时间测量(续):进行纳秒级测量                    start=TimeStamp.GetCycleCount();
12dotnet下时间测量(续):进行纳秒级测量                    end=TimeStamp.GetCycleCount();
13dotnet下时间测量(续):进行纳秒级测量                    sw.WriteLine("{0}",(double)(end-start)/frequency);
14dotnet下时间测量(续):进行纳秒级测量                }

15dotnet下时间测量(续):进行纳秒级测量            }

16dotnet下时间测量(续):进行纳秒级测量        }


作统计图:

dotnet下时间测量(续):进行纳秒级测量

中间有几个10ms的突出点,这应该是测试程序被中断了。

为了避免被外界中断,把当前进程设置为实时进程,所得结果作图如下。

dotnet下时间测量(续):进行纳秒级测量

可以看出,i>50后,结果稳定在0.1646微秒。也就是说,必须至少运行100次TimeStamp.GetCycleCount(),然后每次TimeStamp.GetCycleCount()调用的时间才会降低到0.1646微秒,所得的测试结果才比较精确。为安全起见,建议每次测试之前运行至少500次TimeStamp.GetCycleCount()。

比较实时进程和非实时进程运行的前100次结果,得图:

dotnet下时间测量(续):进行纳秒级测量

可以看出,开始一次TimeStamp.GetCycleCount()调用运行耗时差不多都在2-3微秒,当i在50附近,非实时进程测试中这一值跳到3-4,然后迅速下降到0.164,实时进程跳到>25,然后下降到0.164。这应该和运行时优化有关。

为了得到更精确的时间,需要对测试结果进行修正,扣除一次TimeStamp.GetCycleCount()调用的时间。测量修正效果,代码如下:

 1dotnet下时间测量(续):进行纳秒级测量
 2dotnet下时间测量(续):进行纳秒级测量        static void TestRelativeError()
 3dotnet下时间测量(续):进行纳秒级测量        {
 4dotnet下时间测量(续):进行纳秒级测量            long start,end;
 5dotnet下时间测量(续):进行纳秒级测量            int frequency = TimeStamp.GetFrequency();
 6dotnet下时间测量(续):进行纳秒级测量            long check=0;
 7dotnet下时间测量(续):进行纳秒级测量
 8dotnet下时间测量(续):进行纳秒级测量            using (StreamWriter sw = new StreamWriter("Results.txt")) 
 9dotnet下时间测量(续):进行纳秒级测量            {
10dotnet下时间测量(续):进行纳秒级测量                for(int i=0;i<500;i++)
11dotnet下时间测量(续):进行纳秒级测量                {
12dotnet下时间测量(续):进行纳秒级测量                    TimeStamp.GetCycleCount();
13dotnet下时间测量(续):进行纳秒级测量                }

14dotnet下时间测量(续):进行纳秒级测量
15dotnet下时间测量(续):进行纳秒级测量                for(int i=0;i<500;i++)
16dotnet下时间测量(续):进行纳秒级测量                {
17dotnet下时间测量(续):进行纳秒级测量                    start=TimeStamp.GetCycleCount();
18dotnet下时间测量(续):进行纳秒级测量                    end=TimeStamp.GetCycleCount();
19dotnet下时间测量(续):进行纳秒级测量                    check+=end-start;
20dotnet下时间测量(续):进行纳秒级测量                }

21dotnet下时间测量(续):进行纳秒级测量                check = check/500;
22dotnet下时间测量(续):进行纳秒级测量
23dotnet下时间测量(续):进行纳秒级测量                Console.WriteLine(check);
24dotnet下时间测量(续):进行纳秒级测量
25dotnet下时间测量(续):进行纳秒级测量                for(int i=0;i<10000;i++)
26dotnet下时间测量(续):进行纳秒级测量                {
27dotnet下时间测量(续):进行纳秒级测量                    start=TimeStamp.GetCycleCount();
28dotnet下时间测量(续):进行纳秒级测量                    end=TimeStamp.GetCycleCount();
29dotnet下时间测量(续):进行纳秒级测量                    sw.WriteLine("{0}",(double)(end-start-check)/frequency);
30dotnet下时间测量(续):进行纳秒级测量                }

31dotnet下时间测量(续):进行纳秒级测量            }

32dotnet下时间测量(续):进行纳秒级测量        }


得到下图:

dotnet下时间测量(续):进行纳秒级测量

可以看出,大多数情况下,修正后时间测量误差是-0.0006微秒,也就是-0.6纳秒,一个时钟周期左右!!!!够精确了!!!!!!!!!考虑极少数异常现象,误差也在20纳秒之内,但这些异常点可以通过多次测试去除最高最低点排除在外。

总结:

通过调用rdtsc,并经过时间校正,可以得到纳米级的时间测量工具。

具体步骤如下:

(1)设置运行进程为实时进程;

(2)运行至少100次(最好500次)TimeStamp.GetCycleCount()热热身;

(3)运行500次左右TimeStamp.GetCycleCount()得到TimeStamp.GetCycleCount()调用平均周期数;

(4)进行时间测量,测量结果减去第(3)步中的调用时间

(5)重复第4步n次(n>5)既可,去掉一个最高的,去掉一个最低的,剩下的取平均值,这样测量误差应该小于1纳秒。

本文转自xiaotie博客园博客,原文链接http://www.cnblogs.com/xiaotie/archive/2005/08/17/216876.html如需转载请自行联系原作者


xiaotie 集异璧实验室(GEBLAB)