In this paper, we explore a unified way to construct smoothing functions for solving the absolute value equation associated with second-order cone (SOCAVE). Numerical comparisons are presented, which illustrate what kinds of smoothing functions work well along with the smoothing Newton algorithm. In particular, the numerical experiments show that the well known loss function widely used in engineering community is the worst one among the constructed smoothing functions, which indicates that the other proposed smoothing functions can be employed for solving engineering problems.

The main purpose of this paper is to give stability analysis and error estimates of the local discontinuous Galerkin (LDG) methods coupled with three specific implicit-explicit (IMEX) Runge-Kutta time discretization methods up to third order accuracy, for solving one-dimensional time-dependent linear fourth order partial differential equations. In the time discretization, all the lower order derivative terms are treated explicitly and the fourth order derivative term is treated implicitly. By the aid of energy analysis, we show that the IMEX-LDG schemes are unconditionally energy stable, in the sense that the time step $\dt$ is only required to be upper-bounded by a constant which is independent of the mesh size $h$. The optimal error estimate is also derived by the aid of the elliptic projection and the adjoint argument. Numerical experiments are given to verify that the corresponding IMEX-LDG schemes can achieve optimal error accuracy.

Unmixing problems in many areas such as hyperspectral imaging and differential optical absorption spectroscopy (DOAS) often require finding sparse nonnegative linear combinations of dictionary elements that match observed data. We show how aspects of these problems, such as misalignment of DOAS references and uncertainty in hyperspectral endmembers, can be modeled by expanding the dictionary with grouped elements and imposing a structured sparsity assumption that the combinations within each group should be sparse or even 1-sparse. If the dictionary is highly coherent, it is difficult to obtain good solutions using convex or greedy methods, such as nonnegative least squares (NNLS) or orthogonal matching pursuit. We use penalties related to the Hoyer measure, which is the ratio of the l_1 and l_1 norms, as sparsity penalties to be added to the objective in NNLS-type models. For solving the

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.

We introduce a new iterative regularization procedure for inverse problems based on the use of Bregman distances, with particular focus on problems arising in image processing. We are motivated by the problem of restoring noisy and blurry images via variational methods by using total variation regularization. We obtain rigorous convergence results and effective stopping criteria for the general procedure. The numerical results for denoising appear to give significant improvement over standard models, and preliminary results for deblurring/denoising are very encouraging.