The continuous
scaling for fabrication technologies of
electronic circuits demands the design of new
and improved simulation techniques for
integrated circuits. Therefore, this work shows
how the hypersphere technique can be adapted and
applied to trace a multiparameter homotopy.
Besides, we present a path-following technique
based on circles (evolved from hypersphere),
which is faster, and simpler to be implemented
than hypersphere technique. Last, a comparative
analysis between both techniques applied to
simulation of circuits with bipolar transistors
will be shown.
1. Introduction
The increment of the complexity of circuits influence the scientific progress in the simulation techniques area for integrated circuits. Also, homotopy techniques have been introduced as a useful tool in the area of operating point solution for circuits [1–5], due to the Newton-Raphson (NR) method (widely used) which shows convergence problems [6] like oscillation and divergence.
2. Multiparameter Homotopy
The first step to formulate a homotopy is to establish the equilibrium equation to be solved; it is formulated from Kirchhoff laws, being defined asf(x)=0,wheref:∈Rn⟶Rn,
where x represents the electrical variables of the circuit and n is the number of electrical variables.
Multiparameter homotopies [7–9] are characterized by adding more than one extra homotopy parameter to the equilibrium equation. When homotopy parameters are adjusted to zero, the solution for H(·) becomes trivial, and when parameters reach value of one, then the operating point is located. The multiparameter homotopy function can be represented asH(f(x),λ1,λ2,…,λk)=0,
where homotopy parameters are λ1,λ2,…,λk∈[0,1] and k is the number of homotopy parameters.
Multiparameter homotopy [7] has been proposed in order to avoid fork bifurcations, singularities, among other problems that can be encountered with homotopy paths. Besides, as for the uniparametric [2] and multiparameter homotopies, the tracing technique [10, 11] is a fundamental tool capable of affecting the convergence, speed, and number of solutions located. Therefore, it is proposed to apply two tracing techniques for multiparameter homotopy, both will be described in the following sections.
3. Tracing Techniques
In order to apply tracing techniques described in this paper, a biparametric homotopy based in Newton’s homotopy method will be used as an example:H(f(x),λ1,λ2)=f(x,λ2)-(1-λ1)f(xi,0),whereH:∈Rn+1×R⟶Rn.
With the existence of two parameters (λ1 and λ2), two simultaneous deformations or transformations are produced: one in function f and another in function H. When [x,λ1,λ2]=[xi,0,0], thenH(f(x),λ1,λ2)=f(xi,0)-f(xi,0)=0.
Hence, homotopy function is satisfied. Besides, when [λ1,λ2]=[1,1] becomesH(f(x),λ1,λ2)=f(x),
so the found solution of H is the solution of the equilibrium equation. Nonetheless, as function H has two extra variables, it is necessary to add two equations to the system H in order to be solved using more conventional techniques like NR.
(1) Equation n+1
One equation is added to define path λ1-λ2, which will be named parametric function M(λ1,λ2). This equation traverses three points [λ1,λ2]: p1=[0,0], p2=[A,B], and p3=[1,1]. The proposed equation is
M(λ1,λ2)=-λ1+(λ2+(B(-1+A)/(AB+1-2A)))(-((-1+2A-B)λ2/(AB+1-2A))+2(B(-1+A)/(AB+1-2A))),
where p2 is defined by user, as shown in Figure 1(a). The range of values for A and B is [0,1].
Hypersphere equation is added [12]:
S(⋅)=(x1-c1)2+(x2-c2)2+⋯+(λ1-cn+1)2+(λ2-cn+2)2-r2,
where c is the center of the hypersphere (which adjusts its value each iteration) and r≪1 is the hypersphere radius (step size).
The summary of the procedures consists in the following steps [12] (see Figure 1(b)).
The first sphere is established S0 with center located at t0=[xi,p1] and the equation system is solved (3.1), (3.4), and (3.5) using the NR method (setting t0 as an initial point), locating point t1.
A new hypersphere S1 is created with center at t1.
Using points t0 and t1, it is possible to create a prediction, which touches hypersphere S1 at point k1; it is used as initial point for the NR method, until locating point t2 on the homotopy path.
Steps 2 and 3 are successively repeated (updating hypersphere’s center after each iteration) until crossing point p3.
Points before and after p3 are used to perform an interpolation [13]. The type of interpolation employed in this paper is linear multidimensional interpolation, which produces an approximation xa of solution xs for the equilibrium equation.
Finally, using as initial point xa in the NR method, the precision for the operating point xs is improved.
It is possible to replace (3.5) for the circle equation, in function of the homotopy parameters:C(⋅)=(λ1-cn+1)2+(λ2-cn+2)2-r2,
where r≪1. The rest of the steps to implement the numerical continuation are the same as the hypersphere technique already described.
4. Study Case: Circuit with Bipolar Transistors and a Diode
The following circuit [14] (see Figure 2) contains nine solutions and has become the reference circuit for the homotopy applied to circuit analysis. Using the system reported by [14], equilibrium equation is augmentedf(v1,v2,v3,v4,λ2)={f1=6.103168Is(e40v1-1)λ2+4.36634v2+2.863168Is(e40v2-1)-12,f2=5.4v1+3.58Is(e40v1-1)λ2+6.62Is(e40v2-1)+v3+0.7Is(e40v3-1)+0.5Is(e40v4-1)-22,f3=6.103168Is(e40v3-1)+2.863168Is(e40v4-1)λ2+4.36634v4-12,f4=v1+0.7Is(e40v1-1)λ2+0.5Is(e40v2-1)+5.4v3+3.58Is(e40v3-1)+6.62Is(e40v4-1)λ2-20,
where Is=10-6, [v1,v2,v3,v4] are the voltage drop between base-emitter terminals for each transistor in the circuit (see Figure 2) and λ2 is the second homotopy parameter. The complete homotopy formulation can be established by using the following equations:
the augmented equilibrium equation (4.1),
the homotopy function equation (3.1),
the parametric function (3.4),
the hypersphere function (3.5) or circle function (3.6), which depends on the selected technique.
Chua’s circuit.
Table 1 presents in summary the results of performing the tracing of four paths with different initial points (xi1, xi2, xi3, and xi4), each one. This process was repeated for both tracing techniques, showing the results in Figures 3(a) and 3(b). Comparing both techniques, there are two interesting conclusions to be highlighted: first, homotopy paths traced from the same initial point lead to the same solution; in fact, comparing paths point by point, it can be observed that it is the same path; second, despite that paths are identical, the circle technique required a fixed number of iterations (48); those are much less required than using the hypersphere method. In fact, from Table 1, it can be concluded that, at best (choosing initial point at xi1), the tracing technique for circles required 8.95 times less CPU time (by using MAPLE 15 software in an Intel Quad Core i7 processor at 2.6 GHz) than the hypersphere technique. Both tracing techniques employed radius of r=0.03 and a parametric function M with p2=[0.2,0.3].
Relevant points for homotopy simulations.
Hypersphere—Init. point where [λ1,λ2]=[0,0]
No. Iter
Time (Sec)
Operating point [v1,v2,v3,v4]where [λ1,λ2]=[1,1]
xi1=[-5,-5,-5,-5]
519
7.70
xs1=[0.3830,-3.5446,0.3851,-4.0990]
xi2=[-1,-2,-1,0]
202
3.59
xs2=[0.3869,-4.6321,-0.8002,0.3775]
xi3=[-5,-0.5,-5,0]
216
4.06
xs3=[-0.5136,0.3775,-0.9682,0.3775]
xi4=[-1,0,0,0]
168
3.13
xs4=[-1.0510,0.3775,0.3845,-3.9542]
Circle—Init. point where [λ1,λ2]=[0,0]
No. Iter
Time (Sec)
Operating point [v1,v2,v3,v4]where [λ1,λ2]=[1,1]
xi1=[-5,-5,-5,-5]
48
0.86
xs1=[0.3830,-3.5446,0.3851,-4.0990]
xi2=[-1,-2,-1,0]
48
0.89
xs2=[0.3869,-4.6321,-0.8002,0.3775]
xi3=[-5,-0.5,-5,0]
48
0.74
xs4=[-0.5136,0.3775,-0.9682,0.3775]
xi4=[-1,0,0,0]
48
0.82
xs4=[-1.0510,0.3775,0.3845,-3.9542]
Homotopy trajectories v2-λ1.
Hypersphere technique
Circle technique
Circles technique can be modified changing one of the homotopy parameters by an electrical variable of interest. For instance, the simulation was repeated from initial point xi1, only changing the circle from (3.6) byC(v1,λ1)=v12+λ12-r2,
where v1,r∈ℜ. The result was that the homotopy path already known was traced (see Figure 3(b)) with a total of 191 iterations (locating the same solution at xs1). Also, it is possible to use one of the two homotopy parameters with more than one electrical variable, to implement a reduced hypersphere. Therefore, in a forthcoming work the study of circles technique will be expanded and a possible application to simulate VLSI circuits will also be discussed.
5. Conclusion
This work showed that it is possible to use the hypersphere technique to trace multiparameter homotopies. Besides, a tracing technique derived from hypersphere (circles) was introduced, which is simpler to program and faster than the hypersphere technique. These results make the circles technique an attractive tool to trace multiparameter homotopies.
GoldgeisserL. B.GreenM. M.A method for automatically finding multiple operating points in nonlinear circuits200552477678410.1109/TCSI.2005.8443592132258MelvilleR. C.TrajkovićL.Artificial parameter homotopy methods for the DC operating point problem1997126861877TrajkovićL.MelvilleR. C.FangS. C.Passivity and no-gain properties establish global convergence of a homotopy method for DC operating pointsProceedings of the IEEE International Symposium on Circuits and Systems19902-s2.0-0025665877Vázquez-LealH.Hernández-MartínezL.Sarmiento-ReyesA.Double-bounded homotopy for analysing nonlinear resistive circuitsProceedings of the IEEE International Symposium on Circuits and SystemsMay 20052-s2.0-6764911155310.1109/ISCAS.2005.1465309RoychowdhuryJ. S.MelvilleR. C.Homotopy techniques for obtaining a DC solution of large-scale MOS circuitsProceedings of the 33rd Annual Design Automation Conference19962-s2.0-0029710304KelleyC. T.2003Philadelphia, Pa, USASIAMFundamentals of Algorithms1998383WolfD. M.SandersS. R.Multiparameter homotopy methods for finding DC operating points of nonlinear circuits1996431082483810.1109/81.5389891416903RoychowdhuryJ.MelvilleR.Delivering global DC convergence for large mixed-signal circuits via homotopy/continuation methods200625166782-s2.0-3134445191610.1109/TCAD.2005.852461Vázquez-LealH.Hernández-MartínezL.Sarmiento-ReyesA.Murphy-ArteagaR. S.Improving multi-parameter homotopy via symbolic analysis techniques for circuit simulation2Proceedings of the European Conference on Circuit Theory and Design2003402405Vázquez-LealH.Hernández-MartínezL.Sarmiento-ReyesA.Numerical continuation scheme for tracing the double bounded homotopy for analysing nonlinear circuitsProceedings of the International Conference on Communications, Circuits and Systems20052-s2.0-29844434602AllgowerE. L.GeorgK.1994YamamuraK.Simple algorithms for tracing solution curves19934085375412-s2.0-002764617610.1109/81.242328ZBL0802.65150SosonkinaM.WatsonL. T.StewartD. E.Note on the end game in homotopy zero curve tracking199622328128710.1145/232826.2328431417312ZBL0884.65044UshidaA.ChuaL. O.Tracing solution curves of nonlinear equations with sharp turning points198412112174112010.1002/cta.4490120102