A time-dependent Hermite-Galerkin spectral method (THGSM) is investigated in this paper for the nonlinear convection-diffusion equations in the unbounded domains. The time-dependent scaling factor and translating factor are introduced in the definition of the generalized Hermite functions (GHF). As a consequence, the THGSM based on these GHF has many advantages, not only in theoretical proofs, but also in numerical implementations. The stability and spectral convergence of our proposed method have been established in this paper. The Korteweg-de Vries-Burgers (KdVB) equation and its special cases, including the heat equation and the Burgers' equation, as the examples, have been numerically solved by our method. The numerical results are presented, and it surpasses the existing methods in accuracy. Our theoretical proof of the spectral convergence has been supported by the numerical results.

No abstract uploaded!

Surface based 3D shape analysis plays a fundamental role in computer vision and medical imaging. This work proposes to use optimal mass transport map for shape matching and comparison, focusing on two important applications including surface registration and shape space. The computation of the optimal mass transport map is based on Monge-Brenier theory, in comparison to the conventional method based on Monge-Kantorovich theory, this method significantly improves the efficiency by reducing computational complexity from O(n ² ) to O(n). For surface registration problem, one commonly used approach is to use conformal map to convert the shapes into some canonical space. Although conformal mappings have small angle distortions, they may introduce large area distortions which are likely to cause numerical instability thus resulting failures of shape analysis. This work proposes to compose the

We study the minimization problem of a non-convex sparsity promoting penalty function, the transformed l_1 (TL1), and its application in compressed sensing (CS). The TL1 penalty interpolates l_1 and l_1 norms through a nonnegative parameter l_1 , similar to l_1 with l_1 , and is known to satisfy unbiasedness, sparsity and Lipschitz continuity properties. We first consider the constrained minimization problem and discuss the exact recovery of l_1 norm minimal solution based on the null space property (NSP). We then prove the stable recovery of l_1 norm minimal solution if the sensing matrix l_1 satisfies a restricted isometry property (RIP). Next, we present difference of convex algorithms for TL1 (DCATL1) in computing TL1-regularized constrained and unconstrained problems in CS. The inner loop concerns an l_1 minimization problem on which we employ the Alternating Direction Method of Multipliers (ADMM). For the unconstrained problem, we prove convergence of DCATL1 to a stationary point satisfying the first order optimality condition. In numerical experiments, we identify the optimal value l_1 , and compare DCATL1 with other CS algorithms on two classes of sensing matrices: Gaussian random matrices and over-sampled discrete cosine transform matrices (DCT). We find that for both classes of sensing matrices, the performance of DCATL1 algorithm (initiated with l_1 minimization) always ranks near the top (if not the top), and is the most robust choice insensitive to the conditioning of the sensing matrix l_1 . DCATL1 is also competitive in comparison with DCA on other non-convex penalty functions commonly used in statistics with two

What is the shape of a uniformly massive object that generates a gravitational potential equivalent to that of two equal point-masses? If the weight of each point-mass is sufficiently small compared to the distance between the points then the answer is a pair of balls of equal radius, one centered at each of the two points, but otherwise it is a certain domain of revolution about the axis passing through the two points. The existence and uniqueness of such a domain is known, but an explicit parameterization is known only in the plane where the region is referred to as a Neumann oval. We construct a four-dimensional “Neumann ovaloid”, solving explicitly this inverse potential problem.