Introduction to Model Order Reduction

1. Introduction to Model Order Reduction Luca Daniel Massachusetts Institute of Technology luca@mit.edu http://onigo.mit.edu/~dluca/2006PisaMOR www.rle.mit.edu/cpg

2. Model Order Reduction of Linear Systems • via Modal Analysis • via Rational function fitting (point matching) • via Quasi Convex Optimization • via Pade approximation (AWE) • Projection Framework • SVD, PCA, LVD, POD • Krylov Subspace Moment Matching Projection Methods • Arnoldi • PVL • PRIMA • Truncated Balance Realization (TBR) • Positive Real TBR • Distributed Systems (with Frequency Dependent Matrices)

3. Distributed Linear Systems Examples: • ODE’s with delays (e.g. full-waveintegral equation solvers) • frequency-dependent basis functions • frequency dependent discretizations • solvers using layered-media Green functions(e.g. for handling substrate or dielectrics) • NOTE: Distributed systems may have infinite order (e.g. delay)!!

4. Polynomial interpolation [Phillips96] • Polynomial approximation e.g. Taylor expansion, or a polynomial interpolation for A(s) • Convert to non-distributed model reduction problem • Performance: Fast and accurate in the frequency band of interest • Problem: Can not be used in a time domain circuit simulator because does not guarantee stability and passivity

5. Passivity condition on transfer function • For systems with immittance matrix representation, passivity is equivalent to positive-realness of the transfer function (no unstable poles) (impulse response is real) (no negative resistors) It means its real part is a positive for any frequency. Note: it is a global property!!!! FOR ANY FREQUENCY

6. Positive real transfer functionin the complex plane for different frequencies original system Active region Passive region

7. Why does polynomial interpolation failwhen applied to the Laplace parameter ‘s’? • Although accurate in the frequency band of interest • Polynomial interpolation is unlikely to preserve GLOBAL properties such as positive realness because it is GLOBALLY not well-behaved original system Active region Passive region

8. Observation: practical systems have some loss at any frequency • Most systems are non-ideal i.e. contain some small loss at any frequency i.e. can be described by strictly positive real functions original system Active region Passive region

9. reduced system • Proof: just choose accuracy of interpolation smaller than minimum distance from imaginary axis Using global uniformly convergent interpolants • If A(s) is strictly positive real, a GLOBALLY and UNIFORMLY convergent interpolant will eventually get close enough (for a large enough order M of the interpolant) and be positive-real as well. original system Active region Passive region

10. A good example of uniformly convergent interpolants: the Laguerre basis functions • Consider the family of basis functions: • They form a complete, rational, orthonormal basis over the imaginary axis which gives a uniformly convergent interpolant • no poles in RHP (stable) • (real time-domain representation)

11. Calculation of interpolation coefficients • Note: it is a bilinear transform that maps the Laguerre basis to Fourier series on the unit circle. Im{z} Im{s} Re{z} Re{s} • Hence in practice one can use FFT to calculate the interpolation coefficients: very efficient! • Note: FFT coefficients typically drop quickly and the series can be truncated to the first few M coefficients because field solver matrices A(s) are often smooth.

12. An implementation example:Two wires on a MCM package [D. DAC02] • Discretize Maxwell equations in integral form using PEEC • NOTE: system matrices are frequency dependent because the substrate is handled by layered Green functions package s x=- A(s)x+b u

13. Step 2: Interpolation (example) FFT coefficients of A (s) 10-4 10-5 10-6 10-7 10-8 0 10 20 30 40 50 60

14. Step 2: Interpolation (example) A (s) reconstructed fromfirst 5 out of 64FFT coefficients and compared to original A (s) Imaginary part Real part nH 5 4 3 2 1 0 0 10 20 30 40 50 60

15. Reduction procedure [D. and Phillips DAC01]

16. Reduction procedure [D. and Phillips DAC01]

17. Step 3: Realization (Multichip Module MCM example) • Real part of frequency response • Inductive part of frequency response nH x104Ohm original system 3000 distrib 5 5 reduced system 30 4 4 3 3 2 original system 2 reduced system 108 1010 106 107 108 109 1010 106 107 109 frequency frequency

18. Open issues for distributed systems • Guaranteeing positive realness relies on accuracy of the uniform interpolant. Hence if the matrices are NOT smooth, we might need a large order of the interpolant. • working on internal matrices might give smoother matrices • Laguerre basis functions are efficient since they use FFT. However equally spaced points on the unit circle correspond to non-equally spaced points on the imaginary axis accumulating around a reference center frequency.

19. Model Order Reduction of Linear Systems • via Modal Analysis • via Rational function fitting (point matching) • via Quasi Convex Optimization • via Pade approximation (AWE) • Projection Framework • SVD, PCA, LVD, POD • Krylov Subspace Moment Matching Projection Methods • Arnoldi • PVL • PRIMA • Truncated Balance Realization (TBR) • Positive Real TBR • Laguerre interpolation for Distributed Systems

20. Model Order Reduction of Linear Systems • via Modal Analysis • via Rational function fitting (point matching) • via Quasi Convex Optimization (use this for LARGE and passive systems, or for model construction from measurements, or for distributed systems) • via Pade approximation (AWE) • Projection Framework • SVD, PCA, LVD, POD • Krylov Subspace Moment Matching Projection Methods • Arnoldi • PVL (use this for HUGE systems if passivity is not an issue) • PRIMA (use this for HUGE and passive systems) • Truncated Balance Realization (TBR) (use this for SMALL systems) • Positive Real TBR (use this for SMALL and passive systems) • Laguerre interpolation for Distributed Systems (use this for LARGE systems with frequency dependent matrices, e.g. delays)