Related papers: Evaluating the Evans function: Order reduction in …
A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Stokes system with Flory-Huggins energy functional. A convex splitting is applied to the chemical potential, which in turns leads to the implicit treatment…
We consider numerical approximation to the solution of non-autonomous evolution equations. The order of convergence of the simplest possible Magnus method will be investigated.
Monitoring the integrity of elastic structures using ultrasonic waves requires the efficient identification of material parameters from measured surface displacements. The displacement field is governed by Cauchy's equation of motion, i.e.,…
The monotonicity of the Mittag-Leffler function $E_{\alpha}$ with respect to the parameter $\alpha$ is investigated, via some convex ordering properties for related random variables. In particular, it is shown that the mapping…
In applications involving ordinal predictors, common approaches to reduce dimensionality are either extensions of unsupervised techniques such as principal component analysis, or variable selection procedures that rely on modeling the…
This paper introduces Magnus-based methods for solving stochastic delay-differential equations (SDDEs). We construct Magnus--Euler--Maruyama (MEM) and Magnus--Milstein (MM) schemes by combining stochastic Magnus integrators with Taylor…
We study the Langevin equation with stationary-increment Gaussian noise. We show the strong consistency and the asymptotic normality with Berry--Esseen bound of the so-called alternative estimator of the mean reversion parameter. The…
In this paper, we study the affine phase retrieval problem, which aims to recover signals from the magnitudes of affine measurements. We develop second-order optimization methods based on Newton and Gauss-Newton iterations and establish…
The ranking problem is to order a collection of units by some unobserved parameter, based on observations from the associated distribution. This problem arises naturally in a number of contexts, such as business, where we may want to rank…
This paper concerns a spectral estimation problem in which we want to find a spectral density function that is consistent with estimated second-order statistics. It is an inverse problem admitting multiple solutions, and selection of a…
We investigate the (functional) convex order of for various continuous martingale processes, either with respect to their diffusions coefficients for L\'evy-driven SDEs or their integrands for stochastic integrals. Main results are bordered…
Electrodynamics in curved spacetime can be studied in the Eastwood--Singer gauge, which has the advantage of respecting the invariance under conformal rescalings of the Maxwell equations. Such a construction is here studied in Einstein…
The Gaussian sequence model is a canonical model in nonparametric estimation. In this study, we introduce a multivariate version of the Gaussian sequence model and investigate adaptive estimation over the multivariate Sobolev ellipsoids,…
Welch's method provides an estimator of the power spectral density that is statistically consistent. This is achieved by averaging over periodograms calculated from overlapping segments of a time series. For a finite length time series,…
In a recent paper by Kamrani et al. (2024), exponential Euler method for stiff stochastic differential equations with additive fractional Brownian noise was discussed, and the convergence order close to the Hurst parameter H was proved.…
The Bayesian methods for linear inverse problems is studied using hierarchical Gaussian models. The problems are considered with different discretizations, and we analyze the phenomena which appear when the discretization becomes finer. A…
To address the magnetization dynamics in ferromagnetic materials described by the Landau-Lifshitz-Gilbert equation under large damping parameters, a third-order accurate numerical scheme is developed by building upon a second-order method…
We develop two new stochastic Gauss-Newton algorithms for solving a class of non-convex stochastic compositional optimization problems frequently arising in practice. We consider both the expectation and finite-sum settings under standard…
Resonance based numerical schemes are those in which cancellations in the oscillatory components of the equation are taken advantage of in order to reduce the regularity required of the initial data to achieve a particular order of error…
This is the second paper in a series describing a numerical implementation of the conformal Einstein equation. This paper deals with the technical details of the numerical code used to perform numerical time evolutions from a "minimal" set…