在C/C++算法設計中使用任意位寬

開發定點(fixed-point)算法時,通常需要在設計功能性、數字精度建模、及驗證(仿真)速度之間取得一個平衡。現在,一種新的數據類可使此過程簡單化,由此得到更簡單精確的建模精度、更好的數字求精、及更快的驗證周期,而ANSI C/C++正是開發這種數字求精算法的最佳語言。

某此算法天生就適用于操作整數,或那些理想中的實數(如數字濾波器的系數),它們也可能會使用浮點或定點類型。一般而言,在算法開發的早期階段,會經常用到C語言的float或double浮點類型,因爲它們可提供一個非常大的動態數據範圍,且對大多數程序來說都是適用的。見圖1:

在C/C++算法設計中使用任意位寬

使用C內置的float類型來建模一個FIR濾波器

算法可進行數字求精,以便使用定點算術來降低最終硬件或軟件實現的複雜性。在硬件方面,將整數或定點算術限制爲最小位寬,可在本質上滿足性能、空間、能耗的需要;假如實現中用到了DSP處理器,那麽把算法限制爲整數或定點算術,就可爲特定程序使用盡可能便宜的處理器。

定點算術的建模可通過C語言內置的浮點或整數類型來完成,這做的話,需要顯式編碼並受限于C中浮點數及整數可表示的最大數:64位整數或53位尾數;這些都會給操作數的位寬帶來更多的限制,例如,2個33位的數相乘,會超過64位C整數可表示的範圍。圖2演示了一個FIR濾波器的例子,但temp變量限制爲15位的定點精度,其中10位用于整數位。在這個實現中,LSB的右部位被舍棄(量化模型的截斷),而MSB的左部位也被舍棄(包裝的溢出模型),應該意識到,使用float(或double)的模型在精度上是受限的,且不能再次合成(synthesis)。同樣,由于有取整模型的嚴格位精度定義有先,又由于內置浮點類型的取整將會先被應用,所以對除法這樣的操作來說,就非常難實現了。

在C/C++算法設計中使用任意位寬

使用float建模定點行爲

當許多算法都能依靠本地C數據類型的精度來編寫時,對支持任意長度的整數及定點算法,大家就會抱有極大的期望,而硬件描述語言(HDL)如VHDL,走的也是同一條路。隨著C/C++越來越多地被用于高級合成與驗證工具(High-Level Synthesis and Verification tools),也證實了這種語言本質上有一個足以滿足當前及未來程序需要的數據類型庫。任意長度類型的支持,也可使數據類型的行爲有一個統一的定義,而統一的語義則避免了人工實現上的一些限制。

算法C數據類型

算法C數據類型是一種基于類的C++庫,其實現了任意長度的整數及定點類型,而這些可自由訪問的類型有一系列好處,包括統一及良好定義的語義,還有媲美C/C++內置數據類型的運行時速度,對比SystemC中相應的類型,其運行速度也超過10倍以上。這些數據類型能用于任何符合C++或SystemC規範標准的程序中,並擁有高度可合成的語義。

語義

語義的統一性與一致性是避免在算法中,發生功能性錯誤的要害,以下的例子,也說明了這點:

在C/C++算法設計中使用任意位寬

衆所周知,變量ActLength的範圍爲1至255,萬一編譯器的合成不知道其範圍,就不能進行相應的優化,它的聲明就會從int變爲更嚴格的sc_uint<8>類型;雖然合成會得到更好的結果,但設計就仿真得不正確了。在經過一番調試之後,找到了問題的源頭:在比較表達式k >= ActLength中,兩個操作數變成了一個signed int與一個unsigned long long(爲64位無符號整數,其是sc_uint類型的基類型)之間的比較。對此的解釋是:C/C++整數提升規則指定了在進行比較之前,會把操作數int提升爲一個unsigned long long,例如,假如k的值爲 -1,在提升爲unsigned long long之後,它會變成2^64 – 1。

像這樣語義中的問題一般會非常難以察覺,且是與位寬相關的,例如,可能有人想擴大某個現有算法的位寬,只有看到結果時,才知道是行不通的;這個問題也可能是與特定平台相關的,例如,對1 << 32(兩個操作數都是int類型,結果也是int類型)大家期望返回0,但在大多數平台(或編譯器)上,它都會返回1(沒有移位,只因爲第二個操作數較低的五位被計算進來了);當第一個操作數是一個64位整數時,平台依靠性會表現得更加明顯及頻繁。主要的問題是C/C++標准沒有指定在32位整數情況下移位值(第二個操作數)超出0至31範圍、或在64位整數情況下移位值超出0至63範圍時的行爲。不幸的是,像sc_int、sc_uint這樣的數據類型也不能爲用戶避免這類平台依靠性的問題。

算法C數據類型被設計用于提供統一且一致的語義,因此,它們是可猜測的,例如,對有符號數混用一個無符號操作數仍會産生期望的結果;這些類型的長度不受限制,所以就不存在所謂的精度問題。所有的操作——包括移位和除法——都有完整且一致性的定義,混合不同的類型也能得到期望的結果,如,當x爲一個C內置類型,而y是一個算法C類型時,表達式x+y和y+x均能返回相同的結果。

運行時間

我們的目的是爲了在支持任意長度類型及避免用戶碰到前述語義問題的前提下,得到使用內置類型(位寬不超過64位)手工C/C++編碼優化過的運行時間。算法C類型是爲快速執行及易于合成的語義而設計及實現的,所有操作的位寬由C++編譯器靜態確定,這就避免了動態內存分配,減少了運行時間,也使得語義更加易于合成。另外,實現也爲速度進行了優化,因此可能會調用更多的專用及高效代碼,充分利用了當今編譯器的優化特性。

在C/C++算法設計中使用任意位寬

表1:規格化爲ac_fixed的運行時間比較

表1是當定點算術用算法C定點類型ac_fixed來建模時,各種不同的運行時比較;float的實現在圖2中,sc_fixed_fast爲SystemC中精度受限的定點數據類型,sc_fixed爲任意精度的定點類型。實際中對FIR濾波器進行10^8次調用,TRN/WRAP的運行時間爲6.5秒,RND/SAT爲29秒。其他類型的運行時間也能從這張表中依次推出,如sc_fixed TRN/WRAP將花費6.5s × 227 = 1476s(將近25分鍾)。作爲參考,圖1中的算法(使用無定點建模的float)花費時間爲3.5s(比起使用定點建模的ac_fixed,慢了近兩倍)。

上述的運行時間數據,均由GCC 4.1.1測量得來,而在之前版本的GCC或Visual C++ 2005中得到的數據大致接近。

另外,運行時間也能通過整型數據類型或位操作進一步縮短。表2爲一個DCT算法的相應結果,它由一個每次讀寫2位的移位操作得來,與此對比的運行時間爲SystemC精度受限的sc_int與任意長度的sc_bigint。

在C/C++算法設計中使用任意位寬

表2:規格化爲ac_int的運行時間比較

結論

基于通用標准ANSI C++,這種新的整數與定點算法C類型答應算法及系統設計者指定任意位寬,從而提供比傳統數據類型高200倍的仿真效率。這些新數據類型可成爲C-to-RTL設計鏈中非常有價值的一環,及在整個實現流程中保證了任意位寬的精度。

· 湖北宜昌三峽壩區水面驚現神秘動物

近日,湖北宜昌,一段視頻在當地熱傳:有網友在三峽壩區拍到神秘動物,體型碩大數米長...

· 什麽是語段?語段的類型以及和句群、段落的區別與聯系是什麽?

句群是最高級的語言單位。 段落(自然段)是章法單位...

· 十八部好看的賭石類小說

以下是十八部(排名不分先後)好看的賭石類小說的簡介,喜歡的朋友可以去搜索書名閱讀...

 
任意規模指派問題的C++類實現
一.指派問題在生活中經常遇見這樣的問題,有n項任務要求n個人完成,這n個人完成各項任務的效率(或所需時間)不同,于是産生指派哪個人去完成哪項任務的問題,這類問題稱爲指派問題或分派問題。1.指派問題的數學模...查看完整版>>任意規模指派問題的C++類實現
 
閑談C++算法封裝:窮舉法
  將算法獨立抽象出來,在C++中算不上新鮮:STL中就封裝了不少高效、健壯、靈活的泛型組件及對應的基礎算法,工藝之高、適用性之強,非平常我輩所輕易能及。這裏不打算(也暫沒有能力打算)以STL這樣的工業級要求來...查看完整版>>閑談C++算法封裝:窮舉法
 
遺傳算法的一個例子(C/C++)
遺傳算法的一個例子(C/C++)作者:newsuppy摘要: m個工件分配給m架機床的效益最優化問題,使用一種遺傳算法解決。 一,已知效益矩陣EeijM1M2M3M4M5J15648 3J264985J343254J472453J536455 Jn爲工件,Mn爲機床,矩陣對應...查看完整版>>遺傳算法的一個例子(C/C++)
 
C++ Primer 泛型算法閱讀有感
這個月做VC++界面又有些精進,經驗值也有所提高,但對于GDI操作還是接觸不多,抽點時間又看了十幾頁C++ Primer,模板的推演、泛型算法部分也是朦胧:主要是沒實際應用過,但是不能說是白看了,但對于思維的延伸還是有...查看完整版>>C++ Primer 泛型算法閱讀有感
 
閑談C++算法封裝:窮舉法
將算法獨立抽象出來,在C++中算不上新鮮:STL中就封裝了不少高效、健壯、靈活的泛型組件及對應的基礎算法,工藝之高、適用性之強,非尋常我輩所輕易能及。這裏不打算(也暫沒有能力打算)以STL這樣的工業級要求來...查看完整版>>閑談C++算法封裝:窮舉法