We study analytical and numerical properties of the L 1 - L 2 minimization problem for sparse representation of a signal over a highly coherent dictionary. Though the L 1 - L 2 metric is non-convex, it is Lipschitz continuous. The difference of convex algorithm (DCA) is readily applicable for computing the sparse representation coefficients. The L 1 - L 2 minimization appears as an initialization step of DCA. We further integrate DCA with a non-standard simulated annealing methodology to approximate globally sparse solutions. Non-Gaussian random perturbations are more effective than standard Gaussian perturbations for improving sparsity of solutions. In numerical experiments, we conduct an extensive comparison among sparse penalties such as L 1 - L 2 for L 1 - L 2 based on data from three specific applications (over-sampled discreet cosine basis, differential absorption optical spectroscopy, and image denoising) where

In this paper, we generalize the point integral method to solve the general elliptic PDEs with variable coefficients and corresponding eigenvalue problems with Neumann, Robin and Dirichlet boundary conditions on point cloud. The main idea is using integral equations to approximate the original PDEs. The integral equations are easy to discretize on the point cloud. The truncation error of the integral approximation is analyzed. Numerical examples are presented to demonstrate that PIM is an effective method to solve the elliptic PDEs with smooth coefficients on point cloud.

No abstract uploaded!

We study a robust and efficient eigensolver for computing a few smallest positive eigenvalues of the three-dimensional Maxwell’s transmission eigenvalue problem. The discretized governing equations by the Nedelec edge element result in a large-scale quadratic eigenvalue problem (QEP) for which the spectrum contains many zero eigenvalues and the coefficient matrices consist of patterns in the matrix form XY^{−1}Z, both of which prevent existing eigenvalue solvers from being efficient. To remedy these difficulties, we rewrite the QEP as a particular nonlinear eigenvalue problem and develop a secant-type iteration, together with an indefinite locally optimal block preconditioned conjugate gradient method (LOBPCG), to sequentially compute the desired positive eigenvalues. Furthermore, we propose a novel method to solve the linear systems in each iteration of LOBPCG. Intensive numerical experiments show that our proposed method is robust, although the desired real eigenvalues are surrounded by complex eigenvalues.

An important property for finite difference schemes designed on curvilinear meshes is the exact preservation of free-stream solutions. This property is difficult to fulfill for high order conservative essentially non-oscillatory (WENO) finite difference schemes. In this paper we explore an alternative flux formulation for such finite difference schemes which can preserve free-stream solutions, based on the numerical technique for the metric terms, which can be applied to this alternative flux formulation but is difficult to be applied to the standard finite difference formulation. Free-stream and vortex preservation properties are investigated, and comparison with standard finite difference WENO schemes is made. Theoretical derivation and numerical results show that the finite difference WENO schemes based on the alternative flux formulation can preserve free-stream and vortex solutions on both stationary and dynamically generalized coordinate systems, hence giving much better performance than the standard finite difference WENO schemes for such problems.