10,717 research outputs found
On the Generation of Large Passive Macromodels for Complex Interconnect Structures
This paper addresses some issues related to the passivity of interconnect macromodels computed from measured or simulated port responses. The generation of such macromodels is usually performed via suitable least squares fitting algorithms. When the number of ports and the dynamic order of the macromodel is large, the inclusion of passivity constraints in the fitting process is cumbersome and results in excessive computational and storage requirements. Therefore, we consider in this work a post-processing approach for passivity enforcement, aimed at the detection and compensation of passivity violations without compromising the model accuracy. Two complementary issues are addressed. First, we consider the enforcement of asymptotic passivity at high frequencies based on the perturbation of the direct coupling term in the transfer matrix. We show how potential problems may arise when off-band poles are present in the model. Second, the enforcement of uniform passivity throughout the entire frequency axis is performed via an iterative perturbation scheme on the purely imaginary eigenvalues of associated Hamiltonian matrices. A special formulation of this spectral perturbation using possibly large but sparse matrices allows the passivity compensation to be performed at a cost which scales only linearly with the order of the system. This formulation involves a restarted Arnoldi iteration combined with a complex frequency hopping algorithm for the selective computation of the imaginary eigenvalues to be perturbed. Some examples of interconnect models are used to illustrate the performance of the proposed technique
A bibliography on parallel and vector numerical algorithms
This is a bibliography of numerical methods. It also includes a number of other references on machine architecture, programming language, and other topics of interest to scientific computing. Certain conference proceedings and anthologies which have been published in book form are listed also
High-performance Parallel Solver for Integral Equations of Electromagnetics Based on Galerkin Method
A new parallel solver for the volumetric integral equations (IE) of
electrodynamics is presented. The solver is based on the Galerkin method which
ensures the convergent numerical solution. The main features include: (i) the
memory usage is 8 times lower, compared to analogous IE based algorithms,
without additional restriction on the background media; (ii) accurate and
stable method to compute matrix coefficients corresponding to the IE; (iii)
high degree of parallelism. The solver's computational efficiency is shown on a
problem of magnetotelluric sounding of the high conductivity contrast media. A
good agreement with the results obtained with the second order finite element
method is demonstrated. Due to effective approach to parallelization and
distributed data storage the program exhibits perfect scalability on different
hardware platforms.Comment: The main results of this paper were presented at IAMG 2015 conference
Frieberg, Germany. 28 pages, 11 figure
Hermite matrix in Lagrange basis for scaling static output feedback polynomial matrix inequalities
Using Hermite's formulation of polynomial stability conditions, static output
feedback (SOF) controller design can be formulated as a polynomial matrix
inequality (PMI), a (generally nonconvex) nonlinear semidefinite programming
problem that can be solved (locally) with PENNON, an implementation of a
penalty method. Typically, Hermite SOF PMI problems are badly scaled and
experiments reveal that this has a negative impact on the overall performance
of the solver. In this note we recall the algebraic interpretation of Hermite's
quadratic form as a particular Bezoutian and we use results on polynomial
interpolation to express the Hermite PMI in a Lagrange polynomial basis, as an
alternative to the conventional power basis. Numerical experiments on benchmark
problem instances show the substantial improvement brought by the approach, in
terms of problem scaling, number of iterations and convergence behavior of
PENNON
Passivity Enforcement via Perturbation of Hamiltonian Matrices
This paper presents a new technique for the passivity enforcement of linear time-invariant multiport systems in statespace form. This technique is based on a study of the spectral properties of related Hamiltonian matrices. The formulation is applicable in case the system input-output transfer function is in admittance, impedance, hybrid, or scattering form. A standard test for passivity is first performed by checking the existence of imaginary eigenvalues of the associated Hamiltonian matrix. In the presence of imaginary eigenvalues the system is not passive. In such a case, a new result based on first-order perturbation theory is presented for the precise characterization of the frequency bands where passivity violations occur. This characterization is then used for the design of an iterative perturbation scheme of the state matrices, aimed at the displacement of the imaginary eigenvalues of the Hamiltonian matrix. The result is an effective algorithm leading to the compensation of the passivity violations. This procedure is very efficient when the passivity violations are small, so that first-order perturbation is applicable. Several examples illustrate and validate the procedure
GPU-resident sparse direct linear solvers for alternating current optimal power flow analysis
Integrating renewable resources within the transmission grid at a wide scale poses significant challenges for economic dispatch as it requires analysis with more optimization parameters, constraints, and sources of uncertainty. This motivates the investigation of more efficient computational methods, especially those for solving the underlying linear systems, which typically take more than half of the overall computation time. In this paper, we present our work on sparse linear solvers that take advantage of hardware accelerators, such as graphical processing units (GPUs), and improve the overall performance when used within economic dispatch computations. We treat the problems as sparse, which allows for faster execution but also makes the implementation of numerical methods more challenging. We present the first GPU-native sparse direct solver that can execute on both AMD and NVIDIA GPUs. We demonstrate significant performance improvements when using high-performance linear solvers within alternating current optimal power flow (ACOPF) analysis. Furthermore, we demonstrate the feasibility of getting significant performance improvements by executing the entire computation on GPU-based hardware. Finally, we identify outstanding research issues and opportunities for even better utilization of heterogeneous systems, including those equipped with GPUs
- âŠ