作者：Sharad Sinha
博士生
新加坡南洋理工大學(xué)
[email protected]

Vivado HLS配合C語言等高級語言能幫助您在FPGA上快速實現(xiàn)算法,。

為了幫助用戶了解Vivado HLS如何工作,，我們不妨以矩陣乘法設(shè)計為例逐步剖析從設(shè)計描述（C/C++/SystemC）到FPGA實現(xiàn)整個端對端綜合流程,。矩陣乘法在許多應(yīng)用中都很常見,，并廣泛用于圖像和視頻處理、科學(xué)計算和數(shù)字通信,。本項目中的所有結(jié)果均使用Vivado HLS 2012.4生成,，搭配使用賽靈思 ISE®軟件（14.4版）進行物理綜合和布局布線。此外,，這一流程還采用了ModelSim和GCC-4.2.1-mingw32vc9進行RTL協(xié)同仿真,。

圖1顯示了簡單的綜合流程，從C/C++/SystemC設(shè)計開始,。C/C++/SystemC testbench用于驗證設(shè)計功能的正確性,，同時還可用于RTL和C的協(xié)同仿真。協(xié)同仿真包括驗證生成的RTL設(shè)計（.v或.vhd）功能,，這要使用C/C++/SystemC測試平臺而不是RTL測試平臺或者采用e或Vera驗證語言編寫的測試平臺,。時鐘周期約束設(shè)置了設(shè)計應(yīng)該運行的目標時鐘周期。設(shè)計將被映射到目標FPGA器件——賽靈思FPGA上,。

C語言的矩陣乘法
為了充分利用我們的矩陣乘法實例,，我們將探索矩陣乘法C語言實現(xiàn)方案的各種修訂版本，從而展示它們對綜合設(shè)計的影響,。這一過程將凸顯您在使用HLS進行原型設(shè)計和實際設(shè)計時需要注意的重要問題,。我們將跳過創(chuàng)建工程的有關(guān)步驟，因為您能很方便地在工具文檔中找到相關(guān)參考材料,。我們將重點介紹設(shè)計和實現(xiàn)等方面,。

在典型的Vivado HLS流程中，我們需要三個C/C++文件：源文件（包括待綜合的C函數(shù)）,、頭文件和通過main()函數(shù)調(diào)用描述testbench的文件,。

頭文件不僅包括源文件中使用的函數(shù)的聲明，也包括支持具有特定位寬的用戶定義數(shù)據(jù)類型的指令,。這也使得設(shè)計人員能夠采用與C/C++所定義標準位寬不同的位寬,。舉例來說，整形數(shù)據(jù)類型（int）在C語言中通常為32位長,，但是在Vivado HLS中您可指定用戶定義的數(shù)據(jù)類型,，例如只使用16位的“data”。

圖2顯示了用于矩陣乘法的簡單C函數(shù),。兩個矩陣mat1和mat2進行乘法,。為了簡單起見，兩個矩陣大小一樣,，都是兩行兩列。

在HLS流程中執(zhí)行的步驟如下：
• 第一步：創(chuàng)建工程
• 第二步：測試功能
• 第三步：綜合
• 第四步：RTL協(xié)同仿真
• 第五步：導(dǎo)出RTL / RTL實現(xiàn)

第一步編譯工程并在不同的設(shè)計文件中測試語法錯誤等,。第二步測試待實現(xiàn)的函數(shù)（在源文件中）功能是否正確,。在這一步驟中您將使用testbench執(zhí)行函數(shù)調(diào)用,，驗證其功能是否正確。如果功能驗證失敗,，您就需要返回來修改設(shè)計文件,。

第三步進行綜合，Vivado HLS綜合源文件中定義的函數(shù),。這一步的輸出包括C函數(shù)的Verilog和VHDL代碼（RTL設(shè)計）,，也包括目標FPGA的資源利用率估算和時鐘周期估算。此外,，Vivado HLS還可生成latency估算和回路相關(guān)的度量指標等,。

第四步是使用C testbench仿真生成的RTL。這一步叫做RTL協(xié)同仿真,，因為工具采用的就是之前用于驗證C源代碼的testbench,，現(xiàn)在則測試RTL的功能正確性。要成功完成這一步,，您系統(tǒng)（Windows或Linux）中的PATH環(huán)境變量應(yīng)包含ModelSim安裝的路徑,。此外，您還應(yīng)在ModelSim安裝文件夾中包含GCC-4.2.1-mingw32vc9套件,。

最后,，第五步就要將RTL導(dǎo)出為IP模塊，用于更大的設(shè)計中,，并由其它賽靈思工具進行處理,。您可將RTL導(dǎo)出為IP-XACT格式的IP模塊，也可導(dǎo)出為System Generator IP模塊或pcore格式的IP模塊,，進而用于賽靈思嵌入式設(shè)計套件,。導(dǎo)出Vivado生成的RTL時，您可以選擇工具的“評估”選項來評估布局布線后的性能并且運行RTL實現(xiàn),。在此情況下,， Vivado HLS工具會調(diào)用賽靈思ISE工具。要實現(xiàn)這一目的,，您的系統(tǒng)PATH環(huán)境變量需包括ISE安裝路徑,，Vivado HLS將會搜索ISE安裝。

當然,，您也不一定非要將Vivado生成的RTL導(dǎo)出為以上三種格式之一的IP模塊不可,。導(dǎo)出的格式文件可放在三個不同路徑下：
//impl/
或project_directory>//impl/或
//impl/。此外,，您也可在較大設(shè)計中使用Vivado生成的RTL,，或者將其本身用作頂層設(shè)計。當較大設(shè)計中例化導(dǎo)出的RTL時,，您應(yīng)注意相關(guān)接口要求,。

當綜合圖2中的C函數(shù)時,，您將獲得如圖3所示的RTL級實現(xiàn)方案。您會發(fā)現(xiàn),，實現(xiàn)方案中的矩陣1和矩陣2的元素被讀取到函數(shù),，并且積矩陣的元素被寫出。這樣,，實現(xiàn)方案假定“矩陣乘法”實體以外的存儲器能用來存儲矩陣1,、矩陣2和乘積矩陣的元素。表1介紹了信號描述,，表2則介紹了設(shè)計度量指標,。

表2：用于圖3所示設(shè)計的設(shè)計度量指標

表1：面向圖3中設(shè)計的信號描述

在表1中,，start,、done和idle信號與設(shè)計中控制數(shù)據(jù)路徑的有限狀態(tài)機（FSM）有關(guān)。您會發(fā)現(xiàn),，Vivado HLS生成的Verilog假定運算始于start信號,，并且輸出數(shù)據(jù)在ap_done信號從低變高開始有效。Vivado HLS生成的Verilog/VHDL將始終保持至少三個基本信號：ap_start,、ap_done和ap_idle,，此外還有ap_clk信號。這意味著不管您使用Vivado HLS實現(xiàn)什么設(shè)計,，設(shè)計latency都會限制您的流吞吐量,。圖2中的設(shè)計latency為69個時鐘周期，目標時鐘周期為3納秒,。這意味著在此特定案例中,，所有積矩陣元素需要69個時鐘周期可輸出。這樣,，您在至少69個時鐘周期前不能為設(shè)計提供新一組的輸入矩陣,。

圖3中所示的實現(xiàn)方案現(xiàn)在可能并不是您在FPGA上實現(xiàn)矩陣乘法時所預(yù)想的結(jié)果。您或許希望一款實現(xiàn)方案能讓您輸入矩陣,，并在內(nèi)部進行存儲和計算,，隨后讀取積矩陣元素。這顯然是圖2所示實現(xiàn)方案無法做到的,。該實現(xiàn)方案需要外部存儲器提供矩陣數(shù)據(jù)的輸入和輸出,。

調(diào)整代碼
圖4中的代碼能夠滿足您的需求，它是源文件的一部分,，應(yīng)該屬于C++文件而非此前的C文件,。您應(yīng)在頭文件matrixmultiply.h中包含另外兩個相關(guān)頭文件：hls_stream.h和ap_int.h。請注意，在圖2中,，當源文件為C文件時,，頭文件包含了ap_cint.h。頭文件ap_int.h和ap_cint.h有助于分別為C++和C源文件定義用戶定義的任意位寬的數(shù)據(jù)類型,。需要頭文件hls_stream.h來充分利用流接口，并且只有在源文件為C++語言時才能使用,。

圖4:用于矩陣乘法的重組源代碼

為了讓設(shè)計只接收輸入矩陣流,，并輸出積矩陣流，您應(yīng)在代碼中實現(xiàn)讀和寫數(shù)據(jù)流,。流接口就像FIFO,。默認情況下這個FIFO的深度為1。

表3 - 圖5中設(shè)計的信號描述

表4：圖5所示設(shè)計的設(shè)計度量指標

為了讓設(shè)計只接受輸入矩陣流,，并輸出積矩陣流，您應(yīng)在代碼中實現(xiàn)讀和寫數(shù)據(jù)流,。代碼hls::stream<> stream_name用于為讀和寫數(shù)據(jù)流命名,。這樣，d_mat1和d_mat2為讀取流而d_product為寫入流,。流接口就像FIFO那樣工作,。默認情況下，F(xiàn)IFO的深度為1,。您應(yīng)在Vivado HLS指令面板中通過選擇定義的數(shù)據(jù)流設(shè)置深度,。對于圖4中的代碼而言，每個數(shù)據(jù)流的深度都為4個數(shù)據(jù)單元,。請注意,，這里的（i,j）回路在（p,q）回路之前執(zhí)行，這是C++代碼的順序特性使然,。因此,，d_mat2數(shù)據(jù)流會在d_mat1數(shù)據(jù)流之后填滿。

完成數(shù)據(jù)流接口后,，您可應(yīng)用指令RESOURCE并通過指令面板選擇一個核,，從而選擇將矩陣映射到BRAM。否則將用觸發(fā)器和查找表（LUT）實現(xiàn)矩陣,。請注意,，指令面板只有當源文件在綜合視圖中保持有效時才是有效的。

圖5顯示了圖4中代碼的設(shè)計實現(xiàn)情況,。表3介紹了設(shè)計接口上可用的信號情況,。在表3中，d_product_V_full_n是低有效信號，當需要通知內(nèi)核積矩陣已滿時該信號為低,。但在實現(xiàn)方案中通常不需要這樣,。

表4顯示了3納秒時鐘周期約束下布局布線后的不同設(shè)計度量指標，包含了矩陣陣列映射到BRAM或分布式RAM的情況和未映射的情況,。您從表4中可以看到,，矩陣映射到單端口BRAM時，設(shè)計無法滿足3納秒的時序約束,。表中專門包含了這個結(jié)果,，說明您可用這種方法生成具有不同面積—時序參數(shù)的各種設(shè)計。此外,，您也可從表1看出,，雖然圖2中代碼的時延為69個時鐘周期，低于圖4中調(diào)整后的代碼的設(shè)計方案,，但這種設(shè)計需要矩陣乘法實體以外的存儲器,，這一點我們在上面已經(jīng)解釋過了。

實現(xiàn)方案的精度
就這里顯示的結(jié)果而言,，我將“data”這種數(shù)據(jù)類型定義為16位寬,。因此，所有矩陣元素（左,、右和積矩陣）都為16位寬,。矩陣乘法和加法運算不能實現(xiàn)全精度。您可選擇在頭文件中定義另一種32位寬的數(shù)據(jù)類型data_t1,，積矩陣的所有元素都采用這種數(shù)據(jù)類型,，因為16位數(shù)（左側(cè)矩陣元素）乘以另一個16位數(shù)（右側(cè)矩陣元素）最多得到32位寬。這樣,，資源利用率和時序結(jié)果將不同于表1和表4中的結(jié)果,。

調(diào)整后的源代碼顯示出同樣的源文件會帶來多種不同設(shè)計解決方案。在本例中,，一個設(shè)計解決方案采用BRAM,，而另一個沒有采用。在每個Vivado HLS工程目錄中,，您會看到Vivado HLS為不同的解決方案生成了不同的目錄,。在每個解決方案目錄中都有一個名叫impl（也就是implementation，實現(xiàn)方案）的子目錄,。在這個子目錄中,，您會看到名為Verilog或VHDL的目錄，具體取決于RTL實現(xiàn)階段使用什么樣的源代碼,。這個子目錄中也包含賽靈思ISE工程文件（文件擴展名為.xise）,。如果Vivado HLS生成的設(shè)計是您的頂層設(shè)計，那么您可以雙擊這個文件來啟動賽靈思ISE運行這個解決方案，并生成用于門級時序和功能仿真的布局布線后模型,。但您在Vivado HLS中不能做這種仿真,。

在ISE中啟動解決方案后，您應(yīng)給設(shè)計分配I/O引腳,。隨后您可在ISE Project Navigator中選擇“Generating Programming File”以生成比特流,。

在這一練習中，我們一步步完成了Vivado HLS一個實際的端對端流程,，并在FPGA上實現(xiàn)算法,。對于Vivado HLS中的許多高級特性而言，您應(yīng)了解您需要什么樣的硬件架構(gòu),，從而進行源代碼的調(diào)整。如需了解更多詳情,，《Vivado高層次綜合教程》（UG871; http://www.xilinx.com/support/documentation/sw_manuals/xilinx2012_2/ug87... ）和《Vivado設(shè)計套件用戶指南》（UG002; http://www.xilinx.com/support/documentation/sw_manuals/xil-inx2012_2/ug9... ）這兩個技術(shù)文檔對您大有裨益,。