数值分析

Standard gradient-based iteration algorithms for optimization, such as gradient descent and its various proximal-based extensions to nonsmooth problems, are known to converge slowly for ill-conditioned problems, sometimes requiring many…

This paper presents a class of efficient manifold optimization algorithms for computing the ground state solutions of a semilinear elliptic system, which are unstable saddle points of the variational functional. Variational arguments show…

In this paper, we investigate the asymptotic distribution of the normalized error for the Mittag--Leffler Euler (MLE) method applied to a class of multidimensional fractional stochastic differential equations. These equations are…

Neural networks have demonstrated significant potential in solving partial differential equations (PDEs). While global approaches such as Physics-Informed Neural Networks (PINNs) offer promising capabilities, they often lack inherent…

We prove strong convergence of an upwind-type finite volume method to a weak solution of the Navier-Stokes-Fourier system with the Dirichlet boundary conditions. The limit solution satisfies a weak form of the mass and momentum equations,…

BoostConv has been introduced in earlier works as an effective acceleration technique for nonlinear iterative processes and has been successfully employed in a variety of applications to enhance convergence rates or to compute unstable…

In this paper, we develop a Localized Orthogonal Decomposition (LOD) method for the two-dimensional time-dependent nonlinear Schr\"{o}dinger equation with a wave operator. We prove that our method preserves conservation laws and admits a…

We introduce a novel residual-based a posteriori error estimator for the conforming $C^1$ Virtual Element Method (VEM) applied to the buckling eigenvalue problem, incorporating nonlinear plane stress effects in both two and three…

In this paper, we apply acceleration to the inverse-free preconditioned Krylov subspace method introduced by Golub and Ye, which solves the symmetric generalized eigenvalue problem for the algebraically smallest eigenvalue. As the method is…

We propose a generalized version of the bisection method where the cutting point between the two subintervals is chosen at random following an arbitrary distribution. We compute expected convergence rates with respect to any arbitrary a…

We present a mimetic finite-difference approach for solving Maxwell's equations in one and two spatial dimensions. After introducing the governing equations and the classical Finite-Difference Time-Domain (FDTD) method, we describe mimetic…

We extend the Weak Adversarial Neural Pushforward Method (WANPM) to the McKean--Vlasov mean-field Fokker--Planck equation, covering both the stationary and time-dependent cases. The key observation is that the mean-field nonlinearity -- an…

Two-stage orthogonalization is essential in numerical algorithms such as Krylov subspace methods. For this task we need to orthogonalize a matrix $A$ against another matrix $V$ with orthonormal columns. A common approach is to employ the…

We present FastLSQ, a framework for PDE solving and inverse problems built on trigonometric random Fourier features with exact analytical derivatives. Trigonometric features admit closed-form derivatives of any order in $\mathcal{O}(1)$,…

In this paper, we propose and analyze a mixed virtual element method for the approximation of the eigenvalues and eigenfunctions of the two-dimensional elasticity eigenvalue problem. Under standard assumptions on polygonal meshes, we prove…

In this work, we consider wave propagation in materials characterized by nonlinear properties or damage. To accelerate the simulations of the resulting high-dimensional problems, we apply model order reduction methods. Depending on the…

We consider the numerical solution of high-frequency scattering problems modeled by the Helmholtz equation with a bounded obstacle. Although the analysis of this problem dates back at least 50 years, over the past decade or so, tools and…

Tomographic techniques are vital in modern medicine, allowing doctors to observe patients' interior features. Individual steps in the measurement process are modeled by `single projection operators' $p$. These are line integral operators…

This work extends the discrete compactness results of Walkington (SIAM J. Numer. Anal., 47(6):4680--4710, 2010) for high-order discontinuous Galerkin time discretizations of parabolic problems to more general function space settings. In…

This paper describes a fast direct boundary element method for elastodynamic transmission problems in two dimensions, which can be used for analyzing elastic wave scattering by an inclusion. We develop an efficient solver based on a…