No abstract uploaded!

No abstract uploaded!

We consider an elliptic pseudo-differential equation with a highly oscillating linear potential modeled as a stationary ergodic random field. The random field is a function composed with a centered long-range correlated Gaussian process. In the limiting of vanishing correlation length, the heterogeneous solution converges to a deterministic solution obtained by averaging the random potential. We characterize the deterministic and stochastic correctors. With proper rescaling, the mean-zero stochastic corrector converges to a Gaussian random process in probability and weakly in the spatial variables. In addition, for two prototype equations involving the Laplacian and the fractional Laplacian operators, we prove that the limit holds in distribution in some Hilbert spaces. We also determine the size of the deterministic corrector when it is larger than the stochastic corrector. Depending on the correlation structure of the random field and on the singularities of the Green’s function, we show that either the deterministic or the random part of the corrector dominates.

In this paper, we consider the harmonic extension problem, which is widely used in many applications of machine learning. We formulate the harmonic extension as solving a Laplace-Beltrami equation with Dirichlet boundary condition. We use the point integral method (PIM) to solve the Laplace-Beltrami equation. The basic idea of the PIM method is to approximate the Laplace equation using an integral equation, which is easy to be discretized from points. Based on the integral equation, we found that traditional graph Laplacian method (GLM) may fail to approximate the harmonic functions in the classical sense. For the Laplace-Beltrami equation with Dirichlet boundary, we can prove the convergence of the point integral method. The point integral method is also very easy to implement, which only requires a minor modification of the graph Laplacian. One important application of the harmonic extension in machine learning is semi-supervised learning. We run a popular semi-supervised learning algorithm by Zhu et al. over a couple of well-known datasets and compare the performance of the aforementioned approaches. Our experiments show the PIM performs the best. We also apply PIM to an image recovery problem and show it outperforms GLM. Finally, on a model problem of Laplace-Beltrami equation with Dirichlet boundary, we prove the convergence of the point integral method.

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.