23,776 research outputs found
MATEX: A Distributed Framework for Transient Simulation of Power Distribution Networks
We proposed MATEX, a distributed framework for transient simulation of power
distribution networks (PDNs). MATEX utilizes matrix exponential kernel with
Krylov subspace approximations to solve differential equations of linear
circuit. First, the whole simulation task is divided into subtasks based on
decompositions of current sources, in order to reduce the computational
overheads. Then these subtasks are distributed to different computing nodes and
processed in parallel. Within each node, after the matrix factorization at the
beginning of simulation, the adaptive time stepping solver is performed without
extra matrix re-factorizations. MATEX overcomes the stiff-ness hinder of
previous matrix exponential-based circuit simulator by rational Krylov subspace
method, which leads to larger step sizes with smaller dimensions of Krylov
subspace bases and highly accelerates the whole computation. MATEX outperforms
both traditional fixed and adaptive time stepping methods, e.g., achieving
around 13X over the trapezoidal framework with fixed time step for the IBM
power grid benchmarks.Comment: ACM/IEEE DAC 2014. arXiv admin note: substantial text overlap with
arXiv:1505.0669
A Parallel SSOR Preconditioner for Lattice QCD
A parallelizable SSOR preconditioning scheme for Krylov subspace iterative
solvers in lattice QCD applications involving Wilson fermions is presented. In
actual Hybrid Monte Carlo and quark propagator calculations it helps to reduce
the number of iterations by a factor of 2 compared to conventional odd-even
preconditioning. This corresponds to a gain in cpu-time of 30\% - 70\% over
odd-even preconditioning.Comment: Talk presented at LATTICE96(algorithms), 3 pages, LaTeX file, 3
epsf-files include
Parallel solution of high-order numerical schemes for solving incompressible flows
A new parallel numerical scheme for solving incompressible steady-state flows is presented. The algorithm uses a finite-difference approach to solving the Navier-Stokes equations. The algorithms are scalable and expandable. They may be used with only two processors or with as many processors as are available. The code is general and expandable. Any size grid may be used. Four processors of the NASA LeRC Hypercluster were used to solve for steady-state flow in a driven square cavity. The Hypercluster was configured in a distributed-memory, hypercube-like architecture. By using a 50-by-50 finite-difference solution grid, an efficiency of 74 percent (a speedup of 2.96) was obtained
Solution of partial differential equations on vector and parallel computers
The present status of numerical methods for partial differential equations on vector and parallel computers was reviewed. The relevant aspects of these computers are discussed and a brief review of their development is included, with particular attention paid to those characteristics that influence algorithm selection. Both direct and iterative methods are given for elliptic equations as well as explicit and implicit methods for initial boundary value problems. The intent is to point out attractive methods as well as areas where this class of computer architecture cannot be fully utilized because of either hardware restrictions or the lack of adequate algorithms. Application areas utilizing these computers are briefly discussed
Automated Netlist Generation for 3D Electrothermal and Electromagnetic Field Problems
We present a method for the automatic generation of netlists describing
general three-dimensional electrothermal and electromagnetic field problems.
Using a pair of structured orthogonal grids as spatial discretisation, a
one-to-one correspondence between grid objects and circuit elements is obtained
by employing the finite integration technique. The resulting circuit can then
be solved with any standard available circuit simulator, alleviating the need
for the implementation of a custom time integrator. Additionally, the approach
straightforwardly allows for field-circuit coupling simulations by
appropriately stamping the circuit description of lumped devices. As the
computational domain in wave propagation problems must be finite, stamps
representing absorbing boundary conditions are developed as well.
Representative numerical examples are used to validate the approach. The
results obtained by circuit simulation on the generated netlists are compared
with appropriate reference solutions.Comment: This is a pre-print of an article published in the Journal of
Computational Electronics. The final authenticated version is available
online at: https://dx.doi.org/10.1007/s10825-019-01368-6. All numerical
results can be reproduced by the Matlab code openly available at
https://github.com/tc88/ANTHE
- …