9,216 research outputs found
An Algorithmic Framework for Efficient Large-Scale Circuit Simulation Using Exponential Integrators
We propose an efficient algorithmic framework for time domain circuit
simulation using exponential integrator. This work addresses several critical
issues exposed by previous matrix exponential based circuit simulation
research, and makes it capable of simulating stiff nonlinear circuit system at
a large scale. In this framework, the system's nonlinearity is treated with
exponential Rosenbrock-Euler formulation. The matrix exponential and vector
product is computed using invert Krylov subspace method. Our proposed method
has several distinguished advantages over conventional formulations (e.g., the
well-known backward Euler with Newton-Raphson method). The matrix factorization
is performed only for the conductance/resistance matrix G, without being
performed for the combinations of the capacitance/inductance matrix C and
matrix G, which are used in traditional implicit formulations. Furthermore, due
to the explicit nature of our formulation, we do not need to repeat LU
decompositions when adjusting the length of time steps for error controls. Our
algorithm is better suited to solving tightly coupled post-layout circuits in
the pursuit for full-chip simulation. Our experimental results validate the
advantages of our framework.Comment: 6 pages; ACM/IEEE DAC 201
Time-domain analysis of large-scale circuits by matrix exponential method with adaptive control
We propose an explicit numerical integration method based on matrix exponential operator for transient analysis of large-scale circuits. Solving the differential equation analytically, the limiting factor of maximum time step changes largely from the stability and Taylor truncation error to the error in computing the matrix exponential operator. We utilize Krylov subspace projection to reduce the computation complexity of matrix exponential operator. We also devise a prediction-correction scheme tailored for the matrix exponential approach to dynamically adjust the step size and the order of Krylov subspace approximation. Numerical experiments show the advantages of the proposed method compared with the implicit trapezoidal method. © 1982-2012 IEEE.published_or_final_versio
Circuit simulation via matrix exponential method for stiffness handling and parallel processing
We propose an advanced matrix exponential method (MEXP) to handle the transient simulation of stiff circuits and enable parallel simulation. We analyze the rapid decaying of fast transition elements in Krylov subspace approximation of matrix exponential and leverage such scaling effect to leap larger steps in the later stage of time marching. Moreover, matrix-vector multiplication and restarting scheme in our method provide better scalability and parallelizability than implicit methods. The performance of ordinary MEXP can be improved up to 4.8 times for stiff cases, and the parallel implementation leads to another 11 times speedup. Our approach is demonstrated to be a viable tool for ultra-large circuit simulations (with 1.6M ∼ 12M nodes) that are not feasible with existing implicit methods. © 2012 ACM.published_or_final_versio
Globally stable, highly parallelizable fast transient circuit simulation via faber series
Time-domain circuit simulation based on matrix exponential has attracted renewed interested, owing to its explicit nature and global stability that enable millionth-order circuit simulation. The matrix exponential is commonly computed by Krylov subspace methods, which become inefficient when the circuit is stiff, namely, when the time constants of the circuit differ by several orders. In this paper, we utilize the truncated Faber Series for accurate evaluation of the matrix exponential even under a highly stiff system matrix arising from practical circuits. Experiments have shown that the proposed approach is globally stable, highly accurate and parallelizable, and avoids excessive memory storage demanded by Krylov subspace methods. © 2012 IEEE.published_or_final_versio
Differential-Algebraic Equations
Differential-Algebraic Equations (DAE) are today an independent field of research, which is gaining in importance and becoming of increasing interest for applications and mathematics itself. This workshop has drawn the balance after about 25 years investigations of DAEs and the research aims of the future were intensively discussed
Physics-based passivity-preserving parameterized model order reduction for PEEC circuit analysis
The decrease of integrated circuit feature size and the increase of operating frequencies require 3-D electromagnetic methods, such as the partial element equivalent circuit (PEEC) method, for the analysis and design of high-speed circuits. Very large systems of equations are often produced by 3-D electromagnetic methods, and model order reduction (MOR) methods have proven to be very effective in combating such high complexity. During the circuit synthesis of large-scale digital or analog applications, it is important to predict the response of the circuit under study as a function of design parameters such as geometrical and substrate features. Traditional MOR techniques perform order reduction only with respect to frequency, and therefore the computation of a new electromagnetic model and the corresponding reduced model are needed each time a design parameter is modified, reducing the CPU efficiency. Parameterized model order reduction (PMOR) methods become necessary to reduce large systems of equations with respect to frequency and other design parameters of the circuit, such as geometrical layout or substrate characteristics. We propose a novel PMOR technique applicable to PEEC analysis which is based on a parameterization process of matrices generated by the PEEC method and the projection subspace generated by a passivity-preserving MOR method. The proposed PMOR technique guarantees overall stability and passivity of parameterized reduced order models over a user-defined range of design parameter values. Pertinent numerical examples validate the proposed PMOR approach
Exponential integrators for a Markov chain model of the fast sodium channel of cardiomyocytes
The modern Markov chain models of ionic channels in excitable membranes are
numerically stiff. The popular numerical methods for these models require very
small time steps to ensure stability. Our objective is to formulate and test
two methods addressing this issue, so that the timestep can be chosen based on
accuracy rather than stability.
Both proposed methods extend Rush-Larsen technique, which was originally
developed to Hogdkin-Huxley type gate models. One method, "Matrix Rush-Larsen"
(MRL) uses a matrix reformulation of the Rush-Larsen scheme, where the matrix
exponentials are calculated using precomputed tables of eigenvalues and
eigenvectors. The other, "hybrid operator splitting" (HOS) method exploits
asymptotic properties of a particular Markov chain model, allowing explicit
analytical expressions for the substeps.
We test both methods on the Clancy and Rudy (2002) INa Markov chain model.
With precomputed tables for functions of the transmembrane voltage, both
methods are comparable to the forward Euler method in accuracy and
computational cost, but allow longer time steps without numerical instability.
We conclude that both methods are of practical interest. MRL requires more
computations than HOS, but is formulated in general terms which can be readily
extended to other Markov Chain channel models, whereas the utility of HOS
depends on the asymptotic properties of a particular model.
The significance of the methods is that they allow a considerable speed-up of
large-scale computations of cardiac excitation models by increasing the time
step, while maintaining acceptable accuracy and preserving numerical stability.Comment: 9 pages, 5 figures main text + 14 pages, 1 figure appendix, as
submitted in final form to IEEE TBME 2014/11/11. Copyright IEEE (2014
- …