30 [P]. 2012-4-3.
[0018] (10) Gang Peter Fang, et al. Method for generating and evaluating a table model for circuit simulation:!]. S. Patents, 341, 566[P]. 2012-12-25.
[0019] (11) Shima T, Sugawara T, Moriyama S, et al. Three-dimensional table look-up MOSFET model for precise circuitsimulation[J]. Solid-State Circuits,IEEE Journal of, 1982, 17 (3): 449-454.
[0020] (12) Jiying X, Tao L, Zhiping Y. Accurate and fast table look-up models for leakage current analysis in 65nm CMOS technology[J]. Journal of Semiconductors, 2009, 30 (2): 024004.
[0021] (13) Fang G P, Yeh D C, et al. Fast, accurate MOS table model for circuit simulation using an unstructured grid and preserving monotonicity [C]// Proceedings of the2005Asia and South Pacific Design Automation Conference. ACM, 2005:1102-1106.
[0022] (14) Kiran Kumar Gullapalli, et al. Circuit Simulation acceleration using model caching:!]. S. Patent2013/0054217[P]. 2013-2-28.
[0023](巧）Subramaniam P. Modeling MOS VLSI circuits for transient analysis[J]. Solid-State Circuits, IE邸Journal of, 1986,21 (2) :276_285.
[0024] (1巧Lewis D M. Device model approximation using 2^ trees [J].Computer-Aided Design of Integrated Circuits and Systems，IEEE Transactionson，1990,9(1):30-38. [002引（17) Gondro E, Kowarik 0, Knobl inger G, et al. When do we need non-quasistaticCMOS民F-models?[C]//Custom Integrated Circuits, 2001，IEEE Conferenceon. IEEE, 2001:377-380.
[002引（18) Ng A F L, Ko P K, ChanM. Determining the onset frequency of nonquasistatic effects of the MOSFET in AC simulation[J]. Electron Device Letters, IE邸，2002, 23(1) : 37-39.
[0027] (19)Schrom G, Stach A, Se化erherr S. An inteipolation based MOSFET model for low-voltage applications[J]. Microelectronics journal, 1998, 29(8):529-534

【發(fā)明內(nèi)容】

[002引本發(fā)明的目的是解決當前基于"樹狀"結(jié)構(gòu)模型近似方法（TBMA)構(gòu)建的非均勻網(wǎng) 格表格模型中網(wǎng)格單元查找速度慢的問題，提供一種用于集成電路仿真的晶體管（MOSFET) 精確近似表格查找模型的建模和估值方法。本發(fā)明能實現(xiàn)在非均勻網(wǎng)格模型中不受劃分規(guī) 模限制的快速表格查找方法且繼承了TBMA自適應(yīng)劃分的優(yōu)勢，同時保證所建立晶體管表 格模型估值計算的連續(xù)性和插值的單調(diào)性，提出了相應(yīng)的表格模型建立和估值方法。
[0029] 為了達到上述目的，本發(fā)明的技術(shù)內(nèi)容是：
[0030] 本方法通過建立非線性電路中晶體管的非均勻網(wǎng)格表格模型，借助簡單的哈希映 射和輔助的查找表實現(xiàn)仿真過程中待估點的快速查找及其對應(yīng)物理參數(shù)的估值。
[0031] 本發(fā)明中由二元空間分割樹根據(jù)預(yù)先設(shè)定的精度要求自適應(yīng)地構(gòu)建晶體管物理 參數(shù)（電流、電荷）的非均勻網(wǎng)格模型，加細劃分將樹狀的網(wǎng)格單元存儲結(jié)構(gòu)轉(zhuǎn)化為多維數(shù) 組的網(wǎng)格單元存儲結(jié)構(gòu)，得到精細化的非均勻網(wǎng)格w保證表格模型估值計算的連續(xù)性和單 調(diào)性；
[0032] 在非均勻網(wǎng)格中對物理參數(shù)估值所需網(wǎng)格單元進行查找的關(guān)鍵思想是在非均勻 網(wǎng)格的每個維度上由均勻標尺構(gòu)建輔助的查找表；將查找表項映射為相應(yīng)維度上的網(wǎng)格單 元下標，便可通個待估點在該維度上坐標分量的簡單哈希映射確定其在該維度上的網(wǎng)格下 標；相比于TBMA，本方法能顯著提高非均勻網(wǎng)格表格模型中網(wǎng)格單元的查找速度，對待估 點查找的時間復(fù)雜度由O(logN)下降為0(C)，并且其查找速度與網(wǎng)格規(guī)模無關(guān)，從而可W 在不犧牲查找速度的前提下，提高非均勻網(wǎng)格的劃分粒度W獲得精度更高的表格模型。
[0033] 本發(fā)明提出的方法理論上可W用于任意N維（N> 2)解析模型的近似，除了本發(fā) 明實現(xiàn)的BSIM3中的二維和H維表格模型，還可W為其它器件巧日H維器件，雙極型晶體管 等）建立所需維度的表格模型。
[0034] 具體的，本發(fā)明的晶體管精確近似表格查找模型的建模和估值方法其特征在于， 其包括步驟：
[00巧]步驟一：給定一個器件物理參量函數(shù)的N維封閉定義域空間，在該定義域空間上 構(gòu)建二元空間分割樹（BSPtree)，根據(jù)預(yù)先設(shè)定的誤差闊值和解析的物理參量函數(shù)值對N 維定義域空間進行劃分并記錄下每個維度上的劃分點；
[0036] 步驟二：精細化步驟一得到的劃分空間，設(shè)N個維度上劃分點的數(shù)目分別為 叫，叫，k,由步驟一記錄的劃分點構(gòu)建規(guī)模為（ni-l)X(n廠1)的的N維非均勻網(wǎng)格 (non-uniformgrid)，網(wǎng)格的每一個單元（Cell)代表精細化后的劃分空間；
[0037] 步驟H;為N維非均勻網(wǎng)格在N個維度上的下標分別構(gòu)建N個查找表(一維數(shù)組)， 每一個查找表對應(yīng)一個間隔均勻的標尺（間隔遠小于所有網(wǎng)格單元在該維度上的最小邊界 長度)，查找表項數(shù)由預(yù)先設(shè)定的粒度值決定；
[0038] 步驟四：對N維定義域空間中的一個坐標點進行估值時，首先將該坐標在N個維度 上的分量映射到步驟H構(gòu)建的N個查找表中，得到包含該坐標點的劃分空間在N個維度上 的網(wǎng)格下標，再由網(wǎng)格下標訪問到相應(yīng)網(wǎng)格單元，最終在該網(wǎng)格單元中對待求坐標點進行 H次插值得到其函數(shù)值W及其在N個維度上的近似偏微分值。
[0039] 本發(fā)明中，步驟一，二元空間分割樹中的每一個結(jié)點（BSPnode)與N維定義域空間 中的每一個劃分區(qū)域（凸空間）一一對應(yīng)，此數(shù)據(jù)結(jié)構(gòu)用于存儲相應(yīng)劃分區(qū)域的相關(guān)頂點的 解析函數(shù)值、邊界W及個凸邊中點的近似偏微分值（有限差分近似）。
[0040] 本發(fā)明步驟二中，N維網(wǎng)格的單元與二元空間分割樹的結(jié)點具有相同的數(shù)據(jù)結(jié)構(gòu)； 對于精細化前后邊界不變的劃分區(qū)域所對應(yīng)的網(wǎng)格單元可直接復(fù)制二元分割樹中的相應(yīng) 結(jié)點，若邊界發(fā)生變化，則網(wǎng)格單元中的數(shù)據(jù)成員需要重新計算。
[0041] 本發(fā)明步驟H中，對于在各維度上網(wǎng)格單元邊界沒有和其對應(yīng)標尺"刻度"對齊的 地方，可移動網(wǎng)格邊界來對齊標尺W保證在網(wǎng)格邊界附近的正確查找訪問。
[0042] 本發(fā)明步驟四中，查找表的下標是坐標分量的哈希映射，該哈希函數(shù)由標尺間隔 簡單確定，將給定的坐標分量帶入該哈希函數(shù)便可得到查找表下標，該下標對應(yīng)的查找表 項的值即是相應(yīng)維度上的網(wǎng)格下標。
[0043] 本發(fā)明中，當給定的N維物理參量函數(shù)在某一個維度上變化較慢時，可W固定函