[发明专利]用去耦合电容实现集成电路供电网络噪声优化的快速方法

摘要

用去耦合电容对集成电路供电网络进行噪声优化的方法属于VLSI物理设计领域，尤其是布局布线领域中RLC电源线/地线网络噪声优化的技术范畴。其特点在于：它是一种利用计算机针对基于标准单元布图模式的专用集成电路(ASIC)的布图结构特点及其供电网络的结构特点，对供电网络的优化问题进行建模，并采用非线性规划方法对问题进行有效求解的过程。发明中，新颖之处是将去耦电容的面积作为优化的目标，使优化后得到的去耦电容面积尽可能小，为后面阶段的工作，例如插buffer，留下了充足的空白面积。求解过程中我们利用了以前工作中的优秀成果——等效电路的方法，使求解规模大大降低，很大提高了求解速度并且节省了大量的计算机内存。实验证明，本发明可以快速有效的对供电网络进行优化，并具有优化大规模集成电路供电网络的能力。

主权项

1.用去耦合电容实现集成电路供电网络噪声优化的快速方法，含有用放置去耦合电容于布图空白面积中以提高供电网络可靠性的步骤，其特征在于，它是一种利用计算机针对基于标准单元布图模式的专用集成电路(ASIC)的布图结构特点及其供电网络的结构特点，将用以保持节点电压平稳的去耦合电容的面积作为优化的目标，采用非线性规划的方法对供电网络噪声优化问题进行有效求解的方法，它把罚函数法、等效电路、共轭梯度法和扩展特勒根伴随网络的方法形成一个整体共同作用于电源/地线网络优化上，取得很好的优化效果；具体而言，它依次含有以下步骤：(1)计算机读入包含供电网络信息和单元布局信息的文件：供电网络信息包括供电网络节点之间的关联结构；节点之间的电阻值、电感值、电感初始电压和电流值、片上固有去耦合电容值(单元模块静止时表现出的电容特性)、电容初始电压和电流值以及各个节点连接的单元模块随时间变化的吸纳电流波形(利用PWL表示)，据此在计算机内建立电路的信息，并且标记所有中间点(Middle Node)和交汇点(Cross Node)；单元布局信息包括每个单元行中的单元数目、单元坐标以及单元之间的相对位置关系，根据这些信息，我们可以得到单元之间的空白面积的信息和各单元行上空白面积的分布情况，附加的去耦合电容就添加在这些空白区域上；(2)利用输入的参数信息对供电网络进行基于等效电路的瞬态分析，得到每个节点的电压波形；(3)若各节点的电压没有违反约束的情况出现，则不需对供电网络做优化；否则执行以下步骤：(4)将各节点的电压波形分段线性纪录，用于优化过程中梯度的计算；(5)确定优化的目标函数为添加的去耦电容的面积，在对应的约束条件下求最小值：minA=Σj∈M(wj×H)--(1),]]>M是允许添加去耦合电容的节点集合，M＝{1，...，m}，MN；N为供电网络中所有节点的集合，N＝{1，...，n}；H是标准单元的高，w_j是添加在节点j上的去耦电容c_j的宽度；约束条件如下：(5.1)电压降约束：即各节点相对于电源电压的电压降，用能有效衡量电压降噪声的量值，即每个节点上低于电压阈值的电压部分的积分值S_i表示：si=∫0Tmax(Vmin-vi(t),0)dt=∫t1t2(Vmin-vi(t))dt--(2)]]>[t₁，t₂]为i节点违反电压降约束的时间区间，V_min为电压闽值，v_i(t)为i节点的电压值；(5.2)电迁移约束：用[t’₁，t’₂]时间内，供电网络分支(p，q)两端节点电压差的绝对值和工艺允许的最大电流密度为σ时在分支长度l_pq上的电压降之差的积分值来表示：up,q=∫0Tmax[(|vp(t)-vq(t)|-ρlp,qσ),0]dt=∫t′1t′2(|vp(t)-vq(t)|-ρlp,qσ)dt--(3)]]>其中，v_p和v_q分别是电源分枝(p，q)两端节点的电压，v_p-v_q即该分枝上的电压差；ρ是方块电阻；l_pq是分枝(p，q)的长度；(5.3)去耦电容面积约束：即每单元行中，添加的去耦电容面积小于等于本单元行中的空白面积，而单元行的高度H是定值，因此可直接用宽度表示面积大小：dwr=max(Σy∈ND(r)wr,y-rwr,0),r∈NR--(4)]]>NR是单元行的集合；ND(r)是第r行中去耦电容的位置集合；w_r，y是指位于第r行y位置的去耦电容的宽度；rw_r是第r行上的空白面积的总宽度；dw_r是第r行上违反去耦电容面积约束的惩罚量，将在后面的罚函数法中应用；(5.4)去耦电容最大宽度约束：在某单元行中添加的任意一个去耦电容的宽度小于或等于本行中最大空白面积的宽度；ew_r，y＝max(w_r，y-rw_r，0)，r∈NR，y∈ND(r)，——(5)，ew_r，y是第r行第y个位置上违反去耦电容最大宽度约束的惩罚量，将在后面的罚函数法中应用；(6)构造罚函数f，将非线性约束优化问题转化为无约束优化问题：minf=A+pt]]>=A+α·(Σi∈Nsi2+Σ(p,q)∈Bupq2+Σr∈NRdwr2+Σr∈NR,y∈ND(r)ewr,y2),--(6),]]>其中，A是去耦电容总面积，B是电源网络中所有分枝的集合，P_t是罚项，α为罚因子；(7)设置罚因子α的初始值；每个可添加去耦电容的位置上的去耦电容初始值都设为允许的最小电容宽度，得到电容宽度向量作为下面无约束优化的初始解；另外给出误差限ε₁＞0；(8)求解当前无约束优化问题，得到当前最优解电容宽度向量W⁽¹⁾；(8.1)设当前无约束优化问题的初始解为电容宽度向量W⁽⁰⁾，误差限为ε₂＞0；(8.2)初始优化方向设为罚函数关于电容宽度向量W⁽⁰⁾的负梯度方向：设初始优化方向为P⁽⁰⁾，梯度方向为f(W⁽⁰⁾)，得P⁽⁰⁾＝-f(W⁽⁰⁾)    —(7)，(8.3)用扩展的伴随网络法求罚函数相对于电容宽度向量的梯度：▿f(W)=[∂f∂w1,∂f∂w2,···,∂f∂wj,···,∂f∂wm]T,j∈M--(8),]]>宽度为w_j的电容c_j连接在供电网络的节点j上，假设节点j位于第rj行的yj位置，则罚函数f相对于w_j的偏导数可以表示为：∂f∂wj=H+2α(dwrj+ewrj,yj)+2α[Σi∈Nsi∂si∂wj+Σ(p,q)∈Bupq∂upq∂wj]--(9),]]>rj∈NR，yj∈ND(rj)根据电容大小和宽度的关系：c=ϵoxTox×w×H---(10),]]>ε_ox和T_ox分别是门氧化层的介电常数和厚度；得∂f∂wj=H+2α(dwrj+ewrj,yj)+2αϵoxHTox[Σi∈Nsi∂si∂cj+Σ(p,q)∈Bupq∂upq∂cj]--(11),]]>rj∈NR，yj∈ND(rj)通过对原电路网络建立对应电路中各节点的一系列伴随网络，并作电路灵敏度分析后，得到：∂si∂cj=∫0Tvj,Ei′(t)×v.j(T-t)dt--(12),]]>∂upq∂cj=∫0Tvj,tpq′(t)×v.j(T-t)dt--(13),]]>=∫0T[vj,Ep′(t)-vj,Eq′(t)]×v.j(T-t)dt]]>其中，v.j(T-t)=dvj(T-t)dt]]>是原电路网络中j节点的电压对时间t的导数；v′_j，Ei(t)是对应于节点i建立的伴随网络中j节点的电压(v′_j，Ep(t)，v′_j，Eq(t)同理)，它可以表示如下：vj,Ei′(t)=Z(j)vEi′(t)=Z(j)G-1[Ei(t)+ECi(t)+ELi(t)]--(14),]]>其中，Z(j)是一个位置选择向量，该向量除了j位置上的元素取值为1，其它元素全部为0，即Z(j)＝[0，0，…，0，1，0，…，0]；G是节点电压方程组G*V＝I中的关系矩阵；E_i(t)是电流源向量，它的值是除节点i对应的元素值之外，其余的元素值全部为0，节点i对应的元素值当该节点违反电压降约束时为-1，当该节点不违反约束时为0；EC_i(t)和EL_i(t)是对应于节点i建立的伴随网络中的电容和电感元素做等效后得到的电流源向量；于是得到：Σi∈Nsi∂si∂cj+Σ(p,q)∈Bupq∂upq∂cj]]>=Σi∈Nsi(∫0Tvj,Ei′(t)×v.j(T-t)dt)+Σ(p,q)∈Bupq(∫0T[vj,Ep′(t)-vj,Eq′(t)]×v.j(T-t)dt)--(15),]]>=∫0T{v.j(T-t)×{Σi∈N(si×vj,Ei′(t))+Σ(p,q)∈B{upq×[vj,Ep′(t)-vj,Eq′(t)]}}}dt]]>依据式(14)，(15)得到：Σi∈Nsi∂si∂cj+Σ(p,q)∈Bupq∂upq∂cj=Z(j)G-1Itotal(t)--(16),]]>其中I_total＝I_new(t)+EC(t)+EL(t)            ——(17)，Inew(t)=Σi∈N[si+Ei(t)]+Σ(p,q)∈B{upq×[Ep(t)-Eq(t)]}--(18),]]>其中，EC(t)和EL(t)是将I_new作为伴随网络的电流激励时，电容和电感元素做等效后得到的电流源向量；从而得到：∂f∂wj=H+2α(dwrj+ewrj,yj)+2αϵoxHTox∫0T{[Z(j)×G-1×Itotal(t)]×v.j(T-t)}dt--(19),]]>rj∈NR，yj∈ND(rj)，j∈M(8.4)设迭代已经进行到了第k步，当前电容宽度向量和优化方向分别是为W^(k)和P^(k)，沿P^(k)进行一维搜索，找到非负一维搜索因子λ_k，使罚函数在P^(k)方向上达到最优，即：f(W^(k)+λ_kP^(k))＝minf(W^(k)+λP^(k))         ——(20)，(8.5)根据下式更新电容宽度向量W^(k+1)：W^(k+1)＝W^(k)+λ_kP^(k)    ——(21)，(8.6)如果‖f(W^(k+1))‖＜ε₂，本次无约束优化结束退出，转步骤9；否则设电容宽度向量的维数是m，如果k＝m，则W⁽⁰⁾＝W^(k+1)，f(G⁽⁰⁾)＝f(G^(k+1))，转向步骤8.2；如果k＜m，根据下式更新优化方向，回到步骤8.4继续优化；(9)根据优化结果计算罚项的值，如果p_t＜ε₁，当前电容宽度向量即是所求最优解，算法结束退出；否则更新罚因子α，设l＝l+1，回到步骤8继续优化；(10)由于实际可以添加的去耦电容值不是连续的，具有最小单位数值的限制，我们对最终得到的优化结果进行调整；若某位置添加的去耦电容值小于最小数值则将其合并到邻近的节点上，以形成较大的电容，并将这些电容调整为最小单位数值的倍数。