In this paper, we introduce some new accelerated hybrid algorithms for solving a pseudomonotone equilibrium problem with a Lipschitz-type condition in a Hilbert space. The algorithms are constructed around the extragradient method, the inertial technique, the hybrid (or outer approximation) method, and the shrinking projection method. The algorithms are designed to work either with or without the prior

In this paper, we propose a novel overlapping domain decomposition method that can be applied to various problems in variational imaging such as total variation minimization. Most of recent domain decomposition methods for total variation minimization adopt the Fenchel–Rockafellar duality, whereas the proposed method is based on the primal formulation. Thus, the proposed method can be applied not only

We present an hp-version of the C0-continuous Petrov-Galerkin time stepping method for nonlinear second-order initial value problems. We derive a priori error bound in the H1-norm that is fully explicit with respect to the local time steps and the local approximation orders. Moreover, we prove that the hp-version of the C0-continuous Petrov-Galerkin time stepping method based on geometrically refined

In this work, we explore the advantages of end-to-end learning of multilayer maps offered by feedforward neural networks (FFNNs) for learning and predicting dynamics from transient flow data. While data-driven learning (and machine learning) in general depends on data quality and quantity relative to the underlying dynamics of the system, it is important for a given data-driven learning architecture

We investigate an inverse and ill-posed problem for the two-dimensional inhomogeneous heat equation in the presence of a general source term. The goal here consists of recovering not only the temperature distribution but also the thermal flux from the measure data. With the appearance of the general source term, this model gets far worse than its homogeneous counterpart. Based on an analysis of the

In this paper, we are concerned with multivariate Gegenbauer approximation of functions defined in the d-dimensional hypercube. Two new and sharper bounds for the coefficients of multivariate Gegenbauer expansion of analytic functions are presented based on two different extensions of the Bernstein ellipse. We then establish an explicit error bound for the multivariate Gegenbauer approximation associated

In this paper, we consider the matrices approximated in \({\mathscr{H}}^{2}\) 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 \({\mathscr{H}}^{2}\) matrix that allows to substitute direct \({\mathscr{H}}^{2}\) solution with the solution of the system with an equivalent sparse

One of the most popular time-frequency representations is certainly the Wigner distribution. Its quadratic nature is, however, at the origin of unwanted interferences or artefacts. The desire to suppress these artefacts is the reason why engineers, mathematicians and physicists have been looking for related time-frequency distributions, many of them being members of the Cohen class. Among these, the

This work is concerned with adaptive hybridizable discontinuous Galerkin methods of nonstationary convection diffusion problems. We address first the spatially semidiscrete case and then move to the fully discrete scheme by introducing a backward Euler discretization in time. More specifically, the computable a posteriori error estimator for the time-dependent problem is obtained by using the idea

We focus on a two-dimensional time-space diffusion equation with fractional derivatives in space. The use of Crank-Nicolson in time and finite differences in space leads to dense Toeplitz-like linear systems. Multigrid strategies that exploit such structure are particularly effective when the fractional orders are both close to 2. We seek to investigate how structure-based multigrid approaches can

In this paper, we consider numerical approximations for solving the Cahn-Hilliard phase-field model with the Flory-Huggins-de Gennes free energy for homopolymer blends. We develop an efficient, second-order accurate, and unconditionally energy stable scheme that combines the SAV approach with the stabilization technique, in which the H1 norm is split from the total free energy and two extra linear

We consider fully discrete schemes based on the scalar auxiliary variable (SAV) approach and stabilized SAV approach in time and the Fourier-spectral method in space for the phase field crystal (PFC) equation. Unconditionally, energy stability is established for both first- and second-order fully discrete schemes. In addition to the stability, we also provide a rigorous error estimate which shows that

The present work deals with rational model order reduction methods based on the single-point Least-Square (LS) Padé approximation techniques introduced in Bonizzoni et al. (ESAIM Math. Model. Numer. Anal., 52(4), 1261–1284 2018, Math. Comput. 89, 1229–1257 2020). Algorithmical aspects concerning the construction of rational LS-Padé approximants are described. In particular, we show that the computation

In computational geometry, different ways of space partitioning have been developed, including the Voronoi diagram of points and the power diagram of balls. In this article, a generalized Voronoi partition of overlapping d-dimensional balls, called the boundary-partition-based diagram, is proposed. The definition, properties, and applications of this diagram are presented. Compared to the power diagram

A high-order cut finite element method is formulated for solving the elastic wave equation. Both a single domain problem and an interface problem are treated. The boundary or interface is allowed to cut through the background mesh. To avoid problems with small cuts, stabilizing terms are added to the bilinear forms corresponding to the mass and stiffness matrix. The stabilizing terms penalize jumps

In this paper we introduce concepts from uncertainty quantification (UQ) and numerical analysis for the efficient evaluation of stochastic high dimensional Newton iterates. In particular, we develop complex analytic regularity theory of the solution with respect to the random variables. This justifies the application of sparse grids for the computation of statistical measures. Convergence rates are

In this article, we are concerned with the analysis on the numerical reconstruction of the spatial component in the source term of a time-fractional diffusion equation. This ill-posed problem is solved through a stabilized nonlinear minimization system by an appropriately selected Tikhonov regularization. The existence and the stability of the optimization system are demonstrated. The nonlinear optimization

Hankel-norm approximation is a model reduction method for linear time-invariant systems, which provides the best approximation in the Hankel semi-norm. In this paper, the computation of the optimal Hankel-norm approximation is generalized to the case of linear time-invariant continuous-time descriptor systems. A new algebraic characterization of all-pass descriptor systems is developed and used to

The Trefftz discontinuous Galerkin (TDG) method provides natural well-balanced (WB)and asymptotic preserving (AP) discretization, since exact solutions are used locally in the basis functions. However, one difficult point may be the construction of such solutions, which is a necessary first step in order to apply the TDG scheme. This work deals with the construction of solutions to Friedrichs systems

We establish a construction of hybrid multiscale systems on a bounded domain \({\Omega } \subset \mathbb {R}^{2}\) consisting of shearlets and boundary-adapted wavelets, which satisfy several properties advantageous for applications to imaging science and the numerical analysis of partial differential equations. More precisely, we construct hybrid shearlet-wavelet systems that form frames for the Sobolev

In 2012, Găvruţa introduced the notions of K-frame and of atomic system for a linear bounded operator K in a Hilbert space \({\mathcal{H}}\), in order to decompose its range \(\mathcal {R}(K)\) with a frame-like expansion. In this article, we revisit these concepts for an unbounded and densely defined operator \(A:\mathcal {D}(A)\to {\mathcal{H}}\) in two different ways. In one case, we consider a

Two fast L1 time-stepping methods, including the backward Euler and stabilized semi-implicit schemes, are suggested for the time-fractional Allen-Cahn equation with Caputo’s derivative. The time mesh is refined near the initial time to resolve the intrinsically initial singularity of solution, and unequal time steps are always incorporated into our approaches so that a adaptive time-stepping strategy

The study of non-convex penalties has recently received considerable attentions in sparse approximation. The existing non-convex penalties are proposed on the principle of seeking for a continuous alternative to the ℓ0-norm penalty. In this paper, we come up with a cubic spline penalty (CSP) which is also continuous but closer to ℓ0-norm penalty compared to the existing ones. As a result, it produces

Recently, the structured backward errors for the generalized saddle point problems with some different structures have been studied by some authors, but their results involve some Kronecker products, the vec-permutation matrices, and the orthogonal projection of a large block matrix which make them very expensive to compute when utilized for testing the stability of a practical algorithm or as an effective

A second-order surface reconstruction (SR) method for the shallow water equations with a discontinuous bottom topography and a Manning friction source term is presented. We redefine the water surface level at the cell interface by using the minimum difference between the bottom level and the original water surface level. The reconstructed water surface level is used to define the intermediate bottom

In general, it needs to take about nearly 10 grid points inside a wall boundary layer for low accuracy order methods to get satisfactory wall-normal gradient related results, such as friction coefficient and wall heat flux. If there exist extreme points inside the boundary layer, this situation becomes even worse. In this work, with the help of the analytic solution of the one-dimensional steady compressible

In this article, we propose a new method for low-rank completion of a large sparse matrix, subject to non-negativity constraint. As a challenging prototype of this problem, we have in mind the well-known Netflix problem. Our method is based on the derivation of a constrained gradient system and its numerical integration. The methods we propose are based on the constrained minimization of a functional

Quantifying the error that is induced by numerical approximation techniques is an important task in many fields of applied mathematics. Two characteristic properties of error bounds that are desirable are reliability and efficiency. In this article, we present an error estimation procedure for general nonlinear problems and, in particular, for parameter-dependent problems. With the presented auxiliary

Given the scalar sequence \(\{f_{m}\}^{\infty }_{m=0}\) that satisfies $$ f_{m} = \sum\limits_{i=1}^{k} {a_{i}}{\zeta}_{i}^{m},\quad m=0,1,\ldots, $$ where \(a_{i}, \zeta _{i}\in \mathbb {C}\) and ζi are distinct, the algorithm of Prony concerns the determination of the ai and the ζi from a finite number of the fm. This algorithm is also related to Padé approximants from the infinite power series \({\sum

This paper considers Lévy noise driven nonlinear stochastic Volterra integral equations with doubly weakly singular kernels, whose singular points include both s = 0 and s = t. The existence and uniqueness theorem of the true solution as well as the strong convergence rate of the Euler–Maruyama (EM) method are developed via establishing some fine estimates. Compared with the corresponding results by

This paper deals with Pythagorean hodograph curves of the spaces of algebraic trigonometric functions and trigonometric polynomials, \(\text {span}\left \{1,t,\left \{{\cos \limits } \left ({kt}\right ), {\sin \limits } \left ({kt}\right )\right \}_{k=1}^{m}\right \} \) and \( \text {span}\left \{1,\left \{{\cos \limits } \left ({kt}\right ), {\sin \limits } \left ({kt}\right )\right \}_{k=1}^{m}\right

The paper deals with numerical discretizations of separable nonlinear Hamiltonian systems with additive noise. For such problems, the expected value of the total energy, along the exact solution, drifts linearly with time. We present and analyze a time integrator having the same property for all times. Furthermore, strong and weak convergence of the numerical scheme along with efficient multilevel

In this paper, we propose a discontinuous finite volume element method to solve a phase field model for two immiscible incompressible fluids. In this finite volume element scheme, discontinuous linear finite element basis functions are used to approximate the velocity, phase function, and chemical potential while piecewise constants are used to approximate the pressure. This numerical method is efficient

In this paper, a sine pseudo-spectral-difference scheme that preserves the discrete mass and energy is presented and analyzed for the coupled Gross–Pitaevskii equations with Dirichlet boundary conditions in several spatial dimensions. The Crank–Nicolson finite difference method is employed for approximating the time derivative, and the second-order sine spectral differentiation matrix is deduced and

In this paper, we discuss a priori error estimates of two-grid mixed finite element methods for a class of nonlinear parabolic equations. The lowest order Raviart-Thomas mixed finite element and Crank-Nicolson scheme are used for the spatial and temporal discretization. First, we derive the optimal a priori error estimates for all variables. Second, we present a two-grid scheme and analyze its convergence

In this paper, we propose a linearly implicit Fourier pseudo-spectral scheme, which preserves the total mass and energy conservation laws for the damped nonlinear Schrödinger equation in three dimensions. With the aid of the semi-norm equivalence between the Fourier pseudo-spectral method and the finite difference method, an optimal L2-error estimate for the proposed method without any restriction

Let \(\mathbf d=(d_{j})_{j\in \mathbb {I}_{m}}\in \mathbb {N}^{m}\) be a finite sequence (of dimensions) and \(\alpha =(\alpha _{i})_{i\in \mathbb {I}_{n}}\) be a sequence of positive numbers (of weights), where \(\mathbb {I}_{k}=\{1,\ldots ,k\}\) for \(k\in \mathbb {N}\). We introduce the (α, d)-designs, i.e., m-tuples \({\Phi }=(\mathcal F_{j})_{j\in \mathbb {I}_{m}}\) such that \(\mathcal F_{j}=\{f_{ij}\}_{i\in

Abstract We study worst-case optimal approximation of positive linear functionals in reproducing kernel Hilbert spaces induced by increasingly flat Gaussian kernels. This provides a new perspective and some generalisations to the problem of interpolation with increasingly flat radial basis functions. When the evaluation points are fixed and unisolvent, we show that the worst-case optimal method converges

Abstract The eigenvalues and eigenfunctions of the Stokes operator have been the subject of intense analytical investigation and have applications in the study and simulation of the Navier–Stokes equations. As the Stokes operator is second order and has the divergence-free constraint, computing these eigenvalues and the corresponding eigenfunctions is a challenging task, particularly in complex geometries

Abstract We introduce stochastic variants of the classical Bernstein polynomials associated with a continuous function f, built up from a general triangular array of random variables. We discuss the uniform convergence in probability of the approximation process that they represent, providing at the same time rates of convergence. In the particular case in which the triangular array of random variables

Abstract We present the nonconforming virtual element method for the fourth-order singular perturbation problem. The virtual element proposed in this paper is a variant of the C0-continuous nonconforming virtual element presented in our previous work and allows to compute two different projection operators that are used for the construction of the discrete scheme. We show the optimal convergence in

Abstract Two fundamental difficulties are encountered in the numerical evaluation of time-dependent layer potentials. One is the quadratic cost of history dependence, which has been successfully addressed by splitting the potentials into two parts—a local part that contains the most recent contributions and a history part that contains the contributions from all earlier times. The history part is smooth

Abstract This article presents a new boundary integral approach for finding optimal shapes of peristaltic pumps that transport a viscous fluid. Formulas for computing the shape derivatives of the standard cost functionals and constraints are derived. They involve evaluating physical variables (traction, pressure, etc.) on the boundary only. By employing these formulas in conjunction with a boundary

Abstract A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium. A mixed finite element method is employed to approximate the pressure and the Darcy velocity, and a characteristic finite element method is used to approximate the concentration. Twice Newton iteration is applied on the fine grid to linearize the fully discrete problem using the coarse-grid

Abstract In this paper, we provide the rate of convergence for approximation of an eigenvalue problem describing vibrations of axisymmetric revolution elastic shells. The eigenvalue problem approximation is constructed by the method of finite elements. We study a system of differential equations with variable coefficients, some of which are non-integrable improperly near 0. To this aim, certain special

Abstract In this paper, we introduce the multi-region discontinuous Galerkin composite finite element method (MRDGCFEM) with hp-adaptivity for the discretization of second-order elliptic partial differential equations with discontinuous coefficients. This method allows for the approximation of problems posed on computational domains where the jumps in the diffusion coefficient form a micro-structure

Abstract We propose and analyze a family of decoupled and positivity-preserving discrete duality finite volume (DDFV) schemes for diffusion problems on general polygonal meshes. These schemes are further extensions of the scheme in Su et al. (J. Comput. Phys.372, 773–798, 2018) and share some benefits. First, the two sets of finite volume (FV) equations are constructed for the vertex-centered and cell-centered

Abstract Model-order reduction methods tackle the following general approximation problem: find an “easily-computable” but accurate approximation \(\hat {\boldsymbol {h}}\) of some target solution h⋆. In order to achieve this goal, standard methodologies combine two main ingredients: (i) a set of partial observations of h⋆ and (ii) some “simple” prior model on the set of target solutions. The most

In a similar fashion to estimates shown for Harmonic, Wachspress, and Sibson coordinates in [Gillette et al., AiCM, to appear], we prove interpolation error estimates for the mean value coordinates on convex polygons suitable for standard finite element analysis. Our analysis is based on providing a uniform bound on the gradient of the mean value functions for all convex polygons of diameter one satisfying

We prove the optimal convergence estimate for first order interpolants used in finite element methods based on three major approaches for generalizing barycentric interpolation functions to convex planar polygonal domains. The Wachspress approach explicitly constructs rational functions, the Sibson approach uses Voronoi diagrams on the vertices of the polygon to define the functions, and the Harmonic