亚洲男人的天堂2018av,欧美草比,久久久久久免费视频精选,国色天香在线看免费,久久久久亚洲av成人片仓井空

This paper proposes a regularization of the Monge-Amp\`ere equation in planar convex domains through uniformly elliptic Hamilton-Jacobi-Bellman equations. The regularized problem possesses a unique strong solution $u_\varepsilon$ and is accessible to the discretization with finite elements. This work establishes locally uniform convergence of $u_\varepsilon$ to the convex Alexandrov solution $u$ to the Monge-Amp\`ere equation as the regularization parameter $\varepsilon$ approaches $0$. A mixed finite element method for the approximation of $u_\varepsilon$ is proposed, and the regularized finite element scheme is shown to be locally uniformly convergent. Numerical experiments provide empirical evidence for the efficient approximation of singular solutions $u$.

When expressed in Lagrangian variables, the equations of motion for compressible (barotropic) fluids have the structure of a classical Hamiltonian system in which the potential energy is given by the internal energy of the fluid. The dissipative counterpart of such a system coincides with the porous medium equation, which can be cast in the form of a gradient flow for the same internal energy. Motivated by these related variational structures, we propose a particle method for both problems in which the internal energy is replaced by its Moreau-Yosida regularization in the L2 sense, which can be efficiently computed as a semi-discrete optimal transport problem. Using a modulated energy argument which exploits the convexity of the problem in Eulerian variables, we prove quantitative convergence estimates towards smooth solutions. We verify such estimates by means of several numerical tests.

This paper focuses on the construction and analysis of explicit numerical methods of high dimensional stochastic nonlinear Schrodinger equations (SNLSEs). We first prove that the classical explicit numerical methods are unstable and suffer from the numerical divergence phenomenon. Then we propose a kind of explicit splitting numerical methods and prove that the structure-preserving splitting strategy is able to enhance the numerical stability. Furthermore, we establish the regularity analysis and strong convergence analysis of the proposed schemes for SNLSEs based on two key ingredients. One ingredient is proving new regularity estimates of SNLSEs by constructing a logarithmic auxiliary functional and exploiting the Bourgain space. Another one is providing a dedicated error decomposition formula and a novel truncated stochastic Gronwall's lemma, which relies on the tail estimates of underlying stochastic processes. In particular, our result answers the strong convergence problem of numerical methods for 2D SNLSEs emerged from [C. Chen, J. Hong and A. Prohl, Stoch. Partial Differ. Equ. Anal. Comput. 4 (2016), no. 2, 274-318] and [J. Cui and J. Hong, SIAM J. Numer. Anal. 56 (2018), no. 4, 2045-2069].

In this paper, we extend the work of Brenner and Sung [Math. Comp. 59, 321--338 (1992)] and present a regularity estimate for the elastic equations in concave domains. Based on the regularity estimate we prove that the constants in the error estimates of the nonconforming Crouzeix-Raviart element approximations for the elastic equations/eigenvalue problem are independent of the Lame constant, which means the nonconforming Crouzeix-Raviart element approximations are locking-free. We also establish two kinds of two-grid discretization schemes for the elastic eigenvalue problem and analyze that when the mesh sizes of the coarse grid and fine grid satisfy some relationship, the resulting solutions can achieve optimal accuracy. Numerical examples are provided to show the efficiency of two-grid schemes for the elastic eigenvalue problem.

A recent paper [J. A. Evans, D. Kamensky, Y. Bazilevs, "Variational multiscale modeling with discretely divergence-free subscales", Computers & Mathematics with Applications, 80 (2020) 2517-2537] introduced a novel stabilized finite element method for the incompressible Navier-Stokes equations, which combined residual-based stabilization of advection, energetic stability, and satisfaction of a discrete incompressibility condition. However, the convergence analysis and numerical tests of the cited work were subject to the restrictive assumption of a divergence-conforming choice of velocity and pressure spaces, where the pressure space must contain the divergence of every velocity function. The present work extends the convergence analysis to arbitrary inf-sup-stable velocity-pressure pairs (while maintaining robustness in the advection-dominated regime) and demonstrates the convergence of the method numerically, using both the traditional and isogeometric Taylor-Hood elements.

Stochastic variance reduced gradient (SVRG) is a popular variance reduction technique for accelerating stochastic gradient descent (SGD). We provide a first analysis of the method for solving a class of linear inverse problems in the lens of the classical regularization theory. We prove that for a suitable constant step size schedule, the method can achieve an optimal convergence rate in terms of the noise level (under suitable regularity condition) and the variance of the SVRG iterate error is smaller than that by SGD. These theoretical findings are corroborated by a set of numerical experiments.

In neuroscience, the distribution of a decision time is modelled by means of a one-dimensional Fokker--Planck equation with time-dependent boundaries and space-time-dependent drift. Efficient approximation of the solution to this equation is required, e.g., for model evaluation and parameter fitting. However, the prescribed boundary conditions lead to a strong singularity and thus to slow convergence of numerical approximations. In this article we demonstrate that the solution can be related to the solution of a parabolic PDE on a rectangular space-time domain with homogeneous initial and boundary conditions by transformation and subtraction of a known function. We verify that the solution of the new PDE is indeed more regular than the solution of the original PDE and proceed to discretize the new PDE using a space-time minimal residual method. We also demonstrate that the solution depends analytically on the parameters determining the boundaries as well as the drift. This justifies the use of a sparse tensor product interpolation method to approximate the PDE solution for various parameter ranges. The predicted convergence rates of the minimal residual method and that of the interpolation method are supported by numerical simulations.

Additional grid points are often introduced for the higher-order polynomial of a numerical solution with curvilinear elements. However, those points are likely to be located slightly outside the domain, even when the vertices of the curvilinear elements lie within the curved domain. This misallocation of grid points generates a mesh error, called geometric approximation error. This error is smaller than the discretization error but large enough to significantly degrade a long-time integration. Moreover, this mesh error is considered to be the leading cause of conservation error. Two novel schemes are proposed to improve conservation error and/or discretization error for long-time integration caused by geometric approximation error: The first scheme retrieves the original divergence of the original domain; the second scheme reconstructs the original path of differentiation, called connection, thus retrieving the original connection. The increased accuracies of the proposed schemes are demonstrated by the conservation error for various partial differential equations with moving frames on the sphere.

In this paper, we study an adaptive finite element method for the elliptic equation with line Dirac delta functions as a source term.We investigate the regularity of the solution and the corresponding transmission problem to obtain the jump of normal derivative of the solution on line fractures. To handle the singularity of the solution, we adopt the meshes that conform to line fractures, and propose a novel a posteriori error estimator, in which the edge jump residual essentially use the jump of the normal derivative of the solution on line fractures. The error estimator is proven to be both reliable and efficient, finally an adaptive finite element algorithm is proposed based on the error estimator and the bisection refinement method. Numerical tests are presented to justify the theoretical findings.

We consider the task of learning the parameters of a {\em single} component of a mixture model, for the case when we are given {\em side information} about that component, we call this the "search problem" in mixture models. We would like to solve this with computational and sample complexity lower than solving the overall original problem, where one learns parameters of all components. Our main contributions are the development of a simple but general model for the notion of side information, and a corresponding simple matrix-based algorithm for solving the search problem in this general setting. We then specialize this model and algorithm to four common scenarios: Gaussian mixture models, LDA topic models, subspace clustering, and mixed linear regression. For each one of these we show that if (and only if) the side information is informative, we obtain parameter estimates with greater accuracy, and also improved computation complexity than existing moment based mixture model algorithms (e.g. tensor methods). We also illustrate several natural ways one can obtain such side information, for specific problem instances. Our experiments on real data sets (NY Times, Yelp, BSDS500) further demonstrate the practicality of our algorithms showing significant improvement in runtime and accuracy.