Related papers: Multidimensional Generalized Riemann Problem Solve…
The goal in {\em reconfiguration problems} is to compute a {\em gradual transformation} between two feasible solutions of a problem such that all intermediate solutions are also feasible. In the {\em Matching Reconfiguration Problem} (MRP),…
Recovering a low-CP-rank tensor from noisy linear measurements is a central challenge in high-dimensional data analysis, with applications spanning tensor PCA, tensor regression, and beyond. We exploit the intrinsic geometry of rank-one…
We further research on the accelerated optimization phenomenon on Riemannian manifolds by introducing accelerated global first-order methods for the optimization of $L$-smooth and geodesically convex (g-convex) or $\mu$-strongly g-convex…
We present a five-wave Riemann solver for the equations of ideal relativistic magnetohydrodynamics. Our solver can be regarded as a relativistic extension of the five-wave HLLD Riemann solver initially developed by Miyoshi and Kusano for…
In this paper, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic…
We present a time-explicit discontinuous Galerkin (DG) solver for the time-domain acoustic wave equation on hybrid meshes containing vertex-mapped hexahedral, wedge, pyramidal and tetrahedral elements. Discretely energy-stable formulations…
We present a general and practical procedure to solve the general relativistic hydrodynamic equations by using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows.…
The Discontinuous Galerkin time-domain method is well suited for adaptive algorithms to solve the time-domain Maxwell's equations and depends on robust and economically computable drivers. Adaptive algorithms utilize local indicators to…
We develop a multi-dimensional hybrid discontinuous Galerkin method for multi-scale kinetic equations. This method is based on moment realizability matrices, a concept introduced by D. Levermore, W. Morokoff and B. Nadiga for one…
Equilibrium computation on Riemannian manifolds provides a unifying framework for numerous problems in machine learning and data analytics. One of the simplest yet most fundamental methods is Riemannian gradient descent (RGD). While its…
This work is concerned with constructing a second-order, invariant-domain preserving approximation of the compressible multi-species Euler equations where each species is modeled by an ideal gas equation of state. We give the full solution…
In the present work, a high order finite element type residual distribution scheme is designed in the framework of multidimensional compressible Euler equations of gas dynamics. The strengths of the proposed approximation rely on the…
This paper develops and analyzes an interior penalty discontinuous Galerkin (IPDG) method using piecewise linear polynomials for the indefinite time harmonic Maxwell equations with the impedance boundary condition in the three dimensional…
In this paper, error analysis is established for Runge-Kutta discontinuous Galerkin (RKDG) methods to solve the Vlasov-Maxwell system. This nonlinear hyperbolic system describes the time evolution of collisionless plasma particles of a…
Markov decision processes (MDPs) with large number of states are of high practical interest. However, conventional algorithms to solve MDP are computationally infeasible in this scenario. Approximate dynamic programming (ADP) methods tackle…
This paper proposes a Riemannian Multiobjective Proximal Gradient Method (RMPGM) for composite optimization problems on manifolds. Unlike scalarization-based approaches, the proposed framework directly handles vector-valued objectives and…
The present paper is devoted to existence results for time-periodic solutions of generalized nonlinear wave equations in a closed Riemannian manifold M. Our main focus lies on the doubly degenerate setting where the associated generalized…
Low-rank tensor completion aims to recover a tensor from partially observed entries, and it is widely applicable in fields such as quantum computing and image processing. Due to the significant advantages of the tensor train (TT) format in…
The Riemann problem is fundamental in the computational modeling of hyperbolic partial differential equations, enabling the development of stable and accurate upwind schemes. While exact solvers provide robust upwinding fluxes, their high…
With the advance of supercomputers we can now afford simulations with very large ranges of scales. In astrophysical applications, e.g. simulating Solar, stellar and planetary atmospheres, interstellar medium, etc; physical quantities, like…