Mathematics
Physics-informed neural networks (PINNs) formulate the solution of partial differential equations as residual minimization problems over neural network parameterizations. Although highly flexible, optimization of PINNs using modern variants…
Given a matrix $A$, a matrix nearness problem seeks an $X$ that most closely approximates $A$ in the sense of minimizing $\lVert A - X\rVert$ under a variety of constraints on $X$. A generalized matrix nearness problem seeks the same but…
We consider the numerical solution of partial differential equations with coefficients that are strongly heterogeneous in space. We provide an overview of higher-order localized orthogonal decomposition (LOD) methods for the elliptic…
A novel mixed spectral-Galerkin method based on generalized ball polynomials is proposed for solving the biharmonic equation on a unit ball. By introducing an auxiliary variable to decouple the biharmonic equation into a system of…
Developing high-order numerical schemes for two-phase flow in porous media that preserve key physical properties remains a significant challenge in numerical analysis. In this article, we propose a general framework to construct fully…
We prove undecidability for every positive relevant logic extending the system axiomatized by hypothetical syllogism, prefixing, and suffixing and contained in the logic of the semilattice frame $(P_{\mathrm{fin}}(\mathbb{N}), \cup,…
We present a high-order implicit-explicit discontinuous Galerkin (IMEX-DG) solver for the compressible Euler equations to account for rotational effects within a fully compressible atmospheric framework. Time integration follows a…
In this paper, we present a quantum implicit-explicit (IMEX) scheme for multiscale ordinary and partial differential equations whose discretization parameters are independent of the scaling parameter $\varepsilon$. A key ingredient of our…
We propose a high-order numerical methodology for computing the ground state and time evolution of the two-dimensional Gross-Pitaevskii equation with harmonic trapping potential. The ground state is obtained by combining normalized gradient…
In this short note we give a negative answer to the following open question: \emph{Let $X$ be a $\sigma$-compact paratopological group. Does there exist a continuous isomorphism of $X$ onto a topological group $G$?} Specifically, we…
This work presents a weighted quadrature (WQ) method to fast assemble Galerkin matrices based on unstructured spline surfaces. The method is developed upon a particular variant of unstructured splines, namely the bicubic analysis-suitable…
In this paper, we further investigate and refine the subspace-constrained preconditioning technique to enhance the theoretical and numerical convergence properties of randomized iterative methods for solving linear systems. In particular,…
This paper aims to employ the weak Galerkin method to solve a class of nonlinear eigenvalue problems. We proved the weak Galerkin scheme produces lower bound for the energy. Moreover, by the post-processing technique, we obtain lower bound…
This paper introduces a Laws of Form version of the Quaternions. We call this the Q-Calculus, a 16-valued extension of Laws of Form (LoF) which is closely related to the BF Calculus (where we have a single square root of the mark) and the…
This paper introduces a multifidelity formulation that reduces the computational cost of the proper orthogonal decomposition (POD) of a high-fidelity model by leveraging data from cheaper, lower-fidelity models. POD is a prevalent technique…
Centered finite-difference discretizations of convection--diffusion equations may oscillate when convection dominates at the mesh scale. For homogeneous Dirichlet problems with constant coefficients on uniform Cartesian grids, we derive…
The simulation of extreme Mach astrophysical flows is traditionally viewed through the lens of deterministic positivity-preserving schemes. However, due to Kelvin--Helmholtz instabilities and shock anomalies, the multi-dimensional Euler…
We introduce a new class of fractional backward orthogonal functions designed for the spectral approximation of weakly singular adjoint Volterra integral equations. These basis functions generate an approximation space that naturally…
Conventional mode tracking operates in the dark: it traces dispersion branches on the non-Hermitian eigenvalue manifold using only local continuity, unaware of the global Riemann-sheet topology. When exceptional points (EPs) lie close to…
Conservation laws are conventionally discretized through floating-point flux evaluation, with invariants obtained by cancellation of approximate interface contributions and admissible weak solutions selected by reconstruction and Riemann…