数值分析
Randomized sketching is currently introduced into every area of numerical linear algebra. In Krylov subspace methods, it allows runtime savings at the cost of small accuracy reductions. This work offers a different view on sketching in…
We propose a computational framework for computing low-rank approximations to the ensemble of solutions of a parametrized system of the form $A(\xi)x(\xi)+g(x(\xi))=b(\xi)$ for multiple parameter values. The central idea is to reinterpret…
Spherical centroidal Voronoi tessellations (SCVTs), currently used in numerical weather forecasting models such as the Model for Prediction Across Scales (MPAS), are a type of spherical grid that is highly flexible, allowing the…
We present and analyze a discontinuous Galerkin method for the numerical solution of a class of second-order linear mixed-type partial differential equations, i.e. equations that change their nature from elliptic to hyperbolic through the…
The paper studies a time-nonlocal multiphysics finite element method with Crank-Nicolson scheme for poroelasticity model with secondary consolidation. For the case where the physical parameters $\lambda,\lambda^*$ and $c_0$ are all finite…
Shape optimisation of thin-shell structures requires a flexible, differentiable geometric representation suitable for gradient-based optimisation. We propose a neural parametric representation (NRep) for the shell mid-surface based on a…
We develop a numerical scheme for solving the advection equation of $\mathbb{S}^2$-valued functions of real variables, which models the time-evolution of a $\mathbb{S}^2$-valued mapping on the real line by a known velocity field. The idea…
In this paper we analyze Nitsche's method for the stationary Boussinesq system with Navier's slip and a nonlinear boundary condition. Our analysis of the formulation establishes the robustness of a finite elements scheme in arbitrarily…
In this work we study the theoretical Lipschitz stability and propose a low-rank-assisted numerical method for the inverse medium scattering beyond the Born region. The proposed low-rank structure is based on the disk prolate spheroidal…
In this paper, we propose two monolithic fully discrete finite element methods for fluid-structure interaction (FSI) based on a novel Piola-type Arbitrary Lagrangian-Eulerian (ALE) mapping. For the temporal discretization, we apply the…
We propose a novel fractal based interpolation scheme termed Rational Cubic Trigonometric Zipper Fractal Interpolation Functions (RCTZFIFs) designed to model and preserve the inherent geometric property, positivity, in given datasets. The…
We construct a structure preserving non-conforming finite element approximation scheme for the bi-harmonic wave maps into spheres equation. It satisfies a discrete energy law and preserves the non-convex sphere constraint of the continuous…
This paper is concerned with the construction of discrete counterparts of the Laplace-Beltrami operator on Riemannian manifolds that can be effectively used in the numerical solution of partial differential equations. Since existing…
This paper is concerned with the inverse problem of retrieving the initial value of a time-fractional fourth order parabolic equation from source and final time observation. The considered problem is an {\it ill-posed problem.} We obtain…
The computational cost of concurrent multiscale finite element methods is dominated by the repeated solution of microscopic representative volume element (RVE) problems at macroscopic quadrature points. In this work, we introduce a…
We propose a neural-enhanced weak Galerkin (WG) finite element method for second-order elliptic problems with low-regularity solutions. The method augments the classical WG approximation space with neural network functions constructed via a…
In this paper, we study the computational question of whether the Steklov spectrum of smooth simply connected planar domains can be approximated accurately by a boundary-only formulation based on harmonic conjugation. For the unit disk, the…
Electroencephalography (EEG) source imaging aims to reconstruct the spatial distribution of neural activity within the brain from non-invasive scalp measurements. This inverse problem is severely ill-posed due to the low spatial resolution…
The generalized Gearhart-Koshy acceleration is a recent exact affine search technique designed for the method of cyclic projections onto hyperplanes, i.e., the Kaczmarz method. However, its convergence properties, particularly the linear…
In this paper, for solving inconsistent matrix equations we propose a dual-space residual-based randomized extended Kaczmarz method and its version with Nesterov momentum. Without the full column rank assumptions on coefficient matrices, we…