Variational models with \ell_1-norm based regularization, in particular total variation (TV) and its variants, have long been known to offer superior image restoration quality, but processing speed remained a bottleneck, preventing their widespread use in the practice of color image processing. In this paper, by extending the grayscale image deblurring algorithm proposed in [Y. Wang, J. Yang, W. Yin, and Y. Zhang, <i>SIAM J. Imaging Sci.</i>, 1 (2008), pp. 248272], we construct a simple and efficient algorithm for multichannel image deblurring and denoising, applicable to both within-channel and cross-channel blurs in the presence of additive Gaussian noise. The algorithm restores an image by minimizing an energy function consisting of an \ell_1-norm fidelity term and a regularization term that can be either TV, weighted TV, or regularization functions based on higher-order derivatives. Specifically, we use a multichannel

No abstract uploaded!

Solutions for many problems of interest exhibit singular behaviors at domain corners or points where boundary condition changes type. For this type of problems, direct spectral methods with usual polynomial basis functions do not lead to a satisfactory convergence rate. We develop in this paper a M\"untz-Galerkin method which is based on specially tuned M\"untz polynomials to deal with the singular behaviors of the underlying problems. By exploring the relations between Jacobi polynomials and M\"untz polynomials, we develop efficient implementation procedures for the M\"untz-Galerkin method, and provide optimal error estimates. As examples of applications, we consider the Poisson equation with mixed Dirichlet-Neumann boundary conditions, whose solution behaves like $O(r^{1/2})$ near the singular point, and demonstrate that the M\"untz-Galerkin method greatly improves the rates of convergence of the usual spectral method.

High order path-conservative schemes have been developed for solving nonconservative hyperbolic systems. Recently, it has been observed that this approach may have some computational issues and shortcomings. In this paper, a modification to the high order path-conservative scheme in is proposed to improve its computational performance and to overcome some of the shortcomings. This modification is based on the high order finite volume WENO scheme with subcell resolution and it uses an exact Riemann solver to catch the right paths at the discontinuities. An application to one-dimensional compressible two-medium flows of nonconservative or primitive Euler equations is carried out to show the effectiveness of this new approach.

In this paper, we compute the stationary states of the multicomponent phase-field crystal model by formulating it as a block constrained minimization problem. The original infinite-dimensional non-convex minimization problem is approximated by a finite-dimensional constrained non-convex minimization problem after an appropriate spatial discretization. To efficiently solve the above optimization problem, we propose a so-called adaptive block Bregman proximal gradient (AB-BPG) algorithm that fully exploits the problem's block structure. The proposed method updates each order parameter alternatively, and the update order of blocks can be chosen in a deterministic or random manner. Besides, we choose the step size by developing a practical linear search approach such that the generated sequence either keeps energy dissipation or has a controllable subsequence with energy dissipation. The convergence property of the proposed method is established without the requirement of global Lipschitz continuity of the derivative of the bulk energy part by using the Bregman divergence. The numerical results on computing stationary ordered structures in binary, ternary, and quinary component coupled-mode Swift-Hohenberg models have shown a significant acceleration over many existing methods.