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

We propose a third-order numerical integrator based on the Neumann series and the Filon quadrature, designed mainly for highly oscillatory partial differential equations. The method can be applied to equations that exhibit small or moderate oscillations; however, counter-intuitively, large oscillations increase the accuracy of the scheme. With the proposed approach, the convergence order of the method can be easily improved. Error analysis of the method is also performed. We consider linear evolution equations involving first- and second-time derivatives that feature elliptic differential operators, such as the heat equation or the wave equation. Numerical experiments consider the case in which the space dimension is greater than one and confirm the theoretical study.

The fractional differential equation $L^\beta u = f$ posed on a compact metric graph is considered, where $\beta>0$ and $L = \kappa^2 - \nabla(a\nabla)$ is a second-order elliptic operator equipped with certain vertex conditions and sufficiently smooth and positive coefficients $\kappa, a$. We demonstrate the existence of a unique solution for a general class of vertex conditions and derive the regularity of the solution in the specific case of Kirchhoff vertex conditions. These results are extended to the stochastic setting when $f$ is replaced by Gaussian white noise. For the deterministic and stochastic settings under generalized Kirchhoff vertex conditions, we propose a numerical solution based on a finite element approximation combined with a rational approximation of the fractional power $L^{-\beta}$. For the resulting approximation, the strong error is analyzed in the deterministic case, and the strong mean squared error as well as the $L_2(\Gamma\times \Gamma)$-error of the covariance function of the solution are analyzed in the stochastic setting. Explicit rates of convergences are derived for all cases. Numerical experiments for ${L = \kappa^2 - \Delta, \kappa>0}$ are performed to illustrate the results.

Learning nonparametric systems of Ordinary Differential Equations (ODEs) dot x = f(t,x) from noisy data is an emerging machine learning topic. We use the well-developed theory of Reproducing Kernel Hilbert Spaces (RKHS) to define candidates for f for which the solution of the ODE exists and is unique. Learning f consists of solving a constrained optimization problem in an RKHS. We propose a penalty method that iteratively uses the Representer theorem and Euler approximations to provide a numerical solution. We prove a generalization bound for the L2 distance between x and its estimator and provide experimental comparisons with the state-of-the-art.

Perturbation theory plays a crucial role in sensitivity analysis, which is extensively used to assess the robustness of numerical techniques. To quantify the relative sensitivity of any problem, it becomes essential to investigate structured condition numbers (CNs) via componentwise perturbation theory. This paper addresses and analyzes structured mixed condition number (MCN) and componentwise condition number (CCN) for the Moore-Penrose (M-P) inverse and the minimum norm least squares (MNLS) solution involving rank-structured matrices, which include the Cauchy-Vandermonde (CV) matrices and $\{1,1\}$-quasiseparable (QS) matrices. A general framework has been developed to compute the upper bounds for MCN and CCN of rank deficient parameterized matrices. This framework leads to faster computation of upper bounds of structured CNs for CV and $\{1,1\}$-QS matrices. Furthermore, comparisons of obtained upper bounds are investigated theoretically and experimentally. In addition, the structured effective CNs for the M-P inverse and the MNLS solution of $\{1,1\}$-QS matrices are presented. Numerical tests reveal the reliability of the proposed upper bounds as well as demonstrate that the structured effective CNs are computationally less expensive and can be substantially smaller compared to the unstructured CNs.

This work develops an energy-based discontinuous Galerkin (EDG) method for the nonlinear Schr\"odinger equation with the wave operator. The focus of the study is on the energy-conserving or energy-dissipating behavior of the method with some simple mesh-independent numerical fluxes we designed. We establish error estimates in the energy norm that require careful selection of a test function for the auxiliary equation involving the time derivative of the displacement variable. A critical part of the convergence analysis is to establish the L2 error bounds for the time derivative of the approximation error in the displacement variable by using the equation that determines its mean value. Using a specially chosen test function, we show that one can create a linear system for the time evolution of the unknowns even when dealing with nonlinear properties in the original problem. Extensive numerical experiments are provided to demonstrate the optimal convergence of the scheme in the L2 norm with our choices of the numerical flux.

This paper presents an efficient method for obtaining the least squares Hermitian solutions of the reduced biquaternion matrix equation $(AXB, CXD) = (E, F )$. The method leverages the real representation of reduced biquaternion matrices. Furthermore, we establish the necessary and sufficient conditions for the existence and uniqueness of the Hermitian solution, along with a general expression for it. Notably, this approach differs from the one previously developed by Yuan et al. $(2020)$, which relied on the complex representation of reduced biquaternion matrices. In contrast, our method exclusively employs real matrices and utilizes real arithmetic operations, resulting in enhanced efficiency. We also apply our developed framework to find the Hermitian solutions for the complex matrix equation $(AXB, CXD) = (E, F )$, expanding its utility in addressing inverse problems. Specifically, we investigate its effectiveness in addressing partially described inverse eigenvalue problems. Finally, we provide numerical examples to demonstrate the effectiveness of our method and its superiority over the existing approach.

This paper studies two hybrid discontinuous Galerkin (HDG) discretizations for the velocity-density formulation of the compressible Stokes equations with respect to several desired structural properties, namely provable convergence, the preservation of non-negativity and mass constraints for the density, and gradient-robustness. The later property dramatically enhances the accuracy in well-balanced situations, such as the hydrostatic balance where the pressure gradient balances the gravity force. One of the studied schemes employs an H(div)-conforming velocity ansatz space which ensures all mentioned properties, while a fully discontinuous method is shown to satisfy all properties but the gradient-robustness. Also higher-order schemes for both variants are presented and compared in three numerical benchmark problems. The final example shows the importance also for non-hydrostatic well-balanced states for the compressible Navier-Stokes equations.

We present a space-time ultra-weak discontinuous Galerkin discretization of the linear Schr\"odinger equation with variable potential. The proposed method is well-posed and quasi-optimal in mesh-dependent norms for very general discrete spaces. Optimal $h$-convergence error estimates are derived for the method when test and trial spaces are chosen either as piecewise polynomials, or as a novel quasi-Trefftz polynomial space. The latter allows for a substantial reduction of the number of degrees of freedom and admits piecewise-smooth potentials. Several numerical experiments validate the accuracy and advantages of the proposed method.

Solutions to many important partial differential equations satisfy bounds constraints, but approximations computed by finite element or finite difference methods typically fail to respect the same conditions. Chang and Nakshatrala enforce such bounds in finite element methods through the solution of variational inequalities rather than linear variational problems. Here, we provide a theoretical justification for this method, including higher-order discretizations. We prove an abstract best approximation result for the linear variational inequality and estimates showing that bounds-constrained polynomials provide comparable approximation power to standard spaces. For any unconstrained approximation to a function, there exists a constrained approximation which is comparable in the $W^{1,p}$ norm. In practice, one cannot efficiently represent and manipulate the entire family of bounds-constrained polynomials, but applying bounds constraints to the coefficients of a polynomial in the Bernstein basis guarantees those constraints on the polynomial. Although our theoretical results do not guaruntee high accuracy for this subset of bounds-constrained polynomials, numerical results indicate optimal orders of accuracy for smooth solutions and sharp resolution of features in convection-diffusion problems, all subject to bounds constraints.

A numerical method is proposed for simulation of composite open quantum systems. It is based on Lindblad master equations and adiabatic elimination. Each subsystem is assumed to converge exponentially towards a stationary subspace, slightly impacted by some decoherence channels and weakly coupled to the other subsystems. This numerical method is based on a perturbation analysis with an asymptotic expansion. It exploits the formulation of the slow dynamics with reduced dimension. It relies on the invariant operators of the local and nominal dissipative dynamics attached to each subsystem. Second-order expansion can be computed only with local numerical calculations. It avoids computations on the tensor-product Hilbert space attached to the full system. This numerical method is particularly well suited for autonomous quantum error correction schemes. Simulations of such reduced models agree with complete full model simulations for typical gates acting on one and two cat-qubits (Z, ZZ and CNOT) when the mean photon number of each cat-qubit is less than 8. For larger mean photon numbers and gates with three cat-qubits (ZZZ and CCNOT), full model simulations are almost impossible whereas reduced model simulations remain accessible. In particular, they capture both the dominant phase-flip error-rate and the very small bit-flip error-rate with its exponential suppression versus the mean photon number.