Related papers: How to avoid mass matrix for linear hyperbolic pro…
In the paper, we introduce a matrix method to constructively determine spaces of polynomial solutions (in general, multiplied by exponentials) to a system of constant coefficient linear PDE's with polynomial (multiplied by exponentials)…
We present a sparse matrix permutation from graph theory that gives stable incomplete Lower-Upper (LU) preconditioners necessary for iterative solutions to the steady state density matrix for quantum optomechanical systems. This reordering…
We combine the parameterization method for invariant manifolds with the finite element method for elliptic PDEs,to obtain a new computational framework for high order approximation of invariant manifolds attached to unstable equilibrium…
Linear algebraic expressions are the essence of many computationally intensive problems, including scientific simulations and machine learning applications. However, translating high-level formulations of these expressions to efficient…
In this paper, we study diagonalizable hyperbolic systems in one space dimension. Based on a new gradient entropy estimate, we prove the global existence of a continuous solution, for large and nondecreasing initial data. Moreover, we show…
A version of the Dynamical Systems Method (DSM) for solving ill-conditioned linear algebraic systems is studied in this paper. An {\it a priori} and {\it a posteriori} stopping rules are justified. An iterative scheme is constructed for…
In this paper we will consider Rohde's general form of {1}-inverse of a matrix A. The necessary and sufficient condition for consistency of a linear system Ax=c will be represented. We will also be concerned with the minimal number of free…
Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier…
The effectiveness of the hyperbolic relaxation method for solving the Einstein constraint equations numerically is studied here on a variety of compact orientable three-manifolds. Convergent numerical solutions are found using this method…
An air pollution model is generally described by a system of PDEs on unbounded domain. Transformation of the independent variable is used to convert the problem for nonlinear air pollution on finite computational domain. We investigate the…
In this paper, we consider the matrices approximated in H2 format. The direct solution, as well as the preconditioning, of systems with such matrices is a challenging problem. We propose a non-extensive sparse factorization of the H2 matrix…
The parametrization theorem is derived in a flat nD pseudo-complex affine space. The pseudo-complex hyperbolic space accomodates n-number of uncompactified time-like extra dimensions with sugnature (s,r), where s and r are the numbers of…
We semiclassically derive the leading off-diagonal correction to the spectral form factor of quantum systems with a chaotic classical counterpart. To this end we present a phase space generalization of a recent approach for uniformly…
In this paper, we develop a class of robust numerical methods for solving dynamical systems with multiple time scales. We first represent the solution of a multiscale dynamical system as a transformation of a slowly varying solution. Then,…
Models with dominant advection always posed a difficult challenge for projection-based reduced order modelling. Many methodologies that have recently been proposed are based on the pre-processing of the full-order solutions to accelerate…
This paper deals with the index reduction problem for the class of quasi-regular DAE systems. It is shown that any of these systems can be transformed to a generically equivalent first order DAE system consisting of a single purely…
We consider systems of nonlinear magnetostatics and quasistatics that typically arise in the modeling and simulation of electric machines. The nonlinear problems, eventually obtained after time discretization, are usually solved by…
A hyperbolic system approach is proposed for robust computation of anisotropic diffusion equations that appear in quasineutral plasmas. Though the approach exhibits merits of high extensibility and accurate flux computation, the…
Snapshot matrices built from solutions to hyperbolic partial differential equations exhibit slow decay in singular values, whereas fast decay is crucial for the success of projection- based model reduction methods. To overcome this problem,…
When high-dimensional non-uniformly hyperbolic chaotic systems undergo dynamical perturbations, their long-time statistics are generally observed to respond differentiably with respect to the perturbation. Although important in…