Preconditioning optimal in-domain control of navier-stokes equation using multilevel sequentially semiseparable matrix computations

In this manuscript, we study preconditioning techniques for optimal in-domain control of the Navier-Stokes equation, where the control only acts on a few parts of the domain. Optimization and linearization of the optimal in-domain control problem results in a generalized linear saddle-point system. The Schur complement for the generalized saddle...

Induced Dimension Reduction method for solving linear matrix equations

This paper discusses the solution of large-scale linear matrix equations using the Induced Dimension reduction method (IDR(s)). IDR(s) was originally presented to solve system of linear equations, and is based on the IDR(s) theorem. We generalize the IDR(s) theorem to solve linear problems in any finite-dimensional space. This generalization...

Preconditioning Optimal In-Domain Control of Navier-Stokes Equation Using Multilevel Sequentially Semiseparable Matrix Computations

In this manuscript, we study preconditioning techniques for optimal in-domain control of the Navier-Stokes equation, where the control only acts on a few parts of the domain. Optimization and linearization of the optimal in-domain control problem results in a generalized linear saddle-point system. The Schur complement for the generalized saddle...

Efficient Preconditioners for PDE-Constrained Optimization Problems with a Multi-level Sequentially Semi-Separable Matrix Structure

In this paper, we consider preconditioning for PDE-constrained optimization problems. The underlying problems yield a linear saddle-point system. We study a class of preconditioners based on multilevel sequentially semiseparable (MSSS) matrix computations. The novel global preconditioner is to make use of the global structure of the saddle-point...

Evaluation of Multilevel Sequentially Semiseparable Preconditioners on CFD Benchmark Problems Using IFISS

This paper studies a new preconditioning technique for sparse systems arising from discretized partial differential equations (PDEs) in computational fluid dynamics (CFD), which exploit the multilevel sequentially semiseparable (MSSS) structure of the system matrix. MSSS matrix computations give a data-sparse way to approximate the LU...

A class of efficient preconditioners with multilevel sequentially semiseparable matrix structure

This paper presents a class of preconditioners for sparse systems arising from discretized partial differential equations (PDEs). In this class of preconditioners, we exploit the multilevel sequentially semiseparable (MSSS) structure of the system matrix. The off-diagonal blocks of MSSS matrices are of low-rank, which enables fast computations...

On the use of rigid body modes in the deflated preconditioned conjugate gradient method

Large discontinuities in material properties, such as those encountered in composite materials, lead to ill-conditioned systems of linear equations. These discontinuities give rise to small eigenvalues that may negatively affect the convergence of iterative solution methods such as the preconditioned conjugate gradient method. This paper...

Comparison of the deflated preconditioned conjugate gradient method and algebraic multigrid for composite materials

Many applications in computational science and engineering concern composite materials, which are characterized by large discontinuities in the material properties. Such applications require fine-scale finite-element meshes, which lead to large linear systems that are challenging to solve with current direct and iterative solutions algorithms....

Comparison of the deflated preconditioned conjugate gradient method and parallel direct solver for composite materials

The demand for large FE meshes increases as parallel computing becomes the standard in FE simulations. Direct and iterative solution methods are used to solve the resulting linear systems. Many applications concern composite materials, which are characterized by large discontinuities in the material properties. An example of such a material is...

On the use of rigid body modes in the deflated preconditioned conjugate gradient method

Large discontinuities in material properties, such as encountered in composite materials, lead to ill-conditioned systems of linear equations. These discontinuities give rise to small eigenvalues that may negatively affect the convergence of iterative solution methods such as the Preconditioned Conjugate Gradient (PCG) method. This paper...

The Deflated Relaxed Incomplete Cholesky CG method for use in a real-time ship simulator

Ship simulators are used for training purposes and therefore have to calculate realistic wave patterns around the moving ship in real time. We consider a wave model that is based on the variational Boussinesq formulation, which results in a set of partial differential equations. Discretization of these equations gives a large system of linear...

Convergence bounds for preconditioned GMRES using element-by-element estimates of the field of values

By combining element-by-element estimates for the field of values of a preconditioned matrix with GMRES-convergence estimates it is possible to derive an easily computable upper bound on the GMRES-residual norm. This method can be applied to general finite element systems, but the preconditioner has to be Hermitian and positive definite. The...

Convergence bounds for preconditioned GMRES using element-by-element estimates of the field of values

By combining element-by-element estimates for the field of values of a preconditioned matrix with GMRES-convergence estimates it is possible to derive an easily computable upper bound on the GMRES-residual norm. This method can be applied to general finite element systems, but the preconditioner has to be Hermitian and positive definite. The...