相关论文: Tensor-structured PCG for finite difference solver…
A metric tensor for Riemann manifold Monte Carlo particularly suited for non-linear Bayesian hierarchical models is proposed. The metric tensor is built from symmetric positive semidefinite log-density gradient covariance (LGC) matrices,…
The conjugate gradient (CG) method, a standard and vital way of minimizing the energy of a variational state, is applied to solve several problems in Skyrmion physics. The single-Skyrmion profile optimizing the energy of a two-dimensional…
The paper considers grad-div stabilized equal-order finite elements (FE) methods for the linearized Navier-Stokes equations. A block triangular preconditioner for the resulting system of algebraic equations is proposed which is closely…
In this paper, we study the biharmonic equation with the Navier boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the 4th-order problem into a system of Poisson equations. Different…
Discontinuous Petrov-Galerkin (DPG) methods are new discontinuous Galerkin methods with interesting properties. In this article we consider a domain decomposition preconditioner for a DPG method for the Poisson problem.
Several preconditioned AOR methods have been proposed to solve system of linear equations $Ax=b$, where $A \in \mathbb{R}^{n \times n}$ is a unit Z-matrix. The aim of this paper is to give a comparison result for a class of preconditioners…
Domain decomposition (DD) methods are widely used as preconditioner techniques. Their effectiveness relies on the choice of a locally constructed coarse space. Thus far, this construction was mostly achieved using non-assembled matrices…
We derive a Kronecker product approximation for the micromagnetic long range interactions in a collocation framework by means of separable sinc quadrature. Evaluation of this operator for structured tensors (Canonical format, Tucker format,…
The convergence of the conjugate gradient method for solving large-scale and sparse linear equation systems depends on the spectral properties of the system matrix, which can be improved by preconditioning. In this paper, we develop a…
We propose Riemannian preconditioned algorithms for the tensor completion problem via tensor ring decomposition. A new Riemannian metric is developed on the product space of the mode-2 unfolding matrices of the core tensors in tensor ring…
The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled…
The projected gradient descent (PGD) method has shown to be effective in recovering compressed signals described in a data-driven way by a generative model, i.e., a generator which has learned the data distribution. Further reconstruction…
We present optimal preconditioners for a recently introduced hybridized discontinuous Galerkin finite element discretization of the Stokes equations. Typical of hybridized discontinuous Galerkin methods, the method has degrees-of-freedom…
We construct new Crouzeix-Raviart (CR) spaces of even degree $p$ that are spanned by basis functions mimicking those for the odd degree case. Compared to the standard CR gospel, the present construction allows for the use of nested bases of…
We present two novel classes of fully discrete energy-preserving algorithms for the sine-Gordon equation subject to Neumann boundary conditions. The cosine pseudo-spectral method is first used to develop structure-preserving spatial…
This work introduces the high-order Boris-SDC method for integrating the equations of motion for electrically charged particles in an electric and magnetic field. Boris-SDC relies on a combination of the Boris-integrator with spectral…
In this paper we propose an efficiently preconditioned Newton method for the computation of the leftmost eigenpairs of large and sparse symmetric positive definite matrices. A sequence of preconditioners based on the BFGS update formula is…
We revisit the Hierarchical Poincar\'e-Steklov (HPS) method in a preconditioned iterative setting for variable-coefficient Helmholtz problems with impedance boundary conditions. HPS is commonly presented as a direct solver based on nested…
The aim of this paper is to solve large-and-sparse linear Semidefinite Programs (SDPs) with low-rank solutions. We propose to use a preconditioned conjugate gradient method within second-order SDP algorithms and introduce a new efficient…
We present a direct Poisson solver for massively parallel simulations on three-dimensional Cartesian grids with non-uniform spacing. The method uses a tensor-based formulation in which the operator is diagonalized numerically along two…