We use a positive parameter to develop a dierentiable perturbed reconstruction model to solve the L1–L2 magnetic resonance image (MRI) reconstruction problem. We use Bregman iterative formulation to solve the dierentiable perturbed L1–L2 model, and lagged diusivity xed-point iteration to solve the minimization problem in the Bregman iteration. Two Fourier transforms and an inverse Fourier transform are used to accelerate L1–L2 MRI reconstruction. Real MR images are used to test the method in numerical experiments. The results demonstrate that the proposed method is very ecient for L1–L2MRI reconstruction.

We construct a local Lax-Friedrichs type positivity-preserving flux for compressible Navier-Stokes equations, which can be easily extended to high dimensions for generic forms of equations of state, shear stress tensor and heat flux. With this positivity-preserving flux, any finite volume type schemes including discontinuous Galerkin (DG) schemes with strong stability preserving Runge-Kutta time discretizations satisfy a weak positivity property. With a simple and effi- cient positivity-preserving limiter, high order explicit Runge-Kutta DG schemes are rendered preserving the positivity of density and internal energy without los- ing local conservation or high order accuracy. Numerical tests suggest that the positivity-preserving flux and the positivity-preserving limiter do not induce excessive artificial viscosity, and the high order positivity-preserving DG schemes without other limiters can produce satisfying non-oscillatory solutions when the nonlinear diffusion in compressible Navier-Stokes equations is accurately resolved.

In this paper we study the stochastic Galerkin approximation for the linear transport equation with random inputs and diffusive scal- ing. We first establish uniform (in the Knudsen number) stability results in the random space for the transport equation with uncertain scattering coefficients, and then prove the uniform spectral conver- gence (and consequently the sharp stochastic Asymptotic-Preserving property) of the stochastic Galerkin method. A micro-macro decom- position based fully discrete scheme is adopted for the problem and proved to have a uniform stability. Numerical experiments are conducted to demonstrate the stability and asymptotic properties of the method.

We establish theoretical recovery guarantees of a family of Riemannian optimization algorithms for low rank matrix recovery, which is about recovering an mimes n rank mimes n matrix from mimes n number of linear measurements. The algorithms are first interpreted as iterative hard thresholding algorithms with subspace projections. Based on this connection, we show that provided the restricted isometry constant mimes n of the sensing operator is less than mimes n, the Riemannian gradient descent algorithm and a restarted variant of the Riemannian conjugate gradient algorithm are guaranteed to converge linearly to the underlying rank mimes n matrix if they are initialized by one step hard thresholding. Empirical evaluation shows that the algorithms are able to recover a low rank matrix from nearly the minimum number of measurements necessary.

In this paper, we study the rate of convergence of the cyclic projection algorithm applied to finitely many basic semialgebraic convex sets. We establish an explicit convergence rate estimate which relies on the maximum degree of the polynomials that generate the basic semialgebraic convex sets and the dimension of the underlying space. We achieve our results by exploiting the algebraic structure of the basic semialgebraic convex sets.