284 research outputs found
On the numerical solution of a class of systems of linear matrix equations
We consider the solution of systems of linear matrix equations in two or three unknown matrices. For dense problems we derive algorithms that determine the numerical solution by only involving matrices of the same size as those in the original problem, thus requiring low computational resources. For large and structured systems we show how the problem properties can be exploited to design effective algorithms with low memory and operation requirements. Numerical experiments illustrate the performance of the new methods
Incomplete Augmented Lagrangian Preconditioner for Steady Incompressible Navier-Stokes Equations
An incomplete augmented Lagrangian preconditioner, for the steady incompressible Navier-Stokes equations discretized by stable finite elements, is proposed. The eigenvalues of the preconditioned matrix are analyzed. Numerical experiments show that the incomplete augmented Lagrangian-based preconditioner proposed is very robust and performs quite well by the Picard linearization or the Newton linearization over a wide range of values of the viscosity on both uniform and stretched grids
Graph coarsening: From scientific computing to machine learning
The general method of graph coarsening or graph reduction has been a
remarkably useful and ubiquitous tool in scientific computing and it is now
just starting to have a similar impact in machine learning. The goal of this
paper is to take a broad look into coarsening techniques that have been
successfully deployed in scientific computing and see how similar principles
are finding their way in more recent applications related to machine learning.
In scientific computing, coarsening plays a central role in algebraic multigrid
methods as well as the related class of multilevel incomplete LU
factorizations. In machine learning, graph coarsening goes under various names,
e.g., graph downsampling or graph reduction. Its goal in most cases is to
replace some original graph by one which has fewer nodes, but whose structure
and characteristics are similar to those of the original graph. As will be
seen, a common strategy in these methods is to rely on spectral properties to
define the coarse graph
Interior-point methods for PDE-constrained optimization
In applied sciences PDEs model an extensive variety of phenomena. Typically the final goal of simulations is a system which is optimal in a certain sense. For instance optimal control problems identify a control to steer a system towards a desired state. Inverse problems seek PDE parameters which are most consistent with measurements. In these optimization problems PDEs appear as equality constraints. PDE-constrained optimization problems are large-scale and often nonconvex. Their numerical solution leads to large ill-conditioned linear systems. In many practical problems inequality constraints implement technical limitations or prior knowledge.
In this thesis interior-point (IP) methods are considered to solve nonconvex large-scale PDE-constrained optimization problems with inequality constraints. To cope with enormous fill-in of direct linear solvers, inexact search directions are allowed in an inexact interior-point (IIP) method. This thesis builds upon the IIP method proposed in [Curtis, Schenk, Wächter, SIAM Journal on Scientific Computing, 2010]. SMART tests cope with the lack of inertia information to control Hessian modification and also specify termination tests for the iterative linear solver.
The original IIP method needs to solve two sparse large-scale linear systems in each optimization step. This is improved to only a single linear system solution in most optimization steps. Within this improved IIP framework, two iterative linear solvers are evaluated: A general purpose algebraic multilevel incomplete L D L^T preconditioned SQMR method is applied to PDE-constrained optimization problems for optimal server room cooling in three space dimensions and to compute an ambient temperature for optimal cooling. The results show robustness and efficiency of the IIP method when compared with the exact IP method.
These advantages are even more evident for a reduced-space preconditioned (RSP) GMRES solver which takes advantage of the linear system's structure. This RSP-IIP method is studied on the basis of distributed and boundary control problems originating from superconductivity and from two-dimensional and three-dimensional parameter estimation problems in groundwater modeling. The numerical results exhibit the improved efficiency especially for multiple PDE constraints.
An inverse medium problem for the Helmholtz equation with pointwise box constraints is solved by IP methods. The ill-posedness of the problem is explored numerically and different regularization strategies are compared. The impact of box constraints and the importance of Hessian modification on the optimization algorithm is demonstrated. A real world seismic imaging problem is solved successfully by the RSP-IIP method
BATTPOWER Toolbox: Memory-Efficient and High-Performance Multi-Period AC Optimal Power Flow Solver
With the introduction of massive renewable energy sources and storage
devices, the traditional process of grid operation must be improved in order to
be safe, reliable, fast responsive and cost efficient, and in this regard power
flow solvers are indispensable. In this paper, we introduce an Interior
Point-based (IP) Multi-Period AC Optimal Power Flow (MPOPF) solver for the
integration of Stationary Energy Storage Systems (SESS) and Electric Vehicles
(EV). The primary methodology is based on: 1) analytic and exact calculation of
partial differential equations of the Lagrangian sub-problem, and 2) exploiting
the sparse structure and pattern of the coefficient matrix of Newton-Raphson
approach in the IP algorithm. Extensive results of the application of proposed
methods on several benchmark test systems are presented and elaborated, where
the advantages and disadvantages of different existing algorithms for the
solution of MPOPF, from the standpoint of computational efficiency, are brought
forward. We compare the computational performance of the proposed
Schur-Complement algorithm with a direct sparse LU solver. The comparison is
performed for two different applicational purposes: SESS and EV. The results
suggest the substantial computational performance of Schur-Complement algorithm
in comparison with that of a direct LU solver when the number of storage
devices and optimisation horizon increase for both cases of SESS and EV. Also,
some situations where computational performance is inferior are discussed.Comment: 24 pages, 15 figures, Accepted for publication in IEEE Transactions
on Power System
- …