In this paper, we propose a novel low dimensional manifold model (LDMM) and apply it to some image processing problems. LDMM is based on the fact that the patch manifolds of many natural images have low dimensional structure. Based on this fact, the dimension of the patch manifold is used as a regularization to recover the image. The key step in LDMM is to solve a Laplace-Beltrami equation over a point cloud which is solved by the point integral method. The point integral method enforces the sample point constraints correctly and gives better results than the standard graph Laplacian. Numerical simulations in image denoising, inpainting and super-resolution problems show that LDMM is a powerful method in image processing.

In this work a simple method to enforce the positivity-preserving property for general high-order conservative schemes is proposed. The method detects critical numerical fluxes which may lead to negative density and pressure, and then imposes a simple flux limiter combining the high-order numerical flux with the first-order Lax-Friedrichs flux to satisfy a sufficient condition for preserving positivity. Though an extra time-step size condition is required to maintain the formal order of accuracy, it is less restrictive than those in previous works. A number of numerical examples suggest that this method, when applied on an essentially non-oscillatory scheme, can be used to prevent positivity failure when the flow involves vacuum or near vacuum and very strong discontinuities.

The positivity-preserving property is an important and challenging issue for the numerical solution of radiative transfer equations. In the past few decades, different numerical techniques have been proposed to guarantee positivity of the radiative intensity in several schemes, however it is difficult to maintain both high order accuracy and positivity. The discontinuous Galerkin (DG) finite element method is a high order numerical method which is widely used to solve the neutron/photon transfer equations, due to its distinguished advantages such as high order accuracy, geometric flexibility, suitability for $h$- and $p$-adaptivity, parallel efficiency, and a good theoretical foundation for stability and error estimates. In this paper, we construct arbitrarily high order accurate DG schemes which preserve positivity of the radiative intensity in the simulation of both steady and unsteady radiative transfer equations in one- and two-dimensional geometry by using a combined technique of the scaling positivity-preserving limiter in \cite{ZS2} and a new rotational positivity-preserving limiter. This combined limiter is simple to implement and we prove the properties of positivity-preserving and high order accuracy rigorously. One- and two-dimensional numerical results are provided to verify the good properties of the positivity-preserving DG schemes.

Motivated by big data applications, first-order methods have been extremely popular in recent years. However, naive gradient methods generally converge slowly. Hence, much efforts have been made to accelerate various first-order methods. This paper proposes two accelerated methods towards solving structured linearly constrained convex programming, for which we assume composite convex objective that is the sum of a differentiable function and a possibly nondifferentiable one. The first method is the accelerated linearized augmented Lagrangian method (LALM). At each update to the primal variable, it allows linearization to the differentiable function and also the augmented term, and thus it enables easy subproblems. Assuming merely convexity, we show that LALM owns O(1/t) convergence if parameters are kept fixed during all the iterations and can be accelerated to O(1/t^2) if the parameters are adapted, where t is the number of total iterations. The second method is the accelerated linearized alternating direction method of multipliers (LADMM). In addition to the composite convexity, it further assumes two-block structure on the objective. Different from classic ADMM, our method allows linearization to the objective and also augmented term to make the update simple. Assuming strong convexity on one block variable, we show that LADMM also enjoys O(1=t2) convergence with adaptive parameters. This result is a significant improvement over that in [Goldstein et. al, SIIMS'14], which requires strong convexity on both block variables and no linearization to the objective or augmented term. Numerical experiments are performed on quadratic programming, image denoising, and support vector machine. The proposed accelerated methods are compared to nonaccelerated ones and also existing accelerated methods. The results demonstrate the validness of acceleration and superior performance of the proposed methods over existing ones.

In this paper, we develop, analyze and test the Fourier spectral methods for solving the Degasperis-Procesi(DP) equation which contains nonlinear high order derivatives, and possibly discontinuous or sharp transition solutions. The $L^2$ stability is obtained for general numerical solutions of the Fourier Galerkin method and Fourier collocation (pseudospectral) method.By applying the Gegenbauer reconstruction technique as a post-processing method to the Fourier spectral solution, we reduce the oscillations arising from the discontinuity successfully. The numerical simulation results for different types of solutions of the nonlinear Degasperis-Procesi equation are provided to illustrate the accuracy and capability of the methods.